Polygenic risk for prostate cancer: Decreasing relative risk with age but little impact on absolute risk.

Polygenic risk scores (PRSs) for a variety of diseases have recently been shown to have relative risks that depend on age, and genetic relative risks decrease with increasing age. A refined understanding of the age dependency of PRSs for a disease is important for personalized risk predictions and risk stratification. To further evaluate how the PRS relative risk for prostate cancer depends on age, we refined analyses for a validated PRS for prostate cancer by using 64,274 prostate cancer cases and 46,432 controls of diverse ancestry (82.8% European, 9.8% African American, 3.8% Latino, 2.8% Asian, and 0.8% Ghanaian). Our strategy applied a novel weighted proportional hazards model to case-control data to fully utilize age to refine how the relative risk decreased with age. We found significantly greater relative risks for younger men (age 30-55 years) compared with older men (70-88 years) for both relative risk per standard deviation of the PRS and dichotomized according to the upper 90th percentile of the PRS distribution. For the largest European ancestral group that could provide reliable resolution, the log-relative risk decreased approximately linearly from age 50 to age 75. Despite strong evidence of age-dependent genetic relative risk, our results suggest that absolute risk predictions differed little from predictions that assumed a constant relative risk over ages, from short-term to long-term predictions, simplifying implementation of risk discussions into clinical practice.

Penalized mediation models for multivariate data.

Statistical methods to integrate multiple layers of data, from exposures to intermediate traits to outcome variables, are needed to guide interpretation of complex data sets for which variables are likely contributing in a causal pathway from exposure to outcome. Statistical mediation analysis based on structural equation models provide a general modeling framework, yet they can be difficult to apply to high-dimensional data and they are not automated to select the best fitting model. To overcome these limitations, we developed novel algorithms and software to simultaneously evaluate multiple exposure variables, multiple intermediate traits, and multiple outcome variables. Our penalized mediation models are computationally efficient and simulations demonstrate that they produce reliable results for large data sets. Application of our methods to a study of vascular disease demonstrates their utility to identify novel direct effects of single-nucleotide polymorphisms (SNPs) on coronary heart disease and peripheral artery disease, while disentangling the effects of SNPs on the intermediate risk factors including lipids, cigarette smoking, systolic blood pressure, and type 2 diabetes.

Integration of multiomics data shows down regulation of mismatch repair and tubulin pathways in triple-negative chemotherapy-resistant breast tumors.

BACKGROUND: Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype. Patients with TNBC are primarily treated with neoadjuvant chemotherapy (NAC). The response to NAC is prognostic, with reductions in overall survival and disease-free survival rates in those patients who do not achieve a pathological complete response (pCR). Based on this premise, we hypothesized that paired analysis of primary and residual TNBC tumors following NAC could identify unique biomarkers associated with post-NAC recurrence. METHODS AND RESULTS: We investigated 24 samples from 12 non-LAR TNBC patients with paired pre- and post-NAC data, including four patients with recurrence shortly after surgery (< 24 months) and eight who remained recurrence-free (> 48 months). These tumors were collected from a prospective NAC breast cancer study (BEAUTY) conducted at the Mayo Clinic. Differential expression analysis of pre-NAC biopsies showed minimal gene expression differences between early recurrent and nonrecurrent TNBC tumors; however, post-NAC samples demonstrated significant alterations in expression patterns in response to intervention. Topological-level differences associated with early recurrence were implicated in 251 gene sets, and an independent assessment of microarray gene expression data from the 9 paired non-LAR samples available in the NAC I-SPY1 trial confirmed 56 gene sets. Within these 56 gene sets, 113 genes were observed to be differentially expressed in the I-SPY1 and BEAUTY post-NAC studies. An independent (n = 392) breast cancer dataset with relapse-free survival (RFS) data was used to refine our gene list to a 17-gene signature. A threefold cross-validation analysis of the gene signature with the combined BEAUTY and I-SPY1 data yielded an average AUC of 0.88 for six machine-learning models. Due to the limited number of studies with pre- and post-NAC TNBC tumor data, further validation of the signature is needed. CONCLUSION: Analysis of multiomics data from post-NAC TNBC chemoresistant tumors showed down regulation of mismatch repair and tubulin pathways. Additionally, we identified a 17-gene signature in TNBC associated with post-NAC recurrence enriched with down-regulated immune genes.

