Whole exome sequence analysis in 51 624 participants identifies novel genes and variants associated with refractive error and myopia.

Refractive errors are associated with a range of pathological conditions, such as myopic maculopathy and glaucoma, and are highly heritable. Studies of missense and putative loss of function (pLOF) variants identified via whole exome sequencing (WES) offer the prospect of directly implicating potentially causative disease genes. We performed a genome-wide association study for refractive error in 51 624 unrelated adults, of European ancestry, aged 40-69 years from the UK and genotyped using WES. After testing 29 179 pLOF and 495 263 missense variants, 1 pLOF and 18 missense variants in 14 distinct genomic regions were taken forward for fine-mapping analysis. This yielded 19 putative causal variants of which 18 had a posterior inclusion probability >0.5. Of the 19 putative causal variants, 12 were novel discoveries. Specific variants were associated with a more myopic refractive error, while others were associated with a more hyperopic refractive error. Association with age of onset of spectacle wear (AOSW) was examined in an independent validation sample (38 100 early AOSW cases and 74 243 controls). Of 11 novel variants that could be tested, 8 (73%) showed evidence of association with AOSW status. This work identified COL4A4 and ATM as novel candidate genes associated with refractive error. In addition, novel putative causal variants were identified in the genes RASGEF1, ARMS2, BMP4, SIX6, GSDMA, GNGT2, ZNF652 and CRX. Despite these successes, the study also highlighted the limitations of community-based WES studies compared with high myopia case-control WES studies.

Transmission disequilibrium analysis of whole genome data in childhood-onset systemic lupus erythematosus.

Childhood-onset systemic lupus erythematosus (cSLE) patients are unique, with hallmarks of Mendelian disorders (early-onset and severe disease) and thus are an ideal population for genetic investigation of SLE. In this study, we use the transmission disequilibrium test (TDT), a family-based genetic association analysis that employs robust methodology, to analyze whole genome sequencing data. We aim to identify novel genetic associations in an ancestrally diverse, international cSLE cohort. Forty-two cSLE patients and 84 unaffected parents from 3 countries underwent whole genome sequencing. First, we performed TDT with single nucleotide variant (SNV)-based (common variants) using PLINK 1.9, and gene-based (rare variants) analyses using Efficient and Parallelizable Association Container Toolbox (EPACTS) and rare variant TDT (rvTDT), which applies multiple gene-based burden tests adapted for TDT, including the burden of rare variants test. Applying the GWAS standard threshold (5.0 × 10-8) to common variants, our SNV-based analysis did not return any genome-wide significant SNVs. The rare variant gene-based TDT analysis identified many novel genes significantly enriched in cSLE patients, including HNRNPUL2, a DNA repair protein, and DNAH11, a ciliary movement protein, among others. Our approach identifies several novel SLE susceptibility genes in an ancestrally diverse childhood-onset lupus cohort.

A rare FGF5 candidate variant (rs112475347) for predisposition to nonsquamous, nonsmall-cell lung cancer.

A unique approach with rare resources was used to identify candidate variants predisposing to familial nonsquamous nonsmall-cell lung cancers (NSNSCLC). We analyzed sequence data from NSNSCLC-affected cousin pairs belonging to high-risk lung cancer pedigrees identified in a genealogy of Utah linked to statewide cancer records to identify rare, shared candidate predisposition variants. Variants were tested for association with lung cancer risk in UK Biobank. Evidence for linkage with lung cancer was also reviewed in families from the Genetic Epidemiology of Lung Cancer Consortium. Protein prediction modeling compared the mutation with reference. We sequenced NSNSCLC-affected cousin pairs from eight high-risk lung cancer pedigrees and identified 66 rare candidate variants shared in the cousin pairs. One variant in the FGF5 gene also showed significant association with lung cancer in UKBiobank. This variant was observed in 3/163 additional sampled Utah lung cancer cases, 2 of whom were related in another independent pedigree. Modeling of the predicted protein predicted a second binding site for SO4 that may indicate binding differences. This unique study identified multiple candidate predisposition variants for NSNSCLC, including a rare variant in FGF5 that was significantly associated with lung cancer risk and that segregated with lung cancer in the two pedigrees in which it was observed. FGF5 is an oncogenic factor in several human cancers, and the mutation found here (W81C) changes the binding ability of heparan sulfate to FGF5, which might lead to its deregulation. These results support FGF5 as a potential NSNSCLC predisposition gene and present additional candidate predisposition variants.

