The utility of human two plus one small pronucleated zygotes (2.1PN) based on clinical outcomes and the focused ploidy analysis.

PURPOSE: Are human embryos arising from two plus one small pronucleated zygotes, called 2.1 pronuclei (PN), clinically useful? METHODS: In a retrospective embryo cohort study and prospective experimental study, a total of 287 cycles in which at least one 2.1PN was identified in the fertilization check were included. Embryonic development and clinical outcome were compared for the 1395 2PN zygotes and 304 2.1PN zygotes that were siblings. All embryos were individually cultured in time-lapse systems. Twenty-five 2.1PN-derived blastocysts, donated for research, were used in focused single-nucleotide variant ploidy analysis to identify the distribution pattern of heterozygosity. RESULTS: The average diameter of PN was 24.9 ± 2.4 µm for large PN and 10.2 ± 2.4 µm for small PN; 79.9% of small PN was derived from female pronuclei. Blastocyst formation rate and good-quality blastocyst rate were significantly lower with 2.1PN embryos than with 2PN embryos (40.0% vs. 57.7%, 21.4% vs. 33.5%, respectively). A total of 13 embryos derived from 2.1PN were transferred, and three healthy babies were born. In ploidy constitutions of trophectoderm (TE), 2.1PN-derived blastocyst TE was shown to be mostly diploid (95.8%, 23/24), and only one blastocyst showed triploid. CONCLUSIONS: It was suggested that 2.1PN embryos have lower embryonic developmental potential than 2PN embryos, but most of the 2.1PN were diploid, indicating that they are likely to be clinically usable. It is recommended to perform embryo transfer following a combination of PGT-A and ploidy analysis.

Vaginal Microbiota and Pregnancy Outcomes of Patients with Conization Histories.

Background: One of the major risks of preterm birth is a history of conization. However, the risk of infection due to this procedure is still not well known. Using next-generation sequencing, we aimed to reveal the influence of conization on vaginal microbiota in the following pregnancy, and their relationship between spontaneous preterm birth (sPTB). Methods: We conducted a prospective cohort study including 133 pregnant patients, of whom 25 had conization histories and 108 did not. Vaginal microbiome samples were collected using swabs by an obstetrician upon inclusion in the first trimester and during delivery. V1-V2 of the 16S rRNA gene were amplified and analyzed to identify the bacteria. Results: The conization group had a significantly lower delivery week (34 weeks vs. 36 weeks, p = 0.003) and higher sPTB rate (64% vs. 8.3%, p ≤ 0.001) than the control group. In the conization group, alpha (Chao 1, p = 0.02; phylogenetic diversity whole tree, p = 0.04) and beta diversity (permutational multivariate analysis of variance test, p = 0.04) of the vaginal microbiota was significantly higher during delivery in patients who delivered preterm than in those who delivered term. Community-state type IV in the first trimester was significantly associated with sPTB (overall odds ratio 3.80, 95% confidence interval 1.33-10.8, p = 0.01). Conclusions: Conization is a risk factor for sPTB. Increased risk of sPTB in patients after conization may belong to the vulnerable defense mechanism, due to the shortened cervix and decreased cervical mucus.

Endometrial and vaginal microbiomes influence assisted reproductive technology outcomes.