Nine-gene pharmacogenomics profile service: The Mayo Clinic experience.

PURPOSE: The Pharmacogenomics (PGx) Profile Service was a proof-of-concept project to implement PGx in patient care at Mayo Clinic. METHODS: Eighty-two healthy individuals aged 18 and older underwent genotyping of CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, SLCO1B1, HLA-B*58:01, and VKORC1. A PGx pharmacist was involved in ordering, meeting with patients, interpreting, reviewing, and documenting results. RESULTS: Ninety three percent were CYP1A2 rapid metabolizers, 92% CYP3A4 normal metabolizers, and 88% CYP3A5 poor metabolizers; phenotype frequencies for CYP2C19 and CYP2D6 varied. Seventy-three percent had normal functioning SLCO1B1 transporter, 4% carried the HLA-B*58:01 risk variant, and 35% carried VKORC1 and CYP2C9 variants that increased warfarin sensitivity. CONCLUSION: Pre-emptive PGx testing offered medication improvement opportunity in 56% of participants for commonly used medications. A collaborative approach involving a PGx pharmacist integrated within a clinical practice with regards to utility of PGx results allowed for implementation of the PGx Profile Service. KEY POINTS: The Mayo Clinic PGx (PGx) Profile Service was a proof-of-concept project to utilize PGx testing as another clinical tool to enhance medication selection and decrease serious adverse reactions or medication failures. Over one-half of participants in the pilot using the PGx Profile Service were predicted to benefit from pre-emptive PGx testing to guide pharmacotherapy. PGx pharmacists played a crucial role in the PGx Profile Service by educating participants, identifying medication-gene interactions, and providing evidence-based (CPIC and DPWG) PGx recommendations for past, current, and future medication us.

Penalized models for analysis of multiple mediators.

Mediation analysis attempts to determine whether the relationship between an independent variable (e.g., exposure) and an outcome variable can be explained, at least partially, by an intermediate variable, called a mediator. Most methods for mediation analysis focus on one mediator at a time, although multiple mediators can be jointly analyzed by structural equation models (SEMs) that account for correlations among the mediators. We extend the use of SEMs for the analysis of multiple mediators by creating a sparse group lasso penalized model such that the penalty considers the natural groupings of parameters that determine mediation, as well as encourages sparseness of the model parameters. This provides a way to simultaneously evaluate many mediators and select those that have the most impact, a feature of modern penalized models. Simulations are used to illustrate the benefits and limitations of our approach, and application to a study of DNA methylation and reactive cortisol stress following childhood trauma discovered two novel methylation loci that mediate the association of childhood trauma scores with reactive cortisol stress levels. Our new methods are incorporated into R software called regmed.

Penalized variance components for association of multiple genes with traits.

Variance component models have gained popularity for genetic analyses, driven by their flexibility to simultaneously analyze multiple genetic variants in a gene by kernel statistics, and their ability to account for population stratification via genomic relationship matrices. For exploratory analyses with modest sample sizes and a potentially large number of variance components, it can be challenging to use standard maximum-likelihood or restricted maximum-likelihood methods to estimate variance components, because these iterative methods often fail to converge when likelihood surfaces are fairly flat, and standard-likelihood ratio statistical tests are not adequate. To overcome these limitations, we developed a penalized-likelihood model, whereby the penalty function follows the popular elastic-net approach, applying both L1 and L2 penalties to the variance components. By simulations, we demonstrate the potential gain in power by using both L1 and L2 penalties, and results from our simulations suggest that assigning 80% of the penalty parameter to the L1 penalty and 20% to the L2 penalty provides a reasonable balance between false-positive and false-negative results. Larger sample size improves the properties of our methods, at the cost of longer computation time. Application of our methods to a study of the influence of DNA methylation on levels of cortisol in reaction to stress testing shows how our method can be used to prioritize findings for further functional studies.