What makes a good prediction? Feature importance and beginning to open the black box of machine learning in genetics.

Genetic data have become increasingly complex within the past decade, leading researchers to pursue increasingly complex questions, such as those involving epistatic interactions and protein prediction. Traditional methods are ill-suited to answer these questions, but machine learning (ML) techniques offer an alternative solution. ML algorithms are commonly used in genetics to predict or classify subjects, but some methods evaluate which features (variables) are responsible for creating a good prediction; this is called feature importance. This is critical in genetics, as researchers are often interested in which features (e.g., SNP genotype or environmental exposure) are responsible for a good prediction. This allows for the deeper analysis beyond simple prediction, including the determination of risk factors associated with a given phenotype. Feature importance further permits the researcher to peer inside the black box of many ML algorithms to see how they work and which features are critical in informing a good prediction. This review focuses on ML methods that provide feature importance metrics for the analysis of genetic data. Five major categories of ML algorithms: k nearest neighbors, artificial neural networks, deep learning, support vector machines, and random forests are described. The review ends with a discussion of how to choose the best machine for a data set. This review will be particularly useful for genetic researchers looking to use ML methods to answer questions beyond basic prediction and classification.

Linear mixed models for association analysis of quantitative traits with next-generation sequencing data.

We develop linear mixed models (LMMs) and functional linear mixed models (FLMMs) for gene-based tests of association between a quantitative trait and genetic variants on pedigrees. The effects of a major gene are modeled as a fixed effect, the contributions of polygenes are modeled as a random effect, and the correlations of pedigree members are modeled via inbreeding/kinship coefficients. F -statistics and χ 2 likelihood ratio test (LRT) statistics based on the LMMs and FLMMs are constructed to test for association. We show empirically that the F -distributed statistics provide a good control of the type I error rate. The F -test statistics of the LMMs have similar or higher power than the FLMMs, kernel-based famSKAT (family-based sequence kernel association test), and burden test famBT (family-based burden test). The F -statistics of the FLMMs perform well when analyzing a combination of rare and common variants. For small samples, the LRT statistics of the FLMMs control the type I error rate well at the nominal levels α = 0.01 and 0.05 . For moderate/large samples, the LRT statistics of the FLMMs control the type I error rates well. The LRT statistics of the LMMs can lead to inflated type I error rates. The proposed models are useful in whole genome and whole exome association studies of complex traits.

Genome-wide scans of myopia in Pennsylvania Amish families reveal significant linkage to 12q15, 8q21.3 and 5p15.33.

Myopia is one of the most common ocular disorders in the world, yet the genetic etiology of the disease remains poorly understood. Specialized founder populations, such as the Pennsylvania Amish, provide the opportunity to utilize exclusive genomic architecture, like unique haplotypes, to better understand the genetic causes of myopia. We perform genetic linkage analysis on Pennsylvania Amish families that have a strong familial history of myopia to map any potential causal variants and genes for the disease. 293 individuals from 25 extended families were genotyped on the Illumina ExomePlus array and merged with previous microsatellite data. We coded myopia affection as a binary phenotype; myopia was defined as having a mean spherical equivalent (MSE) of less than or equal to - 1 D (diopters). Two-point and multipoint parametric linkage analyses were performed under an autosomal dominant model. When allowing for locus heterogeneity, we identified two novel genome-wide significantly linked variants at 12q15 (heterogeneity LOD, HLOD = 3.77) in PTPRB and at 8q21.3 (HLOD = 3.35) in CNGB3. We identified further three genome-wide significant variants within a single family. These three variants were located in exons of SLC6A18 at 5p15.33 (LODs ranged from 3.51 to 3.37). Multipoint analysis confirmed the significant signal at 5p15.33 with six genome-wide significant variants (LODs ranged from 3.6 to 3.3). Further suggestive evidence of linkage was observed in several other regions of the genome. All three novel linked regions contain strong candidate genes, especially CNGB3 on 8q21.3, which has been shown to affect photoreceptors and cause complete color blindness. Whole genome sequencing on these regions is planned to conclusively elucidate the causal variants.

Exome genotyping and linkage analysis identifies two novel linked regions and replicates two others for myopia in Ashkenazi Jewish families.

