Effect of preexamination conditions in a centralized-testing model of non-invasive prenatal screening.

OBJECTIVES: Non-invasive prenatal testing requires the presence of fetal DNA in maternal plasma. Understanding how preexamination conditions affect the integrity of cell-free DNA (cfDNA) and fetal fraction (FF) are a prerequisite for test implementation. Therefore, we examined the adjusted effect that EDTA and Streck tubes have on the cfDNA quantity and FF. METHODS: A total of 3,568 maternal blood samples across Canada were collected in either EDTA, or Streck tubes, and processing metrics, maternal body mass index (BMI), gestational age and fetal karyotype and sex were recorded. Plasma samples were sequenced using two different sequencing platforms in separate laboratories. Sequencing data were processed with SeqFF to estimate FF. Linear regression and multivariate imputation by chained equations were used to estimate the adjusted effect of tube type on cfDNA and FF. RESULTS: We found a positive association between cfDNA quantity and blood shipment time in EDTA tubes, which is significantly reduced with the use of Streck tubes. Furthermore, we show the storage of plasma at -80 °C is associated with a 4.4% annual relative decrease in cfDNA levels. FF was not associated with collection tube type when controlling for confounding variables. However, FF was positively associated with gestational age and trisomy 21, while negatively associated with BMI, male fetus, trisomy 18, Turners syndrome and triploidy. CONCLUSIONS: Preexamination, maternal and fetal variables are associated with cfDNA quantity and FF. The consideration of these variables in future studies may help to reduce the number of pregnant women with inconclusive tests as a result of low FF.

Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium.

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n = 14 260), velocity of sound (VOS; n = 15 514) and BMD (n = 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n = 11 452) and new genotyping in 15 cohorts (de novo n = 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 × 10(-8)) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 × 10(-14)). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 × 10(-6) also had the expected direction of association with any fracture (P < 0.05), including three SNPs with P < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, this GWA study reveals the effect of several genes common to central DXA-derived BMD and heel ultrasound/DXA measures and points to a new genetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Improving Fetal Fraction of Noninvasive Prenatal Screening Samples Collected in EDTA-Gel Tubes Using Gel Size Selection: A Head-To-Head Comparison of Methods.

The aim of this study was to compare the use of EDTA-gel blood collection tubes with and without size selection to cell-stabilizing collection tubes for remote blood sampling for noninvasive prenatal screening (NIPS). Sixty-one pregnant women at 10 to 14 weeks' gestation undergoing NIPS were recruited. Participants were phlebotomized with Streck and EDTA-gel tubes. EDTA-gel tubes were centrifuged before shipping. Libraries prepared from cell-free DNA (cfDNA) extracted from both types of tubes were sequenced on Illumina NextSeq 500, and fetal fraction was estimated using SeqFF. EDTA-gel tube libraries were size selected on agarose gel to eliminate cfDNA fragments >160 bp and resequenced. The main outcome measure was fetal fraction expressed as percentage of total cfDNA sequenced, calculated from sequence read counts (SeqFF). Streck tube samples showed an average 1% higher fetal fraction than centrifuged EDTA-gel tubes without size selection. This difference increased with temperature. When EDTA-gel samples' libraries were size selected, the mean fetal fraction increased from 7% to 13%, with no sample having fetal fraction <4%. Using EDTA-gel tubes reduces NIPS sampling cost and tube processing time in the laboratory. Also, using EDTA-gel tubes does not lead to cfDNA degradation. Size selection increases fetal fraction, reduces the number of test failures, increases NIPS clinical performance, and may be helpful in situations asking for a higher fetal fraction, such as twin pregnancies or screening for sub-chromosomal imbalances.

Validation of a New Protocol to Collect and Isolate Plasma from Pregnant Women for Noninvasive Prenatal Testing (NIPT).

