A novel ARNSHI (autosomal recessive non-syndromic hearing impairment) locus, DFNB71, was localized to 8p22-21.3. To map the locus, a whole-genome scan was carried out using DNA samples from a consanguineous seven-generational Pakistani family with profound prelingual ARNSHI. A maximum multipoint LOD score of 4.2 occurred at marker D8S261. The DFNB71 region of homozygosity and 3-unit support interval is flanked by markers D8S1130 and D8S1786. This region has a genetic distance of 19.1 cM and contains 10.6 Mb of sequence.

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

We report two consanguineous families with probands that exhibit intellectual disability, developmental delay, short stature, aphasia, and hypotonia in which homozygous non-synonymous variants were identified in IQSEC1 (GenBank: NM_001134382.3). In a Pakistani family, the IQSEC1 segregating variant is c.1028C>T (p.Thr343Met), while in a Saudi Arabian family the variant is c.962G>A (p.Arg321Gln). IQSEC1-3 encode guanine nucleotide exchange factors for the small GTPase ARF6 and their loss affects a variety of actin-dependent cellular processes, including AMPA receptor trafficking at synapses. The ortholog of IQSECs in the fly is schizo and its loss affects growth cone guidance at the midline in the CNS, also an actin-dependent process. Overexpression of the reference IQSEC1 cDNA in wild-type flies is lethal, but overexpression of the two variant IQSEC1 cDNAs did not affect viability. Loss of schizo caused embryonic lethality that could be rescued to 2nd instar larvae by moderate expression of the human reference cDNA. However, the p.Arg321Gln and p.Thr343Met variants failed to rescue embryonic lethality. These data indicate that the variants behave as loss-of-function mutations. We also show that schizo in photoreceptors is required for phototransduction. Finally, mice with a conditional Iqsec1 deletion in cortical neurons exhibited an increased density of dendritic spines with an immature morphology. The phenotypic similarity of the affecteds and the functional experiments in flies and mice indicate that IQSEC1 variants are the cause of a recessive disease with intellectual disability, developmental delay, and short stature, and that axonal guidance and dendritic projection defects as well as dendritic spine dysgenesis may underlie disease pathogenesis.

Cerebellar malformations are diverse congenital anomalies frequently associated with developmental disability. Although genetic and prenatal non-genetic causes have been described, no systematic analysis has been performed. Here, we present a large-exome sequencing study of Dandy-Walker malformation (DWM) and cerebellar hypoplasia (CBLH). We performed exome sequencing in 282 individuals from 100 families with DWM or CBLH, and we established a molecular diagnosis in 36 of 100 families, with a significantly higher yield for CBLH (51%) than for DWM (16%). The 41 variants impact 27 neurodevelopmental-disorder-associated genes, thus demonstrating that CBLH and DWM are often features of monogenic neurodevelopmental disorders. Though only seven monogenic causes (19%) were identified in more than one individual, neuroimaging review of 131 additional individuals confirmed cerebellar abnormalities in 23 of 27 genetic disorders (85%). Prenatal risk factors were frequently found among individuals without a genetic diagnosis (30 of 64 individuals [47%]). Single-cell RNA sequencing of prenatal human cerebellar tissue revealed gene enrichment in neuronal and vascular cell types; this suggests that defective vasculogenesis may disrupt cerebellar development. Further, de novo gain-of-function variants in PDGFRB, a tyrosine kinase receptor essential for vascular progenitor signaling, were associated with CBLH, and this discovery links genetic and non-genetic etiologies. Our results suggest that genetic defects impact specific cerebellar cell types and implicate abnormal vascular development as a mechanism for cerebellar malformations. We also confirmed a major contribution for non-genetic prenatal factors in individuals with cerebellar abnormalities, substantially influencing diagnostic evaluation and counseling regarding recurrence risk and prognosis.