OBJECTIVE: The role of Lactobacillus-dominant microbiota in the endometrium in reproductive function is unclear. We therefore aimed to explore the impact of the balance of Lactobacillus and pathological bacteria in the endometrial and vaginal microbiomes on the pregnancy outcomes of women treated with assisted reproductive technology (ART). METHODS: This study included 35 women with infertility submitted to good-quality embryo transfers. The cutoff values for abundance of Lactobacillus species (spp.) and pathological bacteria in the endometrium and vagina were calculated. Women with Lactobacillus spp. and pathological bacteria abundance above the cutoff values were categorized in the high-abundance group, whereas those with abundance below cutoff values were categorized in the low abundance group. We divided the patients into four groups based on the combination of high/low abundance of Lactobacillus spp. and pathological bacteria. RESULTS: The 35 cases of good-quality embryo transfer resulted in 21 pregnancies. Pregnant women were present in significantly higher proportions in the high Lactobacillus spp. abundance and low pathological bacteria abundance group, whereas the opposite combination (i.e., low Lactobacillus spp. abundance and high pathological bacteria abundance) saw a significantly higher proportion of nonpregnant women (p=0.022). CONCLUSIONS: The balance between Lactobacillus and pathological bacterial abundance in the endometrial and vaginal microbiomes is associated with pregnancy from ART.

Microbiome analysis in women with endometriosis: Does a microbiome exist in peritoneal fluid and ovarian cystic fluid?

Purpose: To investigate the relationship between the microbiome of the female genital tract and endometriosis. Methods: This prospective cohort study included 36 women who underwent laparoscopic surgery for ovarian tumor from July 2019 to April 2020. Of them, 18 had endometriosis, and 18 did not have endometriosis. Vaginal secretions, endometrial fluid, peritoneal fluid, and ovarian cystic fluid were collected during surgery. Next-generation sequencing of bacterial 16S rRNA was performed to characterize the microbiome. Results: Specific microbiomes were not detected in either peritoneal fluid or ovarian cystic fluid regardless of the presence or absence of endometriosis and the type of cyst. When the cutoff value of infectious bacterial abundance in the vagina was set as 64.3%, there were many cases more than a cutoff value in the endometriosis group significantly (p = 0.01). When the cutoff value of infectious bacterial abundance in the endometrium was set as 18.6%, there were many cases more than a cutoff level in the endometriosis cases significantly (p = 0.02). Conclusion: Peritoneal fluid and ovarian cystic fluid are almost sterile, although dysbiosis may occur in the vaginal and endometrial microbiome in women with endometriosis.

Next-generation sequencing of 16S rRNA for identification of invasive bacterial pathogens in a formalin-fixed paraffin-embedded placental specimen: a case report of perinatal fulminant Streptococcus pyogenes infection.

Intrauterine infection is one of the most important causes of maternal death. In perinatal emergency, we often miss an opportunity to obtain culture specimens. In this study, we tried to examine whether we investigated whether bacteria causing infection can be detected from a formalin-fixed paraffin-embedded (FFPE) placental specimen. We examined the placenta from a maternal invasive infection that resulted in infectious abortion at 18 weeks of gestation. The case was diagnosed by acute fever and abdominal pain, and the patient was cured after 3 weeks of intensive antimicrobial treatment. Four Streptococcus pyogenes strains were isolated from vaginal fluid and blood cultures of the patient. All of the strain types were emm1/ST28. We amplified the V1-V2 region of 16S rRNA from an FFPE placental specimen and sequencing was performed using a next-generation sequencer (NGS). Taxonomic analysis was then performed for sequenced data. We succeeded in detecting causative pathogens from the FFPE placenta: 69.1% of the predominantly identified bacteria were S. pyogenes and other small populations of bacteria were detected. Our results revealed the utility of NGS for 16S rRNA analysis of an FFPE placenta. This method may reveal previous perinatal invasive infections of unknown origin retrospectively.

Genome-wide association study of epilepsy in a Japanese population identified an associated region at chromosome 12q24.

