An integrative variant analysis suite for whole exome next-generation sequencing data

Genomic Balancing Act: deciphering DNA rearrangements in the complex chromosomal aberration involving 5p15.2, 2q31.1, and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and a terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long-read genome sequencing unveiled the presence of four distinct junctions in the father, who is a carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndromes and to better understand genetic disease mechanisms.

A novel variant in the FLNB gene associated with spondylocarpotarsal synostosis syndrome

OBJECTIVES

Genetic disorders involved in skeleton system arise due to the disturbance in skeletal development, growth and homeostasis. Filamin B is an actin binding protein which is large dimeric protein which cross link actin cytoskeleton filaments into dynamic structure. A single nucleotide changes in the FLNB gene causes spondylocarpotarsal synostosis syndrome, a rare bone disorder due to which the fusion of carpels and tarsals synostosis occurred along with fused vertebrae. In the current study we investigated a family residing in north-western areas of Pakistan.

METHODS

The whole exome sequencing of proband was performed followed by Sanger sequencing of all family members of the subject to validate the variant segregation within the family. Bioinformatics tools were utilized to assess the pathogenicity of the variant.

RESULTS

Whole Exome Sequencing revealed a novel variant (NM_001457: c.209C>T and p.Pro70Leu) in the FLNB gene which was homozygous missense mutation in the FLNB gene. The variant was further validated and visualized by Sanger sequencing and protein structure studies respectively as mentioned before.

CONCLUSIONS

The findings have highlighted the importance of the molecular diagnosis in SCT (spondylocarpotarsal synostosis syndrome) for genetic risk counselling in consanguineous families.

Germline mutations in cancer predisposition genes among pediatric patients with cancer and congenital anomalies

BACKGROUND

Childhood cancer has a poorly known etiology, and investigating the underlying genetic background may provide novel insights. A recognized association exists between non-chromosomal birth defects and childhood cancer susceptibility.

METHODS

We performed whole-exome sequencing and chromosomal microarray analysis in a cohort of childhood cancer (22 individuals, 50% with congenital anomalies) to unravel deleterious germline variants.

RESULTS

A diagnostic yield of 14% was found, encompassing heterozygous variants in bona fide dominant Cancer Predisposition Genes (CPGs). Considering candidate and recessive CPGs harboring monoallelic variants, which were also deemed to play a role in the phenotype, the yield escalated to 45%. Most of the deleterious variants were mapped in genes not conventionally linked to the patient's tumor type. Relevant findings were detected in 55% of the syndromic individuals, mostly variants potentially underlying both phenotypes.

CONCLUSION

We uncovered a remarkable prevalence of germline deleterious CPG variants, highlighting the significance of a comprehensive genetic analysis in pediatric cancer, especially when coupled with additional clinical signs. Moreover, our findings emphasized the potential for oligogenic inheritance, wherein multiple genes synergistically increase cancer risk. Lastly, our investigation unveiled potentially novel genotype-phenotype associations, such as SETD5 in neuroblastoma, KAT6A in gliomas, JAG1 in hepatoblastomas, and TNFRSF13B in Langerhans cell histiocytosis.

IMPACT

Novel gene-phenotype associations and candidate genes for pediatric cancer were unraveled, such as KAT6A in gliomas, SETD5 in neuroblastoma, JAG1 in hepatoblastomas, and TNFRSF13B in Langerhans cell histiocytosis. Our analysis revealed a high frequency of deleterious germline variants, particularly in cases accompanied by additional clinical signs, highlighting the importance of a comprehensive genetic evaluation in childhood cancer. Our findings also underscored the potential for oligogenic inheritance in pediatric cancer risk. Understanding the cancer etiology is crucial for genetic counseling, often influencing therapeutic decisions and offering valuable insights into molecular targets for the development of oncological therapies.

Deep sequencing of candidate genes identified 14 variants associated with smoking abstinence in an ethnically diverse sample

Despite the large public health toll of smoking, genetic studies of smoking cessation have been limited with few discoveries of risk or protective loci. We investigated common and rare variant associations with success in quitting smoking using a cohort from 8 randomized controlled trials involving 2231 participants and a total of 10,020 common and 24,147 rare variants. We identified 14 novel markers including 6 mapping to genes previously related to psychiatric and substance use disorders, 4 of which were protective (CYP2B6 (rs1175607105), HTR3B (rs1413172952; rs1204720503), rs80210037 on chr15), and 2 of which were associated with reduced cessation (PARP15 (rs2173763), SCL18A2 (rs363222)). The others mapped to areas associated with cancer including FOXP1 (rs1288980) and ZEB1 (rs7349). Network analysis identified significant canonical pathways for the serotonin receptor signaling pathway, nicotine and bupropion metabolism, and several related to tumor suppression. Two novel markers (rs6749438; rs6718083) on chr2 are flanked by genes associated with regulation of bodyweight. The identification of novel loci in this study can provide new targets of pharmacotherapy and inform efforts to develop personalized treatments based on genetic profiles.

Genomic Balancing Act: Deciphering DNA rearrangements in the Complex Chromosomal Aberration involving 5p15.2, 2q31.1 and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long read genome sequencing unveiled the presence of four distinct junctions in the father, who is carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndrome and to better understand genetic disease mechanisms.

Human whole-exome genotype data for Alzheimer's disease

The heterogeneity of the whole-exome sequencing (WES) data generation methods present a challenge to a joint analysis. Here we present a bioinformatics strategy for joint-calling 20,504 WES samples collected across nine studies and sequenced using ten capture kits in fourteen sequencing centers in the Alzheimer's Disease Sequencing Project. The joint-genotype called variant-called format (VCF) file contains only positions within the union of capture kits. The VCF was then processed specifically to account for the batch effects arising from the use of different capture kits from different studies. We identified 8.2 million autosomal variants. 96.82% of the variants are high-quality, and are located in 28,579 Ensembl transcripts. 41% of the variants are intronic and 1.8% of the variants are with CADD > 30, indicating they are of high predicted pathogenicity. Here we show our new strategy can generate high-quality data from processing these diversely generated WES samples. The improved ability to combine data sequenced in different batches benefits the whole genomics research community.

Dominant negative variants in KIF5B cause osteogenesis imperfecta via down regulation of mTOR signaling

BACKGROUND

Kinesin motor proteins transport intracellular cargo, including mRNA, proteins, and organelles. Pathogenic variants in kinesin-related genes have been implicated in neurodevelopmental disorders and skeletal dysplasias. We identified de novo, heterozygous variants in KIF5B, encoding a kinesin-1 subunit, in four individuals with osteogenesis imperfecta. The variants cluster within the highly conserved kinesin motor domain and are predicted to interfere with nucleotide binding, although the mechanistic consequences on cell signaling and function are unknown.

METHODS

To understand the in vivo genetic mechanism of KIF5B variants, we modeled the p.Thr87Ile variant that was found in two patients in the C. elegans ortholog, unc-116, at the corresponding position (Thr90Ile) by CRISPR/Cas9 editing and performed functional analysis. Next, we studied the cellular and molecular consequences of the recurrent p.Thr87Ile variant by microscopy, RNA and protein analysis in NIH3T3 cells, primary human fibroblasts and bone biopsy.

RESULTS

C. elegans heterozygous for the unc-116 Thr90Ile variant displayed abnormal body length and motility phenotypes that were suppressed by additional copies of the wild type allele, consistent with a dominant negative mechanism. Time-lapse imaging of GFP-tagged mitochondria showed defective mitochondria transport in unc-116 Thr90Ile neurons providing strong evidence for disrupted kinesin motor function. Microscopy studies in human cells showed dilated endoplasmic reticulum, multiple intracellular vacuoles, and abnormal distribution of the Golgi complex, supporting an intracellular trafficking defect. RNA sequencing, proteomic analysis, and bone immunohistochemistry demonstrated down regulation of the mTOR signaling pathway that was partially rescued with leucine supplementation in patient cells.

CONCLUSION

We report dominant negative variants in the KIF5B kinesin motor domain in individuals with osteogenesis imperfecta. This study expands the spectrum of kinesin-related disorders and identifies dysregulated signaling targets for KIF5B in skeletal development.

PLS3 missense variants affecting the actin-binding domains cause X-linked congenital diaphragmatic hernia and body-wall defects

Congenital diaphragmatic hernia (CDH) is a relatively common and genetically heterogeneous structural birth defect associated with high mortality and morbidity. We describe eight unrelated families with an X-linked condition characterized by diaphragm defects, variable anterior body-wall anomalies, and/or facial dysmorphism. Using linkage analysis and exome or genome sequencing, we found that missense variants in plastin 3 (PLS3), a gene encoding an actin bundling protein, co-segregate with disease in all families. Loss-of-function variants in PLS3 have been previously associated with X-linked osteoporosis (MIM: 300910), so we used in silico protein modeling and a mouse model to address these seemingly disparate clinical phenotypes. The missense variants in individuals with CDH are located within the actin-binding domains of the protein but are not predicted to affect protein structure, whereas the variants in individuals with osteoporosis are predicted to result in loss of function. A mouse knockin model of a variant identified in one of the CDH-affected families, c.1497G>C (p.Trp499Cys), shows partial perinatal lethality and recapitulates the key findings of the human phenotype, including diaphragm and abdominal-wall defects. Both the mouse model and one adult human male with a CDH-associated PLS3 variant were observed to have increased rather than decreased bone mineral density. Together, these clinical and functional data in humans and mice reveal that specific missense variants affecting the actin-binding domains of PLS3 might have a gain-of-function effect and cause a Mendelian congenital disorder.

Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