BACKGROUND: Myopia is one of most common eye diseases in the world and affects 1 in 4 Americans. It is a complex disease caused by both environmental and genetics effects; the genetics effects are still not well understood. In this study, we performed genetic linkage analyses on Ashkenazi Jewish families with a strong familial history of myopia to elucidate any potential causal genes. METHODS: Sixty-four extended Ashkenazi Jewish families were previously collected from New Jersey. Genotypes from the Illumina ExomePlus array were merged with prior microsatellite linkage data from these families. Additional custom markers were added for candidate regions reported in literature for myopia or refractive error. Myopia was defined as mean spherical equivalent (MSE) of -1D or worse and parametric two-point linkage analyses (using TwoPointLods) and multi-point linkage analyses (using SimWalk2) were performed as well as collapsed haplotype pattern (CHP) analysis in SEQLinkage and association analyses performed with FBAT and rv-TDT. RESULTS: Strongest evidence of linkage was on 1p36(two-point LOD = 4.47) a region previously linked to refractive error (MYP14) but not myopia. Another genome-wide significant locus was found on 8q24.22 with a maximum two-point LOD score of 3.75. CHP analysis also detected the signal on 1p36, localized to the LINC00339 gene with a maximum HLOD of 3.47, as well as genome-wide significant signals on 7q36.1 and 11p15, which overlaps with the MYP7 locus. CONCLUSIONS: We identified 2 novel linkage peaks for myopia on chromosomes 7 and 8 in these Ashkenazi Jewish families and replicated 2 more loci on chromosomes 1 and 11, one previously reported in refractive error but not myopia in these families and the other locus previously reported in the literature. Strong candidate genes have been identified within these linkage peaks in our families. Targeted sequencing in these regions will be necessary to definitively identify causal variants under these linkage peaks.

Myopia in Chinese families shows linkage to 10q26.13.

Purpose: To determine genetic linkage between myopia and Han Chinese patients with a family history of the disease. Methods: One hundred seventy-six Han Chinese patients from 34 extended families were given eye examinations, and mean spherical equivalent (MSE) in diopters (D) was calculated by adding the spherical component of the refraction to one-half the cylindrical component and taking the average of both eyes. The MSE was converted to a binary phenotype, where all patients with an MSE of -1.00 D or less were coded as affected. Unaffected individuals had an MSE greater than 0.00 D (ages 21 years and up), +1.50 (ages 11-20), or +2.00 D (ages 6-10 years). Individuals between the given upper threshold and -1.00 were coded as unknown. Patients were genotyped on an exome chip. Three types of linkage analyses were performed: single-variant two-point, multipoint, and collapsed haplotype pattern (CHP) variant two-point. Results: The CHP variant two-point results identified a significant peak (heterogeneity logarithm of the odds [HLOD] = 3.73) at 10q26.13 in TACC2. The single-variant two-point and multipoint analyses showed highly suggestive linkage to the same region. The single-variant two-point results identified 25 suggestive variants at HTRA1, also at 10q26.13. Conclusions: We report a significant genetic linkage between myopia and Han Chinese patients at 10q26.13. 10q26.13 contains several good candidate genes, such as TACC2 and the known age-related macular degeneration gene HTRA1. Targeted sequencing of the region is planned to identify the causal variant(s).

Parametric Linkage Analysis Identifies Five Novel Genome-Wide Significant Loci for Familial Lung Cancer.

OBJECTIVE: One of four American cancer patients dies of lung cancer. Environmental factors such as tobacco smoking are known to affect lung cancer risk. However, there is a genetic factor to lung cancer risk as well. Here, we perform parametric linkage analysis on family-based genotype data in an effort to find genetic loci linked to the disease. METHODS: 197 individuals from families with a high-risk history of lung cancer were recruited and genotyped using an Illumina array. Parametric linkage analyses were performed using an affected-only phenotype model with an autosomal dominant inheritance using a disease allele frequency of 0.01. Three types of analyses were performed: single variant two-point, collapsed haplotype pattern variant two-point, and multipoint analysis. RESULTS: Five novel genome-wide significant loci were identified at 18p11.23, 2p22.2, 14q13.1, 16p13, and 20q13.11. The families most informative for linkage were also determined. CONCLUSIONS: The 5 novel signals are good candidate regions, containing genes that have been implicated as having somatic changes in lung cancer or other cancers (though not in germ line cells). Targeted sequencing on the significant loci is planned to determine the causal variants at these loci.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape.

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