OBJECTIVE: Although a number of genes responsible for epilepsy have been identified through Mendelian genetic approaches, and genome-wide association studies (GWASs) have implicated several susceptibility loci, the role of ethnic-specific markers remains to be fully explored. We aimed to identify novel genetic associations with epilepsy in a Japanese population. METHODS: We conducted a GWAS on 1825 patients with a variety of epilepsies and 7975 control individuals. Expression quantitative trait locus (eQTL) analysis of epilepsy-associated single nucleotide polymorphisms (SNPs) was performed using Japanese eQTL data. RESULTS: We identified a novel region, which is ~2 Mb (lead SNP rs149212747, p = 8.57 × 10-10 ), at chromosome 12q24 as a risk for epilepsy. Most of these loci were polymorphic in East Asian populations including Japanese, but monomorphic in the European population. This region harbors 24 transcripts including genes expressed in the brain such as CUX2, ATXN2, BRAP, ALDH2, ERP29, TRAFD1, HECTD4, RPL6, PTPN11, and RPH3A. The eQTL analysis revealed that the associated SNPs are also correlated to differential expression of genes at 12q24. SIGNIFICANCE: These findings suggest that a gene or genes in the CUX2-RPH3A ~2-Mb region contribute to the pathology of epilepsy in the Japanese population.

16S rRNA long-read nanopore sequencing is feasible and reliable for endometrial microbiome analysis.

RESEARCH QUESTION: Full-length 16S rRNA gene sequencing using nanopore technology is a fast alternative to conventional short-read 16S rRNA gene sequencing with low initial investment costs that has been used for various microbiome studies but has not yet been investigated as an alternative approach for endometrial microbiome analysis. Is in-situ 16S rRNA gene long-read sequencing using portable nanopore sequencing technology feasible and reliable for endometrial microbiome analysis? DESIGN: A prospective experimental study based on 33 patients seeking infertility treatment between January and October 2019. A 16S rRNA gene long-read nanopore sequencing protocol for analysing endometrial microbiome samples was established, including negative controls for contamination evaluation and positive controls for bias evaluation. Contamination caused by kit and exterior sources was identified and excluded from the analysis. Endometrial samples from 33 infertile patients were sequenced using the optimized long-read nanopore sequencing protocol and compared with conventional short-read sequencing carried out by external laboratories. RESULTS: Of the 33 endometrial patient samples, 23 successfully amplified (69.7%) and their microbiome was assessed using nanopore sequencing. Of those 23 samples, 14 (60.9%) were Lactobacillus-dominated (>80% of reads mapping to Lactobacillus), with 10 samples resulting in more than 90% Lactobacillus reads. Our long-read nanopore sequencing revealed results similar to two conventional short-read sequencing approaches and to long-read sequencing validation carried out in external laboratories. CONCLUSION: In this pilot study, 16S rRNA gene long-read nanopore sequencing was established to analyse the endometrial microbiome in situ that could be widely applied owing to its cost efficiency and portable character.

Prevalence and Spectrum of Pathogenic Germline Variants in Japanese Patients With Early-Onset Colorectal, Breast, and Prostate Cancer.

PUPOSE: We investigated the prevalence and spectrum of pathogenic germline variants in patients with early-onset colorectal cancer (CRC), breast cancer (BC), and prostate cancer (PCA) in the Japanese population. We also identified pathogenic variants in other cancer risk genes, giving consideration to future multigene testing panels for this population. METHODS: We performed whole-genome sequencing for 1,037 Japanese individuals, including patients with early-onset CRC (n = 196), BC (n = 237), and PCA (n = 215) and controls (n = 389). We screened for pathogenic variants, including single nucleotide variants and copy number variants, among well-established first-tier cancer genes for each cancer type and examined an expended second-tier panel including cancer-predisposing genes from the Cancer Gene Census. RESULTS: Proportions of patients with germline pathogenic variants differed by cancer subgroup, with the highest in BC (14.8%), followed by CRC (9.2%), and PCA (3.7%). In contrast, 2 of 389 control subjects (0.5%) carried a germline pathogenic variant. In comparison with controls, the proportion of patients with pathogenic variants in the second-tier panel was increased significantly for PCA (3.7% to 11.6%, P = 2.96 × 10-4), but not for CRC or BC, after multitesting adjustment. In patients with PCA, DNA repair pathway genes in the extended panel often contained pathogenic variants (P = .011). CONCLUSION: Our analyses support the clinical usefulness of established cancer gene panels in the Japanese population for 3 major cancer types. Additional genes, especially those involved in DNA repair, might be considered for developing multipanel testing in Japanese patients with early-onset PCA.