Cost-effective Whole Exome Sequencing discovers pathogenic variant causing Neurofibromatosis type 1 in a family from Jammu and Kashmir, India

Neurofibromatosis type 1 (NF1) is a multisystemic hereditary disorder associated with an increased risk of benign and malignant tumor formation predominantly on the skin, bone, and peripheral nervous system. It has been reported that out of all the NF1 cases, more than 95% cases develop the disease due to heterozygous loss-of-function variants in Neurofibromin (NF1) gene. However, identification of NF1 causative variants by presently recommended method of gene-targeted Sanger sequencing is challenging and cost-intensive due to the large size of the NF1gene with 60 exons spanning about 350 kb. Further, conducting the genetic studies is difficult in low resource regions and among families with the limited financial capabilities, restricting them from availing diagnostic as well as proper disease management measures. Here, we studied a three-generation family from Jammu and Kashmir state in India, with multiple affected family members showing clinical indications of NF1. We combinedly used two applications, Whole Exome Sequencing (WES) and Sanger sequencing, for this study and discovered a nonsense variant NM_000267.3:c.2041C>T (NP_000258.1:p.Arg681Ter*) in exon 18 of NF1 gene in a cost effective manner. In silico analyses further substantiated the pathogenicity of this novel variant. The study also emphasized on the role of Next Generation Sequencing (NGS) as a cost-effective method for the discovery of pathogenic variants in disorders with known phenotypes found in large sized candidate genes. The current study is the first study based on the genetic characterization of NF1 from Jammu and Kashmir-India, highlighting the importance of the described methodology adopted for the identification and understanding of the disease in low resource region. The early diagnosis of genetic disorders would open the door to appropriate genetic counseling, reducing the disease burden in the affected families and the general population at large.

Pharmacogenomic and Statistical Analysis

Genetic variants can alter response to drugs and other therapeutic interventions. The study of this phenomenon, called pharmacogenomics, is similar in many ways to other types of genetic studies but has distinct methodological and statistical considerations. Genetic variants involved in the processing of exogenous compounds exhibit great diversity and complexity, and the phenotypes studied in pharmacogenomics are also more complex than typical genetic studies. In this chapter, we review basic concepts in pharmacogenomic study designs, data generation techniques, statistical analysis approaches, and commonly used methods and briefly discuss the ultimate translation of findings to clinical care.

Personalized neoantigen vaccine NEO-PV-01 with chemotherapy and anti-PD-1 as first-line treatment for non-squamous non-small cell lung cancer

Neoantigens arising from mutations in tumor DNA provide targets for immune-based therapy. Here, we report the clinical and immune data from a Phase Ib clinical trial of a personalized neoantigen-vaccine NEO-PV-01 in combination with pemetrexed, carboplatin, and pembrolizumab as first-line therapy for advanced non-squamous non-small cell lung cancer (NSCLC). This analysis of 38 patients treated with the regimen demonstrated no treatment-related serious adverse events. Multiple parameters including baseline tumor immune infiltration and on-treatment circulating tumor DNA levels were highly correlated with clinical response. De novo neoantigen-specific CD4⁺ and CD8⁺ T cell responses were observed post-vaccination. Epitope spread to non-vaccinating neoantigens, including responses to KRAS G12C and G12V mutations, were detected post-vaccination. Neoantigen-specific CD4⁺ T cells generated post-vaccination revealed effector and cytotoxic phenotypes with increased CD4⁺ T cell infiltration in the post-vaccine tumor biopsy. Collectively, these data support the safety and immunogenicity of this regimen in advanced non-squamous NSCLC.

Clinical diagnosis of metabolic disorders using untargeted metabolomic profiling and disease-specific networks learned from profiling data

Untargeted metabolomics is a global molecular profiling technology that can be used to screen for inborn errors of metabolism (IEMs). Metabolite perturbations are evaluated based on current knowledge of specific metabolic pathway deficiencies, a manual diagnostic process that is qualitative, has limited scalability, and is not equipped to learn from accumulating clinical data. Our purpose was to improve upon manual diagnosis of IEMs in the clinic by developing novel computational methods for analyzing untargeted metabolomics data. We employed CTD, an automated computational diagnostic method that "connects the dots" between metabolite perturbations observed in individual metabolomics profiling data and modules identified in disease-specific metabolite co-perturbation networks learned from prior profiling data. We also extended CTD to calculate distances between any two individuals (CTDncd) and between an individual and a disease state (CTDdm), to provide additional network-quantified predictors for use in diagnosis. We show that across 539 plasma samples, CTD-based network-quantified measures can reproduce accurate diagnosis of 16 different IEMs, including adenylosuccinase deficiency, argininemia, argininosuccinic aciduria, aromatic L-amino acid decarboxylase deficiency, cerebral creatine deficiency syndrome type 2, citrullinemia, cobalamin biosynthesis defect, GABA-transaminase deficiency, glutaric acidemia type 1, maple syrup urine disease, methylmalonic aciduria, ornithine transcarbamylase deficiency, phenylketonuria, propionic acidemia, rhizomelic chondrodysplasia punctata, and the Zellweger spectrum disorders. Our approach can be used to supplement information from biochemical pathways and has the potential to significantly enhance the interpretation of variants of uncertain significance uncovered by exome sequencing. CTD, CTDdm, and CTDncd can serve as an essential toolset for biological interpretation of untargeted metabolomics data that overcomes limitations associated with manual diagnosis to assist diagnosticians in clinical decision-making. By automating and quantifying the interpretation of perturbation patterns, CTD can improve the speed and confidence by which clinical laboratory directors make diagnostic and treatment decisions, while automatically improving performance with new case data.

Unraveling the Genetic Architecture of Hepatoblastoma Risk: Birth Defects and Increased Burden of Germline Damaging Variants in Gastrointestinal/Renal Cancer Predisposition and DNA Repair Genes

The ultrarare hepatoblastoma (HB) is the most common pediatric liver cancer. HB risk is related to a few rare syndromes, and the molecular bases remain elusive for most cases. We investigated the burden of rare damaging germline variants in 30 Brazilian patients with HB and the presence of additional clinical signs. A high frequency of prematurity (20%) and birth defects (37%), especially craniofacial (17%, including craniosynostosis) and kidney (7%) anomalies, was observed. Putative pathogenic or likely pathogenic monoallelic germline variants mapped to 10 cancer predisposition genes (CPGs: APC, CHEK2, DROSHA, ERCC5, FAH, MSH2, MUTYH, RPS19, TGFBR2 and VHL) were detected in 33% of the patients, only 40% of them with a family history of cancer. These findings showed a predominance of CPGs with a known link to gastrointestinal/colorectal and renal cancer risk. A remarkable feature was an enrichment of rare damaging variants affecting different classes of DNA repair genes, particularly those known as Fanconi anemia genes. Moreover, several potentially deleterious variants mapped to genes impacting liver functions were disclosed. To our knowledge, this is the largest assessment of rare germline variants in HB patients to date, contributing to elucidate the genetic architecture of HB risk.

Quantitative dissection of multilocus pathogenic variation in an Egyptian infant with severe neurodevelopmental disorder resulting from multiple molecular diagnoses

