Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.

Circulating blood cell counts and indices are important indicators of hematopoietic function and a number of clinical parameters, such as blood oxygen-carrying capacity, inflammation, and hemostasis. By performing whole-exome sequence association analyses of hematologic quantitative traits in 15,459 community-dwelling individuals, followed by in silico replication in up to 52,024 independent samples, we identified two previously undescribed coding variants associated with lower platelet count: a common missense variant in CPS1 (rs1047891, MAF = 0.33, discovery + replication p = 6.38 × 10(-10)) and a rare synonymous variant in GFI1B (rs150813342, MAF = 0.009, discovery + replication p = 1.79 × 10(-27)). By performing CRISPR/Cas9 genome editing in hematopoietic cell lines and follow-up targeted knockdown experiments in primary human hematopoietic stem and progenitor cells, we demonstrate an alternative splicing mechanism by which the GFI1B rs150813342 variant suppresses formation of a GFI1B isoform that preferentially promotes megakaryocyte differentiation and platelet production. These results demonstrate how unbiased studies of natural variation in blood cell traits can provide insight into the regulation of human hematopoiesis.

Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

Genetic variations in promoter sequences that alter gene expression play a prominent role in increasing susceptibility to complex diseases. Also, expression levels of APP are essentially regulated by its core promoter and 5' upstream regulatory region and correlate with amyloid beta levels in Alzheimer disease (AD) brains. Here, we systematically sequenced the proximal promoter (-766/+204) and two functional distal regions (-2634/-2159 and -2096/-1563) of APP in two independent AD series with onset ages < or =70 years (Belgian sample, n=180; Dutch sample, n=111) and identified eight novel sequence variants. Three mutations (-118C-->A, -369C-->G, and -534G-->A) identified only in patients with AD showed, in vitro, a nearly twofold neuron-specific increase in APP transcriptional activity, similar to what is expected from triplication of APP in Down syndrome. These mutations either abolished (AP-2 and HES-1) or created (Oct1) transcription-factor binding sites involved in the development and differentiation of neuronal systems. Also, two of these clustered in the 200-bp region (-540/-340) of the APP promoter that showed the highest degree of species conservation. The present study provides evidence that APP-promoter mutations that significantly increase APP expression levels are associated with AD.

RORα, the RAR-related orphan nuclear receptor alpha, is essential for cerebellar development. The spontaneous mutant mouse staggerer, with an ataxic gait caused by neurodegeneration of cerebellar Purkinje cells, was discovered two decades ago to result from homozygous intragenic Rora deletions. However, RORA mutations were hitherto undocumented in humans. Through a multi-centric collaboration, we identified three copy-number variant deletions (two de novo and one dominantly inherited in three generations), one de novo disrupting duplication, and nine de novo point mutations (three truncating, one canonical splice site, and five missense mutations) involving RORA in 16 individuals from 13 families with variable neurodevelopmental delay and intellectual disability (ID)-associated autistic features, cerebellar ataxia, and epilepsy. Consistent with the human and mouse data, disruption of the D. rerio ortholog, roraa, causes significant reduction in the size of the developing cerebellum. Systematic in vivo complementation studies showed that, whereas wild-type human RORA mRNA could complement the cerebellar pathology, missense variants had two distinct pathogenic mechanisms of either haploinsufficiency or a dominant toxic effect according to their localization in the ligand-binding or DNA-binding domains, respectively. This dichotomous direction of effect is likely relevant to the phenotype in humans: individuals with loss-of-function variants leading to haploinsufficiency show ID with autistic features, while individuals with de novo dominant toxic variants present with ID, ataxia, and cerebellar atrophy. Our combined genetic and functional data highlight the complex mutational landscape at the human RORA locus and suggest that dual mutational effects likely determine phenotypic outcome.

The human RNA helicase DDX6 is an essential component of membrane-less organelles called processing bodies (PBs). PBs are involved in mRNA metabolic processes including translational repression via coordinated storage of mRNAs. Previous studies in human cell lines have implicated altered DDX6 in molecular and cellular dysfunction, but clinical consequences and pathogenesis in humans have yet to be described. Here, we report the identification of five rare de novo missense variants in DDX6 in probands presenting with intellectual disability, developmental delay, and similar dysmorphic features including telecanthus, epicanthus, arched eyebrows, and low-set ears. All five missense variants (p.His372Arg, p.Arg373Gln, p.Cys390Arg, p.Thr391Ile, and p.Thr391Pro) are located in two conserved motifs of the RecA-2 domain of DDX6 involved in RNA binding, helicase activity, and protein-partner binding. We use functional studies to demonstrate that the first variants identified (p.Arg373Gln and p.Cys390Arg) cause significant defects in PB assembly in primary fibroblast and model human cell lines. These variants' interactions with several protein partners were also disrupted in immunoprecipitation assays. Further investigation via complementation assays included the additional variants p.Thr391Ile and p.Thr391Pro, both of which, similarly to p.Arg373Gln and p.Cys390Arg, demonstrated significant defects in P-body assembly. Complementing these molecular findings, modeling of the variants on solved protein structures showed distinct spatial clustering near known protein binding regions. Collectively, our clinical and molecular data describe a neurodevelopmental syndrome associated with pathogenic missense variants in DDX6. Additionally, we suggest DDX6 join the DExD/H-box genes DDX3X and DHX30 in an emerging class of neurodevelopmental disorders involving RNA helicases.

