Longitudinal Changes in Bone Mineral Density in Adults with Cystic Fibrosis.

BACKGROUND: Improved survival in people with cystic fibrosis (pwCF) presents new complexities of care, including CF-related bone disease, a common complication in older pwCF. The trajectory of bone loss with age in this population remains unclear. The objective of this study was to estimate the average rate of change in bone mineral density (BMD) in adults with CF. METHODS: This retrospective study included adults with CF, aged 25-48 years, followed between January 2000 and December 2021. Subjects with at least one dual-energy X-ray absorptiometry (DXA) scan were included. Scans obtained post-transplantation, after the initiation of bisphosphonates or cystic fibrosis transmembrane conductance regulator (CFTR) modulator therapy were excluded. The primary outcome was BMD (g/cm2) at the lumbar spine (LS) and femoral neck (FN). A linear mixed-effects model with both random intercept and random slope terms was used to estimate the average annual change in BMD. RESULTS: A total of 1502 DXA scans in 500 adults (average age 28.4y) were included. There was a statistically significant annual decline in BMD of -0.008 gm/cm2/year (95% CI -0.009, -0.007) at the FN and -0.006 gm/cm2/year (95% CI -0.007, -0.004) at the LS. Relative to BMD at age 25, there was a -18.8% decline at the FN by age 48 years and a -11% decline at the LS. Pancreatic insufficient (PI) subjects had a faster rate of decline in BMD compared to pancreatic sufficient (PS) subjects. After adjusting for markers of disease severity, the annual rate of decline remained significant. CONCLUSIONS: Individuals with CF experience bone loss at an age when it is not anticipated, thereby entering early adulthood, where further bone loss is inevitable especially with the decrease in estrogen during menopause, with suboptimal BMD. As the CF population ages, it will become very important to consider interventions to maximize bone health.

Improved survival in people with cystic fibrosis (pwCF) presents new complexities of care, including CF-related bone disease, a common complication in older pwCF. The objective of this study was to estimate how bone mineral density (BMD) changes over time in adults with CF. The study included 500 adults with CF, aged 25-48 years, followed between January 2000 and December 2021 who underwent a scan to assess BMD. On average, BMD decreased over time both at the spine and the hip. People who have pancreatic insufficiency had a faster rate of decline in BMD than people with pancreatic sufficiency. Individuals with CF experience bone loss at an age when bone density is expected to be stable. As the CF population ages, interventions to maximize bone health should be emphasized to prevent fractures in the future.

A Genomic Counseling Model for Population-Based Sequencing: A Pre-Post Intervention Study.

PURPOSE: Novel uses of genome sequencing (GS) present an opportunity for return of results to healthy individuals, prompting the need for scalable genetic counseling strategies. We evaluate the effectiveness of a genomic counseling model (GCM) and explore preferences for GS findings in the general population. METHODS: Participants (N = 466) completed GS and our GCM (digital genomics platform and group-based webinar) and indicated results preferences. Surveys were administered before (T0) and after (T1) GCM. Change in knowledge and decisional conflict (DC) were evaluated using paired-sample T and Wilcoxon tests. Factors influencing knowledge and results preferences were evaluated using linear and logistic regression models. RESULTS: Participants were 56% female, 58% white, and 53% ≥40 years of age. Mean knowledge scores increased (Limitations: 3.73 to 5.63; Benefits: 4.34 to 5.48, P < .0001), and DC decreased (-21.9, P < .0001) at T1 versus T0. Eighty-six percent of participants wished to learn all GS findings at T1 vs 78% at T0 (P < .0001). Older age, negative/mixed attitudes toward genetics and greater DC were associated with change in preferences after intervention. CONCLUSION: In a population-based cohort undergoing GS interested in learning GS findings, our GCM increased knowledge and reduced DC, illustrating the GCM's potential effectiveness for GS counseling in the general population.

Genome-Wide Association Study of Susceptibility to Pseudomonas aeruginosa Infection in Cystic Fibrosis.