Genomic sequencing and clinical genomics have demonstrated that substantial subsets of atypical and/or severe disease presentations result from multilocus pathogenic variation (MPV) causing blended phenotypes. In an infant with a severe neurodevelopmental disorder, four distinct molecular diagnoses were found by exome sequencing (ES). The blended phenotype that includes brain malformation, dysmorphism, and hypotonia was dissected using the Human Phenotype Ontology (HPO). ES revealed variants in CAPN3 (c.259C > G:p.L87V), MUSK (c.1781C > T:p.A594V), NAV2 (c.1996G > A:p.G666R), and ZC4H2 (c.595A > C:p.N199H). CAPN3, MUSK, and ZC4H2 are established disease genes linked to limb-girdle muscular dystrophy (OMIM# 253600), congenital myasthenia (OMIM# 616325), and Wieacker-Wolff syndrome (WWS; OMIM# 314580), respectively. NAV2 is a retinoic-acid responsive novel disease gene candidate with biological roles in neurite outgrowth and cerebellar dysgenesis in mouse models. Using semantic similarity, we show that no gene identified by ES individually explains the proband phenotype, but rather the totality of the clinically observed disease is explained by the combination of disease-contributing effects of the identified genes. These data reveal that multilocus pathogenic variation can result in a blended phenotype with each gene affecting a different part of the nervous system and nervous system-muscle connection. We provide evidence from this n = 1 study that in patients with MPV and complex blended phenotypes resulting from multiple molecular diagnoses, quantitative HPO analysis can allow for dissection of phenotypic contribution of both established disease genes and novel disease gene candidates not yet proven to cause human disease.

Genetic errors of immunity distinguish pediatric nonmalignant lymphoproliferative disorders

BACKGROUND

Pediatric nonmalignant lymphoproliferative disorders (PLPDs) are clinically and genetically heterogeneous. Long-standing immune dysregulation and lymphoproliferation in children may be life-threatening, and a paucity of data exists to guide evaluation and treatment of children with PLPD.

OBJECTIVE

The primary objective of this study was to ascertain the spectrum of genomic immunologic defects in PLPD. Secondary objectives included characterization of clinical outcomes and associations between genetic diagnoses and those outcomes.

METHODS

PLPD was defined by persistent lymphadenopathy, lymph organ involvement, or lymphocytic infiltration for more than 3 months, with or without chronic or significant Epstein-Barr virus (EBV) infection. Fifty-one subjects from 47 different families with PLPD were analyzed using whole exome sequencing.

RESULTS

Whole exome sequencing identified likely genetic errors of immunity in 51% to 62% of families (53% to 65% of affected children). Presence of a genetic etiology was associated with younger age and hemophagocytic lymphohistiocytosis. Ten-year survival for the cohort was 72.4%, and patients with viable genetic diagnoses had a higher survival rate (82%) compared to children without a genetic explanation (48%, P = .03). Survival outcomes for individuals with EBV-associated disease and no genetic explanation were particularly worse than outcomes for subjects with EBV-associated disease and a genetic explanation (17% vs 90%; P = .002). Ascertainment of a molecular diagnosis provided targetable treatment options for up to 18 individuals and led to active management changes for 12 patients.

CONCLUSIONS

PLPD defines children at high risk for mortality, and whole exome sequencing informs clinical risks and therapeutic opportunities for this diagnosis.

Novel pathogenic variants and quantitative phenotypic analyses of Robinow syndrome: WNT signaling perturbation and phenotypic variability

Robinow syndrome (RS) is a genetically heterogeneous disorder with six genes that converge on the WNT/planar cell polarity (PCP) signaling pathway implicated (DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A). RS is characterized by skeletal dysplasia and distinctive facial and physical characteristics. To further explore the genetic heterogeneity, paralog contribution, and phenotypic variability of RS, we investigated a cohort of 22 individuals clinically diagnosed with RS from 18 unrelated families. Pathogenic or likely pathogenic variants in genes associated with RS or RS phenocopies were identified in all 22 individuals, including the first variant to be reported in DVL2. We retrospectively collected medical records of 16 individuals from this cohort and extracted clinical descriptions from 52 previously published cases. We performed Human Phenotype Ontology (HPO) based quantitative phenotypic analyses to dissect allele-specific phenotypic differences. Individuals with FZD2 variants clustered into two groups with demonstrable phenotypic differences between those with missense and truncating alleles. Probands with biallelic NXN variants clustered together with the majority of probands carrying DVL1, DVL2, and DVL3 variants, demonstrating no phenotypic distinction between the NXN-autosomal recessive and dominant forms of RS. While phenotypically similar diseases on the RS differential matched through HPO analysis, clustering using phenotype similarity score placed RS-associated phenotypes in a unique cluster containing WNT5A, FZD2, and ROR2 apart from non-RS-associated paralogs. Through human phenotype analyses of this RS cohort and OMIM clinical synopses of Mendelian disease, this study begins to tease apart specific biologic roles for non-canonical WNT-pathway proteins.

COPB2 loss of function causes a coatopathy with osteoporosis and developmental delay

Coatomer complexes function in the sorting and trafficking of proteins between subcellular organelles. Pathogenic variants in coatomer subunits or associated factors have been reported in multi-systemic disorders, i.e., coatopathies, that can affect the skeletal and central nervous systems. We have identified loss-of-function variants in COPB2, a component of the coatomer complex I (COPI), in individuals presenting with osteoporosis, fractures, and developmental delay of variable severity. Electron microscopy of COPB2-deficient subjects' fibroblasts showed dilated endoplasmic reticulum (ER) with granular material, prominent rough ER, and vacuoles, consistent with an intracellular trafficking defect. We studied the effect of COPB2 deficiency on collagen trafficking because of the critical role of collagen secretion in bone biology. COPB2 siRNA-treated fibroblasts showed delayed collagen secretion with retention of type I collagen in the ER and Golgi and altered distribution of Golgi markers. copb2-null zebrafish embryos showed retention of type II collagen, disorganization of the ER and Golgi, and early larval lethality. Copb2+/- mice exhibited low bone mass, and consistent with the findings in human cells and zebrafish, studies in Copb2+/- mouse fibroblasts suggest ER stress and a Golgi defect. Interestingly, ascorbic acid treatment partially rescued the zebrafish developmental phenotype and the cellular phenotype in Copb2+/- mouse fibroblasts. This work identifies a form of coatopathy due to COPB2 haploinsufficiency, explores a potential therapeutic approach for this disorder, and highlights the role of the COPI complex as a regulator of skeletal homeostasis.

Exome sequence association study of levels and longitudinal change of cardiovascular risk factor phenotypes in European Americans and African Americans from the Atherosclerosis Risk in Communities Study

Cardiovascular disease (CVD) is responsible for 31% of all deaths worldwide. Among CVD risk factors are age, race, increased systolic blood pressure (BP), and dyslipidemia. Both BP and blood lipids levels change with age, with a dose-dependent relationship between the cumulative exposure to hyperlipidemia and the risk of CVD. We performed an exome sequence association study using longitudinal data with up to 7805 European Americans (EAs) and 3171 African Americans (AAs) from the Atherosclerosis Risk in Communities (ARIC) study. We assessed associations of common (minor allele frequency > 5%) nonsynonymous and splice-site variants and gene-based sets of rare variants with levels and with longitudinal change of seven CVD risk factor phenotypes (BP traits: systolic BP, diastolic BP, pulse pressure; lipids traits: triglycerides, total cholesterol, high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C]). Furthermore, we investigated the relationship of the identified variants and genes with select CVD endpoints. We identified two novel genes: DCLK3 associated with the change of HDL-C levels in AAs and RAB7L1 associated with the change of LDL-C levels in EAs. RAB7L1 is further associated with an increased risk of heart failure in ARIC EAs. Investigation of the contribution of genetic factors to the longitudinal change of CVD risk factor phenotypes promotes our understanding of the etiology of CVD outcomes, stressing the importance of incorporating the longitudinal structure of the cohort data in future analyses.

Oligonucleotide capture sequencing of the SARS-CoV-2 genome and subgenomic fragments from COVID-19 individuals

The newly emerged and rapidly spreading SARS-CoV-2 causes coronavirus disease 2019 (COVID-19). To facilitate a deeper understanding of the viral biology we developed a capture sequencing methodology to generate SARS-CoV-2 genomic and transcriptome sequences from infected patients. We utilized an oligonucleotide probe-set representing the full-length genome to obtain both genomic and transcriptome (subgenomic open reading frames [ORFs]) sequences from 45 SARS-CoV-2 clinical samples with varying viral titers. For samples with higher viral loads (cycle threshold value under 33, based on the CDC qPCR assay) complete genomes were generated. Analysis of junction reads revealed regions of differential transcriptional activity among samples. Mixed allelic frequencies along the 20kb ORF1ab gene in one sample, suggested the presence of a defective viral RNA species subpopulation maintained in mixture with functional RNA in one sample. The associated workflow is straightforward, and hybridization-based capture offers an effective and scalable approach for sequencing SARS-CoV-2 from patient samples.

Short stature and combined immunodeficiency associated with mutations in RGS10

We report the clinical and molecular phenotype of three siblings from one family, who presented with short stature and immunodeficiency and carried uncharacterized variants in RGS10 (c.489_491del:p.E163del and c.G511T:p.A171S). This gene encodes regulator of G protein signaling 10 (RGS10), a member of a large family of GTPase-activating proteins (GAPs) that targets heterotrimeric G proteins to constrain the activity of G protein-coupled receptors, including receptors for chemoattractants. The affected individuals exhibited systemic abnormalities directly related to the RGS10 mutations, including recurrent infections, hypergammaglobulinemia, profoundly reduced lymphocyte chemotaxis, abnormal lymph node architecture, and short stature due to growth hormone deficiency. Although the GAP activity of each RGS10 variant was intact, each protein exhibited aberrant patterns of PKA-mediated phosphorylation and increased cytosolic and cell membrane localization and activity compared to the wild-type protein. We propose that the RGS10 p.E163del and p.A171S mutations lead to mislocalization of the RGS10 protein in the cytosol, thereby resulting in attenuated chemokine signaling. This study suggests that RGS10 is critical for both immune competence and normal hormonal metabolism in humans and that rare RGS10 variants may contribute to distinct systemic genetic disorders.