Gene discovery for Mendelian conditions (MCs) offers a direct path to understanding genome function. Approaches based on next-generation sequencing applied at scale have dramatically accelerated gene discovery and transformed genetic medicine. Finding the genetic basis of ∼6,000-13,000 MCs yet to be delineated will require both technical and computational innovation, but will rely to a larger extent on meaningful data sharing.

We performed exome sequencing to evaluate the underlying molecular cause of a patient with bilateral conductive hearing loss due to multiple ossicular abnormalities as well as symphalangism of the fifth digits. This leads to the identification of a novel heterozygous start codon variant in the NOG gene (c.2T>C:p.Met1?) that hinders normal translation of the noggin protein. Variants in NOG lead to a spectrum of otologic, digit, and joint abnormalities, a combination suggested to be referred to as NOG-related-symphalangism spectrum disorder (NOG-SSD). Conductive hearing loss from such variants may stem from stapes footplate ankylosis, fixation of the malleoincudal joint, or fixation of the incus short process. In this case, the constellation of both stapes and incus fixation, an exceptionally tall stapes suprastructure, thickened long process of the incus, and enlarged incus body was encountered, leading to distinct challenges during otologic surgery to improve hearing thresholds. This case highlights multiple abnormalities to the ossicular chain in a patient with a start codon variant in NOG. We provide detailed imaging data on these malformations as well as surgical considerations and outcomes.

Massively parallel sequencing technologies provide great opportunities for discovering rare susceptibility variants involved in complex disease etiology via large-scale imputation and exome and whole-genome sequence-based association studies. Due to modest effect sizes, large sample sizes of tens to hundreds of thousands of individuals are required for adequately powered studies. Current analytical tools are obsolete when it comes to handling these large datasets. To facilitate the analysis of large-scale sequence-based studies, we developed SEQSpark which implements parallel processing based on Spark to increase the speed and efficiency of performing data quality control, annotation, and association analysis. To demonstrate the versatility and speed of SEQSpark, we analyzed whole-genome sequence data from the UK10K, testing for associations with waist-to-hip ratios. The analysis, which was completed in 1.5 hr, included loading data, annotation, principal component analysis, and single variant and rare variant aggregate association analysis of >9 million variants. For rare variant aggregate analysis, an exome-wide significant association (p < 2.5 × 10-6) was observed with CCDC62 (SKAT-O [p = 6.89 × 10-7], combined multivariate collapsing [p = 1.48 × 10-6], and burden of rare variants [p = 1.48 × 10-6]). SEQSpark was also used to analyze 50,000 simulated exomes and it required 1.75 hr for the analysis of a quantitative trait using several rare variant aggregate association methods. Additionally, the performance of SEQSpark was compared to Variant Association Tools and PLINK/SEQ. SEQSpark was always faster and in some situations computation was reduced to a hundredth of the time. SEQSpark will empower large sequence-based epidemiological studies to quickly elucidate genetic variation involved in the etiology of complex traits.

Non-syndromic cleft lip with or without cleft palate (NS-CL/P) is one of the most common human birth defects and is generally considered a complex trait. Despite numerous loci identified by genome-wide association studies, the effect sizes of common variants are relatively small, with much of the presumed genetic contribution remaining elusive. We report exome-sequencing results in 209 people from 72 multi-affected families with pedigree structures consistent with autosomal-dominant inheritance and variable penetrance. Herein, pathogenic variants are described in four genes encoding components of the p120-catenin complex (CTNND1, PLEKHA7, PLEKHA5) and an epithelial splicing regulator (ESRP2), in addition to the known CL/P-associated gene, CDH1, which encodes E-cadherin. The findings were also validated in a second cohort of 497 people with NS-CL/P, comprising small families and singletons with pathogenic variants in these genes identified in 14% of multi-affected families and 2% of the replication cohort of smaller families. Enriched expression of each gene/protein in human and mouse embryonic oro-palatal epithelia, demonstration of functional impact of CTNND1 and ESRP2 variants, and recapitulation of the CL/P spectrum in Ctnnd1 knockout mice support a causative role in CL/P pathogenesis. These data show that primary defects in regulators of epithelial cell adhesion are the most significant contributors to NS-CL/P identified to date and that inherited and de novo single gene variants explain a substantial proportion of NS-CL/P.

