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.

FAVOR: functional annotation of variants online resource and annotator for variation across the human genome.

Large biobank-scale whole genome sequencing (WGS) studies are rapidly identifying a multitude of coding and non-coding variants. They provide an unprecedented resource for illuminating the genetic basis of human diseases. Variant functional annotations play a critical role in WGS analysis, result interpretation, and prioritization of disease- or trait-associated causal variants. Existing functional annotation databases have limited scope to perform online queries and functionally annotate the genotype data of large biobank-scale WGS studies. We develop the Functional Annotation of Variants Online Resources (FAVOR) to meet these pressing needs. FAVOR provides a comprehensive multi-faceted variant functional annotation online portal that summarizes and visualizes findings of all possible nine billion single nucleotide variants (SNVs) across the genome. It allows for rapid variant-, gene- and region-level queries of variant functional annotations. FAVOR integrates variant functional information from multiple sources to describe the functional characteristics of variants and facilitates prioritizing plausible causal variants influencing human phenotypes. Furthermore, we provide a scalable annotation tool, FAVORannotator, to functionally annotate large-scale WGS studies and efficiently store the genotype and their variant functional annotation data in a single file using the annotated Genomic Data Structure (aGDS) format, making downstream analysis more convenient. FAVOR and FAVORannotator are available at https://favor.genohub.org.

The Genetic Landscape of Familial Pulmonary Fibrosis.

Rationale and Objectives: Up to 20% of idiopathic interstitial lung disease is familial, referred to as familial pulmonary fibrosis (FPF). An integrated analysis of FPF genetic risk was performed by comprehensively evaluating for genetic rare variants (RVs) in a large cohort of FPF kindreds. Methods: Whole-exome sequencing and/or candidate gene sequencing from affected individuals in 569 FPF kindreds was performed, followed by cosegregation analysis in large kindreds, gene burden analysis, gene-based risk scoring, cell-type enrichment analysis, and coexpression network construction. Measurements and Main Results: It was found that 14.9-23.4% of genetic risk in kindreds could be explained by RVs in genes previously linked to FPF, predominantly telomere-related genes. New candidate genes were identified in a small number of families-including SYDE1, SERPINB8, GPR87, and NETO1-and tools were developed for evaluation and prioritization of RV-containing genes across kindreds. Several pathways were enriched for RV-containing genes in FPF, including focal adhesion and mitochondrial complex I assembly. By combining single-cell transcriptomics with prioritized candidate genes, expression of RV-containing genes was discovered to be enriched in smooth muscle cells, type II alveolar epithelial cells, and endothelial cells. Conclusions: In the most comprehensive FPF genetic study to date, the prevalence of RVs in known FPF-related genes was defined, and new candidate genes and pathways relevant to FPF were identified. However, new RV-containing genes shared across multiple kindreds were not identified, thereby suggesting that heterogeneous genetic variants involving a variety of genes and pathways mediate genetic risk in most FPF kindreds.

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.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data.

INTRODUCTION: Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci. METHODS: We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants. RESULTS: Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses. DISCUSSION: This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Advancements in APOE and dementia research: Highlights from the 2023 AAIC Advancements: APOE conference.

INTRODUCTION: The apolipoprotein E gene (APOE) is an established central player in the pathogenesis of Alzheimer's disease (AD), with distinct apoE isoforms exerting diverse effects. apoE influences not only amyloid-beta and tau pathologies but also lipid and energy metabolism, neuroinflammation, cerebral vascular health, and sex-dependent disease manifestations. Furthermore, ancestral background may significantly impact the link between APOE and AD, underscoring the need for more inclusive research. METHODS: In 2023, the Alzheimer's Association convened multidisciplinary researchers at the "AAIC Advancements: APOE" conference to discuss various topics, including apoE isoforms and their roles in AD pathogenesis, progress in apoE-targeted therapeutic strategies, updates on disease models and interventions that modulate apoE expression and function. RESULTS: This manuscript presents highlights from the conference and provides an overview of opportunities for further research in the field. DISCUSSION: Understanding apoE's multifaceted roles in AD pathogenesis will help develop targeted interventions for AD and advance the field of AD precision medicine. HIGHLIGHTS: APOE is a central player in the pathogenesis of Alzheimer's disease. APOE exerts a numerous effects throughout the brain on amyloid-beta, tau, and other pathways. The AAIC Advancements: APOE conference encouraged discussions and collaborations on understanding the role of APOE.

Genetic examination of hematological parameters in SARS-CoV-2 infection and COVID-19.