Multivariate generalized linear model for genetic pleiotropy.

When a single gene influences more than one trait, known as pleiotropy, it is important to detect pleiotropy to improve the biological understanding of a gene. This can lead to improved screening, diagnosis, and treatment of diseases. Yet, most current multivariate methods to evaluate pleiotropy test the null hypothesis that none of the traits are associated with a variant; departures from the null could be driven by just one associated trait. A formal test of pleiotropy should assume a null hypothesis that one or fewer traits are associated with a genetic variant. We recently developed statistical methods to analyze pleiotropy for quantitative traits having a multivariate normal distribution. We now extend this approach to traits that can be modeled by generalized linear models, such as analysis of binary, ordinal, or quantitative traits, or a mixture of these types of traits. Based on methods from estimating equations, we developed a new test for pleiotropy. We then extended the testing framework to a sequential approach to test the null hypothesis that $k+1$ traits are associated, given that the null of $k$ associated traits was rejected. This provides a testing framework to determine the number of traits associated with a genetic variant, as well as which traits, while accounting for correlations among the traits. By simulations, we illustrate the Type-I error rate and power of our new methods, describe how they are influenced by sample size, the number of traits, and the trait correlations, and apply the new methods to a genome-wide association study of multivariate traits measuring symptoms of major depression. Our new approach provides a quantitative assessment of pleiotropy, enhancing current analytic practice.

Deep sequencing across germline genome-wide association study signals relating to breast cancer events in women receiving aromatase inhibitors for adjuvant therapy of early breast cancer.

OBJECTIVE: To identify additional genetic variants beyond those observed in a previous genome-wide association study (GWAS) in women treated on the MA.27 clinical trial in which women were randomized to 5 years of adjuvant therapy with anastrozole or exemestane. PATIENTS AND METHODS: We performed a matched case-control study in 234 women who had a recurrence of breast cancer (cases) and 649 women who had not (controls). The analysis was restricted to White women with an estrogen receptor-positive breast cancer. Multiplex PCR-based targeted deep sequencing was performed of the MIR2052HG region on chromosome 8 between positions 75.4 and 75.7, a span of 300 kb, in an attempt to identify additional functional single nucleotide polymorphisms (SNPs). RESULTS: A total of 4677 unique variants were identified that had not been identified in the previous GWAS. Clinical Annotation of Variants analysis revealed 10 variants, including eight SNPs and two insertion-deletion mutations with moderate or high impact. However, none of the common and variant regions was significant after adjustment for the most significant SNP (rs13260300) identified in our previous GWAS. We performed haplotype analysis that revealed two regions in which the haplotypes lost significance when adjusted for this prior GWAS SNP and one region with two significant haplotypes (P = 0.046 and 0.031) after adjusting for the GWAS SNP. CONCLUSION: We were unable to identify common or rare variant regions that added value to the findings from our previous GWAS. We did find two haplotypes that were significant after adjusting for our top GWAS SNP but these were considered to be of marginal value.

Identification of missing variants by combining multiple analytic pipelines.