Wiedemann-Rautenstrauch syndrome (WRS), also known as neonatal progeroid syndrome, is a rare disorder of unknown etiology. It has been proposed to be autosomal-recessive and is characterized by variable clinical features, such as intrauterine growth restriction and poor postnatal weight gain, characteristic facial features (triangular appearance to the face, convex nasal profile or pinched nose, and small mouth), widened fontanelles, pseudohydrocephalus, prominent scalp veins, lipodystrophy, and teeth abnormalities. A previous report described a single WRS patient with bi-allelic truncating and splicing variants in POLR3A. Here we present seven additional infants, children, and adults with WRS and bi-allelic truncating and/or splicing variants in POLR3A. POLR3A, the largest subunit of RNA polymerase III, is a DNA-directed RNA polymerase that transcribes many small noncoding RNAs that regulate transcription, RNA processing, and translation. Bi-allelic missense variants in POLR3A have been associated with phenotypes distinct from WRS: hypogonadotropic hypogonadism and hypomyelinating leukodystrophy with or without oligodontia. Our findings confirm the association of bi-allelic POLR3A variants with WRS, expand the clinical phenotype of WRS, and suggest specific POLR3A genotypes associated with WRS and hypomyelinating leukodystrophy.

The urgent need for massively scaled clinical testing for SARS-CoV-2, along with global shortages of critical reagents and supplies, has necessitated development of streamlined laboratory testing protocols. Conventional nucleic acid testing for SARS-CoV-2 involves collection of a clinical specimen with a nasopharyngeal swab in transport medium, nucleic acid extraction, and quantitative reverse transcription PCR (RT-qPCR) (1). As testing has scaled across the world, the global supply chain has buckled, rendering testing reagents and materials scarce (2). To address shortages, we developed SwabExpress, an end-to-end protocol developed to employ mass produced anterior nares swabs and bypass the requirement for transport media and nucleic acid extraction.

Following its emergence in Wuhan, China, in late November or early December 2019, the SARS-CoV-2 virus has rapidly spread throughout the world. On March 11, 2020, the World Health Organization declared Coronavirus Disease 2019 (COVID-19) a pandemic. Genome sequencing of SARS-CoV-2 strains allows for the reconstruction of transmission history connecting these infections. Here, we analyze 346 SARS-CoV-2 genomes from samples collected between 20 February and 15 March 2020 from infected patients in Washington State, USA. We found that the large majority of SARS-CoV-2 infections sampled during this time frame appeared to have derived from a single introduction event into the state in late January or early February 2020 and subsequent local spread, strongly suggesting cryptic spread of COVID-19 during the months of January and February 2020, before active community surveillance was implemented. We estimate a common ancestor of this outbreak clade as occurring between 18 January and 9 February 2020. From genomic data, we estimate an exponential doubling between 2.4 and 5.1 days. These results highlight the need for large-scale community surveillance for SARS-CoV-2 introductions and spread and the power of pathogen genomics to inform epidemiological understanding.

We identified ten persons in six consanguineous families with distal arthrogryposis (DA) who had congenital contractures, scoliosis, and short stature. Exome sequencing revealed that each affected person was homozygous for one of two different rare variants (c.470G>T [p.Cys157Phe] or c.469T>C [p.Cys157Arg]) affecting the same residue of myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF). In a seventh family, a c.487G>A (p.Gly163Ser) variant in MYLPF arose de novo in a father, who transmitted it to his son. In an eighth family comprised of seven individuals with dominantly inherited DA, a c.98C>T (p.Ala33Val) variant segregated in all four persons tested. Variants in MYLPF underlie both dominant and recessively inherited DA. Mylpf protein models suggest that the residues associated with dominant DA interact with myosin whereas the residues altered in families with recessive DA only indirectly impair this interaction. Pathological and histological exam of a foot amputated from an affected child revealed complete absence of skeletal muscle (i.e., segmental amyoplasia). To investigate the mechanism for this finding, we generated an animal model for partial MYLPF impairment by knocking out zebrafish mylpfa. The mylpfa mutant had reduced trunk contractile force and complete pectoral fin paralysis, demonstrating that mylpf impairment most severely affects limb movement. mylpfa mutant muscle weakness was most pronounced in an appendicular muscle and was explained by reduced myosin activity and fiber degeneration. Collectively, our findings demonstrate that partial loss of MYLPF function can lead to congenital contractures, likely as a result of degeneration of skeletal muscle in the distal limb.