Rare variant analyses across multiethnic cohorts identify novel genes for refractive error.

Refractive error, measured here as mean spherical equivalent (SER), is a complex eye condition caused by both genetic and environmental factors. Individuals with strong positive or negative values of SER require spectacles or other approaches for vision correction. Common genetic risk factors have been identified by genome-wide association studies (GWAS), but a great part of the refractive error heritability is still missing. Some of this heritability may be explained by rare variants (minor allele frequency [MAF] ≤ 0.01.). We performed multiple gene-based association tests of mean Spherical Equivalent with rare variants in exome array data from the Consortium for Refractive Error and Myopia (CREAM). The dataset consisted of over 27,000 total subjects from five cohorts of Indo-European and Eastern Asian ethnicity. We identified 129 unique genes associated with refractive error, many of which were replicated in multiple cohorts. Our best novel candidates included the retina expressed PDCD6IP, the circadian rhythm gene PER3, and P4HTM, which affects eye morphology. Future work will include functional studies and validation. Identification of genes contributing to refractive error and future understanding of their function may lead to better treatment and prevention of refractive errors, which themselves are important risk factors for various blinding conditions.

Genome-wide linkage search for cancer susceptibility loci in a cohort of non BRCA1/2 families in Sri Lanka.

OBJECTIVE: Although linkage studies have been utilized for the identification of variants associated with cancer in the world, little is known about their role in non BRCA1/2 individuals in the Sri Lankans. Hence we performed linkage analysis to identify susceptibility loci related to the inherited risk of cancer in a cohort of Sri Lankans affected with hereditary breast cancer. The Illumina global screening array having 654,027 single nucleotide polymorphism markers was performed in four families, in which at least three individuals within third degree relatives were affected by breast cancer. Two-point parametric linkage analysis was conducted assuming disease allele frequency of 1%. Penetrance was set at 90% for carriers with a 10% phenocopy rate. RESULTS: Thirty-one variants exhibited genome-wide suggestive HLODs. The top overall HLOD score was at rs1856277, an intronic variant in MYO16 on chromosome 13. The two most informative families also suggested several candidate linked loci in genes, including ERAP1, RPRM, WWOX, CDH1, EXOC1, HUS1B, STIM1 and TUSC1. This study provides the first step in identifying germline variants that may be involved in risk of cancer in cancer-aggregated non-BRCA1/2 families from the understudied Sri Lankan population. Several candidate linked regions showed suggestive evidence of linkage to cancer risk.

Targeted Sequencing of Candidate Regions Associated with Sagittal and Metopic Nonsyndromic Craniosynostosis.

Craniosynostosis (CS) is a major birth defect in which one or more skull sutures fuse prematurely. We previously performed a genome-wide association study (GWAS) for sagittal non-syndromic CS (sNCS), identifying associations downstream from BMP2 on 20p12.3 and intronic to BBS9 on 7p14.3; analyses of imputed variants in DLG1 on 3q29 were also genome-wide significant. We followed this work with a GWAS for metopic non-syndromic NCS (mNCS), discovering a significant association intronic to BMP7 on 20q13.31. In the current study, we sequenced the associated regions on 3q29, 7p14.3, and 20p12.3, including two candidate genes (BMP2 and BMPER) near some of these regions in 83 sNCS child-parent trios, and sequenced regions on 7p14.3 and 20q13.2-q13.32 in 80 mNCS child-parent trios. These child-parent trios were selected from the original GWAS cohorts if the probands carried at least one copy of the top associated GWAS variant (rs1884302 C allele for sNCS; rs6127972 T allele for mNCS). Many of the variants sequenced in these targeted regions are strongly predicted to be within binding sites for transcription factors involved in craniofacial development or bone morphogenesis. Variants enriched in more than one trio and predicted to be damaging to gene function are prioritized for functional studies.

Association analyses of rare variants identify two genes associated with refractive error.