BACKGROUND: After decades of identifying risk factors using array-based genome-wide association studies (GWAS), genetic research of complex diseases has shifted to sequencing-based rare variants discovery. This requires large sample sizes for statistical power and has brought up questions about whether the current variant calling practices are adequate for large cohorts. It is well-known that there are discrepancies between variants called by different pipelines, and that using a single pipeline always misses true variants exclusively identifiable by other pipelines. Nonetheless, it is common practice today to call variants by one pipeline due to computational cost and assume that false negative calls are a small percent of total. RESULTS: We analyzed 10,000 exomes from the Alzheimer's Disease Sequencing Project (ADSP) using multiple analytic pipelines consisting of different read aligners and variant calling strategies. We compared variants identified by using two aligners in 50,100, 200, 500, 1000, and 1952 samples; and compared variants identified by adding single-sample genotyping to the default multi-sample joint genotyping in 50,100, 500, 2000, 5000 and 10,000 samples. We found that using a single pipeline missed increasing numbers of high-quality variants correlated with sample sizes. By combining two read aligners and two variant calling strategies, we rescued 30% of pass-QC variants at sample size of 2000, and 56% at 10,000 samples. The rescued variants had higher proportions of low frequency (minor allele frequency [MAF] 1-5%) and rare (MAF < 1%) variants, which are the very type of variants of interest. In 660 Alzheimer's disease cases with earlier onset ages of ≤65, 4 out of 13 (31%) previously-published rare pathogenic and protective mutations in APP, PSEN1, and PSEN2 genes were undetected by the default one-pipeline approach but recovered by the multi-pipeline approach. CONCLUSIONS: Identification of the complete variant set from sequencing data is the prerequisite of genetic association analyses. The current analytic practice of calling genetic variants from sequencing data using a single bioinformatics pipeline is no longer adequate with the increasingly large projects. The number and percentage of quality variants that passed quality filters but are missed by the one-pipeline approach rapidly increased with sample size.

Establishing and characterizing patient-derived xenografts using pre-chemotherapy percutaneous biopsy and post-chemotherapy surgical samples from a prospective neoadjuvant breast cancer study.

BACKGROUND: Patient-derived xenografts (PDXs) are increasingly used in cancer research as a tool to inform cancer biology and drug response. Most available breast cancer PDXs have been generated in the metastatic setting. However, in the setting of operable breast cancer, PDX models both sensitive and resistant to chemotherapy are needed for drug development and prospective data are lacking regarding the clinical and molecular characteristics associated with PDX take rate in this setting. METHODS: The Breast Cancer Genome Guided Therapy Study (BEAUTY) is a prospective neoadjuvant chemotherapy (NAC) trial of stage I-III breast cancer patients treated with neoadjuvant weekly taxane+/-trastuzumab followed by anthracycline-based chemotherapy. Using percutaneous tumor biopsies (PTB), we established and characterized PDXs from both primary (untreated) and residual (treated) tumors. Tumor take rate was defined as percent of patients with the development of at least one stably transplantable (passed at least for four generations) xenograft that was pathologically confirmed as breast cancer. RESULTS: Baseline PTB samples from 113 women were implanted with an overall take rate of 27.4% (31/113). By clinical subtype, the take rate was 51.3% (20/39) in triple negative (TN) breast cancer, 26.5% (9/34) in HER2+, 5.0% (2/40) in luminal B and 0% (0/3) in luminal A. The take rate for those with pCR did not differ from those with residual disease in TN (p = 0.999) and HER2+ (p = 0.2401) tumors. The xenografts from 28 of these 31 patients were such that at least one of the xenografts generated had the same molecular subtype as the patient. Among the 35 patients with residual tumor after NAC adequate for implantation, the take rate was 17.1%. PDX response to paclitaxel mirrored the patients' clinical response in all eight PDX tested. CONCLUSIONS: The generation of PDX models both sensitive and resistant to standard NAC is feasible and these models exhibit similar biological and drug response characteristics as the patients' primary tumors. Taken together, these models may be useful for biomarker discovery and future drug development.

Exome sequencing reveals frequent deleterious germline variants in cancer susceptibility genes in women with invasive breast cancer undergoing neoadjuvant chemotherapy.