Intracranial aneurysms (IAs), abdominal aortic aneurysms (AAAs), and thoracic aortic aneurysms (TAAs) all have a familial predisposition. Given that aneurysm types are known to co-occur, we hypothesized that there may be shared genetic risk factors for IAs, AAAs, and TAAs.

The recent availability of whole-exome sequencing has opened new possibilities for the evaluation of individuals with genetically undiagnosed intellectual disability.

Co-inheritance of α-thalassemia has a significant protective effect on the severity of complications of sickle cell disease (SCD), including stroke. However, little information exists on the association and interactions for the common African ancestral α-thalassemia mutation (-α3.7 deletion) and β-globin traits (HbS trait [SCT] and HbC trait) on important clinical phenotypes such as red blood cell parameters, anemia, and chronic kidney disease (CKD). In a community-based cohort of 2,916 African Americans from the Jackson Heart Study, we confirmed the expected associations between SCT, HbC trait, and the -α3.7 deletion with lower mean corpuscular volume/mean corpuscular hemoglobin and higher red blood cell count and red cell distribution width. In addition to the recently recognized association of SCT with lower estimated glomerular filtration rate and glycated hemoglobin (HbA1c), we observed a novel association of the -α3.7 deletion with higher HbA1c levels. Co-inheritance of each additional copy of the -α3.7 deletion significantly lowered the risk of anemia and chronic kidney disease among individuals with SCT (P-interaction = 0.031 and 0.019, respectively). Furthermore, co-inheritance of a novel α-globin regulatory variant was associated with normalization of red cell parameters in individuals with the -α3.7 deletion and significantly negated the protective effect of α-thalassemia on stroke in 1,139 patients with sickle cell anemia from the Cooperative Study of Sickle Cell Disease (CSSCD) (P-interaction = 0.0049). Functional assays determined that rs11865131, located in the major alpha-globin enhancer MCS-R2, was the most likely causal variant. These findings suggest that common α- and β-globin variants interact to influence hematologic and clinical phenotypes in African Americans, with potential implications for risk-stratification and counseling of individuals with SCD and SCT.

Homeless shelters are a high-risk setting for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission because of crowding and shared hygiene facilities.

Cohen syndrome (CS) is a rare genetic disorder caused by variants of the VPS13B gene. CS patients are affected with a severe form of retinal dystrophy, and in several cases cataracts also develop. The purpose of this study was to investigate the mechanisms and risk factors for cataract in CS, as well as to report on cataract surgeries in CS patients.

Congenital hearing impairment (HI) is genetically heterogeneous making its genetic diagnosis challenging. Investigation of novel HI genes and variants will enhance our understanding of the molecular mechanisms and to aid genetic diagnosis. We performed exome sequencing and analysis using DNA samples from affected members of two large families from Ghana and Pakistan, segregating autosomal-dominant (AD) non-syndromic HI (NSHI). Using in silico approaches, we modeled and evaluated the effect of the likely pathogenic variants on protein structure and function. We identified two likely pathogenic variants in SLC12A2, c.2935G>A:p.(E979K) and c.2939A>T:p.(E980V), which segregate with NSHI in a Ghanaian and Pakistani family, respectively. SLC12A2 encodes an ion transporter crucial in the homeostasis of the inner ear endolymph and has recently been reported to be implicated in syndromic and non-syndromic HI. Both variants were mapped to alternatively spliced exon 21 of the SLC12A2 gene. Exon 21 encodes for 17 residues in the cytoplasmatic tail of SLC12A2, is highly conserved between species, and preferentially expressed in cochlear tissues. A review of previous studies and our current data showed that out of ten families with either AD non-syndromic or syndromic HI, eight (80%) had variants within the 17 amino acid residue region of exon 21 (48 bp), suggesting that this alternate domain is critical to the transporter activity in the inner ear. The genotypic spectrum of SLC12A2 was expanded and the involvement of SLC12A2 in ADNSHI was confirmed. These results also demonstrate the role that SLC12A2 plays in ADNSHI in diverse populations including sub-Saharan Africans.