Hardware Acceleration of Genomics Data Analysis: Challenges and Opportunities

The significant decline in the cost of genome sequencing has dramatically changed the typical bioinformatics pipeline for analysing sequencing data. Where traditionally, the computational challenge of sequencing is now secondary to genomic data analysis. Short read alignment (SRA) is a ubiquitous process within every modern bioinformatics pipeline in the field of genomics and is often regarded as the principal computational bottleneck. Many hardware and software approaches have been provided to solve the challenge of acceleration. However, previous attempts to increase throughput using many-core processing strategies have enjoyed limited success, mainly due to a dependence on global memory for each computational block. The limited scalability and high energy costs of many-core SRA implementations pose a significant constraint in maintaining acceleration. The Networks-On-Chip (NoC) hardware interconnect mechanism has advanced the scalability of many-core computing systems and, more recently, has demonstrated potential in SRA implementations by integrating multiple computational blocks such as pre-alignment filtering and sequence alignment efficiently, while minimising memory latency and global memory access. This paper provides a state of the art review on current hardware acceleration strategies for genomic data analysis, and it establishes the challenges and opportunities of utilising NoCs as a critical building block in next-generation sequencing (NGS) technologies for advancing the speed of analysis.

Novel compound heterozygous missense variants (c.G955A and c.A1822C) of CACNA2D4 likely causing autosomal recessive retinitis pigmentosa in a Chinese patient

Retinitis pigmentosa (RP) is a rare and heterogeneous group of inherited ocular diseases. However, the relationship between CACNA2D4 mutations and RP is not well understood. In this study, a Chinese autosomal recessive retinitis pigmentosa (arRP) pedigree was enrolled and targeted next-generation sequencing was employed for identifying the causative gene in the proband. These steps were followed by confirmatory Sanger sequencing and segregation analysis. RNA-sequencing (RNA-seq) data and semi-quantitative reverse transcription polymerase chain reaction analysis were then applied to examine the expressions in the human and mouse tissues. Novel compound heterozygous, deleterious missense variants of the CACNA2D4 gene, NM_172364.4: c.G955A (p.D319N) and c.A1822C (p.I608L), were identified in the arRP pedigree, co-segregating with the clinical phenotype in the patient. The CACNA2D4 protein is highly conserved among species. The CACNA2D4 mRNA expression showed the highest expression in the retina of humans and in the later four developmental stages/times of retinal tissues in mice, indicating its role in retina/eye functions and developments. This study is the first to identify novel compound heterozygous mutations c.G955A (p.D319N) and c.A1822C (p.I608L) in the CACNA2D4 gene. These might be disease-causing mutations, thereby extending the mutational spectra. The identification of pathogenic CACNA2D4 variants is expected to enhance our understanding of the genotype-phenotype correlations of arRP for disease diagnosis and genetic counseling. The relationship between the CACNA2D4 variants and diseases/phenotypes other than RP has also been reviewed and discussed in this paper.

Next Generation Sequencing Technology in the Clinic and Its Challenges

Simple Summary

Precise identification and annotation of mutations are of utmost importance in clinical oncology. Insights of the DNA sequence can provide meaningful knowledge to unravel the underlying genetics of disease. Hence, tailoring of personalized medicine often relies on specific genomic alteration for treatment efficacy. The aim of this review is to highlight that sequencing harbors much more than just four nucleotides. Moreover, the gradual transition from first to second generation sequencing technologies has led to awareness for choosing the most appropriate bioinformatic analytic tools based on the aim, quality and demand for a specific purpose. Thus, the same raw data can lead to various results reflecting the intrinsic features of different datamining pipelines.

Abstract

Data analysis has become a crucial aspect in clinical oncology to interpret output from next-generation sequencing-based testing. NGS being able to resolve billions of sequencing reactions in a few days has consequently increased the demand for tools to handle and analyze such large data sets. Many tools have been developed since the advent of NGS, featuring their own peculiarities. Increased awareness when interpreting alterations in the genome is therefore of utmost importance, as the same data using different tools can provide diverse outcomes. Hence, it is crucial to evaluate and validate bioinformatic pipelines in clinical settings. Moreover, personalized medicine implies treatment targeting efficacy of biological drugs for specific genomic alterations. Here, we focused on different sequencing technologies, features underlying the genome complexity, and bioinformatic tools that can impact the final annotation. Additionally, we discuss the clinical demand and design for implementing NGS.

Functional Common and Rare ERBB2 Germline Variants Cooperate in Familial and Sporadic Cancer Susceptibility

