Targeted sequencing for hereditary breast and ovarian cancer in BRCA1/2-negative families reveals complex genetic architecture and phenocopies.

Approximately 20% of breast cancer cases are attributed to increased family risk, yet variation in BRCA1/2 can only explain 20%-25% of cases. Historically, only single gene or single variant testing were common in at-risk family members, and further sequencing studies were rarely offered after negative results. In this study, we applied an efficient and inexpensive targeted sequencing approach to provide molecular diagnoses in 245 human samples representing 134 BRCA mutation-negative (BRCAX) hereditary breast and ovarian cancer (HBOC) families recruited from 1973 to 2019 by Dr. Henry Lynch. Sequencing identified 391 variants, which were functionally annotated and ranked based on their predicted clinical impact. Known pathogenic CHEK2 breast cancer variants were identified in five BRCAX families in this study. While BRCAX was an inclusion criterion for this study, we still identified a pathogenic BRCA2 variant (p.Met192ValfsTer13) in one family. A portion of BRCAX families could be explained by other hereditary cancer syndromes that increase HBOC risk: Li-Fraumeni syndrome (gene: TP53) and Lynch syndrome (gene: MSH6). Interestingly, many families carried additional variants of undetermined significance (VOUSs) that may further modify phenotypes of syndromic family members. Ten families carried more than one potential VOUS, suggesting the presence of complex multi-variant families. Overall, nine BRCAX HBOC families in our study may be explained by known likely pathogenic/pathogenic variants, and six families carried potential VOUSs, which require further functional testing. To address this, we developed a functional assay where we successfully re-classified one family's PMS2 VOUS as benign.

In vitro data suggest a role for PMS2 Kozak sequence mutations in Lynch syndrome risk.

Lynch syndrome (LS) is the most common hereditary cancer syndrome. Heterozygous loss-of-function variants in PMS2 are linked to LS. While these variants are not directly cancer causing, reduced PMS2 function results in the accumulation of somatic variants and increased cancer risk over time due to DNA mismatch repair dysfunction. It is reasonable that other types of genetic variation that impact the expression of PMS2 may also contribute to cancer risk. The Kozak sequence is a highly conserved translation initiation motif among higher eukaryotes and is defined as the nine base pairs upstream of the translation start codon through the first four bases of the translated sequence (5'-[GTT]GCATCCATGG-3'; human PMS2: NM_000535.7). While Kozak sequence variants in PMS2 have been reported in ClinVar in patients with suspected hereditary cancer, all variants upstream of the translation start site are currently classified as variants of undetermined significance (VUSs). We hypothesized that variants significantly disrupting the Kozak sequence of PMS2 would decrease PMS2 protein expression, contributing to increased cancer risk over time. Using a dual-luciferase reporter plasmid and site-directed mutagenesis, we generated the wild-type human PMS2 and the ClinVar VUSs within the PMS2 Kozak sequence. Besides the c.1A>C variant, which is already known to be pathogenic, we implicate six additional variants as American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) pathogenic supporting (PP) variants and classify ten as benign supporting (BP). In summary, we present a method developed for the classification of human PMS2 Kozak sequence variants that can contribute to the re-classification of VUSs identified in patients.

Association of Pseudomonas aeruginosa infection stage with lung function trajectory in children with cystic fibrosis.

