Machine learning-based longitudinal prediction for GJB2-related sensorineural hearing loss.

BACKGROUND: Recessive GJB2 variants, the most common genetic cause of hearing loss, may contribute to progressive sensorineural hearing loss (SNHL). The aim of this study is to build a realistic predictive model for GJB2-related SNHL using machine learning to enable personalized medical planning for timely intervention. METHOD: Patients with SNHL with confirmed biallelic GJB2 variants in a nationwide cohort between 2005 and 2022 were included. Different data preprocessing protocols and computational algorithms were combined to construct a prediction model. We randomly divided the dataset into training, validation, and test sets at a ratio of 72:8:20, and repeated this process ten times to obtain an average result. The performance of the models was evaluated using the mean absolute error (MAE), which refers to the discrepancy between the predicted and actual hearing thresholds. RESULTS: We enrolled 449 patients with 2184 audiograms available for deep learning analysis. SNHL progression was identified in all models and was independent of age, sex, and genotype. The average hearing progression rate was 0.61 dB HL per year. The best MAE for linear regression, multilayer perceptron, long short-term memory, and attention model were 4.42, 4.38, 4.34, and 4.76 dB HL, respectively. The long short-term memory model performed best with an average MAE of 4.34 dB HL and acceptable accuracy for up to 4 years. CONCLUSIONS: We have developed a prognostic model that uses machine learning to approximate realistic hearing progression in GJB2-related SNHL, allowing for the design of individualized medical plans, such as recommending the optimal follow-up interval for this population.

Genetics in neovascular age-related macular degeneration susceptibility and treatment response to anti-VEGF intravitreal injection: A case series study.

BACKGROUND: To identify genotypes associated with neovascular age-related macular degeneration (nAMD) and investigate the associations between genotype variations and anti-vascular endothelial growth factor (VEGF) treatment response. METHODS: This observational, retrospective, case series study enrolled patients diagnosed with nAMD who received anti-VEGF treatment in National Taiwan University Hospital with at least one-year follow-up between 2012 and 2020. A genome-wide association study (GWAS) was conducted on enrolled patients and controls. Correlations between the genotypes identified from GWAS and the treatment response of functional/anatomical biomarkers, including visual acuity (VA), presence of intraretinal or subretinal fluid (SRF), serous or fibrovascular pigmented epithelium detachment (PED), and disruption of the ellipsoid zone (EZ), were analysed. RESULTS: In total, 182 patients with nAMD and 1748 controls were enrolled. GWAS revealed 16 single nucleotide polymorphisms (SNPs) as risk loci for nAMD, including seven loci in CFH and ARMS2/HTRA1 and nine novel loci, including rs117517872 and rs79835234(COPB2-DT), rs7525578(RAP1A), rs2123738(LOC105376755), rs1374879(CNTN3), rs3812692(SAR1A), rs117501587(PRKCA), rs9965945(CNDP1), and rs189769231(MATK). Our study revealed rs800292(CFH), rs11200638(HTRA1), and rs2123738(LOC105376755) correlated with poor treatment response in VA (P = 0.005), SRF (P = 0.044), and fibrovascular PED (P = 0.007), respectively. Rs9965945(CNDP1) was correlated with poor response in disruption of EZ (P = 0.046) and serous PED (P = 0.049). CONCLUSIONS: Among the 16 SNPs found in the GWAS, four loci-CFH, ARMS2/HTRA1, and two novel loci-were correlated with the susceptibility of nAMD and anatomical/functional responses after anti-VEGF treatment.

A novel framework for human leukocyte antigen (HLA) genotyping using probe capture-based targeted next-generation sequencing and computational analysis.