Blood pressure (BP) is a heritable risk factor for cardiovascular disease. To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP), we genotyped ~50,000 SNPs in up to 87,736 individuals of European ancestry and combined these in a meta-analysis. We replicated findings in an independent set of 68,368 individuals of European ancestry. Our analyses identified 11 previously undescribed associations in independent loci containing 31 genes including PDE1A, HLA-DQB1, CDK6, PRKAG2, VCL, H19, NUCB2, RELA, HOXC@ complex, FBN1, and NFAT5 at the Bonferroni-corrected array-wide significance threshold (p < 6 × 10(-7)) and confirmed 27 previously reported associations. Bioinformatic analysis of the 11 loci provided support for a putative role in hypertension of several genes, such as CDK6 and NUCB2. Analysis of potential pharmacological targets in databases of small molecules showed that ten of the genes are predicted to be a target for small molecules. In summary, we identified previously unknown loci associated with BP. Our findings extend our understanding of genes involved in BP regulation, which may provide new targets for therapeutic intervention or drug response stratification.

Pontocerebellar hypoplasia (PCH) is a heterogeneous group of rare recessive disorders with prenatal onset, characterized by hypoplasia of pons and cerebellum. Mutations in a small number of genes have been reported to cause PCH, and the vast majority of PCH cases are explained by mutations in TSEN54, which encodes a subunit of the tRNA splicing endonuclease complex. Here we report three families with homozygous truncating mutations in TBC1D23 who display moderate to severe intellectual disability and microcephaly. MRI data from available affected subjects revealed PCH, small normally proportioned cerebellum, and corpus callosum anomalies. Furthermore, through in utero electroporation, we show that downregulation of TBC1D23 affects cortical neuron positioning. TBC1D23 is a member of the Tre2-Bub2-Cdc16 (TBC) domain-containing RAB-specific GTPase-activating proteins (TBC/RABGAPs). Members of this protein family negatively regulate RAB proteins and modulate the signaling between RABs and other small GTPases, some of which have a crucial role in the trafficking of intracellular vesicles and are involved in neurological disorders. Here, we demonstrate that dense core vesicles and lysosomal trafficking dynamics are affected in fibroblasts harboring TBC1D23 mutation. We propose that mutations in TBC1D23 are responsible for a form of PCH with small, normally proportioned cerebellum and should be screened in individuals with syndromic pontocereballar hypoplasia.

Predictive genomic profiling used to produce personalized nutrition and other lifestyle health recommendations is currently offered directly to consumers. By examining previous meta-analyses and HuGE reviews, we assessed the scientific evidence supporting the purported gene-disease associations for genes included in genomic profiles offered online. We identified seven companies that offer predictive genomic profiling. We searched PubMed for meta-analyses and HuGE reviews of studies of gene-disease associations published from 2000 through June 2007 in which the genotypes of people with a disease were compared with those of a healthy or general-population control group. The seven companies tested at least 69 different polymorphisms in 56 genes. Of the 56 genes tested, 24 (43%) were not reviewed in meta-analyses. For the remaining 32 genes, we found 260 meta-analyses that examined 160 unique polymorphism-disease associations, of which only 60 (38%) were found to be statistically significant. Even the 60 significant associations, which involved 29 different polymorphisms and 28 different diseases, were generally modest, with synthetic odds ratios ranging from 0.54 to 0.88 for protective variants and from 1.04 to 3.2 for risk variants. Furthermore, genes in cardiogenomic profiles were more frequently associated with noncardiovascular diseases than with cardiovascular diseases, and though two of the five genes of the osteogenomic profiles did show significant associations with disease, the associations were not with bone diseases. There is insufficient scientific evidence to conclude that genomic profiles are useful in measuring genetic risk for common diseases or in developing personalized diet and lifestyle recommendations for disease prevention.