People hospitalized with COVID-19 often exhibit altered hematological traits associated with disease prognosis (e.g., lower lymphocyte and platelet counts). We investigated whether inter-individual variability in baseline hematological traits influences risk of acute SARS-CoV-2 infection or progression to severe COVID-19. We report inconsistent associations between blood cell traits with incident SARS-CoV-2 infection and severe COVID-19 in UK Biobank and the Vanderbilt University Medical Center Synthetic Derivative (VUMC SD). Since genetically determined blood cell measures better represent cell abundance across the lifecourse, we also assessed the shared genetic architecture of baseline blood cell traits on COVID-19 related outcomes by Mendelian randomization (MR) analyses. We found significant relationships between COVID-19 severity and mean sphered cell volume after adjusting for multiple testing. However, MR results differed significantly across different freezes of COVID-19 summary statistics and genetic correlation between these traits was modest (0.1), decreasing our confidence in these results. We observed overlapping genetic association signals between other hematological and COVID-19 traits at specific loci such as MAPT and TYK2. In conclusion, we did not find convincing evidence of relationships between the genetic architecture of blood cell traits and either SARS-CoV-2 infection or COVID-19 hospitalization, though we do see evidence of shared signals at specific loci.

Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants.

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

Centers for Mendelian Genomics: A decade of facilitating gene discovery.

PURPOSE: Mendelian disease genomic research has undergone a massive transformation over the past decade. With increasing availability of exome and genome sequencing, the role of Mendelian research has expanded beyond data collection, sequencing, and analysis to worldwide data sharing and collaboration. METHODS: Over the past 10 years, the National Institutes of Health-supported Centers for Mendelian Genomics (CMGs) have played a major role in this research and clinical evolution. RESULTS: We highlight the cumulative gene discoveries facilitated by the program, biomedical research leveraged by the approach, and the larger impact on the research community. Beyond generating a list of gene-phenotype relationships and participating in widespread data sharing, the CMGs have created resources, tools, and training for the larger community to foster understanding of genes and genome variation. The CMGs have participated in a wide range of data sharing activities, including deposition of all eligible CMG data into the Analysis, Visualization, and Informatics Lab-space (AnVIL), sharing candidate genes through the Matchmaker Exchange and the CMG website, and sharing variants in Genotypes to Mendelian Phenotypes (Geno2MP) and VariantMatcher. CONCLUSION: The work is far from complete; strengthening communication between research and clinical realms, continued development and sharing of knowledge and tools, and improving access to richly characterized data sets are all required to diagnose the remaining molecularly undiagnosed patients.

Exome sequencing identifies genetic variants in anophthalmia and microphthalmia.

Anophthalmia and microphthalmia (A/M) are rare birth defects affecting up to 2 per 10,000 live births. These conditions are manifested by the absence of an eye or reduced eye volumes within the orbit leading to vision loss. Although clinical case series suggest a strong genetic component in A/M, few systematic investigations have been conducted on potential genetic contributions owing to low population prevalence. To overcome this challenge, we utilized DNA samples and data collected as part of the National Birth Defects Prevention Study (NBDPS). The NBDPS employed multi-center ascertainment of infants affected by A/M. We performed exome sequencing on 67 family trios and identified numerous genes affected by rare deleterious nonsense and missense variants in this cohort, including de novo variants. We identified 9 nonsense changes and 86 missense variants that are absent from the reference human population (Genome Aggregation Database), and we suggest that these are high priority candidate genes for A/M. We also performed literature curation, single cell transcriptome comparisons, and molecular pathway analysis on the candidate genes and performed protein structure modeling to determine the potential pathogenic variant consequences on PAX6 in this disease.

Genetic counseling for early onset and familial dementia: Patient perspectives on exome sequencing.

Genetic testing has become routine for many inherited conditions; however, little is known about the unique issues that arise when offering genetic testing for inherited forms of dementia. To better understand the patient perspective, we surveyed study participants about their experiences as they underwent genetic counseling and genetic testing for dementia. We recruited 50 pairs of subjects. Each pair was comprised of one person with cognitive impairment and a cognitively intact co-participant. Study participants received pre- and post-test genetic counseling and comprehensive genetic testing for dementia. During the study, participant pairs completed four surveys which asked about their experience. Testing began with a 38 gene dementia panel. Participants with negative panel results or variants of uncertain significance (VUS) were reflexed to exome sequencing (ES). Twenty-nine participants (58%) reported that their primary motivation to join the study was for the benefit to their families. Fifty-two percent of participants initially planned to use their test results to make health and wellness changes, but, six months after disclosure, only 31% had done so. Six months after result disclosure, approximately 90% of participant pairs accurately recalled their genetic test results. Overall satisfaction with testing was high, and decision regret was negligible. This observational study describes the experiences of study participants undergoing genetic counseling and genetic testing for dementia and found that most participant pairs accurately recalled their results up to six months following disclosure while also maintaining high levels of satisfaction without decision regret. These findings suggest that, in the context of genetic counseling, genetic testing can be effectively used in this population.