Human leukocyte antigen (HLA) genes play pivotal roles in numerous immunological applications. Given the immense number of polymorphisms, achieving accurate high-throughput HLA typing remains challenging. This study aimed to harness the human pan-genome reference consortium (HPRC) resources as a potential benchmark for HLA reference materials. We meticulously annotated specific four field-resolution alleles for 11 HLA genes (HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, -DRB3, -DRB4 and -DRB5) from 44 high-quality HPRC personal genome assemblies. For sequencing, we crafted HLA-specific probes and conducted capture-based targeted sequencing of the genomic DNA of the HPRC cohort, ensuring focused and comprehensive coverage of the HLA region of interest. We used publicly available short-read whole-genome sequencing (WGS) data from identical samples to offer a comparative perspective. To decipher the vast amount of sequencing data, we employed seven distinct software tools: OptiType, HLA-VBseq, HISAT genotype, SpecHLA, T1K, QzType, and DRAGEN. Each tool offers unique capabilities and algorithms for HLA genotyping, allowing comprehensive analysis and validation of the results. We then compared these results with benchmarks derived from personal genome assemblies. Our findings present a comprehensive four-field-resolution HLA allele annotation for 44 HPRC samples. Significantly, our innovative targeted next-generation sequencing (NGS) approach for HLA genes showed superior accuracy compared with conventional short-read WGS. An integrated analysis involving QzType, T1K, and DRAGEN was developed, achieving 100% accuracy for all 11 HLA genes. In conclusion, our study highlighted the combination of targeted short-read sequencing and astute computational analysis as a robust approach for HLA genotyping. Furthermore, the HPRC cohort has emerged as a valuable assembly-based reference in this realm.

Complete genomic profiles of 1496 Taiwanese reveal curated medical insights.