BACKGROUND: Most laboratories use specialized tubes (e.g., Streck) to recover circulating cell-free DNA (ccfDNA) for noninvasive prenatal testing (NIPT). We validated a low cost, simple procedure for collecting NIPT samples in remote laboratories that avoids highspeed centrifugation. EDTA gel blood sampling tube allows simple separation of plasma from blood cells. Decanted plasma is filtered to remove cell debris. The procedure can be performed within a few minutes after the blood centrifugation step, and ccfDNA-grade plasma can be frozen for transportation. METHODS: We recruited 51 pregnant women and collected blood in one EDTA-gel Greiner tube and two Streck tubes. All tubes were centrifuged at 1600 g x 10 min within 6 h of sample collection. Plasma from EDTA tubes was poured into a syringe cylinder and filtered through a 0.45 µm Millipore filter. Plasma from Streck tubes was recovered with a pipette and one was filtered as above while the second was centrifuged at 16 000 g. The ccfDNA was isolated and NGS sequencing libraries were prepared and sequenced on an Illumina system. Fetal fractions were estimated using SeqFF. This study had a power of 79% to detect a decrease of 1% in fetal fractions with the new method. RESULTS: We did not observe any significant difference between the three procedures for the fetal fraction nor for the quality or quantity of libraries produced. CONCLUSION: EDTA-gel tubes with filtration provide high quality plasma for ccfDNA analysis and can be sent frozen to the NIPT laboratory. This is economical and it frees the laboratory of time-consuming steps.

Investigating the causal role of MRE11A p.E506* in breast and ovarian cancer.

The nuclease MRE11A is often included in genetic test panels for hereditary breast and ovarian cancer (HBOC) due to its BRCA1-related molecular function in the DNA repair pathway. However, whether MRE11A is a true predisposition gene for HBOC is still questionable. We determined to investigate this notion by dissecting the molecular genetics of the c.1516G > T;p.E506* truncating MRE11A variant, that we pinpointed in two unrelated French-Canadian (FC) HBOC patients. We performed a case-control study for the variant in ~ 2500 breast, ovarian, and endometrial cancer patients from the founder FC population of Quebec. Furthermore, we looked for the presence of second somatic alterations in the MRE11A gene in the tumors of the carriers. In summary, these investigations suggested that the identified variant is not associated with an increased risk of developing breast or ovarian cancer. We finally performed a systematic review for all the previously reported MRE11A variants in breast and ovarian cancer. We found that MRE11A germline variants annotated as pathogenic on ClinVar often lacked evidence for such classification, hence misleading the clinical management for affected patients. In summary, our report suggests the lack of clinical utility of MRE11A testing in HBOC, at least in the White/Caucasian populations.

Precision of Fetal DNA Fraction Estimation by Quantitative Polymerase Chain Reaction Quantification of a Differently Methylated Target in Noninvasive Prenatal Testing.

BACKGROUND: The performance of noninvasive prenatal testing (NIPT) assays is critically determined by the proportion of fetal DNA or fetal fraction (FF). Fetomaternal differential methylation of certain genomic regions has been proposed as a universal marker of fetal origin, and previous reports have suggested the use of methylation-sensitive restriction enzyme (MSRE) assays to estimate FF. METHODS: We analyzed the performance of FF estimation using an MSRE assay with duplex quantitative polymerase chain reaction (qPCR). Mixtures of genomic DNA from placental cells and from adult women were digested with 2 MSRE and FF estimates obtained, for a total of 221 pairwise treatment/control comparisons. RESULTS: The coefficient of variance (CV) of the MSRE assays was high, ranging from 24% to 60%. An alternative in silico FF estimation algorithm, SeqFF, displayed slightly lower variability, with a CV of 22%. CONCLUSION: These results cast doubts on the usefulness of the MSRE-based assay of differentially methylated markers for FF estimation. The lack of a universal method capable of precisely estimating FF remains an incompletely solved issue.

Development of Reference Materials for Noninvasive Prenatal Aneuploidy Testing by Massively Parallel Sequencing: A Proof-of-Concept Study.

BACKGROUND: Noninvasive prenatal aneuploidy testing (NIPT) represents the first large-scale clinical application of massively parallel sequencing technology. However, no NIPT reference material (RM) has yet been widely adopted, impeding the development of quality management systems and standardization. Developing an NIPT RM from a biological sample is complicated by the low concentration of cell-free DNA (cfDNA), which implies pooling specimens and frequent resampling. METHODS: We tested the feasibility of using DNA from immortalized cell lines of a woman and her aneuploid offspring to spike an artificial plasma matrix. Enzymatic fragmentation of extracted DNA was optimized to achieve fragment size profiles with a mode of 150 to 200 bp, similar to biological cfDNA. This synthetic material was compared with routine biological samples from pregnant women by a targeted NIPT assay in a multiplex sequencing run on a Proton platform. RESULTS: Sequencing statistics were similar between artificially prepared material and routine biological samples, as well as relative chromosomal representation, and no matrix effects could be detected. Estimate of fetal fraction (FF) was within the range of expected value, and aneuploidy detection statistic (z-score) was also comparable between both types of samples. CONCLUSIONS: Artificial plasma spiked with DNA from cell lines of mother and offspring is a promising strategy for developing NIPT RM. This type of material would offer the advantage of a constant and stable composition, allowing for greater standardization of NIPT assays. Moreover, it preserves the parental relatedness used by targeted assay to estimate FF by identification of paternal alleles in single-nucleotide polymorphisms or other variable regions.