BACKGROUND: Pseudomonas aeruginosa (Pa) infection in cystic fibrosis (CF) is characterized in stages: never (prior to first positive culture) to incident (first positive culture) to chronic. The association of Pa infection stage with lung function trajectory is poorly understood and the impact of age on this association has not been examined. We hypothesized that FEV1 decline would be slowest prior to Pa infection, intermediate after incident infection and greatest after chronic Pa infection. METHODS: Participants in a large US prospective cohort study diagnosed with CF prior to age 3 contributed data through the U.S. CF Patient Registry. Cubic spline linear mixed effects models were used to evaluate the longitudinal association of Pa stage (never, incident, chronic using 4 different definitions) with FEV1 adjusted for relevant covariates. Models contained interaction terms between age and Pa stage. RESULTS: 1,264 subjects born 1992-2006 provided a median 9.5 (IQR 0.25 to 15.75) years of follow up through 2017. 89% developed incident Pa; 39-58% developed chronic Pa depending on the definition. Compared to never Pa, incident Pa infection was associated with greater annual FEV1 decline and chronic Pa infection with the greatest FEV1 decline. The most rapid FEV1 decline and strongest association with Pa infection stage was seen in early adolescence (ages 12-15). CONCLUSIONS: Annual FEV1 decline worsens significantly with each Pa infection stage in children with CF. Our findings suggest that measures to prevent chronic infection, particularly during the high-risk period of early adolescence, could mitigate FEV1 decline and improve survival.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients.

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

Pleiotropic modifiers of age-related diabetes and neonatal intestinal obstruction in cystic fibrosis.

Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting ∼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF.

Accounting for population structure in genetic studies of cystic fibrosis.

CFTR F508del (c.1521_1523delCTT, p.Phe508delPhe) is the most common pathogenic allele underlying cystic fibrosis (CF), and its frequency varies in a geographic cline across Europe. We hypothesized that genetic variation associated with this cline is overrepresented in a large cohort (N > 5,000) of persons with CF who underwent whole-genome sequencing and that this pattern could result in spurious associations between variants correlated with both the F508del genotype and CF-related outcomes. Using principal-component (PC) analyses, we showed that variation in the CFTR region disproportionately contributes to a PC explaining a relatively high proportion of genetic variance. Variation near CFTR was correlated with population structure among persons with CF, and this correlation was driven by a subset of the sample inferred to have European ancestry. We performed genome-wide association studies comparing persons with CF with one versus two copies of the F508del allele; this allowed us to identify genetic variation associated with the F508del allele and to determine that standard PC-adjustment strategies eliminated the significant association signals. Our results suggest that PC adjustment can adequately prevent spurious associations between genetic variants and CF-related traits and are therefore effective tools to control for population structure even when population structure is confounded with disease severity and a common pathogenic variant.

Club cell secretory protein and lung function in children with cystic fibrosis.

BACKGROUND: Club cell secretory protein (CC16) exerts anti-inflammatory functions in lung disease. We sought to determine the relation of serum CC16 deficits and genetic variants that control serum CC16 to lung function among children with cystic fibrosis (CF). METHODS: We used longitudinal data from CF children (EPIC Study) with no positive cultures for Pseudomonas aeruginosa prior to enrollment. Circulating levels of CC16 and an inflammatory score (generated from CRP, SAA, calprotectin, G-CSF) were compared between participants with the lowest and highest FEV1 levels in adolescence (LLF and HLF groups, respectively; N = 130-per-group). Single nucleotide variants (SNVs) in the SCGB1A1, EHF-APIP loci were tested for association with circulating CC16 and with decline of FEV1 and FEV1/FVC% predicted levels between ages 7-16 using mixed models. RESULTS: Compared with the HLF group, the LLF group had lower levels of CC16 (geometric means: 8.2 vs 6.5 ng/ml, respectively; p = 0.0002) and higher levels of the normalized inflammatory score (-0.21 vs 0.21, p = 0.0007). Participants in the lowest CC16 and highest inflammation tertile had the highest odds for having LLF (p<0.0001 for comparison with participants in the highest CC16 and lowest inflammation tertile). Among seven SNVs associated with circulating CC16, the top SNV rs3741240 was associated with decline of FEV1/FVC and, marginally, FEV1 (p = 0.003 and 0.025, respectively; N = 611 participants, 20,801 lung function observations). CONCLUSIONS: Serum CC16 deficits are strongly associated with severity of CF lung disease and their effects are additive with systemic inflammation. The rs3741240 A allele is associated with low circulating CC16 and, possibly, accelerated lung function decline in CF.