INTRODUCTION: The population of Taiwan has a long history of ethno-cultural evolution. The Taiwanese population was isolated from other large populations such as the European, Han Chinese, and Japanese population. The Taiwan Biobank (TWB) project has built a nationwide database, particularly for personal whole-genome sequence (WGS) to facilitate basic and clinical collaboration nationally and internationally, making it one of the most valuable public datasets of the East Asian population. OBJECTIVES: This study provides comprehensive medical genomic findings from TWB WGS data, for better characterization of disease susceptibility and the choice of ideal treatment regimens in Taiwanese population. METHODS: We reanalyzed 1496 WGS using a PrecisionFDA Truth challenge winner method Sentieon DNAscope. Single nucleotide variants (SNV) and small insertions/deletions (INDEL) were benchmarked. We also analyzed pharmacogenomic (PGx) drug-associated alleles, and copy number variants (CNV). Multiple practicing clinicians reviewed and curated the clinically significant variants. Variant annotations can be browsed at TaiwanGenomes (https://genomes.tw). RESULTS: We found that each participant had an average of 6,870.7 globally novel variants and 75.3% (831/1103) of the participants harbored at least one PharmGKB-selected high evidence level human leukocyte antigen (HLA) risk allele. 54 PharmGKB-reported high-level instances of evidence of Cytochrome P450 variant-drug pairs, with a population frequency of over 13.2%. We also identified 23 variants in the ACMG secondary finding V3 gene list from 25 participants, suggesting that 1.67% (25/1496) of the population is harboring at least one medical actionable variant. Our carrier status analyses suggest that one in 25 couples (3.94%) would risk having offspring with at least one pathogenic variant, which is in line with rates found in Japan and Singapore. For pathogenic CNV, we detected 6.88% and 2.02% carrier rates for alpha thalassemia and spinal muscular atrophy, respectively. CONCLUSION: Our study highlights the overall medical insights of a complete Taiwanese genomic profile.

Genetic Underpinnings and Audiological Characteristics in Children With Unilateral Sensorineural Hearing Loss.

OBJECTIVE: Unilateral sensorineural hearing loss (USNHL) is a condition commonly encountered in otolaryngology clinics. However, its molecular pathogenesis remains unclear. This study aimed to investigate the genetic underpinnings of childhood USNHL and analyze the associated audiological features. STUDY DESIGN: Retrospective analysis of a prospectively recruited cohort. SETTING: Tertiary referral center. METHODS: We enrolled 38 children with USNHL between January 1, 2018, and December 31, 2021, and performed physical, audiological, imaging, and congenital cytomegalovirus (cCMV) examinations as well as genetic testing using next-generation sequencing (NGS) targeting 30 deafness genes. The audiological results were compared across different etiologies. RESULTS: Causative genetic variants were identified in 8 (21.1%) patients, including 5 with GJB2 variants, 2 with PAX3 variants, and 1 with the EDNRB variant. GJB2 variants were found to be associated with mild-to-moderate USNHL in various audiogram configurations, whereas PAX3 and EDNRB variants were associated with profound USNHL in flat audiogram configurations. In addition, whole-genome sequencing and extended NGS targeting 213 deafness genes were performed in 2 multiplex families compatible with autosomal recessive inheritance; yet no definite causative variants were identified. Cochlear nerve deficiency and cCMV infection were observed in 9 and 2, respectively, patients without definite genetic diagnoses. CONCLUSION: Genetic underpinnings can contribute to approximately 20% of childhood USNHL, and different genotypes are associated with various audiological features. These findings highlight the utility of genetic examinations in guiding the diagnosis, counseling, and treatment of USNHL in children.

Combining Panel-Based Next-Generation Sequencing and Exome Sequencing for Genetic Liver Diseases.

OBJECTIVES: To determine how advanced genetic analysis methods may help in clinical diagnosis. STUDY DESIGN: We report a combined genetic diagnosis approach for patients with clinical suspicion of genetic liver diseases in a tertiary referral center, using tools either tier 1: Sanger sequencing on SLC2SA13, ATP8B1, ABCB11, ABCB4, and JAG1 genes, tier 2: panel-based next generation sequencing (NGS), or tier 3: whole-exome sequencing (WES) analysis. RESULTS: In a total of 374 patients undergoing genetic analysis, 175 patients received tier 1 Sanger sequencing based on phenotypic suspicion, and pathogenic variants were identified in 38 patients (21.7%). Tier 2 included 216 patients (39 of tier 1-negative patients) who received panel-based NGS, and pathogenic variants were identified in 60 (27.8%). In tier 3, 41 patients received WES analysis, and 20 (48.8%) obtained genetic diagnosis. Pathogenic variants were detected in 6 of 19 (31.6%) who tested negative in tier 2, and a greater detection rate in 14 of 22 (63.6%) patients with deteriorating/multiorgan disease receiving one-step WES (P = .041). The overall disease spectrum is comprised of 35 genetic defects; 90% of genes belong to the functional categories of small molecule metabolism, ciliopathy, bile duct development, and membrane transport. Only 13 (37%) genetic diseases were detected in more than 2 families. A hypothetical approach using a small panel-based NGS can serve as the first tier with diagnostic yield of 27.8% (98/352). CONCLUSIONS: NGS based genetic test using a combined panel-WES approach is efficient for the diagnosis of the highly diverse genetic liver diseases.

Revisiting Genetic Epidemiology with a Refined Targeted Gene Panel for Hereditary Hearing Impairment in the Taiwanese Population.

Hearing impairment is one of the most common sensory disorders in children, and targeted next-generation sequencing (NGS)-based genetic examinations can assist in its prognostication and management. In 2020, we developed a simplified 30-gene NGS panel from the original 214-gene NGS version based on Taiwanese genetic epidemiology data to increase the accessibility of NGS-based examinations. In this study, we evaluated the diagnostic performance of the 30-gene NGS panel and compared it with that of the original 214-gene NGS panel in patient subgroups with different clinical features. Data on the clinical features, genetic etiologies, audiological profiles, and outcomes were collected from 350 patients who underwent NGS-based genetic examinations for idiopathic bilateral sensorineural hearing impairment between 2020 and 2022. The overall diagnostic yield was 52%, with slight differences in genetic etiology between patients with different degrees of hearing impairment and ages of onset. No significant difference was found in the diagnostic yields between the two panels, regardless of clinical features, except for a lower detection rate of the 30-gene panel in the late-onset group. For patients with negative genetic results, where the causative variant is undetectable on current NGS-based methods, part of the negative results may be due to genes not covered by the panel or yet to be identified. In such cases, the hearing prognosis varies and may decline over time, necessitating appropriate follow-up and consultation. In conclusion, genetic etiologies can serve as references for refining targeted NGS panels with satisfactory diagnostic performance.

Identification of nine novel variants across PAX3, SOX10, EDNRB, and MITF genes in Waardenburg syndrome with next-generation sequencing.

BACKGROUND: Waardenburg syndrome (WS) is a hereditary, genetically heterogeneous disorder characterized by variable presentations of sensorineural hearing impairment and pigmentation anomalies. This study aimed to investigate the clinical features of WS in detail and determine the genetic causes of patients with clinically suspected WS. METHODS: A total of 24 patients from 21 Han-Taiwanese families were enrolled and underwent comprehensive physical and audiological examinations. We applied targeted next-generation sequencing (NGS) to investigate the potential causative variants in these patients and further validated the candidate variants through Sanger sequencing. RESULTS: We identified 19 causative variants of WS in our cohort. Of these variants, nine were novel and discovered in PAX3, SOX10, EDNRB, and MITF genes, including missense, nonsense, deletion, and splice site variants. Several patients presented with skeletal deformities, hypotonia, megacolon, and neurological disorders that were rarely seen in WS. CONCLUSION: This study revealed highly phenotypic variability in Taiwanese WS patients and demonstrated that targeted NGS allowed us to clarify the genetic diagnosis and extend the genetic variant spectrum of WS.

Profiling genes encoding the adaptive immune receptor repertoire with gAIRR Suite.

Adaptive immune receptor repertoire (AIRR) is encoded by T cell receptor (TR) and immunoglobulin (IG) genes. Profiling these germline genes encoding AIRR (abbreviated as gAIRR) is important in understanding adaptive immune responses but is challenging due to the high genetic complexity. Our gAIRR Suite comprises three modules. gAIRR-seq, a probe capture-based targeted sequencing pipeline, profiles gAIRR from individual DNA samples. gAIRR-call and gAIRR-annotate call alleles from gAIRR-seq reads and annotate whole-genome assemblies, respectively. We gAIRR-seqed TRV and TRJ of seven Genome in a Bottle (GIAB) DNA samples with 100% accuracy and discovered novel alleles. We also gAIRR-seqed and gAIRR-called the TR and IG genes of a subject from both the peripheral blood mononuclear cells (PBMC) and oral mucosal cells. The calling results from these two cell types have a high concordance (99% for all known gAIRR alleles). We gAIRR-annotated 36 genomes to unearth 325 novel TRV alleles and 29 novel TRJ alleles. We could further profile the flanking sequences, including the recombination signal sequence (RSS). We validated two structural variants for HG002 and uncovered substantial differences of gAIRR genes in references GRCh37 and GRCh38. gAIRR Suite serves as a resource to sequence, analyze, and validate germline TR and IG genes to study various immune-related phenotypes.

Comprehensive Etiologic Analyses in Pediatric Cochlear Implantees and the Clinical Implications.

Cochlear implantation is the treatment of choice for children with profound sensorineural hearing impairment (SNHI), yet the outcomes of cochlear implants (CI) vary significantly across individuals. To investigate the CI outcomes in pediatric patients with SNHI due to various etiologies, we prospectively recruited children who underwent CI surgery at two tertiary referral CI centers from 2010 to 2021. All patients underwent comprehensive history taking, next generation sequencing (NGS)-based genetic examinations, and imaging studies. The CI outcomes were evaluated using Categories of Auditory Performance (CAP) and Speech Intelligibility Rating (SIR) scores. Of the 160 pediatric cochlear implantees (76 females and 84 males) included in this study, comprehensive etiological work-up helped achieve clinical diagnoses in 83.1% (133/160) of the patients, with genetic factors being the leading cause (61.3%). Imaging studies identified certain findings in 31 additional patients (19.3%). Four patients (2.5%) were identified with congenital cytomegalovirus infection (cCMV), and 27 patients (16.9%) remained with unknown etiologies. Pathogenic variants in the four predominant non-syndromic SNHI genes (i.e., SLC26A4, GJB2, MYO15A, and OTOF) were associated with favorable CI outcomes (Chi-square test, p = 0.023), whereas cochlear nerve deficiency (CND) on imaging studies was associated with unfavorable CI outcomes (Chi-square test, p < 0.001). Our results demonstrated a clear correlation between the etiologies and CI outcomes, underscoring the importance of thorough etiological work-up preoperatively in pediatric CI candidates.

Cochlear Implantation Outcomes in Patients with Auditory Neuropathy Spectrum Disorder of Genetic and Non-Genetic Etiologies: A Multicenter Study.

With diverse etiologies and clinical features, the management of pediatric auditory neuropathy spectrum disorder (ANSD) is often challenging, and the outcomes of cochlear implants (CIs) are variable. This study aimed to investigate CI outcomes in pediatric patients with ANSD of different etiologies. Thirty-six children with ANSD who underwent cochlear implantation between 2001 and 2021 were included. Comprehensive etiological analyses were conducted, including a history review, next-generation sequencing-based genetic examinations, and imaging studies using high-resolution computed tomography and magnetic resonance imaging. Serial behavioral and speech audiometry were performed before and after surgery, and the outcomes with CI were evaluated using the Categories of Auditory Performance (CAP) and Speech Intelligibility Rating (SIR) scores. By etiology, 18, 1, 1, and 10 patients had OTOF-related, WFS1-related, OPA1-related, and cochlear nerve deficiency (CND)-related ANSD, respectively. Six patients had no definite etiology. The average CI-aided behavioral threshold was 28.3 ± 7.8 dBHL, and those with CND-related ANSD were significantly worse than OTOF-related ANSD. The patients' median CAP and SIR scores were 6 and 4, respectively. Favorable CI outcomes were observed in patients with certain etiologies of ANSD, particularly those with OTOF (CAP/SIR scores 5-7/2-5), WFS1 (CAP/SIR score 6/5), and OPA1 variants (CAP/SIR score 7/5). Patients with CND had suboptimal CI outcomes (CAP/SIR scores 2-6/1-3). Identifying the etiologies in ANSD patients is crucial before surgery and can aid in predicting prognoses.

Identification of a wide spectrum of ciliary gene mutations in nonsyndromic biliary atresia patients implicates ciliary dysfunction as a novel disease mechanism.

BACKGROUND: Biliary atresia (BA) is the most common obstructive cholangiopathy in neonates, often progressing to end-stage cirrhosis. BA pathogenesis is believed to be multifactorial, but the genetic contribution, especially for nonsyndromic BA (common form: > 85%) remains poorly defined. METHODS: We conducted whole exome sequencing on 89 nonsyndromic BA trios to identify rare variants contributing to BA etiology. Functional evaluation using patients' liver biopsies, human cell and zebrafish models were performed. Clinical impact on respiratory system was assessed with clinical evaluation, nasal nitric oxide (nNO), high speed video analysis and transmission electron microscopy. FINDINGS: We detected rare, deleterious de novo or biallelic variants in liver-expressed ciliary genes in 31.5% (28/89) of the BA patients. Burden test revealed 2.6-fold (odds ratio (OR) [95% confidence intervals (CI)]= 2.58 [1.15-6.07], adjusted p = 0.034) over-representation of rare, deleterious mutations in liver-expressed ciliary gene set in patients compared to controls. Functional analyses further demonstrated absence of cilia in the BA livers with KIF3B and TTC17 mutations, and knockdown of PCNT, KIF3B and TTC17 in human control fibroblasts and cholangiocytes resulted in reduced number of cilia. Additionally, CRISPR/Cas9-engineered zebrafish knockouts of KIF3B, PCNT and TTC17 displayed reduced biliary flow. Abnormally low level of nNO was detected in 80% (8/10) of BA patients carrying deleterious ciliary mutations, implicating the intrinsic ciliary defects. INTERPRETATION: Our findings support strong genetic susceptibility for nonsyndromic BA. Ciliary gene mutations leading to cholangiocyte cilia malformation and dysfunction could be a key biological mechanism in BA pathogenesis. FUNDING: The study is supported by General Research Fund, HMRF Commissioned Paediatric Research at HKCH and Li Ka Shing Faculty of Medicine Enhanced New Staff Start-up Fund.

A joint study of whole exome sequencing and structural MRI analysis in major depressive disorder.

BACKGROUND: Major depressive disorder (MDD) is a leading cause of disability worldwide and influenced by both environmental and genetic factors. Genetic studies of MDD have focused on common variants and have been constrained by the heterogeneity of clinical symptoms. METHODS: We sequenced the exome of 77 cases and 245 controls of Han Chinese ancestry and scanned their brain. Burden tests of rare variants were performed first to explore the association between genes/pathways and MDD. Secondly, parallel Independent Component Analysis was conducted to investigate genetic underpinnings of gray matter volume (GMV) changes of MDD. RESULTS: Two genes (CSMD1, p = 5.32×10-6; CNTNAP5, p = 1.32×10-6) and one pathway (Neuroactive Ligand Receptor Interactive, p = 1.29×10-5) achieved significance in burden test. In addition, we identified one pair of imaging-genetic components of significant correlation (r = 0.38, p = 9.92×10-6). The imaging component reflected decreased GMV in cases and correlated with intelligence quotient (IQ). IQ mediated the effects of GMV on MDD. The genetic component enriched in two gene sets, namely Singling by G-protein coupled receptors [false discovery rate (FDR) q = 3.23×10-4) and Alzheimer Disease Up (FDR q = 6.12×10-4). CONCLUSIONS: Both rare variants analysis and imaging-genetic analysis found evidence corresponding with the neuroinflammation and synaptic plasticity hypotheses of MDD. The mediation of IQ indicates that genetic component may act on MDD through GMV alteration and cognitive impairment.

Cancer gene mutations in congenital pulmonary airway malformation patients.

BACKGROUND: Newborns affected with congenital pulmonary airway malformations (CPAMs) may present with severe respiratory distress or remain asymptomatic. While surgical resection is the definitive treatment for symptomatic CPAMs, prophylactic elective surgery may be recommended for asymptomatic CPAMs owing to the risk of tumour development. However, the implementation of prophylactic surgery is quite controversial on the grounds that more evidence linking CPAMs and cancer is needed. The large gap in knowledge of CPAM pathogenesis results in uncertainties and controversies in disease management. As developmental genes control postnatal cell growth and contribute to cancer development, we hypothesised that CPAMs may be underlain by germline mutations in genes governing airways development. METHODS: Sequencing of the exome of 19 patients and their unaffected parents. RESULTS: A more than expected number of mutations in cancer genes (false discovery rate q-value <5.01×10-5) was observed. The co-occurrence, in the same patient, of damaging variants in genes encoding interacting proteins is intriguing, the most striking being thyroglobulin (TG) and its receptor, megalin (LRP2). Both genes are highly relevant in lung development and cancer. CONCLUSIONS: The overall excess of mutations in cancer genes may account for the reported association of CPAMs with carcinomas and provide some evidence to argue for prophylactic surgery by some surgeons.

Uncovering the genetic lesions underlying the most severe form of Hirschsprung disease by whole-genome sequencing.

Hirschsprung disease (HSCR) is a complex birth defect characterized by the lack of ganglion cells along a variable length of the distal intestine. A large proportion of HSCR patients remain genetically unexplained. We applied whole-genome sequencing (WGS) on 9 trios where the probands are sporadically affected with the most severe form of the disorder and harbor no coding sequence variants affecting the function of known HSCR genes. We found de novo protein-altering variants in three intolerant to change genes-CCT2, VASH1, and CYP26A1-for which a plausible link with the enteric nervous system (ENS) exists. De novo single-nucleotide and indel variants were present in introns and non-coding neighboring regions of ENS-related genes, including NRG1 and ERBB4. Joint analysis with those inherited rare variants found under recessive and/or digenic models revealed both patient-unique and shared genetic features where rare variants were found to be enriched in the extracellular matrix-receptor (ECM-receptor) pathway (p = 3.4 × 10-11). Delineation of the genetic profile of each patient might help finding common grounds that could lead to the discovery of shared molecules that could be used as drug targets for the currently ongoing cell therapy effort which aims at providing an alternative to the surgical treatment.

The role of gene variants in the pathogenesis of neurodegenerative disorders as revealed by next generation sequencing studies: a review.

The clinical diagnosis of neurodegenerative disorders based on phenotype is difficult in heterogeneous conditions with overlapping symptoms. It does not take into account the disease etiology or the highly variable clinical course even amongst patients diagnosed with the same disorder. The advent of next generation sequencing (NGS) has allowed for a system-wide, unbiased approach to identify all gene variants in the genome simultaneously. With the plethora of new genes being identified, genetic rather than phenotype-based classification of Mendelian diseases such as spinocerebellar ataxia (SCA), hereditary spastic paraplegia (HSP) and Charcot-Marie-Tooth disease (CMT) has become widely accepted. It has also become clear that gene variants play a role in common and predominantly sporadic neurodegenerative diseases such as Parkinson's disease (PD) and amyotrophic lateral sclerosis (ALS). The observation of pleiotropy has emerged, with mutations in the same gene giving rise to diverse phenotypes, which further increases the complexity of phenotype-genotype correlation. Possible mechanisms of pleiotropy include different downstream effects of different mutations in the same gene, presence of modifier genes, and oligogenic inheritance. Future directions include development of bioinformatics tools and establishment of more extensive public genotype/phenotype databases to better distinguish deleterious gene variants from benign polymorphisms, translation of genetic findings into pathogenic mechanisms through in-vitro and in-vivo studies, and ultimately finding disease-modifying therapies for neurodegenerative disorders.

Burden of rare variants in ALS genes influences survival in familial and sporadic ALS.

Genetic variants are implicated in the development of amyotrophic lateral sclerosis (ALS), but it is unclear whether the burden of rare variants in ALS genes has an effect on survival. We performed whole genome sequencing on 8 familial ALS (FALS) patients with superoxide dismutase 1 (SOD1) mutation and whole exome sequencing on 46 sporadic ALS (SALS) patients living in Hong Kong and found that 67% had at least 1 rare variant in the exons of 40 ALS genes; 22% had 2 or more. Patients with 2 or more rare variants had lower probability of survival than patients with 0 or 1 variant (p = 0.001). After adjusting for other factors, each additional rare variant increased the risk of respiratory failure or death by 60% (p = 0.0098). The presence of the rare variant was associated with the risk of ALS (Odds ratio 1.91, 95% confidence interval 1.03-3.61, p = 0.03), and ALS patients had higher rare variant burden than controls (MB, p = 0.004). Our findings support an oligogenic basis with the burden of rare variants affecting the development and survival of ALS.

Robust and rapid algorithms facilitate large-scale whole genome sequencing downstream analysis in an integrative framework.

Whole genome sequencing (WGS) is a promising strategy to unravel variants or genes responsible for human diseases and traits. However, there is a lack of robust platforms for a comprehensive downstream analysis. In the present study, we first proposed three novel algorithms, sequence gap-filled gene feature annotation, bit-block encoded genotypes and sectional fast access to text lines to address three fundamental problems. The three algorithms then formed the infrastructure of a robust parallel computing framework, KGGSeq, for integrating downstream analysis functions for whole genome sequencing data. KGGSeq has been equipped with a comprehensive set of analysis functions for quality control, filtration, annotation, pathogenic prediction and statistical tests. In the tests with whole genome sequencing data from 1000 Genomes Project, KGGSeq annotated several thousand more reliable non-synonymous variants than other widely used tools (e.g. ANNOVAR and SNPEff). It took only around half an hour on a small server with 10 CPUs to access genotypes of ∼60 million variants of 2504 subjects, while a popular alternative tool required around one day. KGGSeq's bit-block genotype format used 1.5% or less space to flexibly represent phased or unphased genotypes with multiple alleles and achieved a speed of over 1000 times faster to calculate genotypic correlation.

Inheritance-mode specific pathogenicity prioritization (ISPP) for human protein coding genes.

MOTIVATION: Exome sequencing studies have facilitated the detection of causal genetic variants in yet-unsolved Mendelian diseases. However, the identification of disease causal genes among a list of candidates in an exome sequencing study is still not fully settled, and it is often difficult to prioritize candidate genes for follow-up studies. The inheritance mode provides crucial information for understanding Mendelian diseases, but none of the existing gene prioritization tools fully utilize this information. RESULTS: We examined the characteristics of Mendelian disease genes under different inheritance modes. The results suggest that Mendelian disease genes with autosomal dominant (AD) inheritance mode are more haploinsufficiency and de novo mutation sensitive, whereas those autosomal recessive (AR) genes have significantly more non-synonymous variants and regulatory transcript isoforms. In addition, the X-linked (XL) Mendelian disease genes have fewer non-synonymous and synonymous variants. As a result, we derived a new scoring system for prioritizing candidate genes for Mendelian diseases according to the inheritance mode. Our scoring system assigned to each annotated protein-coding gene (N = 18 859) three pathogenic scores according to the inheritance mode (AD, AR and XL). This inheritance mode-specific framework achieved higher accuracy (area under curve = 0.84) in XL mode. CONCLUSION: The inheritance-mode specific pathogenicity prioritization (ISPP) outperformed other well-known methods including Haploinsufficiency, Recessive, Network centrality, Genic Intolerance, Gene Damage Index and Gene Constraint scores. This systematic study suggests that genes manifesting disease inheritance modes tend to have unique characteristics. AVAILABILITY AND IMPLEMENTATION: ISPP is included in KGGSeq v1.0 (http://grass.cgs.hku.hk/limx/kggseq/), and source code is available from (https://github.com/jacobhsu35/ISPP.git). CONTACT: mxli@hku.hkSupplementary information: Supplementary data are available at Bioinformatics online.