Prospective head-to-head comparison of accuracy of two sequencing platforms for screening for fetal aneuploidy by cell-free DNA: the PEGASUS study.

We compared clinical validity of two non-invasive prenatal screening (NIPS) methods for fetal trisomies 13, 18, 21, and monosomy X. We recruited prospectively 2203 women at high risk of fetal aneuploidy and 1807 at baseline risk. Three-hundred and twenty-nine euploid samples were randomly removed. The remaining 1933 high risk and 1660 baseline-risk plasma aliquots were assigned randomly between four laboratories and tested with two index NIPS tests, blind to maternal variables and pregnancy outcomes. The two index tests used massively parallel shotgun sequencing (semiconductor-based and optical-based). The reference standard for all fetuses was invasive cytogenetic analysis or clinical examination at birth and postnatal follow-up. For each chromosome of interest, chromosomal ratios were calculated (number of reads for chromosome/total number of reads). Euploid samples' mean chromosomal ratio coefficients of variation were 0.48 (T21), 0.34 (T18), and 0.31 (T13). According to the reference standard, there were 155 cases of T21, 49 T18, 8 T13 and 22 45,X. Using a fetal fraction ≥4% to call results and a chromosomal ratio z-score of ≥3 to report a positive result, detection rates (DR), and false positive rates (FPR) were not statistically different between platforms: mean DR 99% (T21), 100%(T18, T13); 79%(45,X); FPR < 0.3% for T21, T18, T13, and <0.6% for 45,X. Both methods' negative predictive values in high-risk pregnancies were >99.8%, except for 45,X(>99.6%). Threshold analysis in high-risk pregnancies with different fetal fractions and z-score cut-offs suggested that a z-score cutoff to 3.5 for positive results improved test accuracy. Both sequencing platforms showed equivalent and excellent clinical validity.

Non-invasive prenatal aneuploidy testing: Critical diagnostic performance parameters predict sample z-score values.

OBJECTIVES: Non-invasive prenatal aneuploidy testing (NIPT) by next-generation sequencing of circulating cell-free DNA in maternal plasma relies on chromosomal ratio (chrratio) measurements to detect aneuploid values that depart from euploid ratios. Diagnostic performances are known to depend on the fraction of fetal DNA (FF) present in maternal plasma, although how this translates into specific quantitative changes in specificity/positive predictive values and which other variables might also be important is not well understood. DESIGN & METHODS: To explore this issue, theoretical relationships between FF and various measures of diagnostic performances were assessed for a range of parameter values. Empirical data from three NIPT assays were then used to validate theoretical calculations. RESULTS: For a given positivity threshold, dramatic changes in specificity and positive predictive values (PPV) as a function of both FF and the coefficient of variation (CV) of the chrratio measurement were observed. Theoretically predicted and observed chrratio z-scores agreed closely, confirming the determinant impact of small changes in both FF and chrratio CV. CONCLUSIONS: Evaluation of NIPT assay performances therefore requires knowledge of the FF distribution in the population in which the test is intended to be used and, in particular, of the precise value of the assay chrratio CV for each chromosome or genomic region of interest. Laboratories offering NIPT testing should carefully measure these parameters to ensure test reliability and clinical usefulness in interpreting individual patients' results.

Identification of a novel truncating PALB2 mutation and analysis of its contribution to early-onset breast cancer in French-Canadian women.