We investigated a Spanish and Catalan family in which multiple cancer types tracked across three generations, but for which no genetic etiology had been identified. Whole-exome sequencing of germline DNA from multiple affected family members was performed to identify candidate variants to explain this occurrence of familial cancer. We discovered in all cancer-affected family members a single rare heterozygous germline variant (I654V, rs1801201) in ERBB2/HER2, which is located in a transmembrane glycine zipper motif critical for ERBB2-mediated signaling and in complete linkage disequilibrium (D' = 1) with a common polymorphism (I655V, rs1136201) previously reported in some populations as associated with cancer risk. Because multiple cancer types occurred in this family, we tested both the I654V and the I655V variants for association with cancer across multiple tumor types in 6,371 cases of Northern European ancestry drawn from The Cancer Genome Atlas and 6,647 controls, and found that the rare variant (I654V) was significantly associated with an increased risk for cancer (OR = 1.40; P = 0.021; 95% confidence interval (CI), 1.05-1.89). Functional assays performed in HEK 293T cells revealed that both the I655V single mutant (SM) and the I654V;I655V double mutant (DM) stabilized ERBB2 protein and activated ERBB2 signaling, with the DM activating ERBB2 significantly more than the SM alone. Thus, our results suggest a model whereby heritable genetic variation in the transmembrane domain activating ERBB2 signaling is associated with both sporadic and familial cancer risk, with increased ERBB2 stabilization and activation associated with increased cancer risk. PREVENTION RELEVANCE: By performing whole-exome sequencing on germline DNA from multiple cancer-affected individuals belonging to a family in which multiple cancer types track across three generations, we identified and then characterized functional common and rare variation in ERBB2 associated with both sporadic and familial cancer. Our results suggest that heritable variation activating ERBB2 signaling is associated with risk for multiple cancer types, with increases in signaling correlated with increases in risk, and modified by ancestry or family history.

Investigating Genetic Factors Contributing to Variable Expressivity of Class I 17p13.3 Microduplication

17p13.3 microduplications are rare copy number variations (CNVs) associated with variable phenotypes, including facial dysmorphism, developmental delay, intellectual disability, and autism. Typically, when a recognized pathogenic CNV is identified, other genetic factors are not considered. We investigated via whole-exome sequencing the presence of additional variants in four carriers of class I 17p13.3 microduplications. A 730 kb 17p13.3 microduplication was identified in two half-brothers with intellectual disability, but not in a third affected half-brother or blood cells from their normal mother (Family A), thus leading to the hypothesis of maternal germline mosaicism. No additional pathogenic variants were detected in Family A. Two affected siblings carried maternally inherited 450 kb 17p13.3 microduplication (Family B); the three carriers of the microduplication exhibited microcephaly and learning disability/speech impairment of variable degrees. Exome analysis revealed a variant of uncertain significance in RORA, a gene already linked to autism, in the autistic boy; his sister was heterozygous for a CYP1B1 pathogenic variant that could be related to her congenital glaucoma. Besides, both siblings carried a loss-of-function variant in DIP2B, a candidate gene for intellectual disability, which was inherited from their father, who also exhibited learning disability in childhood. In conclusion, additional pathogenic variants were revealed in two affected carriers of class I 17p13.3 microduplication (Family B), probably adding to their phenotypes. These results provided new evidence regarding the contribution of RORA and DIP2B to neurocognitive deficits, and highlighted the importance of full genetic investigation in carriers of CNV syndromes with variable expressivity. Finally, we suggest that microcephaly may be a rare clinical feature also related to the presence of the class I 17p13.3 microduplication.

Whole exome sequencing in a large pedigree with DCM identifies a novel mutation in RBM20

Background: Familial dilated cardiomyopathy (DCM) is genetically heterogeneous and is associated with mutations in at least 40 different genes. Apart from TTN encoding the giant protein Titin, none of these genes have an expected diagnostic yield of more than 5% complicating genetic diagnosis. Whole exome sequencing (WES) is a powerful alternative for the identification of the causal gene, however variant interpretation remains challenging. We report on WES in a large family with autosomal dominant DCM complicated by end stage heart failure and non-sustained ventricular arrhythmias in whom no causative mutation was identified using a targeted gene panel including 28 genes.Methods and results: WES was applied on 2 affected cousins. Stringent filtering of the identified genetic variants was performed including population variant frequencies, in silico analysis, orthologous and paralogous conservation. Subsequently Sanger sequencing was performed for 10 potential disease causing variants in order to confirm the presence of the variant and to evaluate co-segregation. Only one variant in exon 9 of the RBM20 gene (c.2714T > A, p.Met950Lys, NM_001334363) showed full co-segregation in the 7 affected family members resulting in a maximum 2-point LOD score of 2.1 and suggesting this as the pathogenic mutation responsible for the phenotype. Recently mutations in RBM20 have been linked to arrhythmogenic dilated cardiomyopathy caused by defective splicing of the giant sarcomere protein titin and abnormal calcium handling.Conclusions: We report the identification of a novel mutation in RBM20 by WES in a large pedigree with DCM.

PHENOTYPIC VARIABILITY OF RECESSIVE RDH12-ASSOCIATED RETINAL DYSTROPHY

PURPOSE

To characterize the phenotypic variability and report the genetic defects in a cohort of Chinese patients with biallelic variants of the retinol dehydrogenase 12 (RDH12) gene.

METHODS

The study included 38 patients from 38 unrelated families with biallelic pathogenic RDH12 variants. Systematic next-generation sequencing data analysis, Sanger sequencing validation, and segregation analysis were used to identify the pathogenic mutations. Detailed ophthalmic examinations, including electroretinogram, fundus photography, fundus autofluorescence and optical coherence tomography, and statistical analysis were performed to evaluate phenotype variability.

RESULTS

Twenty-five different mutations of RDH12 were identified in the 38 families. Six of these variants were novel. Val146Asp was observed at the highest frequency (23.7%), and it was followed by Arg62Ter (14.5%) and Thr49Met (9.2%). Twenty-three probands were diagnosed with early-onset severe retinal dystrophy, 6 with Leber congenital amaurosis, 7 with autosomal recessive retinitis pigmentosa, and 2 with cone-rod dystrophy. Self-reported nyctalopia occurred in about a half of patients (55.3%) and was significantly more common among older patients (P < 0.01). Nyctalopia was not significantly associated with best-corrected visual acuity (P = 0.72), but older patients had significantly greater best-corrected visual acuity loss (P < 0.01). Only 15.8% of the patients had nystagmus, which was significantly more likely to occur among 36.8% of the patients with hyperopia >3D (P < 0.01) and/or in cases of reduced best-corrected visual acuity (P = 0.01), but was not associated with age (P = 0.87).

CONCLUSION

Several high-frequency RDH12 variants were identified in patients with inherited retinal dystrophies, most of which were missense mutations. Variable but characteristic phenotypes of a progressive nature was observed. Overall, the findings indicated that biallelic RDH12 mutations are a common cause of early-onset retinal dystrophy and a rare cause of cone-rod dystrophy.

A Phase Ib Trial of Personalized Neoantigen Therapy Plus Anti-PD-1 in Patients with Advanced Melanoma, Non-small Cell Lung Cancer, or Bladder Cancer

Neoantigens arise from mutations in cancer cells and are important targets of T cell-mediated anti-tumor immunity. Here, we report the first open-label, phase Ib clinical trial of a personalized neoantigen-based vaccine, NEO-PV-01, in combination with PD-1 blockade in patients with advanced melanoma, non-small cell lung cancer, or bladder cancer. This analysis of 82 patients demonstrated that the regimen was safe, with no treatment-related serious adverse events observed. De novo neoantigen-specific CD4⁺ and CD8⁺ T cell responses were observed post-vaccination in all of the patients. The vaccine-induced T cells had a cytotoxic phenotype and were capable of trafficking to the tumor and mediating cell killing. In addition, epitope spread to neoantigens not included in the vaccine was detected post-vaccination. These data support the safety and immunogenicity of this regimen in patients with advanced solid tumors (Clinicaltrials.gov: NCT02897765).

Congenital diaphragmatic hernia as a prominent feature of a SPECC1L-related syndrome

Congenital diaphragmatic hernias (CDH) confer substantial morbidity and mortality. Genetic defects, including chromosomal anomalies, copy number variants, and sequence variants are identified in ~30% of patients with CDH. A genetic etiology is not yet found in 70% of patients, however there is a growing number of genetic syndromes and single gene disorders associated with CDH. While there have been two reported individuals with X-linked Opitz G/BBB syndrome with MID1 mutations who have CDH as an associated feature, CDH appears to be a much more prominent feature of a SPECC1L-related autosomal dominant Opitz G/BBB syndrome. Features unique to autosomal dominant Opitz G/BBB syndrome include branchial fistulae, omphalocele, and a bicornuate uterus. Here we present one new individual and five previously reported individuals with CDH found to have SPECC1L mutations. These cases provide strong evidence that SPECC1L is a bona fide CDH gene. We conclude that a SPECC1L-related Opitz G/BBB syndrome should be considered in any patient with CDH who has additional features of hypertelorism, a prominent forehead, a broad nasal bridge, anteverted nares, cleft lip/palate, branchial fistulae, omphalocele, and/or bicornuate uterus.

A novel missense variant c.G644A (p.G215E) of the RPGR gene in a Chinese family causes X-linked retinitis pigmentosa

The mutations in patients with X-linked retinitis pigmentosa (xlRP) have not been well described in the Chinese population. In the present study, a five-generation Chinese retinitis pigmentosa (RP) family was recruited; targeted next-generation sequencing (TGS) was used to identify causative genes and Sanger sequencing for co-segregation. RNA-seq data analysis and revere transcriptional-polymerase chain reaction (RT-PCR) were applied to investigate gene expression patterns of RP GTPase regulator (RPGR) in human and Rpgr in mouse. A novel, hemizygous, deleterious and missense variant: c.G644A (p.G215E) in the RPGR gene (NM_000328.2) exon 7 of X-chromosome was identified in the proband, which was co-segregated with the clinical phenotypes in this family. RNA-seq data showed that RPGR is ubiquitously expressed in 27 human tissues with testis in highest, but no eye tissues data. Then the expressions for Rpgr mRNA in mice including eye tissues were conducted and showed that Rpgr transcript is ubiquitously expressed very highly in retina and testis, and highly in other eye tissues including lens, sclera, and cornea; and expressed highly in the six different developmental times of retinal tissue. Ubiquitous expression in different tissues from eye and very high expression in the retina indicated that RPGR plays a vital role in eye functions, particularly in retina. In conclusion, our study is the first to indicate that the novel missense variant c.G644A (p.G215E) in the RPGR gene might be the disease-causing mutation in this xlRP family, expanding mutation spectrum. These findings facilitate better understanding of the molecular pathogenesis of this disease; provide new insights for genetic counseling and healthcare.

Clinical Characteristics and Mutation Analyses of Ovarian Sertoli-Leydig Cell Tumors

BACKGROUND

There are limited studies on Sertoli-Leydig cell tumors (SLCTs) and no data in the population of Chinese patients with SLCTs from the genetic level. In addition, previous studies on SLCTs have focused exclusively on mutations in the DICER1 gene and no data exists on the genetic landscape of SLCTs.

METHODS

Patients with moderately or poorly differentiated SLCTs who underwent surgical resection between January 2012 and October 2018 in our institution were recruited. Whole exome sequencing was performed on formalin-fixed, paraffin-embedded tumor tissue and peripheral blood or normal tissue samples.

RESULTS

Seventeen patients were recruited with 19 tumor samples. The rate of tumor-associated germline mutations was 6 of 17 (35.3%), and that of DICER1 germline mutations was 4 of 17 (23.5%). Regarding clinical relapse, patients with germline tumor-associated mutations had significantly poorer prognosis than those without (p = .007), and those with germline DICER1 mutations were relatively more likely to exhibit clinical relapse, although not to a significant degree (p = .069). Regarding somatic mutations, firstly, the subclone evolution analysis demonstrated that the two tumors on the contralateral ovary were primary tumors, respectively. Secondly, somatic mutations were most commonly found in CDC27 (10/19, 52.6%), DICER1 (4/19, 21.1%), and MUC22 (4/19, 21.1%). And the analysis of cancer cell fractions showed that DICER1 mutations were correlated with tumorigenesis of SLCTs. The rates of germline and somatic DICER1 mutations were higher in patients who were younger than 18 years than those in older patients (p = .022 and p = .001, respectively).

CONCLUSION

Our study indicates that genetic testing may have important clinical significance for patients with SLCTs, particularly for younger patients.

IMPLICATIONS FOR PRACTICE

Bilateral ovarian Sertoli-Leydig cell tumors were verified to be primary tumors from the genetic perspective. The rates of germline and somatic DICER1 mutations were 4 of 17 (23.5%) and 4 of 19 (21.1%), respectively. The rates of germline and somatic DICER1 mutations were higher in patients who were younger than 18 years than those in older patients (p = .022 and p = .001, respectively).

A cost-effective barcode system for maize genetic discrimination based on bi-allelic InDel markers

BACKGROUND

Maize is one of the most important cereal crop all over the world with a complex genome of about 2.3 gigabase, and exhibits tremendous phenotypic and molecular diversity among different germplasms. Along with the phenotype identification, molecular markers have been accepted extensively as an alternative tool to discriminate different genotypes.

RESULTS

By using previous re-sequencing data of 205 lines, bi-allelic insertions and deletions (InDels) all over maize genome were screened, and a barcode system was constructed consisting of 37 bi-allelic insertion-deletion markers with high polymorphism information content (PIC) values, large discriminative size among varieties. The barcode system was measured and determined, different maize hybrids and inbreds were clearly discriminated efficiently with these markers, and hybrids responding parents were accurately determined. Compared with microarray data of more than 200 maize lines, the barcode system can discriminate maize varieties with 1.57% of different loci as a threshold. The barcode system can be used in standardized easy and quick operation with very low cost and minimum equipment requirements.

CONCLUSION

A barcode system was constructed for genetic discrimination of maize lines, including 37 InDel markers with high PIC values and user-friendly. The barcode system was measured and determined for efficient identification of maize lines.

Next Generation Sequencing of 134 Children with Autism Spectrum Disorder and Regression

Approximately 30% of individuals with autism spectrum disorder (ASD) experience developmental regression, the etiology of which remains largely unknown. We performed a complete literature search and identified 47 genes that had been implicated in such cases. We sequenced these genes in a preselected cohort of 134 individuals with regressive autism. In total, 16 variants in 12 genes with evidence supportive of pathogenicity were identified. They were classified as variants of uncertain significance based on ACMG standards and guidelines. Among these were recurring variants in GRIN2A and PLXNB2, variants in genes that were linked to syndromic forms of ASD (GRIN2A, MECP2, CDKL5, SCN1A,PCDH19, UBE3A, and SLC9A6), and variants in the form of oligogenic heterozygosity (EHMT1, SLC9A6, and MFSD8).

Low-level parental somatic mosaic SNVs in exomes from a large cohort of trios with diverse suspected Mendelian conditions

Purpose:

The goal of this study was to assess the scale of low-level parental mosaicism in exome sequencing (ES) databases.

Methods:

We analyzed approximately 2000 family trio ES datasets from the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) and Baylor Genetics (BG). Among apparent de novo single nucleotide variants (SNVs) identified in the affected probands, we selected rare unique variants with variant allele fraction (VAF) between 30-70% in the probands and lower than 10% in one of the parents.