QUESTION: Pseudomonas aeruginosa (Pa) is a common pathogen that contributes to progressive lung disease in Cystic Fibrosis (CF). Genetic factors other than CF-causing CFTR variations contribute approximately 85% of the variation in chronic Pa infection age in CF according to twin studies, but the susceptibility loci remain unknown. Our objective is to advance understanding of the genetic basis of host susceptibility to Pa infection. MATERIALS AND METHODS: We conducted a genome-wide association study (GWAS) of chronic Pa infection age in 1037 Canadians with CF. We subsequently assessed the genetic correlation between chronic Pa infection age and lung function through polygenic risk score (PRS) analysis and inferred their causal relationship through bi-directional Mendelian Randomization analysis. RESULTS: Two novel genome-wide significant loci with lead SNPs rs62369766 (chr5p12; p-value= 1.98 ×10-8) and rs927553 (chr13q12.12; p-value= 1.91 × 10-8) were associated with chronic Pa infection age. The rs62369766 locus was validated using an independent French cohort (N=501). Furthermore, PRS constructed from CF lung function-associated SNPs was significantly associated with chronic Pa infection age (p-value=0.002). Finally, our analysis presented evidence for a causal effect of lung function on the chronic Pa infection age (Beta=0.782 years, p-value= 4.24 × 10-4). In the reverse direction, we observed a moderate effect (Beta=0.002, p-value=0.012). CONCLUSIONS: We identified two novel loci that are associated with chronic Pa infection age in individuals with CF. Additionally, we provided evidence of common genetic contributors and a potential causal relationship between Pa infection susceptibility and lung function in CF. Therapeutics targeting these genetic factors may delay the onset of chronic infections which accounts for significant remaining morbidity in CF.

Directional integration and pathway enrichment analysis for multi-omics data.

Omics techniques generate comprehensive profiles of biomolecules in cells and tissues. However, a holistic understanding of underlying systems requires joint analyses of multiple data modalities. We present DPM, a data fusion method for integrating omics datasets using directionality and significance estimates of genes, transcripts, or proteins. DPM allows users to define how the input datasets are expected to interact directionally given the experimental design or biological relationships between the datasets. DPM prioritises genes and pathways that change consistently across the datasets and penalises those with inconsistent directionality. To demonstrate our approach, we characterise gene and pathway regulation in IDH-mutant gliomas by jointly analysing transcriptomic, proteomic, and DNA methylation datasets. Directional integration of survival information in ovarian cancer reveals candidate biomarkers with consistent prognostic signals in transcript and protein expression. DPM is a general and adaptable framework for gene prioritisation and pathway analysis in multi-omics datasets.

Annotation of nuclear lncRNAs based on chromatin interactions.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Canadian COVID-19 host genetics cohort replicates known severity associations.

The HostSeq initiative recruited 10,059 Canadians infected with SARS-CoV-2 between March 2020 and March 2023, obtained clinical information on their disease experience and whole genome sequenced (WGS) their DNA. We analyzed the WGS data for genetic contributors to severe COVID-19 (considering 3,499 hospitalized cases and 4,975 non-hospitalized after quality control). We investigated the evidence for replication of loci reported by the International Host Genetics Initiative (HGI); analyzed the X chromosome; conducted rare variant gene-based analysis and polygenic risk score testing. Population stratification was adjusted for using meta-analysis across ancestry groups. We replicated two loci identified by the HGI for COVID-19 severity: the LZTFL1/SLC6A20 locus on chromosome 3 and the FOXP4 locus on chromosome 6 (the latter with a variant significant at P < 5E-8). We found novel significant associations with MRAS and WDR89 in gene-based analyses, and constructed a polygenic risk score that explained 1.01% of the variance in severe COVID-19. This study provides independent evidence confirming the robustness of previously identified COVID-19 severity loci by the HGI and identifies novel genes for further investigation.

SLCO5A1 and synaptic assembly genes contribute to impulsivity in juvenile myoclonic epilepsy.