BACKGROUND: PALB2 has recently been identified as a breast cancer susceptibility gene. PALB2 mutations are rare causes of hereditary breast cancer but may be important in countries such as Finland where a founder mutation is present. We sought to estimate the contribution of PALB2 mutations to the burden of breast cancer in French Canadians from Quebec. METHODS: We screened all coding exons of PALB2 in a sample of 50 French-Canadian women diagnosed with either early-onset breast cancer or familial breast cancer at a single Montreal hospital. The genetic variants identified in this sample were then studied in 356 additional women with breast cancer diagnosed before age 50 and in 6,448 newborn controls. RESULTS: We identified a single protein-truncating mutation in PALB2 (c.2323 C>T, resulting in Q775X) in 1 of the 50 high-risk women. This variant was present in 2 of 356 breast cancer cases and in none of 6,440 newborn French-Canadian controls (P = 0.003). We also identified two novel new non-synonymous single nucleotide polymorphisms in exon 4 of PALB2 (c.5038 A>G [I76V] and c.5156 G>T [G115V]). G115V was found in 1 of 356 cases and in 15 of 6,442 controls (P = 0.6). The I76V variant was not identified in either the extended case series or the controls. CONCLUSION: We have identified a novel truncating mutation in PALB2. The mutation was found in approximately 0.5% of unselected French-Canadian women with early-onset breast cancer and appears to have a single origin. Although mutations are infrequent, PALB2 can be added to the list of breast cancer susceptibility genes for which founder mutations have been identified in the French-Canadian population.

LRP5 coding polymorphisms influence the variation of peak bone mass in a normal population of French-Canadian women.

INTRODUCTION: Bone mineral density has a strong genetic component but it is also influenced by environmental factors making it a complex trait to study. LRP5 gene was previously shown to be involved in rare diseases affecting bone mass. Mutations associated with gain-of-function were described as well as loss-of-function mutations. Following this discovery, many frequent LRP5 polymorphisms were tested against the variation of BMD in the normal population. MATERIALS AND METHODS: Heel bone parameters (SOS, BUA) were measured by right calcaneal QUS in 5021 healthy French-Canadian women and for 2104 women, BMD evaluated by DXA at two sites was available (femoral neck (FN) and lumbar spine (LS)). Among women with QUS measures and those with DXA measures, 26.5% and 32.8% respectively were premenopausal, 9.2% and 10.7% were perimenopausal and 64.2% and 56.5% were postmenopausal. About a third of the peri- and postmenopausal women never received hormone therapy. Two single nucleotide coding polymorphisms (Val667Met and Ala1330Val) in LRP5 gene were genotyped by allele-specific PCR. All bone measures were tested individually for associations with each polymorphism by analysis of covariance with adjustment for non genetic risk factors. Furthermore, haplotype analysis was performed to take into account the strong linkage disequilibrium between the two polymorphisms. RESULTS AND CONCLUSION: The two LRP5 polymorphisms were found to be associated with all five bone measures (L2L4 and femoral neck DXA as well as heel SOS, BUA and stiffness index) in the whole sample. Premenopausal women drove the association as expected from the proposed role of LRP5 in peak bone mass. Our results suggest that the Val667Met polymorphism is the causative variant but this remains to be functionally proven.

Functionally Null RAD51D Missense Mutation Associates Strongly with Ovarian Carcinoma.

RAD51D is a key player in DNA repair by homologous recombination (HR), and RAD51D truncating variant carriers have an increased risk for ovarian cancer. However, the contribution of nontruncating RAD51D variants to cancer predisposition remains uncertain. Using deep sequencing and case-control genotyping studies, we show that in French Canadians, the missense RAD51D variant c.620C>T;p.S207L is highly prevalent and is associated with a significantly increased risk for ovarian high-grade serous carcinoma (HGSC; 3.8% cases vs. 0.2% controls). The frequency of the p.S207L variant did not significantly differ from that of controls in breast, endometrial, pancreas, or colorectal adenocarcinomas. Functionally, we show that this mutation impairs HR by disrupting the RAD51D-XRCC2 interaction and confers PARP inhibitor sensitivity. These results highlight the importance of a functional RAD51D-XRCC2 interaction to promote HR and prevent the development of HGSC. This study identifies c.620C>T;p.S207L as the first bona fide pathogenic RAD51D missense cancer susceptibility allele and supports the use of targeted PARP-inhibitor therapies in ovarian cancer patients carrying deleterious missense RAD51D variants. Cancer Res; 77(16); 4517-29. ©2017 AACR.

A frequent regulatory variant of the estrogen-related receptor alpha gene associated with BMD in French-Canadian premenopausal women.