The 2-oxoglutarate dehydrogenase-like (OGDHL) protein is a rate-limiting enzyme in the Krebs cycle that plays a pivotal role in mitochondrial metabolism. OGDHL expression is restricted mainly to the brain in humans. Here, we report nine individuals from eight unrelated families carrying bi-allelic variants in OGDHL with a range of neurological and neurodevelopmental phenotypes including epilepsy, hearing loss, visual impairment, gait ataxia, microcephaly, and hypoplastic corpus callosum. The variants include three homozygous missense variants (p.Pro852Ala, p.Arg244Trp, and p.Arg299Gly), three compound heterozygous single-nucleotide variants (p.Arg673Gln/p.Val488Val, p.Phe734Ser/p.Ala327Val, and p.Trp220Cys/p.Asp491Val), one homozygous frameshift variant (p.Cys553Leufs∗16), and one homozygous stop-gain variant (p.Arg440Ter). To support the pathogenicity of the variants, we developed a novel CRISPR-Cas9-mediated tissue-specific knockout with cDNA rescue system for dOgdh, the Drosophila ortholog of human OGDHL. Pan-neuronal knockout of dOgdh led to developmental lethality as well as defects in Krebs cycle metabolism, which was fully rescued by expression of wild-type dOgdh. Studies using the Drosophila system indicate that p.Arg673Gln, p.Phe734Ser, and p.Arg299Gly are severe loss-of-function alleles, leading to developmental lethality, whereas p.Pro852Ala, p.Ala327Val, p.Trp220Cys, p.Asp491Val, and p.Arg244Trp are hypomorphic alleles, causing behavioral defects. Transcript analysis from fibroblasts obtained from the individual carrying the synonymous variant (c.1464T>C [p.Val488Val]) in family 2 showed that the synonymous variant affects splicing of exon 11 in OGDHL. Human neuronal cells with OGDHL knockout exhibited defects in mitochondrial respiration, indicating the essential role of OGDHL in mitochondrial metabolism in humans. Together, our data establish that the bi-allelic variants in OGDHL are pathogenic, leading to a Mendelian neurodevelopmental disease in humans.

Hearing impairment (HI) is a sensory disorder with a prevalence of 0.0055 live births in South Africa. DNA samples from a South African family presenting with progressive, autosomal dominant non-syndromic HI were subjected to whole-exome sequencing, and a novel monoallelic variant in REST [c.1244GC; p.(C415S)], was identified as the putative causative variant. The co-segregation of the variant was confirmed with Sanger Sequencing. The variant is absent from databases, 103 healthy South African controls, and 52 South African probands with isolated HI. In silico analysis indicates that the p.C415S variant in REST substitutes a conserved cysteine and results in changes to the surrounding secondary structure and the disulphide bonds, culminating in alteration of the tertiary structure of REST. Localization studies using ectopically expressed GFP-tagged Wild type (WT) and mutant REST in HEK-293 cells show that WT REST localizes exclusively to the nucleus; however, the mutant protein localizes throughout the cell. Additionally, mutant REST has an impaired ability to repress its known target AF1q. The data demonstrates that the identified mutation compromises the function of REST and support its implication in HI. This study is the second report, worldwide, to implicate REST in HI and suggests that it should be included in diagnostic HI panels.