Tuberculosis infection and lung adenocarcinoma: Mendelian randomization and pathway analysis of genome-wide association study data from never-smoking Asian women.

We investigated whether genetic susceptibility to tuberculosis (TB) influences lung adenocarcinoma development among never-smokers using TB genome-wide association study (GWAS) results within the Female Lung Cancer Consortium in Asia. Pathway analysis with the adaptive rank truncated product method was used to assess the association between a TB-related gene-set and lung adenocarcinoma using GWAS data from 5512 lung adenocarcinoma cases and 6277 controls. The gene-set consisted of 31 genes containing known/suggestive associations with genetic variants from previous TB-GWAS. Subsequently, we followed-up with Mendelian Randomization to evaluate the association between TB and lung adenocarcinoma using three genome-wide significant variants from previous TB-GWAS in East Asians. The TB-related gene-set was associated with lung adenocarcinoma (p = 0.016). Additionally, the Mendelian Randomization showed an association between TB and lung adenocarcinoma (OR = 1.31, 95% CI: 1.03, 1.66, p = 0.027). Our findings support TB as a causal risk factor for lung cancer development among never-smoking Asian women.

Genome-wide association study suggests four variants influencing outcomes with ranibizumab therapy in exudative age-related macular degeneration.

To identify factors associated with ranibizumab responses in patients with exudative age-related macular degeneration (AMD), we performed a genome-wide association study (GWAS) and a replication study using a total of 919 exudative AMD patients treated with intravitreal ranibizumab in a Japanese population. In the combined analysis of GWAS and the replication study, no loci reached genome-wide significant level; however, we found four variants showed suggestive level of associations with visual loss at month three (rs17822656, rs76150532, rs17296444, and rs75165563: Pcombined < 1.0 × 10-5). Of the candidate genes within these loci, three were relevant to VEGF-related pathway (KCNMA1, SOCS2, and OTX2). The proportions of patients who worsened visual acuity were 13.7%, 38.8%, 58.0%, and 80.0% in patients with 0, 1, 2, and 3 or more identified risk variants, respectively. Changes in visual acuity decreased linearly as the number of risk variants increased (P = 1.67 × 10-12). The area under the curve using age, baseline visual acuity, and history of previous treatment was 0.607, and improved significantly to 0.713 in combination with identified variants (P < 0.0001). Although further study is needed to confirm their associations, our results offer candidate variants influencing response to ranibizumab therapy.

Association between GWAS-identified lung adenocarcinoma susceptibility loci and EGFR mutations in never-smoking Asian women, and comparison with findings from Western populations.

To evaluate associations by EGFR mutation status for lung adenocarcinoma risk among never-smoking Asian women, we conducted a meta-analysis of 11 loci previously identified in genome-wide association studies (GWAS). Genotyping in an additional 10,780 never-smoking cases and 10,938 never-smoking controls from Asia confirmed associations with eight known single nucleotide polymorphisms (SNPs). Two new signals were observed at genome-wide significance (P < 5 × 10-8), namely, rs7216064 (17q24.3, BPTF), for overall lung adenocarcinoma risk, and rs3817963 (6p21.3, BTNL2) which is specific to cases with EGFR mutations. In further sub-analyses by EGFR status, rs9387478 (ROS1/DCBLD1) and rs2179920 (HLA-DPB1) showed stronger estimated associations in EGFR-positive compared to EGFR-negative cases. Comparison of the overall associations with published results in Western populations revealed that the majority of these findings were distinct, underscoring the importance of distinct contributing factors for smoking and non-smoking lung cancer. Our results extend the catalogue of regions associated with lung adenocarcinoma in non-smoking Asian women and highlight the importance of how the germline could inform risk for specific tumour mutation patterns, which could have important translational implications.