PURPOSE: Genetic variants identified through population-based genome-wide studies are generally of high frequency, exerting their action in the central part of the refractive error spectrum. However, the power to identify associations with variants of lower minor allele frequency is greatly reduced, requiring considerable sample sizes. Here we aim to assess the impact of rare variants on genetic variation of refractive errors in a very large general population cohort. METHODS: Genetic association analyses of non-cyclopaedic autorefraction calculated as mean spherical equivalent (SPHE) used whole-exome sequence genotypic information from 50,893 unrelated participants in the UK Biobank of European ancestry. Gene-based analyses tested for association with SPHE using an optimised SNP-set kernel association test (SKAT-O) restricted to rare variants (minor allele frequency < 1%) within protein-coding regions of the genome. All models were adjusted for age, sex and common lead variants within the same locus reported by previous genome-wide association studies. Potentially causal markers driving association at significant loci were elucidated using sensitivity analyses by sequentially dropping the most associated variants from gene-based analyses. RESULTS: We found strong statistical evidence for association of SPHE with the SIX6 (p-value = 2.15 x 10-10, or Bonferroni-Corrected p = 4.41x10-06) and the CRX gene (p-value = 6.65 x 10-08, or Bonferroni-Corrected p = 0.001). The SIX6 gene codes for a transcription factor believed to be critical to the eye, retina and optic disc development and morphology, while CRX regulates photoreceptor specification and expression of over 700 genes in the retina. These novel associations suggest an important role of genes involved in eye morphogenesis in refractive error. CONCLUSION: The results of our study support previous research highlighting the importance of rare variants to the genetic risk of refractive error. We explain some of the origins of the genetic signals seen in GWAS but also report for the first time a completely novel association with the CRX gene.

Myopia in African Americans Is Significantly Linked to Chromosome 7p15.2-14.2.

Purpose: The purpose of this study was to perform genetic linkage analysis and association analysis on exome genotyping from highly aggregated African American families with nonpathogenic myopia. African Americans are a particularly understudied population with respect to myopia. Methods: One hundred six African American families from the Philadelphia area with a family history of myopia were genotyped using an Illumina ExomePlus array and merged with previous microsatellite data. Myopia was initially measured in mean spherical equivalent (MSE) and converted to a binary phenotype where individuals were identified as affected, unaffected, or unknown. Parametric linkage analysis was performed on both individual variants (single-nucleotide polymorphisms [SNPs] and microsatellites) as well as gene-based markers. Family-based association analysis and transmission disequilibrium test (TDT) analysis modified for rare variants was also performed. Results: Genetic linkage analysis identified 2 genomewide significant variants at 7p15.2 and 7p14.2 (in the intergenic region between MIR148A and NFE2L3 and in the noncoding RNA LOC401324) and 2 genomewide significant genes (CRHR2 and AVL9) both at 7p14.3. No genomewide results were found in the association analyses. Conclusions: This study identified a significant linkage peak in African American families for myopia at 7p15.2 to 7p14.2, the first potential risk locus for myopia in African Americans. Interesting candidate genes are located in the region, including PDE1C, which is highly expressed in the eyes, and known to be involved in retinal development. Further identification of the causal variants at this linkage peak will help elucidate the genetics of myopia in this understudied population.

Genetic Variation and Recurrent Haplotypes on Chromosome 6q23-25 Risk Locus in Familial Lung Cancer.

Although lung cancer is known to be caused by environmental factors, it has also been shown to have genetic components, and the genetic etiology of lung cancer remains understudied. We previously identified a lung cancer risk locus on 6q23-25 using microsatellite data in families with a history of lung cancer. To further elucidate that signal, we performed targeted sequencing on nine of our most strongly linked families. Two-point linkage analysis of the sequencing data revealed that the signal was heterogeneous and that different families likely had different risk variants. Three specific haplotypes were shared by some of the families: 6q25.3-26 in families 42 and 44, 6q25.2-25.3 in families 47 and 59, and 6q24.2-25.1 in families 30, 33, and 35. Region-based logarithm of the odds scores and expression data identified the likely candidate genes for each haplotype overlap: ARID1B at 6q25.3, MAP3K4 at 6q26, and UTRN (6q24.1) and PHACTR2 (6q24.2). Further annotation was used to zero in on potential risk variants in those genes. All four genes are good candidate genes for lung cancer risk, having been linked to either lung cancer specifically or other cancers. However, this is the first time any of these genes has been implicated in germline risk. Functional analysis of these four genes is planned for future work. SIGNIFICANCE: This study identifies four genes associated with lung cancer risk, which could help guide future lung cancer prevention and treatment approaches.

IMI 2021 Yearly Digest.