Results:

Out of 102 candidate mosaic variants validated using amplicon-based NGS, droplet digital PCR, or blocker displacement amplification, 27 (26.4%) were confirmed to be low- (VAF between 1-10%) or very low- (VAF <1%) level mosaic. Detection precision in parental samples with two or more alternate reads was 63.6% (BHCMG) and 43.6% (BG). In nine investigated individuals, we observed variability of mosaic ratios among blood, saliva, fibroblast, buccal, hair, and urine samples.

Conclusion:

Our computational pipeline enables robust discrimination between true and false positive candidate mosaic variants and efficient detection of low-level mosaicism in ES samples. We confirm that the presence of two or more alternate reads in the parental sample is a reliable predictor of low-level parental somatic mosaicism.

Characterization of Genetic and Allelic Diversity Amongst Cultivated and Wild Lentil Accessions for Germplasm Enhancement

Intensive breeding of cultivated lentil has resulted in a relatively narrow genetic base, which limits the options to increase crop productivity through selection. Assessment of genetic diversity in the wild gene pool of lentil, as well as characterization of useful and novel alleles/genes that can be introgressed into elite germplasm, presents new opportunities and pathways for germplasm enhancement, followed by successful crop improvement. In the current study, a lentil collection consisting of 467 wild and cultivated accessions that originated from 10 diverse geographical regions was assessed, to understand genetic relationships among different lentil species/subspecies. A total of 422,101 high-confidence SNP markers were identified against the reference lentil genome (cv. CDC Redberry). Phylogenetic analysis clustered the germplasm collection into four groups, namely, Lens culinaris/Lens orientalis, Lens lamottei/Lens odemensis, Lens ervoides, and Lens nigricans. A weak correlation was observed between geographical origin and genetic relationship, except for some accessions of L. culinaris and L. ervoides. Genetic distance matrices revealed a comparable level of variation within the gene pools of L. culinaris (Nei’s coefficient 0.01468–0.71163), L. ervoides (Nei’s coefficient 0.01807–0.71877), and L. nigricans (Nei’s coefficient 0.02188–1.2219). In order to understand any genic differences at species/subspecies level, allele frequencies were calculated from a subset of 263 lentil accessions. Among all cultivated and wild lentil species, L. nigricans exhibited the greatest allelic differentiation across the genome compared to all other species/subspecies. Major differences were observed on six genomic regions with the largest being on Chromosome 1 (c. 1 Mbp). These results indicate that L. nigricans is the most distantly related to L. culinaris and additional structural variations are likely to be identified from genome sequencing studies. This would provide further insights into evolutionary relationships between cultivated and wild lentil germplasm, for germplasm improvement and introgression.

Comprehensive fundamental somatic variant calling and quality management strategies for human cancer genomes

Next-generation sequencing (NGS) technology has revolutionised human cancer research, particularly via detection of genomic variants with its ultra-high-throughput sequencing and increasing affordability. However, the inundation of rich cancer genomics data has resulted in significant challenges in its exploration and translation into biological insights. One of the difficulties in cancer genome sequencing is software selection. Currently, multiple tools are widely used to process NGS data in four stages: raw sequence data pre-processing and quality control (QC), sequence alignment, variant calling and annotation and visualisation. However, the differences between these NGS tools, including their installation, merits, drawbacks and application, have not been fully appreciated. Therefore, a systematic review of the functionality and performance of NGS tools is required to provide cancer researchers with guidance on software and strategy selection. Another challenge is the multidimensional QC of sequencing data because QC can not only report varied sequence data characteristics but also reveal deviations in diverse features and is essential for a meaningful and successful study. However, monitoring of QC metrics in specific steps including alignment and variant calling is neglected in certain pipelines such as the 'Best Practices Workflows' in GATK. In this review, we investigated the most widely used software for the fundamental analysis and QC of cancer genome sequencing data and provided instructions for selecting the most appropriate software and pipelines to ensure precise and efficient conclusions. We further discussed the prospects and new research directions for cancer genomics.

Variants in the Kisspeptin-GnRH Pathway Modulate the Hormonal Profile and Reproductive Outcomes

Kisspeptin has been identified as a key regulatory protein in the release of gonadotropin-releasing hormone (GnRH), which subsequently increases gonadotropin secretion during puberty to establish reproductive function and regulate the hypothalamic-pituitary-gonadal axis. The effects of variants in the KISS1, KISS1R, and GNRHR genes and their possible association with assisted reproduction outcomes remain to be elucidated. In this study, we used next-generation sequencing to investigate the associations of the genetic diversity at the candidate loci for KISS1, KISS1R, and GNRHR with the hormonal profiles and reproductive outcomes in 86 women who underwent in vitro fertilization treatments. Variants in the KISS1 and KISS1R genes were associated with luteinizing hormone (rs35431622:T>C), anti-Mullerian hormone (rs71745629delT), follicle-stimulating hormone (rs73507529:C>A), and estradiol (rs73507527:G>A, rs350130:A>G, and rs73507529:C>A) levels, as well as with reproductive outcomes such as the number of oocytes retrieved (s35431622:T>C), metaphasis II oocytes (rs35431622:T>C), and embryos (rs1132506:G>C). Additionally, variants in the GNRHR UTR3' (rs1038426:C>A, rs12508464:A>C, rs13150734:C>A, rs17635850:A>G, rs35683646:G>A, rs35610027:C>G, rs35845954:T>C, rs17635749:C>T, and rs7666201:C>T) were associated with low prolactin levels. A conjoint analysis of clinical, hormonal, and genetic variables using a generalized linear model identified two variants of the KISS1 gene (rs71745629delT and rs1132506:G>C) that were significantly associated with hormonal variations and reproductive outcomes. The findings suggest that variants in KISS1, KISS1R, and GNRHR genes can modulate hormone levels and reproductive outcomes.

Deficiencies in vesicular transport mediated by TRAPPC4 are associated with severe syndromic intellectual disability