Contribution of a Non-classical HLA Gene, HLA-DOA, to the Risk of Rheumatoid Arthritis.

Despite the progress in human leukocyte antigen (HLA) causal variant mapping, independent localization of major histocompatibility complex (MHC) risk from classical HLA genes is challenging. Here, we conducted a large-scale MHC fine-mapping analysis of rheumatoid arthritis (RA) in a Japanese population (6,244 RA cases and 23,731 controls) population by using HLA imputation, followed by a multi-ethnic validation study including east Asian and European populations (n = 7,097 and 23,149, respectively). Our study identified an independent risk of a synonymous mutation at HLA-DOA, a non-classical HLA gene, on anti-citrullinated protein autoantibody (ACPA)-positive RA risk (p = 1.4 × 10(-9)), which demonstrated a cis-expression quantitative trait loci (cis-eQTL) effect on HLA-DOA expression. Trans-ethnic comparison revealed different linkage disequilibrium (LD) patterns in HLA-DOA and HLA-DRB1, explaining the observed HLA-DOA variant risk heterogeneity among ethnicities, which was most evident in the Japanese population. Although previous HLA fine-mapping studies have identified amino acid polymorphisms of the classical HLA genes as driving genetic susceptibility to disease, our study additionally identifies the dosage contribution of a non-classical HLA gene to disease etiology. Our study contributes to the understanding of HLA immunology in human diseases and suggests the value of incorporating additional ancestry in MHC fine-mapping.

Association of variations in HLA class II and other loci with susceptibility to EGFR-mutated lung adenocarcinoma.

Lung adenocarcinoma driven by somatic EGFR mutations is more prevalent in East Asians (30-50%) than in European/Americans (10-20%). Here we investigate genetic factors underlying the risk of this disease by conducting a genome-wide association study, followed by two validation studies, in 3,173 Japanese patients with EGFR mutation-positive lung adenocarcinoma and 15,158 controls. Four loci, 5p15.33 (TERT), 6p21.3 (BTNL2), 3q28 (TP63) and 17q24.2 (BPTF), previously shown to be strongly associated with overall lung adenocarcinoma risk in East Asians, were re-discovered as loci associated with a higher susceptibility to EGFR mutation-positive lung adenocarcinoma. In addition, two additional loci, HLA class II at 6p21.32 (rs2179920; P =5.1 × 10(-17), per-allele OR=1.36) and 6p21.1 (FOXP4) (rs2495239; P=3.9 × 10(-9), per-allele OR=1.19) were newly identified as loci associated with EGFR mutation-positive lung adenocarcinoma. This study indicates that multiple genetic factors underlie the risk of lung adenocarcinomas with EGFR mutations.

Meta-analysis of genome-wide association studies identifies multiple lung cancer susceptibility loci in never-smoking Asian women.

Genome-wide association studies (GWAS) of lung cancer in Asian never-smoking women have previously identified six susceptibility loci associated with lung cancer risk. To further discover new susceptibility loci, we imputed data from four GWAS of Asian non-smoking female lung cancer (6877 cases and 6277 controls) using the 1000 Genomes Project (Phase 1 Release 3) data as the reference and genotyped additional samples (5878 cases and 7046 controls) for possible replication. In our meta-analysis, three new loci achieved genome-wide significance, marked by single nucleotide polymorphism (SNP) rs7741164 at 6p21.1 (per-allele odds ratio (OR) = 1.17; P = 5.8 × 10(-13)), rs72658409 at 9p21.3 (per-allele OR = 0.77; P = 1.41 × 10(-10)) and rs11610143 at 12q13.13 (per-allele OR = 0.89; P = 4.96 × 10(-9)). These findings identified new genetic susceptibility alleles for lung cancer in never-smoking women in Asia and merit follow-up to understand their biological underpinnings.