Elevated impulsivity is a key component of attention-deficit hyperactivity disorder (ADHD), bipolar disorder and juvenile myoclonic epilepsy (JME). We performed a genome-wide association, colocalization, polygenic risk score, and pathway analysis of impulsivity in JME (n = 381). Results were followed up with functional characterisation using a drosophila model. We identified genome-wide associated SNPs at 8q13.3 (P = 7.5 × 10-9) and 10p11.21 (P = 3.6 × 10-8). The 8q13.3 locus colocalizes with SLCO5A1 expression quantitative trait loci in cerebral cortex (P = 9.5 × 10-3). SLCO5A1 codes for an organic anion transporter and upregulates synapse assembly/organisation genes. Pathway analysis demonstrates 12.7-fold enrichment for presynaptic membrane assembly genes (P = 0.0005) and 14.3-fold enrichment for presynaptic organisation genes (P = 0.0005) including NLGN1 and PTPRD. RNAi knockdown of Oatp30B, the Drosophila polypeptide with the highest homology to SLCO5A1, causes over-reactive startling behaviour (P = 8.7 × 10-3) and increased seizure-like events (P = 6.8 × 10-7). Polygenic risk score for ADHD genetically correlates with impulsivity scores in JME (P = 1.60 × 10-3). SLCO5A1 loss-of-function represents an impulsivity and seizure mechanism. Synaptic assembly genes may inform the aetiology of impulsivity in health and disease.

Characterizing Risk Factors for Hospitalization and Clinical Characteristics in a Cohort of COVID-19 Patients Enrolled in the GENCOV Study.

The GENCOV study aims to identify patient factors which affect COVID-19 severity and outcomes. Here, we aimed to evaluate patient characteristics, acute symptoms and their persistence, and associations with hospitalization. Participants were recruited at hospital sites across the Greater Toronto Area in Ontario, Canada. Patient-reported demographics, medical history, and COVID-19 symptoms and complications were collected through an intake survey. Regression analyses were performed to identify associations with outcomes including hospitalization and COVID-19 symptoms. In total, 966 responses were obtained from 1106 eligible participants (87% response rate) between November 2020 and May 2022. Increasing continuous age (aOR: 1.05 [95%CI: 1.01-1.08]) and BMI (aOR: 1.17 [95%CI: 1.10-1.24]), non-White/European ethnicity (aOR: 2.72 [95%CI: 1.22-6.05]), hypertension (aOR: 2.78 [95%CI: 1.22-6.34]), and infection by viral variants (aOR: 5.43 [95%CI: 1.45-20.34]) were identified as risk factors for hospitalization. Several symptoms including shortness of breath and fever were found to be more common among inpatients and tended to persist for longer durations following acute illness. Sex, age, ethnicity, BMI, vaccination status, viral strain, and underlying health conditions were associated with developing and having persistent symptoms. By improving our understanding of risk factors for severe COVID-19, our findings may guide COVID-19 patient management strategies by enabling more efficient clinical decision making.

Variation in prognosis and treatment outcome in juvenile myoclonic epilepsy: a Biology of Juvenile Myoclonic Epilepsy Consortium proposal for a practical definition and stratified medicine classifications.

Reliable definitions, classifications and prognostic models are the cornerstones of stratified medicine, but none of the current classifications systems in epilepsy address prognostic or outcome issues. Although heterogeneity is widely acknowledged within epilepsy syndromes, the significance of variation in electroclinical features, comorbidities and treatment response, as they relate to diagnostic and prognostic purposes, has not been explored. In this paper, we aim to provide an evidence-based definition of juvenile myoclonic epilepsy showing that with a predefined and limited set of mandatory features, variation in juvenile myoclonic epilepsy phenotype can be exploited for prognostic purposes. Our study is based on clinical data collected by the Biology of Juvenile Myoclonic Epilepsy Consortium augmented by literature data. We review prognosis research on mortality and seizure remission, predictors of antiseizure medication resistance and selected adverse drug events to valproate, levetiracetam and lamotrigine. Based on our analysis, a simplified set of diagnostic criteria for juvenile myoclonic epilepsy includes the following: (i) myoclonic jerks as mandatory seizure type; (ii) a circadian timing for myoclonia not mandatory for the diagnosis of juvenile myoclonic epilepsy; (iii) age of onset ranging from 6 to 40 years; (iv) generalized EEG abnormalities; and (v) intelligence conforming to population distribution. We find sufficient evidence to propose a predictive model of antiseizure medication resistance that emphasises (i) absence seizures as the strongest stratifying factor with regard to antiseizure medication resistance or seizure freedom for both sexes and (ii) sex as a major stratifying factor, revealing elevated odds of antiseizure medication resistance that correlates to self-report of catamenial and stress-related factors including sleep deprivation. In women, there are reduced odds of antiseizure medication resistance associated with EEG-measured or self-reported photosensitivity. In conclusion, by applying a simplified set of criteria to define phenotypic variations of juvenile myoclonic epilepsy, our paper proposes an evidence-based definition and prognostic stratification of juvenile myoclonic epilepsy. Further studies in existing data sets of individual patient data would be helpful to replicate our findings, and prospective studies in inception cohorts will contribute to validate them in real-world practice for juvenile myoclonic epilepsy management.