The conserved transport protein particle (TRAPP) complexes regulate key trafficking events and are required for autophagy. TRAPPC4, like its yeast Trs23 orthologue, is a core component of the TRAPP complexes and one of the essential subunits for guanine nucleotide exchange factor activity for Rab1 GTPase. Pathogenic variants in specific TRAPP subunits are associated with neurological disorders. We undertook exome sequencing in three unrelated families of Caucasian, Turkish and French-Canadian ethnicities with seven affected children that showed features of early-onset seizures, developmental delay, microcephaly, sensorineural deafness, spastic quadriparesis and progressive cortical and cerebellar atrophy in an effort to determine the genetic aetiology underlying neurodevelopmental disorders. All seven affected subjects shared the same identical rare, homozygous, potentially pathogenic variant in a non-canonical, well-conserved splice site within TRAPPC4 (hg19:chr11:g.118890966A>G; TRAPPC4: NM_016146.5; c.454+3A>G). Single nucleotide polymorphism array analysis revealed there was no haplotype shared between the tested Turkish and Caucasian families suggestive of a variant hotspot region rather than a founder effect. In silico analysis predicted the variant to cause aberrant splicing. Consistent with this, experimental evidence showed both a reduction in full-length transcript levels and an increase in levels of a shorter transcript missing exon 3, suggestive of an incompletely penetrant splice defect. TRAPPC4 protein levels were significantly reduced whilst levels of other TRAPP complex subunits remained unaffected. Native polyacrylamide gel electrophoresis and size exclusion chromatography demonstrated a defect in TRAPP complex assembly and/or stability. Intracellular trafficking through the Golgi using the marker protein VSVG-GFP-ts045 demonstrated significantly delayed entry into and exit from the Golgi in fibroblasts derived from one of the affected subjects. Lentiviral expression of wild-type TRAPPC4 in these fibroblasts restored trafficking, suggesting that the trafficking defect was due to reduced TRAPPC4 levels. Consistent with the recent association of the TRAPP complex with autophagy, we found that the fibroblasts had a basal autophagy defect and a delay in autophagic flux, possibly due to unsealed autophagosomes. These results were validated using a yeast trs23 temperature sensitive variant that exhibits constitutive and stress-induced autophagic defects at permissive temperature and a secretory defect at restrictive temperature. In summary we provide strong evidence for pathogenicity of this variant in a member of the core TRAPP subunit, TRAPPC4 that associates with vesicular trafficking and autophagy defects. This is the first report of a TRAPPC4 variant, and our findings add to the growing number of TRAPP-associated neurological disorders.

Detecting sequence variants in clinically important protozoan parasites

Second and third generation sequencing methods are crucial for population genetic studies, and variant detection is a popular approach for exploiting this sequence data. While mini- and microsatellites are historically useful markers for studying important Protozoa such as Toxoplasma and Plasmodium spp., detecting non-repetitive variants such as those found in genes can be fundamental to investigating a pathogen's biology. These variants, namely single nucleotide polymorphisms and insertions and deletions, can help elucidate the genetic basis of an organism's pathogenicity, identify selective pressures, and resolve phylogenetic relationships. They also have the added benefit of possessing a comparatively low mutation rate, which contributes to their stability. However, there is a plethora of variant analysis tools with nuanced pipelines and conflicting recommendations for best practise, which can be confounding. This lack of standardisation means that variant analysis requires careful parameter optimisation, an understanding of its limitations, and the availability of high quality data. This review explores the value of variant detection when applied to non-model organisms such as clinically important protozoan pathogens. The limitations of current methods are discussed, including special considerations that require the end-users' attention to ensure that the results generated are reproducible, and the biological conclusions drawn are valid.

Genomic Profiling of Childhood Tumor Patient-Derived Xenograft Models to Enable Rational Clinical Trial Design

Accelerating cures for children with cancer remains an immediate challenge as a result of extensive oncogenic heterogeneity between and within histologies, distinct molecular mechanisms evolving between diagnosis and relapsed disease, and limited therapeutic options. To systematically prioritize and rationally test novel agents in preclinical murine models, researchers within the Pediatric Preclinical Testing Consortium are continuously developing patient-derived xenografts (PDXs)-many of which are refractory to current standard-of-care treatments-from high-risk childhood cancers. Here, we genomically characterize 261 PDX models from 37 unique pediatric cancers; demonstrate faithful recapitulation of histologies and subtypes; and refine our understanding of relapsed disease. In addition, we use expression signatures to classify tumors for TP53 and NF1 pathway inactivation. We anticipate that these data will serve as a resource for pediatric oncology drug development and will guide rational clinical trial design for children with cancer.

Association of rare non-coding SNVs in the lung-specific FOXF1 enhancer with a mitigation of the lethal ACDMPV phenotype

Haploinsufficiency of FOXF1 causes alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a lethal neonatal lung developmental disorder. We describe two similar heterozygous CNV deletions involving the FOXF1 enhancer and re-analyze FOXF1 missense mutation, all associated with an unexpectedly mitigated disease phenotype. In one case, the deletion of the maternal allele of the FOXF1 enhancer caused pulmonary hypertension and histopathologically diagnosed MPV without the typical ACD features. In the second case, the deletion of the paternal enhancer resulted in ACDMPV rather than the expected neonatal lethality. In both cases, FOXF1 expression in lung tissue was higher than usually seen or expected in patients with similar deletions, suggesting an increased activity of the remaining allele of the enhancer. Sequencing of these alleles revealed two rare SNVs, rs150502618-A and rs79301423-T, mapping to the partially overlapping binding sites for TFAP2s and CTCF in the core region of the enhancer. Moreover, in a family with three histopathologically-diagnosed ACDMPV siblings whose missense FOXF1 mutation was inherited from the healthy non-mosaic carrier mother, we have identified a rare SNV rs28571077-A within 2-kb of the above-mentioned non-coding SNVs in the FOXF1 enhancer in the mother, that was absent in the affected newborns and 13 unrelated ACDMPV patients with CNV deletions of this genomic region. Based on the low population frequencies of these three variants, their absence in ACDMPV patients, the results of reporter assay, RNAi and EMSA experiments, and in silico predictions, we propose that the described SNVs might have acted on FOXF1 enhancer as hypermorphs.

Variable expressivity of syndromic BMP4-related eye, brain, and digital anomalies: A review of the literature and description of three new cases

Microphthalmia with brain and digital anomalies (MCOPS6, MIM# 607932) is an autosomal dominant disorder caused by loss-of-function variants or large deletions involving BMP4, which encodes bone morphogenetic protein 4, a member of the TGF-β protein superfamily. BMP4 has a number of roles in embryonic development including neurogenesis, lens induction, development of cartilage and bone, urogenital development, limb and digit patterning, hair follicle regeneration, as well as tooth formation. In addition to syndromic microphthalmia, BMP4 variants have been implicated in non-syndromic cleft lip with or without cleft palate and congenital healed cleft lip indicating different allelic presentations. MCOPS6 subjects may also lack some of the major phenotypic hallmarks of the disorder, including microphthalmia, indicating variable expressivity. As only a handful of individuals with MCOPS6 have been described, we review the clinical findings in previously reported cases with either deletions or loss-of-function variants in BMP4. We describe three new cases, including two subjects with novel deletions and one subject with a likely pathogenic de novo nonsense variant [c.1052C>G, p.(S351*)] in BMP4. One of the subjects had dual molecular diagnoses including a co-occurring microdeletion at 17q21.31 associated with Koolen de Vries syndrome, which has a partially overlapping disease phenotype. None of these individuals had clinically apparent microphthalmia or anopthalmia, which have been reported in a majority of previously described cases. One subject had exophthalmia and strabismus, while another had bilateral Peters anomaly and sclerocornea, thus expanding the phenotype associated with BMP4 loss-of-function variants.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

The Genomics of Arthrogryposis, a Complex Trait: Candidate Genes and Further Evidence for Oligogenic Inheritance

Arthrogryposis is a clinical finding that is present either as a feature of a neuromuscular condition or as part of a systemic disease in over 400 Mendelian conditions. The underlying molecular etiology remains largely unknown because of genetic and phenotypic heterogeneity. We applied exome sequencing (ES) in a cohort of 89 families with the clinical sign of arthrogryposis. Additional molecular techniques including array comparative genomic hybridization (aCGH) and Droplet Digital PCR (ddPCR) were performed on individuals who were found to have pathogenic copy number variants (CNVs) and mosaicism, respectively. A molecular diagnosis was established in 65.2% (58/89) of families. Eleven out of 58 families (19.0%) showed evidence for potential involvement of pathogenic variation at more than one locus, probably driven by absence of heterozygosity (AOH) burden due to identity-by-descent (IBD). RYR3, MYOM2, ERGIC1, SPTBN4, and ABCA7 represent genes, identified in two or more families, for which mutations are probably causative for arthrogryposis. We also provide evidence for the involvement of CNVs in the etiology of arthrogryposis and for the idea that both mono-allelic and bi-allelic variants in the same gene cause either similar or distinct syndromes. We were able to identify the molecular etiology in nine out of 20 families who underwent reanalysis. In summary, our data from family-based ES further delineate the molecular etiology of arthrogryposis, yielded several candidate disease-associated genes, and provide evidence for mutational burden in a biological pathway or network. Our study also highlights the importance of reanalysis of individuals with unsolved diagnoses in conjunction with sequencing extended family members.

Quality control and integration of genotypes from two calling pipelines for whole genome sequence data in the Alzheimer's disease sequencing project

The Alzheimer's Disease Sequencing Project (ADSP) performed whole genome sequencing (WGS) of 584 subjects from 111 multiplex families at three sequencing centers. Genotype calling of single nucleotide variants (SNVs) and insertion-deletion variants (indels) was performed centrally using GATK-HaplotypeCaller and Atlas V2. The ADSP Quality Control (QC) Working Group applied QC protocols to project-level variant call format files (VCFs) from each pipeline, and developed and implemented a novel protocol, termed "consensus calling," to combine genotype calls from both pipelines into a single high-quality set. QC was applied to autosomal bi-allelic SNVs and indels, and included pipeline-recommended QC filters, variant-level QC, and sample-level QC. Low-quality variants or genotypes were excluded, and sample outliers were noted. Quality was assessed by examining Mendelian inconsistencies (MIs) among 67 parent-offspring pairs, and MIs were used to establish additional genotype-specific filters for GATK calls. After QC, 578 subjects remained. Pipeline-specific QC excluded ~12.0% of GATK and 14.5% of Atlas SNVs. Between pipelines, ~91% of SNV genotypes across all QCed variants were concordant; 4.23% and 4.56% of genotypes were exclusive to Atlas or GATK, respectively; the remaining ~0.01% of discordant genotypes were excluded. For indels, variant-level QC excluded ~36.8% of GATK and 35.3% of Atlas indels. Between pipelines, ~55.6% of indel genotypes were concordant; while 10.3% and 28.3% were exclusive to Atlas or GATK, respectively; and ~0.29% of discordant genotypes were. The final WGS consensus dataset contains 27,896,774 SNVs and 3,133,926 indels and is publicly available.

Targeted sequencing identifies a missense variant in the BEST3 gene associated with antihypertensive response to hydrochlorothiazide

Chromosome 12q15 was identified in Genetic Epidemiology of Response Assessment (GERA) and replicated in Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) for its association with blood pressure (BP) response to hydrochlorothiazide (HCTZ). However, the functional variant is unknown and we aimed to identify the likely functional variants through targeted sequencing. The chromosome 12q15 region was sequenced in 397 best and worst responders to HCTZ in PEAR (N=199) and GERA (N=198) hypertensive study participants. Logistic regression was used for the association analysis adjusting for age, sex, race, and principal components 1 and 2. For validation, the significant single nucleotide polymorphism was tested for association with the change in systolic (ΔSBP) and diastolic BP (ΔDBP) post-treatment in the entire PEAR (N=370) and GERA (N=570) cohorts. A novel missense polymorphism (G>A, Pro383Leu) in BEST3, rs61747221, was significantly associated with better HCTZ response (P=0.0021, odds ratio=2.05). It was validated in the entire cohort of PEAR (ΔSBP: P=0.021, β=-1.60, ΔDBP: P=0.023, β=-1.08) and GERA (ΔSBP: P=0.028, β=-1.95, ΔDBP: P=0.032, β=-1.28). BEST3 encodes the calcium sensitive chloride channel in the vascular smooth muscle implicated in the regulation of BP, especially in response to vasoconstrictors like angiotensin II. These results suggest that BEST3 is involved in the chronic BP lowering mechanism of thiazides and highlight its importance as a genetic predictor of the BP response to thiazide diuretics.