The neuromuscular junction (NMJ) is one of the best-studied cholinergic synapses. Inherited defects of peripheral neurotransmission result in congenital myasthenic syndromes (CMSs), a clinically and genetically heterogeneous group of rare diseases with fluctuating fatigable muscle weakness as the clinical hallmark. Whole-exome sequencing and Sanger sequencing in six unrelated families identified compound heterozygous and homozygous mutations in SLC5A7 encoding the presynaptic sodium-dependent high-affinity choline transporter 1 (CHT), which is known to be mutated in one dominant form of distal motor neuronopathy (DHMN7A). We identified 11 recessive mutations in SLC5A7 that were associated with a spectrum of severe muscle weakness ranging from a lethal antenatal form of arthrogryposis and severe hypotonia to a neonatal form of CMS with episodic apnea and a favorable prognosis when well managed at the clinical level. As expected given the critical role of CHT for multisystemic cholinergic neurotransmission, autonomic dysfunctions were reported in the antenatal form and cognitive impairment was noticed in half of the persons with the neonatal form. The missense mutations induced a near complete loss of function of CHT activity in cell models. At the human NMJ, a delay in synaptic maturation and an altered maintenance were observed in the antenatal and neonatal forms, respectively. Increased synaptic expression of butyrylcholinesterase was also observed, exposing the dysfunction of cholinergic metabolism when CHT is deficient in vivo. This work broadens the clinical spectrum of human diseases resulting from reduced CHT activity and highlights the complexity of cholinergic metabolism at the synapse.

Refractive errors are common eye disorders of public health importance worldwide. Ocular axial length (AL) is the major determinant of refraction and thus of myopia and hyperopia. We conducted a meta-analysis of genome-wide association studies for AL, combining 12,531 Europeans and 8,216 Asians. We identified eight genome-wide significant loci for AL (RSPO1, C3orf26, LAMA2, GJD2, ZNRF3, CD55, MIP, and ALPPL2) and confirmed one previously reported AL locus (ZC3H11B). Of the nine loci, five (LAMA2, GJD2, CD55, ALPPL2, and ZC3H11B) were associated with refraction in 18 independent cohorts (n = 23,591). Differential gene expression was observed for these loci in minus-lens-induced myopia mouse experiments and human ocular tissues. Two of the AL genes, RSPO1 and ZNRF3, are involved in Wnt signaling, a pathway playing a major role in the regulation of eyeball size. This study provides evidence of shared genes between AL and refraction, but importantly also suggests that these traits may have unique pathways.

Many disorders are associated with altered serum protein concentrations, including malnutrition, cancer, and cardiovascular, kidney, and inflammatory diseases. Although these protein concentrations are highly heritable, relatively little is known about their underlying genetic determinants. Through transethnic meta-analysis of European-ancestry and Japanese genome-wide association studies, we identified six loci at genome-wide significance (p < 5 × 10(-8)) for serum albumin (HPN-SCN1B, GCKR-FNDC4, SERPINF2-WDR81, TNFRSF11A-ZCCHC2, FRMD5-WDR76, and RPS11-FCGRT, in up to 53,190 European-ancestry and 9,380 Japanese individuals) and three loci for total protein (TNFRS13B, 6q21.3, and ELL2, in up to 25,539 European-ancestry and 10,168 Japanese individuals). We observed little evidence of heterogeneity in allelic effects at these loci between groups of European and Japanese ancestry but obtained substantial improvements in the resolution of fine mapping of potential causal variants by leveraging transethnic differences in the distribution of linkage disequilibrium. We demonstrated a functional role for the most strongly associated serum albumin locus, HPN, for which Hpn knockout mice manifest low plasma albumin concentrations. Other loci associated with serum albumin harbor genes related to ribosome function, protein translation, and proteasomal degradation, whereas those associated with serum total protein include genes related to immune function. Our results highlight the advantages of transethnic meta-analysis for the discovery and fine mapping of complex trait loci and have provided initial insights into the underlying genetic architecture of serum protein concentrations and their association with human disease.

Spinocerebellar ataxias (SCAs) are a heterogeneous group of autosomal-dominant neurodegenerative disorders involving the cerebellum and 23 different genes. We mapped SCA38 to a 56 Mb region on chromosome 6p in a SCA-affected Italian family by whole-genome linkage analysis. Targeted resequencing identified a single missense mutation (c.689G>T [p.Gly230Val]) in ELOVL5. Mutation screening of 456 independent SCA-affected individuals identified the same mutation in two further unrelated Italian families. Haplotyping showed that at least two of the three families shared a common ancestor. One further missense variant (c.214C>G [p.Leu72Val]) was found in a French family. Both missense changes affect conserved amino acids, are predicted to be damaging by multiple bioinformatics tools, and were not identified in ethnically matched controls or within variant databases. ELOVL5 encodes an elongase involved in the synthesis of polyunsaturated fatty acids of the ω3 and ω6 series. Arachidonic acid and docosahexaenoic acid, two final products of the enzyme, were reduced in the serum of affected individuals. Immunohistochemistry on control mice and human brain demonstrated high levels in Purkinje cells. In transfection experiments, subcellular localization of altered ELOVL5 showed a perinuclear distribution with a signal increase in the Golgi compartment, whereas the wild-type showed a widespread signal in the endoplasmic reticulum. SCA38 and SCA34 are examples of SCAs due to mutations in elongase-encoding genes, emphasizing the importance of fatty-acid metabolism in neurological diseases.