UNLABELLED: Genes are important BMD determinants. We studied the association of an ESRRA gene functional variant with BMD in 1335 premenopausal women. The ESRRA genotype was an independent predictor of L2-L4 BMD, with an effect similar to smoking and equivalent to a 10-kg difference in weight. INTRODUCTION: Several genetic polymorphisms have been associated with osteoporosis or osteoporosis fractures, but no functional effect has been shown for most of these gene variants. Because functional studies have implicated estrogen-related receptor alpha (ESRRA) in bone metabolism, we evaluated whether a recently described regulatory variant of the ESRRA gene is associated with lumbar and hip BMD as measured by DXA and with heel bone parameters as measured by quantitative ultrasound (QUS). MATERIALS AND METHODS: Heel bone parameters were measured by right calcaneal QUS in 1335 healthy French-Canadian premenopausal women, and one-half of these women also had their BMD evaluated at two sites: femoral neck and lumbar spine (L2-L4) by DXA. All bone measures were tested separately for association with the ESRRA genotype by analysis of covariance. The significance of the ESRRA contribution to the model was also assessed by two different permutation tests. RESULTS: A statistically significant association between ESRRA genotype and lumbar spine BMD was observed: women carrying the long ESRRA genotype had a 3.9% (0.045 g/cm2) higher lumbar spine BMD than those carrying the short ESRRA genotype (p = 0.004), independently of other risk factors measured. This effect of ESRRA genotype is similar to the effect of smoking and equivalent to a 10-kg difference in weight. This association was confirmed by permutation tests (p = 0.004). The same trend was observed for femoral neck BMD (2.6%, p = 0.07). However, no association was observed between ESRRA and QUS heel bone measures. CONCLUSION: These results support the genetic influence of this ESRRA regulatory variant on BMD.

Risk of Misdiagnosis Due to Allele Dropout and False-Positive PCR Artifacts in Molecular Diagnostics: Analysis of 30,769 Genotypes.

Quality control is a complex issue for clinical molecular diagnostic applications. In the case of genotyping assays, artifacts such as allele dropout represent a risk of misdiagnosis for amplification-based methods. However, its frequency of occurrence in PCR-based diagnostic assays remains unknown. To maximize the likelihood of detecting allele dropout, our clinical genotyping PCR-based assays are designed with two independent assays for each allele (nonoverlapping primers on each DNA strand). To estimate the incidence of allelic dropout, we took advantage of the capacity of our clinical assays to detect such events. We retrospectively studied their occurrence in the initial PCR assay for 30,769 patient reports for mutations involved in four diseases produced over 8 years. Ninety-three allele dropout events were detected and all were solved before reporting. In addition, 42 cases of artifacts caused by amplification of an allele ultimately confirmed to not be part of the genotype (drop-in events) were detected and solved. These artifacts affected 1:227 genotypes, 94% of which were due to nonreproducible PCR failures rather than sequence variants interfering with the assay, suggesting that careful primer design cannot prevent most of these errors. This provides a quantitative estimate for clinical laboratories to take this phenomenon into account in quality management and to favor assay designs that can detect (and minimize) occurrence of these artifacts in routine clinical use.

Germline RECQL mutations are associated with breast cancer susceptibility.

Several moderate- and high-risk breast cancer susceptibility genes have been discovered, but more are likely to exist. To discover new breast cancer susceptibility genes, we used 2 populations (from Poland and Quebec, Canada) and applied whole-exome sequencing in a discovery phase (n = 195), followed by validation. We identified rare recurrent RECQL mutations in each population. In Quebec, 7 of 1,013 higher-risk breast cancer cases and 1 of 7,136 newborns carried the c.634C>T (p.Arg215*) variant (P = 0.00004). In Poland, 30 of 13,136 unselected breast cancer cases and 2 of 4,702 controls carried the c.1667_1667+3delAGTA (p.K555delinsMYKLIHYSFR) variant (P = 0.008). RECQL is implicated in resolving stalled DNA replication forks to prevent double-stranded DNA (dsDNA) breaks. This function is related to that of other known breast cancer susceptibility genes, many of which are involved in repairing dsDNA breaks. We conclude that RECQL is a breast cancer susceptibility gene.

Atypical presentation of an anomalous origin of the right coronary artery with severe compression between the great vessels.

Anomalous aortic origin of the right coronary artery is a rare coronary anomaly which, in a minority of cases, can cause clinical manifestations such as ischemic chest pain or arrhythmic syncope. We describe a case of anomalous aortic origin of the right coronary artery characterized by signs of left heart failure associated with ventricular tachycardia.

Prevalence of selected genomic deletions and duplications in a French-Canadian population-based sample of newborns.

