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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Inositol polyphosphate 4-phosphatase B as a regulator of bone mass in mice and humans.

Osteoporosis is a multifactorial genetic disease characterized by reduction of bone mass due to dysregulation of osteoclast differentiation or maturation. Herein, we identified a regulator of osteoclastogenesis, the murine homolog of inositol polyphosphate 4-phosphatase type IIα (Inpp4bα). Expression of Inpp4bα is detected from early osteoclast differentiation to activation stage. Targeted expression of native Inpp4bα ex vivo repressed whereas phosphatase-inactive Inpp4bα stimulated osteoclast differentiation. Inpp4bα acts on intracellular calcium level that modulates NFATc1 nuclear translocation and activation. In vivo mice deficient in Inpp4b displayed increased osteoclast differentiation rate and potential resulting in decreased bone mass and osteoporosis. Importantly, INPP4B in human was identified as a susceptibility locus for osteoporosis. This study defined Inpp4b as a major modulator of the osteoclast differentiation and as a gene linked to variability of bone mineral density in mice and humans.

High-density polymorphisms analysis of 23 candidate genes for association with bone mineral density.

Osteoporosis is a bone disease characterized by low bone mineral density (BMD), a highly heritable and polygenic trait. Women are more prone than men to develop osteoporosis due to a lower peak bone mass and accelerated bone loss at menopause. Peak bone mass has been convincingly shown to be due to genetic factors with heritability up to 80%. Menopausal bone loss has been shown to have around 38% to 49% heritability depending on the site studied. To have more statistical power to detect small genetic effects we focused on premenopausal women. We studied 23 candidate genes, some involved in calcium and vitamin-D regulation and others because estrogens strongly induced their gene expression in mice where it was correlated with humerus trabecular bone density. High-density polymorphisms were selected to cover the entire gene variability and 231 polymorphisms were genotyped in a first sample of 709 premenopausal women. Positive associations were retested in a second, independent, sample of 673 premenopausal women. Ten polymorphisms remained associated with BMD in the combined samples and one was further associated in a large sample of postmenopausal women (1401 women). This associated polymorphism was located in the gene CSF3R (granulocyte colony stimulating factor receptor) that had never been associated with BMD before. The results reported in this study suggest a role for CSF3R in the determination of bone density in women.

Association with replication between estrogen-related receptor gamma (ESRRgamma) polymorphisms and bone phenotypes in women of European ancestry.

Osteoporosis is a bone disease characterized by low bone mineral density (BMD), a highly heritable polygenic trait. Women are more prone than men to develop osteoporosis owing to a lower peak bone mass and accelerated bone loss at menopause. Lack of estrogen thus is a major risk factor for osteoporosis. In addition to having strong similarity to the estrogen receptor 1 (ESR1), the orphan nuclear estrogen-related receptor gamma (ESRRgamma) is widely expressed and shows overlap with ESR1 expression in tissues where estrogen has important physiologic functions. For these reasons, we have undertaken a study of ESRRgamma sequence variants in association with bone measurements [heel quantitative ultrasound (QUS) by measurements of broadband ultrasound attenuation (BUA), speed of sound (SOS), and stiffness index (SI) and dual-energy X-ray absorptiometry (DXA) at the femoral neck (FN) and lumbar spine (LS)]. A silent variant was found to be associated with multiple bone measurements (LS, BUA, SOS, and SI), the p values ranging from .006 to .04 in a sample of 5144 Quebec women. The region of this variant was analyzed using the HapMap database and the Gabriel method to define a block of 20 kb. Using the Tagger method, eight TagSNPs were identified and genotyped in a sample of 1335 women. Four of these SNPs capture the five major block haplotypes. One SNP (rs2818964) and one haplotype were significantly associated with multiple bone measures. All SNPs involved in the associations were analyzed in two other sample sets with significant results in the same direction. These results suggest involvement of ESRRgamma in the determination of bone density in women.