Human iris color was one of the first traits for which Mendelian segregation was established. To date, the genetics of iris color is still not fully understood and is of interest, particularly in view of forensic applications. In three independent genome-wide association (GWA) studies of a total of 1406 persons and a genome-wide linkage study of 1292 relatives, all from the Netherlands, we found that the 15q13.1 region is the predominant region involved in human iris color. There were no other regions showing consistent genome-wide evidence for association and linkage to iris color. Single nucleotide polymorphisms (SNPs) in the HERC2 gene and, to a lesser extent, in the neighboring OCA2 gene were independently associated to iris color variation. OCA2 has been implicated in iris color previously. A replication study within two populations confirmed that the HERC2 gene is a new and significant determinant of human iris color variation, in addition to OCA2. Furthermore, HERC2 rs916977 showed a clinal allele distribution across 23 European populations, which was significantly correlated to iris color variation. We suggest that genetic variants regulating expression of the OCA2 gene exist in the HERC2 gene or, alternatively, within the 11.7 kb of sequence between OCA2 and HERC2, and that most iris color variation in Europeans is explained by those two genes. Testing markers in the HERC2-OCA2 region may be useful in forensic applications to predict eye color phenotypes of unknown persons of European genetic origin.

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Intracranial aneurysms (IAs) are acquired cerebrovascular abnormalities characterized by localized dilation and wall thinning in intracranial arteries, possibly leading to subarachnoid hemorrhage and severe outcome in case of rupture. Here, we identified one rare nonsense variant (c.1378A>T) in the last exon of ANGPTL6 (Angiopoietin-Like 6)-which encodes a circulating pro-angiogenic factor mainly secreted from the liver-shared by the four tested affected members of a large pedigree with multiple IA-affected case subjects. We showed a 50% reduction of ANGPTL6 serum concentration in individuals heterozygous for the c.1378A>T allele (p.Lys460Ter) compared to relatives homozygous for the normal allele, probably due to the non-secretion of the truncated protein produced by the c.1378A>T transcripts. Sequencing ANGPTL6 in a series of 94 additional index case subjects with familial IA identified three other rare coding variants in five case subjects. Overall, we detected a significant enrichment (p = 0.023) in rare coding variants within this gene among the 95 index case subjects with familial IA, compared to a reference population of 404 individuals with French ancestry. Among the 6 recruited families, 12 out of 13 (92%) individuals carrying IA also carry such variants in ANGPTL6, versus 15 out of 41 (37%) unaffected ones. We observed a higher rate of individuals with a history of high blood pressure among affected versus healthy individuals carrying ANGPTL6 variants, suggesting that ANGPTL6 could trigger cerebrovascular lesions when combined with other risk factors such as hypertension. Altogether, our results indicate that rare coding variants in ANGPTL6 are causally related to familial forms of IA.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

This study establishes PYROXD1 variants as a cause of early-onset myopathy and uses biospecimens and cell lines, yeast, and zebrafish models to elucidate the fundamental role of PYROXD1 in skeletal muscle. Exome sequencing identified recessive variants in PYROXD1 in nine probands from five families. Affected individuals presented in infancy or childhood with slowly progressive proximal and distal weakness, facial weakness, nasal speech, swallowing difficulties, and normal to moderately elevated creatine kinase. Distinctive histopathology showed abundant internalized nuclei, myofibrillar disorganization, desmin-positive inclusions, and thickened Z-bands. PYROXD1 is a nuclear-cytoplasmic pyridine nucleotide-disulphide reductase (PNDR). PNDRs are flavoproteins (FAD-binding) and catalyze pyridine-nucleotide-dependent (NAD/NADH) reduction of thiol residues in other proteins. Complementation experiments in yeast lacking glutathione reductase glr1 show that human PYROXD1 has reductase activity that is strongly impaired by the disease-associated missense mutations. Immunolocalization studies in human muscle and zebrafish myofibers demonstrate that PYROXD1 localizes to the nucleus and to striated sarcomeric compartments. Zebrafish with ryroxD1 knock-down recapitulate features of PYROXD1 myopathy with sarcomeric disorganization, myofibrillar aggregates, and marked swimming defect. We characterize variants in the oxidoreductase PYROXD1 as a cause of early-onset myopathy with distinctive histopathology and introduce altered redox regulation as a primary cause of congenital muscle disease.