Gene copy number variation and pediatric mental health/neurodevelopment in a general population.

We assessed the relationship of gene copy number variation (CNV) in mental health/neurodevelopmental traits and diagnoses, physical health and cognition in a community sample of 7100 unrelated children and youth of European or East Asian ancestry (Spit for Science). Clinically significant or susceptibility CNVs were present in 3.9% of participants and were associated with elevated scores on a continuous measure of attention-deficit/hyperactivity disorder (ADHD) traits (P = 5.0 × 10-3), longer response inhibition (a cognitive deficit found in several mental health and neurodevelopmental disorders; P = 1.0 × 10-2) and increased prevalence of mental health diagnoses (P = 1.9 × 10-6, odds ratio: 3.09), specifically ADHD, autism spectrum disorder anxiety and learning problems/learning disorder (P's < 0.01). There was an increased burden of rare deletions in gene-sets related to brain function or expression in brain associated with more ADHD traits. With the current mental health crisis, our data established a baseline for delineating genetic contributors in pediatric-onset conditions.

Cost-effectiveness analysis of genetic tools to predict treatment response in patients with cystic fibrosis.

BACKGROUND: Cystic fibrosis (CF) transmembrane conductance regulator (CFTR) modulator therapies show variable efficacy for patients with CF. Patient-derived predictive tools may identify individuals likely to respond to CFTRs, but are not in routine use. We aimed to determine the cost-utility of predictive tool-guided treatment with CFTRs as add-on to standard of care (SoC) for individuals with CF. METHODS: This economic evaluation compared two strategies using an individual level simulation: (i) Treat All, where all patients received CFTRs plus SoC and (ii) Test→Treat, where patients who tested positive on predictive tools received CFTRs plus SoC and those who tested negative received SoC only. We simulated 50,000 individuals over their lifetime, and estimated costs (2020 CAD) per quality-adjusted life year (QALY) from the healthcare payer's perspective, discounted at 1.5% annually. The model was populated using Canadian CF registry data and published literature. Probabilistic and deterministic sensitivity were conducted. RESULTS: The Treat All and Test→Treat and strategies yielded 22.41 and 21.36 QALYs, and cost $4.21 M and $3.15 M respectively. Results of probabilistic sensitivity analysis showed that Test→Treat was highly cost-effective compared to Treat All in 100% of simulations at cost-effectiveness thresholds as high as $500,000 per QALY. Test→Treat may save between $931 K to $1.1 M per QALY lost, depending on sensitivity and specificity of predictive tools. CONCLUSION: The use of predictive tools could optimize the health benefits of CFTR modulators while reducing costs. Our findings support the use of pre-treatment predictive testing and may help inform coverage and reimbursement policies for individuals with CF.

Validating organoid-derived human intestinal monolayers for personalized therapy in cystic fibrosis.

Highly effective drugs modulating the defective protein encoded by the CFTR gene have revolutionized cystic fibrosis (CF) therapy. Preclinical drug-testing on human nasal epithelial (HNE) cell cultures and 3-dimensional human intestinal organoids (3D HIO) are used to address patient-specific variation in drug response and to optimize individual treatment for people with CF. This study is the first to report comparable CFTR functional responses to CFTR modulator treatment among patients with different classes of CFTR gene variants using the three methods of 2D HIO, 3D HIO, and HNE. Furthermore, 2D HIO showed good correlation to clinical outcome markers. A larger measurable CFTR functional range and access to the apical membrane were identified as advantages of 2D HIO over HNE and 3D HIO, respectively. Our study thus expands the utility of 2D intestinal monolayers as a preclinical drug testing tool for CF.

RoPE: A robust profile likelihood method for differential gene expression analysis.

Variation in RNA-Seq data creates modeling challenges for differential gene expression (DE) analysis. Statistical approaches address conventional small sample sizes and implement empirical Bayes or non-parametric tests, but frequently produce different conclusions. Increasing sample sizes enable proposal of alternative DE paradigms. Here we develop RoPE, which uses a data-driven adjustment for variation and a robust profile likelihood ratio DE test. Simulation studies show RoPE can have improved performance over existing tools as sample size increases and has the most reliable control of error rates. Application of RoPE demonstrates that an active Pseudomonas aeruginosa infection downregulates the SLC9A3 Cystic Fibrosis modifier gene.

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.

Identification of brain cell types underlying genetic association with word reading and correlated traits.

Neuroimaging studies implicate multiple cortical regions in reading ability/disability. However, the neural cell types integral to the reading process are unknown. To contribute to this gap in knowledge, we integrated genetic results from genome-wide association studies for word reading (n = 5054) with gene expression datasets from adult/fetal human brain. Linkage disequilibrium score regression (LDSC) suggested that variants associated with word reading were enriched in genes expressed in adult excitatory neurons, specifically layer 5 and 6 FEZF2 expressing neurons and intratelencephalic (IT) neurons, which express the marker genes LINC00507, THEMIS, or RORB. Inhibitory neurons (VIP, SST, and PVALB) were also found. This finding was interesting as neurometabolite studies previously implicated excitatory-inhibitory imbalances in the etiology of reading disabilities (RD). We also tested traits that shared genetic etiology with word reading (previously determined by polygenic risk scores): attention-deficit/hyperactivity disorder (ADHD), educational attainment, and cognitive ability. For ADHD, we identified enrichment in L4 IT adult excitatory neurons. For educational attainment and cognitive ability, we confirmed previous studies identifying multiple subclasses of adult cortical excitatory and inhibitory neurons, as well as astrocytes and oligodendrocytes. For educational attainment and cognitive ability, we also identified enrichment in multiple fetal cortical excitatory and inhibitory neurons, intermediate progenitor cells, and radial glial cells. In summary, this study supports a role of excitatory and inhibitory neurons in reading and excitatory neurons in ADHD and contributes new information on fetal cell types enriched in educational attainment and cognitive ability, thereby improving our understanding of the neurobiological basis of reading/correlated traits.

Large scale genome-wide association analyses identify novel genetic loci and mechanisms in hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is an important cause of morbidity and mortality with both monogenic and polygenic components. We here report results from the largest HCM genome-wide association study (GWAS) and multi-trait analysis (MTAG) including 5,900 HCM cases, 68,359 controls, and 36,083 UK Biobank (UKB) participants with cardiac magnetic resonance (CMR) imaging. We identified a total of 70 loci (50 novel) associated with HCM, and 62 loci (32 novel) associated with relevant left ventricular (LV) structural or functional traits. Amongst the common variant HCM loci, we identify a novel HCM disease gene, SVIL, which encodes the actin-binding protein supervillin, showing that rare truncating SVIL variants cause HCM. Mendelian randomization analyses support a causal role of increased LV contractility in both obstructive and non-obstructive forms of HCM, suggesting common disease mechanisms and anticipating shared response to therapy. Taken together, the findings significantly increase our understanding of the genetic basis and molecular mechanisms of HCM, with potential implications for disease management.

Language and reading impairments are associated with increased prevalence of non-right-handedness.

Handedness has been studied for association with language-related disorders because of its link with language hemispheric dominance. No clear pattern has emerged, possibly because of small samples, publication bias, and heterogeneous criteria across studies. Non-right-handedness (NRH) frequency was assessed in N = 2503 cases with reading and/or language impairment and N = 4316 sex-matched controls identified from 10 distinct cohorts (age range 6-19 years old; European ethnicity) using a priori set criteria. A meta-analysis (Ncases  = 1994) showed elevated NRH % in individuals with language/reading impairment compared with controls (OR = 1.21, CI = 1.06-1.39, p = .01). The association between reading/language impairments and NRH could result from shared pathways underlying brain lateralization, handedness, and cognitive functions.

High-quality read-based phasing of cystic fibrosis cohort informs genetic understanding of disease modification.

Phasing of heterozygous alleles is critical for interpretation of cis-effects of disease-relevant variation. We sequenced 477 individuals with cystic fibrosis (CF) using linked-read sequencing, which display an average phase block N50 of 4.39 Mb. We use these samples to construct a graph representation of CFTR haplotypes, demonstrating its utility for understanding complex CF alleles. These are visualized in a Web app, CFTbaRcodes, that enables interactive exploration of CFTR haplotypes present in this cohort. We perform fine-mapping and phasing of the chr7q35 trypsinogen locus associated with CF meconium ileus, an intestinal obstruction at birth associated with more severe CF outcomes and pancreatic disease. A 20-kb deletion polymorphism and a PRSS2 missense variant p.Thr8Ile (rs62473563) are shown to independently contribute to meconium ileus risk (p = 0.0028, p = 0.011, respectively) and are PRSS2 pancreas eQTLs (p = 9.5 × 10-7 and p = 1.4 × 10-4, respectively), suggesting the mechanism by which these polymorphisms contribute to CF. The phase information from linked reads provides a putative causal explanation for variation at a CF-relevant locus, which also has implications for the genetic basis of non-CF pancreatitis, to which this locus has been reported to contribute.

Genome screening, reporting, and genetic counseling for healthy populations.

Rapid advancements of genome sequencing (GS) technologies have enhanced our understanding of the relationship between genes and human disease. To incorporate genomic information into the practice of medicine, new processes for the analysis, reporting, and communication of GS data are needed. Blood samples were collected from adults with a PCR-confirmed SARS-CoV-2 (COVID-19) diagnosis (target N = 1500). GS was performed. Data were filtered and analyzed using custom pipelines and gene panels. We developed unique patient-facing materials, including an online intake survey, group counseling presentation, and consultation letters in addition to a comprehensive GS report. The final report includes results generated from GS data: (1) monogenic disease risks; (2) carrier status; (3) pharmacogenomic variants; (4) polygenic risk scores for common conditions; (5) HLA genotype; (6) genetic ancestry; (7) blood group; and, (8) COVID-19 viral lineage. Participants complete pre-test genetic counseling and confirm preferences for secondary findings before receiving results. Counseling and referrals are initiated for clinically significant findings. We developed a genetic counseling, reporting, and return of results framework that integrates GS information across multiple areas of human health, presenting possibilities for the clinical application of comprehensive GS data in healthy individuals.

Hypothesis-driven genome-wide association studies provide novel insights into genetics of reading disabilities.

Reading Disability (RD) is often characterized by difficulties in the phonology of the language. While the molecular mechanisms underlying it are largely undetermined, loci are being revealed by genome-wide association studies (GWAS). In a previous GWAS for word reading (Price, 2020), we observed that top single-nucleotide polymorphisms (SNPs) were located near to or in genes involved in neuronal migration/axon guidance (NM/AG) or loci implicated in autism spectrum disorder (ASD). A prominent theory of RD etiology posits that it involves disturbed neuronal migration, while potential links between RD-ASD have not been extensively investigated. To improve power to identify associated loci, we up-weighted variants involved in NM/AG or ASD, separately, and performed a new Hypothesis-Driven (HD)-GWAS. The approach was applied to a Toronto RD sample and a meta-analysis of the GenLang Consortium. For the Toronto sample (n = 624), no SNPs reached significance; however, by gene-set analysis, the joint contribution of ASD-related genes passed the threshold (p~1.45 × 10-2, threshold = 2.5 × 10-2). For the GenLang Cohort (n = 26,558), SNPs in DOCK7 and CDH4 showed significant association for the NM/AG hypothesis (sFDR q = 1.02 × 10-2). To make the GenLang dataset more similar to Toronto, we repeated the analysis restricting to samples selected for reading/language deficits (n = 4152). In this GenLang selected subset, we found significant association for a locus intergenic between BTG3-C21orf91 for both hypotheses (sFDR q < 9.00 × 10-4). This study contributes candidate loci to the genetics of word reading. Data also suggest that, although different variants may be involved, alleles implicated in ASD risk may be found in the same genes as those implicated in word reading. This finding is limited to the Toronto sample suggesting that ascertainment influences genetic associations.