When sequencing blood and tumor samples to identify targetable somatic variants for cancer therapy, clinically relevant germline variants may be uncovered. We evaluated the prevalence of deleterious germline variants in cancer susceptibility genes in women with breast cancer referred for neoadjuvant chemotherapy and returned clinically actionable results to patients. Exome sequencing was performed on blood samples from women with invasive breast cancer referred for neoadjuvant chemotherapy. Germline variants within 142 hereditary cancer susceptibility genes were filtered and reviewed for pathogenicity. Return of results was offered to patients with deleterious variants in actionable genes if they were not aware of their result through clinical testing. 124 patients were enrolled (median age 51) with the following subtypes: triple negative (n = 43, 34.7%), HER2+ (n = 37, 29.8%), luminal B (n = 31, 25%), and luminal A (n = 13, 10.5%). Twenty-eight deleterious variants were identified in 26/124 (21.0%) patients in the following genes: ATM (n = 3), BLM (n = 1), BRCA1 (n = 4), BRCA2 (n = 8), CHEK2 (n = 2), FANCA (n = 1), FANCI (n = 1), FANCL (n = 1), FANCM (n = 1), FH (n = 1), MLH3 (n = 1), MUTYH (n = 2), PALB2 (n = 1), and WRN (n = 1). 121/124 (97.6%) patients consented to return of research results. Thirteen (10.5%) had actionable variants, including four that were returned to patients and led to changes in medical management. Deleterious variants in cancer susceptibility genes are highly prevalent in patients with invasive breast cancer referred for neoadjuvant chemotherapy undergoing exome sequencing. Detection of these variants impacts medical management.

The kinship2 R package for pedigree data.

BACKGROUND: The kinship2 package is restructured from the previous kinship package. Existing features are now enhanced and new features added for handling pedigree objects. METHODS: Pedigree plotting features have been updated to display features on complex pedigrees while adhering to pedigree plotting standards. Kinship matrices can now be calculated for the X chromosome. Other methods have been added to subset and trim pedigrees while maintaining the pedigree structure. CONCLUSION: We make the kinship2 package available for R on the Contributed R Archives Network (CRAN), where data management is built-in and other packages can use the pedigree object.

Multiple genetic variant association testing by collapsing and kernel methods with pedigree or population structured data.

Searching for rare genetic variants associated with complex diseases can be facilitated by enriching for diseased carriers of rare variants by sampling cases from pedigrees enriched for disease, possibly with related or unrelated controls. This strategy, however, complicates analyses because of shared genetic ancestry, as well as linkage disequilibrium among genetic markers. To overcome these problems, we developed broad classes of "burden" statistics and kernel statistics, extending commonly used methods for unrelated case-control data to allow for known pedigree relationships, for autosomes and the X chromosome. Furthermore, by replacing pedigree-based genetic correlation matrices with estimates of genetic relationships based on large-scale genomic data, our methods can be used to account for population-structured data. By simulations, we show that the type I error rates of our developed methods are near the asymptotic nominal levels, allowing rapid computation of P-values. Our simulations also show that a linear weighted kernel statistic is generally more powerful than a weighted "burden" statistic. Because the proposed statistics are rapid to compute, they can be readily used for large-scale screening of the association of genomic sequence data with disease status.

Breast Cancer Polygenic-Risk Score Influence on Risk-Reducing Endocrine Therapy Use: Genetic Risk Estimate (GENRE) Trial 1-Year and 2-Year Follow-Up.