Non-coding variants in MYH11, FZD3, and SORCS3 are associated with dementia in women.

INTRODUCTION: Recent studies suggest that both sex-specific genetic risk factors and those shared between dementia and stroke are involved in dementia pathogenesis. METHODS: We performed both single-variant and gene-based genome-wide association studies of >11,000 whole genome sequences from the Women's Health Initiative cohort to discover loci associated with dementia, with adjustment for age, ethnicity, stroke, and venous thromboembolism status. Evidence for prior evidence of association and differential gene expression in dementia-related tissues and samples was gathered for each locus. RESULTS: Our multiethnic studies identified significant associations between variants within APOE, MYH11, FZD3, SORCS3, and GOLGA8B and risk of dementia. Ten genes implicated by these loci, including MYH11, FZD3, SORCS3, and GOLGA8B, were differentially expressed in the context of Alzheimer's disease. DISCUSSION: Our association of MYH11, FZD3, SORCS3, and GOLGA8B with dementia is supported by independent functional studies in human subjects, model systems, and associations with shared risk factors for stroke and dementia.

Mutations in GDF11 and the extracellular antagonist, Follistatin, as a likely cause of Mendelian forms of orofacial clefting in humans.

Cleft lip with or without cleft palate (CL/P) is generally viewed as a complex trait with multiple genetic and environmental contributions. In 70% of cases, CL/P presents as an isolated feature and/or deemed nonsyndromic. In the remaining 30%, CL/P is associated with multisystem phenotypes or clinically recognizable syndromes, many with a monogenic basis. Here we report the identification, via exome sequencing, of likely pathogenic variants in two genes that encode interacting proteins previously only linked to orofacial clefting in mouse models. A variant in GDF11 (encoding growth differentiation factor 11), predicting a p.(Arg298Gln) substitution at the Furin protease cleavage site, was identified in one family that segregated with CL/P and both rib and vertebral hypersegmentation, mirroring that seen in Gdf11 knockout mice. In the second family in which CL/P was the only phenotype, a mutation in FST (encoding the GDF11 antagonist, Follistatin) was identified that is predicted to result in a p.(Cys56Tyr) substitution in the region that binds GDF11. Functional assays demonstrated a significant impact of the specific mutated amino acids on FST and GDF11 function and, together with embryonic expression data, provide strong evidence for the importance of GDF11 and Follistatin in the regulation of human orofacial development.

Functional annotation of genomic variants in studies of late-onset Alzheimer's disease.

Motivation: Annotation of genomic variants is an increasingly important and complex part of the analysis of sequence-based genomic analyses. Computational predictions of variant function are routinely incorporated into gene-based analyses of rare-variants, though to date most studies use limited information for assessing variant function that is often agnostic of the disease being studied. Results: In this work, we outline an annotation process motivated by the Alzheimer's Disease Sequencing Project, illustrate the impact of including tissue-specific transcript sets and sources of gene regulatory information and assess the potential impact of changing genomic builds on the annotation process. While these factors only impact a small proportion of total variant annotations (∼5%), they influence the potential analysis of a large fraction of genes (∼25%). Availability and implementation: Individual variant annotations are available via the NIAGADS GenomicsDB, at https://www.niagads.org/genomics/ tools-and-software/databases/genomics-database. Annotations are also available for bulk download at https://www.niagads.org/datasets. Annotation processing software is available at http://www.icompbio.net/resources/software-and-downloads/. Supplementary information: Supplementary data are available at Bioinformatics online.

Local ancestry at APOE modifies Alzheimer's disease risk in Caribbean Hispanics.

INTRODUCTION: Although the relationship between APOE and Alzheimer's disease (AD) is well established in populations of European descent, the effects of APOE and ancestry on AD risk in diverse populations is not well understood. METHODS: Logistic mixed model regression and survival analyses were performed in a sample of 3067 Caribbean Hispanics and 3028 individuals of European descent to assess the effects of APOE genotype, local ancestry, and genome-wide ancestry on AD risk and age at onset. RESULTS: Among the Caribbean Hispanics, individuals with African-derived ancestry at APOE had 39% lower odds of AD than individuals with European-derived APOE, after adjusting for APOE genotype, age, and genome-wide ancestry. While APOE E2 and E4 effects on AD risk and age at onset were significant in the Caribbean Hispanics, they were substantially attenuated compared with those in European ancestry individuals. DISCUSSION: These results suggest that additional genetic variation in the APOE region influences AD risk beyond APOE E2/E3/E4.