Leveraging TOPMed imputation server and constructing a cohort-specific imputation reference panel to enhance genotype imputation among cystic fibrosis patients.

Cystic fibrosis (CF) is a severe genetic disorder that can cause multiple comorbidities affecting the lungs, the pancreas, the luminal digestive system and beyond. In our previous genome-wide association studies (GWAS), we genotyped approximately 8,000 CF samples using a mixture of different genotyping platforms. More recently, the Cystic Fibrosis Genome Project (CFGP) performed deep (approximately 30×) whole genome sequencing (WGS) of 5,095 samples to better understand the genetic mechanisms underlying clinical heterogeneity among patients with CF. For mixtures of GWAS array and WGS data, genotype imputation has proven effective in increasing effective sample size. Therefore, we first performed imputation for the approximately 8,000 CF samples with GWAS array genotype using the Trans-Omics for Precision Medicine (TOPMed) freeze 8 reference panel. Our results demonstrate that TOPMed can provide high-quality imputation for patients with CF, boosting genomic coverage from approximately 0.3-4.2 million genotyped markers to approximately 11-43 million well-imputed markers, and significantly improving polygenic risk score (PRS) prediction accuracy. Furthermore, we built a CF-specific CFGP reference panel based on WGS data of patients with CF. We demonstrate that despite having approximately 3% the sample size of TOPMed, our CFGP reference panel can still outperform TOPMed when imputing some CF disease-causing variants, likely owing to allele and haplotype differences between patients with CF and general populations. We anticipate our imputed data for 4,656 samples without WGS data will benefit our subsequent genetic association studies, and the CFGP reference panel built from CF WGS samples will benefit other investigators studying CF.

Comparing encounter-based and annualized chronic pseudomonas infection definitions in cystic fibrosis.

Chronic Pseudomonas aeruginosa (Pa) infection is associated with increased morbidity and mortality in people with cystic fibrosis (CF). There is no gold standard definition of chronic Pa infection in CF. We compared chronic Pa definitions using encounter-based versus annualized data in the Early Pseudomonas Infection Control (EPIC) Observational study cohort, and subsequently compared annualized chronic Pa definitions across a range of U.S. cohorts spanning decades of CF care. We found that an annualized chronic Pa definition requiring at least 1 Pa+ culture in 3 of 4 consecutive years ("Green 3/4") resulted in chronic Pa metrics similar to established encounter-based modified Leeds criteria definitions, including a similar age at and proportion who fulfilled chronic Pa criteria, and a similar proportion with sustained Pa infection after meeting the chronic Pa definition. The Green 3/4 chronic Pa definition will be valuable for longitudinal analyses in cohorts with limited culture frequency.

8q24 genetic variation and comprehensive haplotypes altering familial risk of prostate cancer.

The 8q24 genomic locus is tied to the origin of numerous cancers. We investigate its contribution to hereditary prostate cancer (HPC) in independent study populations of the Nashville Familial Prostate Cancer Study and International Consortium for Prostate Cancer Genetics (combined: 2,836 HPC cases, 2,206 controls of European ancestry). Here we report 433 variants concordantly associated with HPC in both study populations, accounting for 9% of heritability and modifying age of diagnosis as well as aggressiveness; 183 reach genome-wide significance. The variants comprehensively distinguish independent risk-altering haplotypes overlapping the 648 kb locus (three protective, and four risk (peak odds ratios: 1.5, 4, 5, and 22)). Sequence of the near-Mendelian haplotype reveals eleven causal mutation candidates. We introduce a linkage disequilibrium-based algorithm discerning eight independent sentinel variants, carrying considerable risk prediction ability (AUC = 0.625) for a single locus. These findings elucidate 8q24 locus structure and correlates for clinical prediction of prostate cancer risk.