Chromosomal microarray analysis has identified many novel microdeletions or microduplications that produce neurodevelopmental disorders with a recognizable clinical phenotype and that are not observed in normal individuals. However, imbalance of other genomic regions is associated with a variable phenotype with intellectual disability (ID) or autism in some individuals but are also observed in completely normal individuals. Several large studies have reported the prevalence of copy number (CN) variants in people with particular features (e.g., ID, autism, schizophrenia, or epilepsy); few studies have investigated the prevalence of genomic CN changes in the general population. We used a high-throughput method to screen 6813 consecutive cord blood samples from a predominantly French-Canadian population to assess genomic CN in five genomic regions: 1p36, 15q11-q13, 16p11.2, 16p11.2-p12.2, and 22q11.2. We identified one deletion and one duplication within 1p36, two deletions of 15q11-q13, eight deletions of 16p11.2-p12.2, two deletions and five duplications of 16p11.2, and six duplications of 22q11.2. This study provides estimates of the frequency of CN variants in an unselected population. Our findings have important implications for genetic counseling.

Meta-analysis of genome-wide studies identifies WNT16 and ESR1 SNPs associated with bone mineral density in premenopausal women.

Previous genome-wide association studies (GWAS) have identified common variants in genes associated with variation in bone mineral density (BMD), although most have been carried out in combined samples of older women and men. Meta-analyses of these results have identified numerous single-nucleotide polymorphisms (SNPs) of modest effect at genome-wide significance levels in genes involved in both bone formation and resorption, as well as other pathways. We performed a meta-analysis restricted to premenopausal white women from four cohorts (n = 4061 women, aged 20 to 45 years) to identify genes influencing peak bone mass at the lumbar spine and femoral neck. After imputation, age- and weight-adjusted bone-mineral density (BMD) values were tested for association with each SNP. Association of an SNP in the WNT16 gene (rs3801387; p = 1.7 × 10(-9) ) and multiple SNPs in the ESR1/C6orf97 region (rs4870044; p = 1.3 × 10(-8) ) achieved genome-wide significance levels for lumbar spine BMD. These SNPs, along with others demonstrating suggestive evidence of association, were then tested for association in seven replication cohorts that included premenopausal women of European, Hispanic-American, and African-American descent (combined n = 5597 for femoral neck; n = 4744 for lumbar spine). When the data from the discovery and replication cohorts were analyzed jointly, the evidence was more significant (WNT16 joint p = 1.3 × 10(-11) ; ESR1/C6orf97 joint p = 1.4 × 10(-10) ). Multiple independent association signals were observed with spine BMD at the ESR1 region after conditioning on the primary signal. Analyses of femoral neck BMD also supported association with SNPs in WNT16 and ESR1/C6orf97 (p < 1 × 10(-5) ). Our results confirm that several of the genes contributing to BMD variation across a broad age range in both sexes have effects of similar magnitude on BMD of the spine in premenopausal women. These data support the hypothesis that variants in these genes of known skeletal function also affect BMD during the premenopausal period.

Assessment of gene-by-sex interaction effect on bone mineral density.

Sexual dimorphism in various bone phenotypes, including bone mineral density (BMD), is widely observed; however, the extent to which genes explain these sex differences is unclear. To identify variants with different effects by sex, we examined gene-by-sex autosomal interactions genome-wide, and performed expression quantitative trait loci (eQTL) analysis and bioinformatics network analysis. We conducted an autosomal genome-wide meta-analysis of gene-by-sex interaction on lumbar spine (LS) and femoral neck (FN) BMD in 25,353 individuals from 8 cohorts. In a second stage, we followed up the 12 top single-nucleotide polymorphisms (SNPs; p < 1 × 10(-5) ) in an additional set of 24,763 individuals. Gene-by-sex interaction and sex-specific effects were examined in these 12 SNPs. We detected one novel genome-wide significant interaction associated with LS-BMD at the Chr3p26.1-p25.1 locus, near the GRM7 gene (male effect = 0.02 and p = 3.0 × 10(-5) ; female effect = -0.007 and p = 3.3 × 10(-2) ), and 11 suggestive loci associated with either FN- or LS-BMD in discovery cohorts. However, there was no evidence for genome-wide significant (p < 5 × 10(-8) ) gene-by-sex interaction in the joint analysis of discovery and replication cohorts. Despite the large collaborative effort, no genome-wide significant evidence for gene-by-sex interaction was found to influence BMD variation in this screen of autosomal markers. If they exist, gene-by-sex interactions for BMD probably have weak effects, accounting for less than 0.08% of the variation in these traits per implicated SNP. © 2012 American Society for Bone and Mineral Research.

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture.

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10(-8)). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10(-4), Bonferroni corrected), of which six reached P < 5 × 10(-8), including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.