Performance comparison of four commercial human whole-exome capture platforms.

Whole exome sequencing (WXS) is widely used to identify causative genetic mutations of diseases. However, not only have several commercial human exome capture platforms been developed, but substantial updates have been released in the past few years. We report a performance comparison for the latest release of four commercial platforms, Roche/NimbleGen's SeqCap EZ Human Exome Library v3.0, Illumina's Nextera Rapid Capture Exome (v1.2), Agilent's SureSelect XT Human All Exon v5 and Agilent's SureSelect QXT, using the same DNA samples. Agilent XT showed the highest target enrichment efficiency and the best SNV and short indel detection sensitivity in coding regions with the least amount of sequencing. Agilent QXT had slightly inferior target enrichment than Agilent XT. Illumina, with additional sequencing, detected SNVs and short indels at the same quality as Agilent XT, and showed the best performance in coverage of medically interesting mutations. NimbleGen detected more SNVs and indels in untranslated regions than the others. We also found that the platforms, which enzymatically fragment the genomic DNA (gDNA), detected more homozygous SNVs than those using sonicated gDNA. We believe that our analysis will help investigators when selecting a suitable exome capture platform for their particular research.

Construction of a population-specific HLA imputation reference panel and its application to Graves' disease risk in Japanese.

To fine map association signals of human leukocyte antigen (HLA) variants in the major histocompatibility complex (MHC) region, we constructed a Japanese population-specific reference panel (n = 908). We conducted trans-ancestry comparisons of linkage disequilibrium (LD) and haplotype structure for HLA variants using an entropy-based LD measurement, ɛ, and a visualization tool to capture high-dimensional variables. Our Japanese reference panel exhibited stronger LD between HLA genes than European or other East Asian populations, characterized by one population-specific common long-range HLA haplotype. We applied HLA imputation to genome-wide association study (GWAS) data for Graves' disease in Japanese (n = 9,003) and found that amino acid polymorphisms of multiple class I and class II HLA genes independently contribute to disease risk (HLA-DPB1, HLA-A, HLA-B and HLA-DRB1; P < 2.3 × 10(-6)), with the strongest impact at HLA-DPB1 (P = 1.6 × 10(-42)). Our study illustrates the value of population-specific HLA reference panels.

A genome-wide association study of late-onset Alzheimer's disease in a Japanese population.

OBJECTIVE: Although a number of genome-wide association studies (GWASs) of late-onset Alzheimer's disease (LOAD) have been carried out, there have been little GWAS data on East Asian populations. DESIGN: To discover the novel susceptibility loci of LOAD, we carried out a GWAS using 816 LOAD cases and 7992 controls with a replication analysis using an independent panel of 1011 LOAD cases and 7212 controls in a Japanese population. In addition, we carried out a stratified analysis by APOE-ε4 status to eliminate the established effect of APOE region. RESULTS: Our data indicated that 18p11.32 (rs1992269, P = 9.77 × 10(-7)), CNTNAP2 (rs802571, P = 1.26 × 10(-6)), and 12q24.23 (rs11613092, P = 6.85 × 10(-6)) were suggestive loci for susceptibility to LOAD. CONCLUSION: We identified three suggestive loci for susceptibility to LOAD in a Japanese population. Among these, rs802571, located at intron 1 of CNTNAP2, was considered to be a plausible candidate locus from a functional perspective.

A genome-wide association study identifies PLCL2 and AP3D1-DOT1L-SF3A2 as new susceptibility loci for myocardial infarction in Japanese.