Interchromosomal template-switching as a novel molecular mechanism for imprinting perturbations associated with Temple syndrome

Background

Intrachromosomal triplications (TRP) can contribute to disease etiology via gene dosage effects, gene disruption, position effects, or fusion gene formation. Recently, post-zygotic de novo triplications adjacent to copy-number neutral genomic intervals with runs of homozygosity (ROH) have been shown to result in uniparental isodisomy (UPD). The genomic structure of these complex genomic rearrangements (CGRs) shows a consistent pattern of an inverted triplication flanked by duplications (DUP-TRP/INV-DUP) formed by an iterative DNA replisome template-switching mechanism during replicative repair of a single-ended, double-stranded DNA (seDNA), the ROH results from an interhomolog or nonsister chromatid template switch. It has been postulated that these CGRs may lead to genetic abnormalities in carriers due to dosage-sensitive genes mapping within the copy-number variant regions, homozygosity for alleles at a locus causing an autosomal recessive (AR) disease trait within the ROH region, or imprinting-associated diseases.

Methods

Here, we report a family wherein the affected subject carries a de novo 2.2-Mb TRP followed by 42.2 Mb of ROH and manifests clinical features overlapping with those observed in association with chromosome 14 maternal UPD (UPD(14)mat). UPD(14)mat can cause clinical phenotypic features enabling a diagnosis of Temple syndrome. This CGR was then molecularly characterized by high-density custom aCGH, genome-wide single-nucleotide polymorphism (SNP) and methylation arrays, exome sequencing (ES), and the Oxford Nanopore long-read sequencing technology.

Results

We confirmed the postulated DUP-TRP/INV-DUP structure by multiple orthogonal genomic technologies in the proband. The methylation status of known differentially methylated regions (DMRs) on chromosome 14 revealed that the subject shows the typical methylation pattern of UPD(14)mat. Consistent with these molecular findings, the clinical features overlap with those observed in Temple syndrome, including speech delay.

Conclusions

These data provide experimental evidence that, in humans, triplication can lead to segmental UPD and imprinting disease. Importantly, genotype/phenotype analyses further reveal how a post-zygotically generated complex structural variant, resulting from a replication-based mutational mechanism, contributes to expanding the clinical phenotype of known genetic syndromes. Mechanistically, such events can distort transmission genetics resulting in homozygosity at a locus for which only one parent is a carrier as well as cause imprinting diseases.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0633-y) contains supplementary material, which is available to authorized users.

IFT52 as a Novel Candidate for Ciliopathies Involving Retinal Degeneration

Purpose

Mutations in the intraflagellar transport protein 52 homolog (IFT52) gene are reported to interrupt ciliary function and cause short-rib thoracic dysplasia (SRTD), a specific form of skeletal ciliopathy. However, the roles of these mutations in retinal ciliopathy are inexplicit. We herein aim to study the impact of IFT52 mutations in retinopathies.

Methods

A patient with syndromic ciliopathy, presenting mild SRTD (skeletal ciliopathy) and Liber congenital amaurosis (LCA; retinal ciliopathy), and nine unaffected family members were recruited. Comprehensive systemic evaluations, including ophthalmic tests, were received by the patient. Whole genome sequencing (WGS) was applied for genetic annotation. An in vitro cell system was employed to study the pathogenicity of the variant.

Results

WGS identified a homozygous missense variation in IFT52, c.556A>G (p.T186A), carried by the patient but absent in both unaffected siblings. In silico analysis supported the pathogenic nature of this highly conserved variant. Structural analysis suggested that this substitution could generate a novel hydrogen bond between the mutated residue 186 and proline at residue 192, thus potentially interrupting the tertiary structure and the stability of the IFT52 protein. In vitro cellular study indicated that this mutation might disturb the stability of encoded IFT52 protein and dramatically disrupt cilia elongation in hTERT-RPE1 cells in a loss-of-function manner.

Conclusions

This report expands ocular phenotypes of IFT52 mutation-caused ciliopathy to include retinal ciliopathy and demonstrates its deleterious nature in interrupting primary ciliary function. Our study hence highlights the need for screening for IFT52 mutations in LCA patients and ophthalmic reviews of patients carrying IFT52 mutations.

An Exome-Wide Sequencing Study of the GOLDN Cohort Reveals Novel Associations of Coding Variants and Fasting Plasma Lipids

Background: Associations of both common and rare genetic variants with fasting blood lipids have been extensively studied. However, most of the rare coding variants associated with lipids are population-specific, and exploration of genetic data from diverse population samples may enhance the identification of novel associations with rare variants. Results: We searched for novel coding genetic variants associated with fasting lipid levels in 894 samples from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) with exome-wide sequencing-based genotype data. In single variant tests, one variant (rs11171663 in ITGA7) was associated with fasting triglyceride levels (P = 7.66E-08), explaining approximately 3.2% of the total trait variance. In gene-based tests, we found statistically significant associations between ITGA7 (P = 1.77E-07) and SLCO2A1 (P = 7.18E-07) and triglycerides, as well as between POT1 (P = 3.00E-07) and low-density lipoprotein cholesterol. In another independent replication cohort consisting of 3,183 African American samples from Hypertension Genetic Epidemiology Network (HyperGEN) and the Genetic Epidemiology Network of Arteriopathy (GENOA), the top genes achieved P-values of 0.04 (ITGA7), 0.08 (SLCO2A1), and 0.02 (POT1). In GOLDN, gene transcript levels of ITGA7 and SLCO2A1 were associated with fasting triglycerides (P = 0.07 and P = 0.02), highlighting functional relevance of our findings. Conclusion: In this study, we present preliminary evidence of novel rare variant determinants of fasting lipids, and reveal potential underlying molecular mechanisms. Moreover, these results were replicated in an independent cohort. Our findings may inform novel biomarkers of disease risk and treatment targets.