Refinement of breast cancer risk estimates with a polygenic-risk score (PRS) may improve uptake of risk-reducing endocrine therapy (ET). A previous clinical trial assessed the influence of adding a PRS to traditional risk estimates on ET use. We stratified participants according to PRS-refined breast cancer risk and evaluated ET use and ET-related quality of life (QOL) at 1-year (previously reported) and 2-year follow-ups. Of 151 participants, 58 (38.4%) initiated ET, and 22 (14.6%) discontinued ET by 2 years; 42 (27.8%) and 36 (23.8%) participants were using ET at 1- and 2-year follow-ups, respectively. At the 2-year follow-up, 39% of participants with a lifetime breast cancer risk of 40.1% to 100.0%, 18% with a 20.1% to 40.0% risk, and 16% with a 0.0% to 20.0% risk were taking ET (overall P = 0.01). Moreover, 40% of participants whose breast cancer risk increased by 10% or greater with addition of the PRS to a traditional breast cancer-risk model were taking ET versus 0% whose risk decreased by 10% or greater (P = 0.004). QOL was similar for participants taking or not taking ET at 1- and 2-year follow-ups, although most who discontinued ET did so because of adverse effects. However, these QOL results may have been skewed by the long interval between QOL surveys and lack of baseline QOL data. PRS-informed breast cancer prevention counseling has a lasting, but waning, effect over time. Additional follow-up studies are needed to address the effect of PRS on ET adherence, ET-related QOL, supplemental breast cancer screening, and other risk-reducing behaviors. PREVENTION RELEVANCE: Risk-reducing medications for breast cancer are considerably underused. Informing women at risk with precise and individualized risk assessment tools may substantially affect the incidence of breast cancer. In our study, a risk assessment tool (IBIS-polygenic-risk score) yielded promising results, with 39% of women at highest risk starting preventive medication.

Automated mitotic spindle hotspot counts are highly associated with clinical outcomes in systemically untreated early-stage triple-negative breast cancer.

Operable triple-negative breast cancer (TNBC) has a higher risk of recurrence and death compared to other subtypes. Tumor size and nodal status are the primary clinical factors used to guide systemic treatment, while biomarkers of proliferation have not demonstrated value. Recent studies suggest that subsets of TNBC have a favorable prognosis, even without systemic therapy. We evaluated the association of fully automated mitotic spindle hotspot (AMSH) counts with recurrence-free (RFS) and overall survival (OS) in two separate cohorts of patients with early-stage TNBC who did not receive systemic therapy. AMSH counts were obtained from areas with the highest mitotic density in digitized whole slide images processed with a convolutional neural network trained to detect mitoses. In 140 patients from the Mayo Clinic TNBC cohort, AMSH counts were significantly associated with RFS and OS in a multivariable model controlling for nodal status, tumor size, and tumor-infiltrating lymphocytes (TILs) (p < 0.0001). For every 10-point increase in AMSH counts, there was a 16% increase in the risk of an RFS event (HR 1.16, 95% CI 1.08-1.25), and a 7% increase in the risk of death (HR 1.07, 95% CI 1.00-1.14). We corroborated these findings in a separate cohort of systemically untreated TNBC patients from Radboud UMC in the Netherlands. Our findings suggest that AMSH counts offer valuable prognostic information in patients with early-stage TNBC who did not receive systemic therapy, independent of tumor size, nodal status, and TILs. If further validated, AMSH counts could help inform future systemic therapy de-escalation strategies.

Using the gene ontology to scan multilevel gene sets for associations in genome wide association studies.

Gene-set analyses have been widely used in gene expression studies, and some of the developed methods have been extended to genome wide association studies (GWAS). Yet, complications due to linkage disequilibrium (LD) among single nucleotide polymorphisms (SNPs), and variable numbers of SNPs per gene and genes per gene-set, have plagued current approaches, often leading to ad hoc "fixes." To overcome some of the current limitations, we developed a general approach to scan GWAS SNP data for both gene-level and gene-set analyses, building on score statistics for generalized linear models, and taking advantage of the directed acyclic graph structure of the gene ontology when creating gene-sets. However, other types of gene-set structures can be used, such as the popular Kyoto Encyclopedia of Genes and Genomes (KEGG). Our approach combines SNPs into genes, and genes into gene-sets, but assures that positive and negative effects of genes on a trait do not cancel. To control for multiple testing of many gene-sets, we use an efficient computational strategy that accounts for LD and provides accurate step-down adjusted P-values for each gene-set. Application of our methods to two different GWAS provide guidance on the potential strengths and weaknesses of our proposed gene-set analyses.