Genetic Variation in Genes Underlying Diverse Dementias May Explain a Small Proportion of Cases in the Alzheimer's Disease Sequencing Project.

BACKGROUND/AIMS: The Alzheimer's Disease Sequencing Project (ADSP) aims to identify novel genes influencing Alzheimer's disease (AD). Variants within genes known to cause dementias other than AD have previously been associated with AD risk. We describe evidence of co-segregation and associations between variants in dementia genes and clinically diagnosed AD within the ADSP. METHODS: We summarize the properties of known pathogenic variants within dementia genes, describe the co-segregation of variants annotated as "pathogenic" in ClinVar and new candidates observed in ADSP families, and test for associations between rare variants in dementia genes in the ADSP case-control study. The participants were clinically evaluated for AD, and they represent European, Caribbean Hispanic, and isolate Dutch populations. RESULTS/CONCLUSIONS: Pathogenic variants in dementia genes were predominantly rare and conserved coding changes. Pathogenic variants within ARSA, CSF1R, and GRN were observed, and candidate variants in GRN and CHMP2B were nominated in ADSP families. An independent case-control study provided evidence of an association between variants in TREM2, APOE, ARSA, CSF1R, PSEN1, and MAPT and risk of AD. Variants in genes which cause dementing disorders may influence the clinical diagnosis of AD in a small proportion of cases within the ADSP.

The power of representation: Statistical analysis of diversity in US Alzheimer's disease genetics data.

INTRODUCTION: Alzheimer's disease (AD) is a complex disease influenced by genetics and environment. More than 75 susceptibility loci have been linked to late-onset AD, but most of these loci were discovered in genome-wide association studies (GWAS) exclusive to non-Hispanic White individuals. There are wide disparities in AD risk across racially stratified groups, and while these disparities are not due to genetic differences, underrepresentation in genetic research can further exacerbate and contribute to their persistence. We investigated the racial/ethnic representation of participants in United States (US)-based AD genetics and the statistical implications of current representation. METHODS: We compared racial/ethnic data of participants from array and sequencing studies in US AD genetics databases, including National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site (NIAGADS) and NIAGADS Data Sharing Service (dssNIAGADS), to AD and related dementia (ADRD) prevalence and mortality. We then simulated the statistical power of these datasets to identify risk variants from non-White populations. RESULTS: There is insufficient statistical power (probability <80%) to detect single nucleotide polymorphisms (SNPs) with low to moderate effect sizes (odds ratio [OR]<1.5) using array data from Black and Hispanic participants; studies of Asian participants are not powered to detect variants OR <= 2. Using available and projected sequencing data from Black and Hispanic participants, risk variants with OR = 1.2 are detectable at high allele frequencies. Sample sizes remain insufficiently powered to detect these variants in Asian populations. DISCUSSION: AD genetics datasets are largely representative of US ADRD burden. However, there is a wide discrepancy between proportional representation and statistically meaningful representation. Most variation identified in GWAS of non-Hispanic White individuals have low to moderate effects. Comparable risk variants in non-White populations are not detectable given current sample sizes, which could lead to disparities in future studies and drug development. We urge AD genetics researchers and institutions to continue investing in recruiting diverse participants and use community-based participatory research practices.

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.

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.

Exome sequencing identifies novel genes underlying primary congenital glaucoma in the National Birth Defects Prevention Study.

BACKGROUND: Primary congenital glaucoma (PCG) affects approximately 1 in 10,000 live born infants in the United States (U.S.). PCG has a autosomal recessive inheritance pattern, and variable expressivity and reduced penetrance have been reported. Likely causal variants in the most commonly mutated gene, CYP1B1, are less prevalent in the U.S., suggesting that alternative genes may contribute to the condition. This study utilized exome sequencing to investigate the genetic architecture of PCG in the U.S. and to identify novel genes and variants. METHODS: We studied 37 family trios where infants had PCG and were part of the National Birth Defects Prevention Study (births 1997-2011), a U.S. multicenter study of birth defects. Samples underwent exome sequencing and sequence reads were aligned to the human reference sample (NCBI build 37/hg19). Variant filtration was conducted under de novo and Mendelian inheritance models using GEMINI. RESULTS: Among candidate variants, CYP1B1 was most represented (five trios, 13.5%). Twelve probands (32%) had potentially pathogenic variants in other genes not previously linked to PCG but important in eye development and/or to underlie Mendelian conditions with potential phenotypic overlap (e.g., CRYBB2, RXRA, GLI2). CONCLUSION: Variation in the genes identified in this population-based study may help to further explain the genetics of PCG.