Despite considerable progress in preventive and therapeutic strategies, myocardial infarction (MI) is one of the leading causes of death throughout the world. A total of 55 susceptibility genes have been identified mostly in European genome-wide association studies (GWAS). Nevertheless, large-scale GWAS from other population could possibly find additional susceptibility loci. To identify as many MI susceptibility loci as possible, we performed a large-scale genomic analysis in Japanese population. To identify MI susceptibility loci in Japanese, we conducted a GWAS using 1666 cases and 3198 controls using the Illumina Human610-Quad BeadChip and HumanHap550v3 Genotyping BeadChip. We performed replication studies using a total of 11,412 cases and 28,397 controls in the Japanese population. Our study identified two novel susceptibility loci for MI: PLCL2 on chromosome 3p24.3 (rs4618210:A>G, P = 2.60 × 10(-9), odds ratio (OR) = 0.91) and AP3D1-DOT1L-SF3A2 on chromosome 19p13.3 (rs3803915:A>C, P = 3.84 × 10(-9), OR = 0.89). Besides, a total of 14 previously reported MI susceptibility loci were replicated in our study. In particular, we validated a strong association on chromosome 12q24 (rs3782886:A>G: P = 1.14 × 10(-14), OR = 1.46). Following pathway analysis using 265 genes related to MI or coronary artery disease, we found that these loci might be involved in the pathogenesis of MI via the promotion of atherosclerosis. In the present large-scale genomic analysis, we identified PLCL2 and AP3D1-DOT1L-SF3A2 as new susceptibility loci for MI in the Japanese population. Our findings will add novel findings for MI susceptibility loci.

Genome-wide association analysis in East Asians identifies breast cancer susceptibility loci at 1q32.1, 5q14.3 and 15q26.1.

In a three-stage genome-wide association study among East Asian women including 22,780 cases and 24,181 controls, we identified 3 genetic loci newly associated with breast cancer risk, including rs4951011 at 1q32.1 (in intron 2 of the ZC3H11A gene; P=8.82×10(-9)), rs10474352 at 5q14.3 (near the ARRDC3 gene; P=1.67×10(-9)) and rs2290203 at 15q26.1 (in intron 14 of the PRC1 gene; P=4.25×10(-8)). We replicated these associations in 16,003 cases and 41,335 controls of European ancestry (P=0.030, 0.004 and 0.010, respectively). Data from the ENCODE Project suggest that variants rs4951011 and rs10474352 might be located in an enhancer region and transcription factor binding sites, respectively. This study provides additional insights into the genetics and biology of breast cancer.

Genome-wide association study of breast cancer in the Japanese population.

Breast cancer is the most common malignancy among women in worldwide including Japan. Several studies have identified common genetic variants to be associated with the risk of breast cancer. Due to the complex linkage disequilibrium structure and various environmental exposures in different populations, it is essential to identify variants associated with breast cancer in each population, which subsequently facilitate the better understanding of mammary carcinogenesis. In this study, we conducted a genome-wide association study (GWAS) as well as whole-genome imputation with 2,642 cases and 2,099 unaffected female controls. We further examined 13 suggestive loci (P<1.0 × 10(-5)) using an independent sample set of 2,885 cases and 3,395 controls and successfully validated two previously-reported loci, rs2981578 (combined P-value of 1.31 × 10(-12), OR = 1.23; 95% CI = 1.16-.30) on chromosome 10q26 (FGFR2), rs3803662 (combined P-value of 2.79 × 10(-11), OR = 1.21; 95% CI = 1.15-.28) and rs12922061 (combined P-value of 3.97 × 10(-10), OR = 1.23; 95% CI = 1.15-.31) on chromosome 16q12 (TOX3-LOC643714). Weighted genetic risk score on the basis of three significantly associated variants and two previously reported breast cancer associated loci in East Asian population revealed that individuals who carry the most risk alleles in category 5 have 2.2 times higher risk of developing breast cancer in the Japanese population than those who carry the least risk alleles in reference category 1. Although we could not identify additional loci associated with breast cancer, our study utilized one of the largest sample sizes reported to date, and provided genetic status that represent the Japanese population. Further local and international collaborative study is essential to identify additional genetic variants that could lead to a better, accurate prediction for breast cancer.