Purpose: The International Myopia Institute (IMI) Yearly Digest highlights new research considered to be of importance since the publication of the first series of IMI white papers. Methods: A literature search was conducted for articles on myopia between 2019 and mid-2020 to inform definitions and classifications, experimental models, genetics, interventions, clinical trials, and clinical management. Conference abstracts from key meetings in the same period were also considered. Results: One thousand articles on myopia have been published between 2019 and mid-2020. Key advances include the use of the definition of premyopia in studies currently under way to test interventions in myopia, new definitions in the field of pathologic myopia, the role of new pharmacologic treatments in experimental models such as intraocular pressure-lowering latanoprost, a large meta-analysis of refractive error identifying 336 new genetic loci, new clinical interventions such as the defocus incorporated multisegment spectacles and combination therapy with low-dose atropine and orthokeratology (OK), normative standards in refractive error, the ethical dilemma of a placebo control group when myopia control treatments are established, reporting the physical metric of myopia reduction versus a percentage reduction, comparison of the risk of pediatric OK wear with risk of vision impairment in myopia, the justification of preventing myopic and axial length increase versus quality of life, and future vision loss. Conclusions: Large amounts of research in myopia have been published since the IMI 2019 white papers were released. The yearly digest serves to highlight the latest research and advances in myopia.

Whole Exome Sequencing of Highly Aggregated Lung Cancer Families Reveals Linked Loci for Increased Cancer Risk on Chromosomes 12q, 7p, and 4q.

BACKGROUND: Lung cancer kills more people than any other cancer in the United States. In addition to environmental factors, lung cancer has genetic risk factors as well, though the genetic etiology is still not well understood. We have performed whole exome sequencing on 262 individuals from 28 extended families with a family history of lung cancer. METHODS: Parametric genetic linkage analysis was performed on these samples using two distinct analyses-the lung cancer only (LCO) analysis, where only patients with lung cancer were coded as affected, and the all aggregated cancers (AAC) analysis, where other cancers seen in the pedigree were coded as affected. RESULTS: The AAC analysis yielded a genome-wide significant result at rs61943670 in POLR3B at 12q23.3. POLR3B has been implicated somatically in lung cancer, but this germline finding is novel and is a significant expression quantitative trait locus in lung tissue. Interesting genome-wide suggestive haplotypes were also found within individual families, particularly near SSPO at 7p36.1 in one family and a large linked haplotype spanning 4q21.3-28.3 in a different family. The 4q haplotype contains potential causal rare variants in DSPP at 4q22.1 and PTPN13 at 4q21.3. CONCLUSIONS: Regions on 12q, 7p, and 4q are linked to increased cancer risk in highly aggregated lung cancer families, 12q across families and 7p and 4q within a single family. POLR3B, SSPO, DSPP, and PTPN13 are currently the best candidate genes. IMPACT: Functional work on these genes is planned for future studies and if confirmed would lead to potential biomarkers for risk in cancer.

Small posterior fossa in Chiari I malformation affected families is significantly linked to 1q43-44 and 12q23-24.11 using whole exome sequencing.

The posterior fossa of the cranium contains the cerebellum and brainstem. Processes that reduce the volume of the posterior fossa squeeze the cerebellum and brainstem caudally, resulting in Chiari I malformation (CM1). CM1 causes neck pain, balance issues, decreased motor skills and headaches in those affected. We have posterior fossa measurements and whole exome sequence data on individuals from 7 extended families from Russia that have a family history of CM1. We performed parametric linkage analyses using an autosomal dominant inheritance model with a disease allele frequency of 0.01 and a penetrance of 0.8 for carriers and 0.0 for non-carriers. Variant-based two-point linkage analysis and gene-based linkage analysis was performed. Our results found a genome-wide significant signal on chromosome 1q43-44 (max HLOD = 3.3) in the variant-based analysis and 12q23 (max HLOD = 4.2) in the gene-based analysis. In both cases, the signal was driven by a single (different) family that contained a long, linked haplotype across the region in question. Using functional annotation, we were able to identify several rare nonsynonymous variants that were enriched in each family. The best candidate genes were rs765865412:G>A in MYBPC1 for the 12q haplotype and rs61749963:A>G in COX20 for the 1q haplotype. Good candidate variants in the 1q haplotype were also identified in CEP170 and AKT. Further laboratory work is planned to verify the causality of these genes.