Detecting genomic clustering of risk variants from sequence data: cases versus controls.

As the ability to measure dense genetic markers approaches the limit of the DNA sequence itself, taking advantage of possible clustering of genetic variants in, and around, a gene would benefit genetic association analyses, and likely provide biological insights. The greatest benefit might be realized when multiple rare variants cluster in a functional region. Several statistical tests have been developed, one of which is based on the popular Kulldorff scan statistic for spatial clustering of disease. We extended another popular spatial clustering method--Tango's statistic--to genomic sequence data. An advantage of Tango's method is that it is rapid to compute, and when single test statistic is computed, its distribution is well approximated by a scaled χ(2) distribution, making computation of p values very rapid. We compared the Type-I error rates and power of several clustering statistics, as well as the omnibus sequence kernel association test. Although our version of Tango's statistic, which we call "Kernel Distance" statistic, took approximately half the time to compute than the Kulldorff scan statistic, it had slightly less power than the scan statistic. Our results showed that the Ionita-Laza version of Kulldorff's scan statistic had the greatest power over a range of clustering scenarios.

Regression modeling of allele frequencies and testing Hardy Weinberg Equilibrium.

BACKGROUND/AIMS: Tests for whether observed genotype proportions fit Hardy Weinberg Equilibrium (HWE) are widely used in population genetics analyses, as well as to evaluate quality of genotype data. To date, all methods testing for HWE require subjects to be classified into discrete categories, yet it is becoming clear that the distribution of allele frequencies tends to be smooth over geographic regions. METHODS: To evaluate the HWE assumption, we develop new approaches to model allele frequencies as functions of covariates, and use these models to test whether there is residual correlation between the two alleles of subjects; lack of residual correlation supports the null hypothesis of HWE, but conditional on how the covariates influence the allele frequencies. RESULTS: By simulations, we illustrate that a simple statistical test of residual correlation of alleles adequately controls the type I error rate, while maintaining power that is comparable to standard tests for HWE. CONCLUSION: Our approach can be implemented in standard software, enabling more flexible and powerful ways to evaluate the association of covariates with allele frequencies and whether these associations 'explain' departures from HWE when the covariates are ignored, opening new strategies to evaluate the quality of genotype data generated by next-generation sequencing assays.

Prevalence of heavy chain MGUS by race and family history risk groups using a high-sensitivity screening method.

Mass-spectrometry (MS) assays detect lower levels of monoclonal proteins and result in earlier detection of monoclonal gammopathy of undetermined significance (MGUS). We examined heavy chain MGUS prevalence using matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) MS among 3 risk groups, ages 50 or older: 327 African Americans (AA) and 1223 European Americans (EA) from a clinical biobank and 1093 unaffected first-degree relatives (FDR) of patients with hematologic disorders. Age- and sex-adjusted prevalence rates were directly standardized to 2010 United States population. Prevalence ratios were estimated for comparisons of AA and FDR to the EA group using the Poisson distribution. Results were also compared with population-based prevalence using conventional gel-based methods. Risk groups had similar sex and age distributions. MALDI-TOF MGUS prevalence was higher in the AA (16.5% [95% confidence interval (CI), 12.2%, 20.8%]) and FDR (18.3% [95% CI, 16.6%, 21.6%]) than in EA (10.8% [95% CI, 8.8%, 12.7%]), translating to prevalence ratios of 1.73 (95% CI, 1.31, 2.29) and 1.90 (95% CI, 1.55, 2.34), respectively. MALDI-TOF EA prevalence was over threefold higher than conventional estimates but showed similar age trends. Thus, the MALDI-TOF assay found greater numbers with MGUS but similar relative differences by race, family history, and age as prior studies.