Provision of a molecularly confirmed diagnosis in a timely manner for children and adults with rare genetic diseases shortens their "diagnostic odyssey," improves disease management, and fosters genetic counseling with respect to recurrence risks while assuring reproductive choices. In a general clinical genetics setting, the current diagnostic rate is approximately 50%, but for those who do not receive a molecular diagnosis after the initial genetics evaluation, that rate is much lower. Diagnostic success for these more challenging affected individuals depends to a large extent on progress in the discovery of genes associated with, and mechanisms underlying, rare diseases. Thus, continued research is required for moving toward a more complete catalog of disease-related genes and variants. The International Rare Diseases Research Consortium (IRDiRC) was established in 2011 to bring together researchers and organizations invested in rare disease research to develop a means of achieving molecular diagnosis for all rare diseases. Here, we review the current and future bottlenecks to gene discovery and suggest strategies for enabling progress in this regard. Each successful discovery will define potential diagnostic, preventive, and therapeutic opportunities for the corresponding rare disease, enabling precision medicine for this patient population.

We performed whole exome sequence (WES) to identify genetic modifiers on 184 individuals with 22q11.2 deletion syndrome (22q11DS), of whom 89 case subjects had severe congenital heart disease (CHD) and 95 control subjects had normal hearts. Three genes including JMJD1C (jumonji domain containing 1C), RREB1 (Ras responsive element binding protein 1), and SEC24C (SEC24 family member C) had rare (MAF < 0.001) predicted deleterious single-nucleotide variations (rdSNVs) in seven case subjects and no control subjects (p = 0.005; Fisher exact and permutation tests). Because JMJD1C and RREB1 are involved in chromatin modification, we investigated other histone modification genes. Eighteen case subjects (20%) had rdSNVs in four genes (JMJD1C, RREB1, MINA, KDM7A) all involved in demethylation of histones (H3K9, H3K27). Overall, rdSNVs were enriched in histone modifier genes that activate transcription (Fisher exact p = 0.0004, permutations, p = 0.0003, OR = 5.16); however, rdSNVs in control subjects were not enriched. This implicates histone modification genes as influencing risk for CHD in presence of the deletion.

Advances in genomics have transformed our ability to identify the genetic causes of rare diseases (RDs), yet we have a limited understanding of the mechanistic roles of most genes in health and disease. When a novel RD gene is first discovered, there is minimal insight into its biological function, the pathogenic mechanisms of disease-causing variants, and how therapy might be approached. To address this gap, the Canadian Rare Diseases Models and Mechanisms (RDMM) Network was established to connect clinicians discovering new disease genes with Canadian scientists able to study equivalent genes and pathways in model organisms (MOs). The Network is built around a registry of more than 500 Canadian MO scientists, representing expertise for over 7,500 human genes. RDMM uses a committee process to identify and evaluate clinician-MO scientist collaborations and approve 25,000 Canadian dollars in catalyst funding. To date, we have made 85 clinician-MO scientist connections and funded 105 projects. These collaborations help confirm variant pathogenicity and unravel the molecular mechanisms of RD, and also test novel therapies and lead to long-term collaborations. To expand the impact and reach of this model, we made the RDMM Registry open-source, portable, and customizable, and we freely share our committee structures and processes. We are currently working with emerging networks in Europe, Australia, and Japan to link international RDMM networks and registries and enable matches across borders. We will continue to create meaningful collaborations, generate knowledge, and advance RD research locally and globally for the benefit of patients and families living with RD.

Organ fibrosis is a shared endpoint of many diseases, yet underlying mechanisms are not well understood. Several pathways governed by the primary cilium, a sensory antenna present on most vertebrate cells, have been linked with fibrosis. Ciliopathies usually start early in life and represent a considerable disease burden. We performed massively parallel sequencing by using cohorts of genetically unsolved individuals with unexplained liver and kidney failure and correlated this with clinical, imaging, and histopathological analyses. Mechanistic studies were conducted with a vertebrate model and primary cells. We detected bi-allelic deleterious variants in TULP3, encoding a critical adaptor protein for ciliary trafficking, in a total of 15 mostly adult individuals, originating from eight unrelated families, with progressive degenerative liver fibrosis, fibrocystic kidney disease, and hypertrophic cardiomyopathy with atypical fibrotic patterns on histopathology. We recapitulated the human phenotype in adult zebrafish and confirmed disruption of critical ciliary cargo composition in several primary cell lines derived from affected individuals. Further, we show interaction between TULP3 and the nuclear deacetylase SIRT1, with roles in DNA damage repair and fibrosis, and report increased DNA damage ex vivo. Transcriptomic studies demonstrated upregulation of profibrotic pathways with gene clusters for hypertrophic cardiomyopathy and WNT and TGF-β signaling. These findings identify variants in TULP3 as a monogenic cause for progressive degenerative disease of major organs in which affected individuals benefit from early detection and improved clinical management. Elucidation of mechanisms crucial for DNA damage repair and tissue maintenance will guide novel therapeutic avenues for this and similar genetic and non-genomic diseases.

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise healthy individuals can provide insights into genes and biologic pathways involved in production, differentiation, or clearance of particular WBC lineages (myeloid, lymphoid) and also potentially inform the genetic basis of autoimmune, allergic, and blood diseases. We performed an exome array-based meta-analysis of total WBC and subtype counts (neutrophils, monocytes, lymphocytes, basophils, and eosinophils) in a multi-ancestry discovery and replication sample of ∼157,622 individuals from 25 studies. We identified 16 common variants (8 of which were coding variants) associated with one or more WBC traits, the majority of which are pleiotropically associated with autoimmune diseases. Based on functional annotation, these loci included genes encoding surface markers of myeloid, lymphoid, or hematopoietic stem cell differentiation (CD69, CD33, CD87), transcription factors regulating lineage specification during hematopoiesis (ASXL1, IRF8, IKZF1, JMJD1C, ETS2-PSMG1), and molecules involved in neutrophil clearance/apoptosis (C10orf54, LTA), adhesion (TNXB), or centrosome and microtubule structure/function (KIF9, TUBD1). Together with recent reports of somatic ASXL1 mutations among individuals with idiopathic cytopenias or clonal hematopoiesis of undetermined significance, the identification of a common regulatory 3' UTR variant of ASXL1 suggests that both germline and somatic ASXL1 mutations contribute to lower blood counts in otherwise asymptomatic individuals. These association results shed light on genetic mechanisms that regulate circulating WBC counts and suggest a prominent shared genetic architecture with inflammatory and autoimmune diseases.

Neurodevelopmental disorders (NDD) are genetically and phenotypically heterogeneous conditions due to defects in genes involved in development and function of the nervous system. Individuals with NDD, in addition to their primary neurodevelopmental phenotype, may also have accompanying syndromic features that can be very helpful diagnostically especially those with recognizable facial appearance. In this study, we describe ten similarly affected individuals from six unrelated families of different ethnic origins having bi-allelic truncating variants in TMEM94, which encodes for an uncharacterized transmembrane nuclear protein that is highly conserved across mammals. The affected individuals manifested with global developmental delay/intellectual disability, and dysmorphic facial features including triangular face, deep set eyes, broad nasal root and tip and anteverted nostrils, thick arched eye brows, hypertrichosis, pointed chin, and hypertelorism. Birthweight in the upper normal range was observed in most, and all but one had congenital heart defects (CHD). Gene expression analysis in available cells from affected individuals showed reduced expression of TMEM94. Global transcriptome profiling using microarray and RNA sequencing revealed several dysregulated genes essential for cell growth, proliferation and survival that are predicted to have an impact on cardiotoxicity hematological system and neurodevelopment. Loss of Tmem94 in mouse model generated by CRISPR/Cas9 was embryonic lethal and led to craniofacial and cardiac abnormalities and abnormal neuronal migration pattern, suggesting that this gene is important in craniofacial, cardiovascular, and nervous system development. Our study suggests the genetic etiology of a recognizable dysmorphic syndrome with NDD and CHD and highlights the role of TMEM94 in early development.

Rare-variant association studies in common, complex diseases are customarily conducted under an additive risk model in both single-variant and burden testing. Here, we describe a method to improve detection of rare recessive variants in complex diseases termed RAFT (recessive-allele-frequency-based test). We found that RAFT outperforms existing approaches when the variant influences disease risk in a recessive manner on simulated data. We then applied our method to 1,791 Finnish individuals with type 2 diabetes (T2D) and 2,657 matched control subjects. In BBS10, we discovered a rare variant (c.1189A>G [p.Ile397Val]; rs202042386) that confers risk of T2D in a recessive state (p = 1.38 × 10(-6)) and would be missed by conventional methods. Testing of this variant in an established in vivo zebrafish model confirmed the variant to be pathogenic. Taken together, these data suggest that RAFT can effectively reveal rare recessive contributions to complex diseases overlooked by conventional association tests.

A challenge in standard genetic studies is maintaining good power to detect associations, especially for low prevalent diseases and rare variants. The traditional methods are most powerful when evaluating the association between variants in balanced study designs. Without accounting for family correlation and unbalanced case-control ratio, these analyses could result in inflated type I error. One cost-effective solution to increase statistical power is exploitation of available family history (FH) that contains valuable information about disease heritability. Here, we develop methods to address the aforementioned type I error issues while providing optimal power to analyze aggregates of rare variants by incorporating additional information from FH. With enhanced power in these methods exploiting FH and accounting for relatedness and unbalanced designs, we successfully detect genes with suggestive associations with Alzheimer disease, dementia, and type 2 diabetes by using the exome chip data from the Framingham Heart Study.

The study of recent natural selection in human populations has important applications to human history and medicine. Positive natural selection drives the increase in beneficial alleles and plays a role in explaining diversity across human populations. By discovering traits subject to positive selection, we can better understand the population level response to environmental pressures including infectious disease. Our study examines unusual population differentiation between three large data sets to detect natural selection. The populations examined, African Americans, Nigerians, and Gambians, are genetically close to one another (F(ST) < 0.01 for all pairs), allowing us to detect selection even with moderate changes in allele frequency. We also develop a tree-based method to pinpoint the population in which selection occurred, incorporating information across populations. Our genome-wide significant results corroborate loci previously reported to be under selection in Africans including HBB and CD36. At the HLA locus on chromosome 6, results suggest the existence of multiple, independent targets of population-specific selective pressure. In addition, we report a genome-wide significant (p = 1.36 × 10(-11)) signal of selection in the prostate stem cell antigen (PSCA) gene. The most significantly differentiated marker in our analysis, rs2920283, is highly differentiated in both Africa and East Asia and has prior genome-wide significant associations to bladder and gastric cancers.

Exome sequencing has readily enabled the discovery of the genetic mutations responsible for a wide range of diseases. This success has been particularly remarkable in the severe epilepsies and other neurodevelopmental diseases for which rare, often de novo, mutations play a significant role in disease risk. Despite significant progress, the high genetic heterogeneity of these disorders often requires large sample sizes to identify a critical mass of individuals with disease-causing mutations in a single gene. By pooling genetic findings across multiple studies, we have identified six individuals with severe developmental delay (6/6), refractory seizures (5/6), and similar dysmorphic features (3/6), each harboring a de novo mutation in PPP3CA. PPP3CA encodes the alpha isoform of a subunit of calcineurin. Calcineurin encodes a calcium- and calmodulin-dependent serine/threonine protein phosphatase that plays a role in a wide range of biological processes, including being a key regulator of synaptic vesicle recycling at nerve terminals. Five individuals with de novo PPP3CA mutations were identified among 4,760 trio probands with neurodevelopmental diseases; this is highly unlikely to occur by chance (p = 1.2 × 10-8) given the size and mutability of the gene. Additionally, a sixth individual with a de novo mutation in PPP3CA was connected to this study through GeneMatcher. Based on these findings, we securely implicate PPP3CA in early-onset refractory epilepsy and further support the emerging role for synaptic dysregulation in epilepsy.

With advances in whole-genome sequencing (WGS) technology, more advanced statistical methods for testing genetic association with rare variants are being developed. Methods in which variants are grouped for analysis are also known as variant-set, gene-based, and aggregate unit tests. The burden test and sequence kernel association test (SKAT) are two widely used variant-set tests, which were originally developed for samples of unrelated individuals and later have been extended to family data with known pedigree structures. However, computationally efficient and powerful variant-set tests are needed to make analyses tractable in large-scale WGS studies with complex study samples. In this paper, we propose the variant-set mixed model association tests (SMMAT) for continuous and binary traits using the generalized linear mixed model framework. These tests can be applied to large-scale WGS studies involving samples with population structure and relatedness, such as in the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program. SMMATs share the same null model for different variant sets, and a virtue of this null model, which includes covariates only, is that it needs to be fit only once for all tests in each genome-wide analysis. Simulation studies show that all the proposed SMMATs correctly control type I error rates for both continuous and binary traits in the presence of population structure and relatedness. We also illustrate our tests in a real data example of analysis of plasma fibrinogen levels in the TOPMed program (n = 23,763), using the Analysis Commons, a cloud-based computing platform.

Facioscapulohumeral muscular dystrophy type 1 (FSHD1) is caused by contraction of the D4Z4 repeat array on chromosome 4 to a size of 1-10 units. The residual number of D4Z4 units inversely correlates with clinical severity, but significant clinical variability exists. Each unit contains a copy of the DUX4 retrogene. Repeat contractions are associated with changes in D4Z4 chromatin structure that increase the likelihood of DUX4 expression in skeletal muscle, but only when the repeat resides in a genetic background that contains a DUX4 polyadenylation signal. Mutations in the structural maintenance of chromosomes flexible hinge domain containing 1 (SMCHD1) gene, encoding a chromatin modifier of D4Z4, also result in the increased likelihood of DUX4 expression in individuals with a rare form of FSHD (FSHD2). Because SMCHD1 directly binds to D4Z4 and suppresses somatic expression of DUX4, we hypothesized that SMCHD1 may act as a genetic modifier in FSHD1. We describe three unrelated individuals with FSHD1 presenting an unusual high clinical severity based on their upper-sized FSHD1 repeat array of nine units. Each of these individuals also carries a mutation in the SMCHD1 gene. Familial carriers of the FSHD1 allele without the SMCHD1 mutation were only mildly affected, suggesting a modifier effect of the SMCHD1 mutation. Knocking down SMCHD1 in FSHD1 myotubes increased DUX4 expression, lending molecular support to a modifier role for SMCHD1 in FSHD1. We conclude that FSHD1 and FSHD2 share a common pathophysiological pathway in which the FSHD2 gene can act as modifier for disease severity in families affected by FSHD1.

Multiple primary tumors (MPTs) affect a substantial proportion of cancer survivors and can result from various causes, including inherited predisposition. Currently, germline genetic testing of MPT-affected individuals for variants in cancer-predisposition genes (CPGs) is mostly targeted by tumor type. We ascertained pre-assessed MPT individuals (with at least two primary tumors by age 60 years or at least three by 70 years) from genetics centers and performed whole-genome sequencing (WGS) on 460 individuals from 440 families. Despite previous negative genetic assessment and molecular investigations, pathogenic variants in moderate- and high-risk CPGs were detected in 67/440 (15.2%) probands. WGS detected variants that would not be (or were not) detected by targeted resequencing strategies, including low-frequency structural variants (6/440 [1.4%] probands). In most individuals with a germline variant assessed as pathogenic or likely pathogenic (P/LP), at least one of their tumor types was characteristic of variants in the relevant CPG. However, in 29 probands (42.2% of those with a P/LP variant), the tumor phenotype appeared discordant. The frequency of individuals with truncating or splice-site CPG variants and at least one discordant tumor type was significantly higher than in a control population (χ2 = 43.642; p ≤ 0.0001). 2/67 (3%) probands with P/LP variants had evidence of multiple inherited neoplasia allele syndrome (MINAS) with deleterious variants in two CPGs. Together with variant detection rates from a previous series of similarly ascertained MPT-affected individuals, the present results suggest that first-line comprehensive CPG analysis in an MPT cohort referred to clinical genetics services would detect a deleterious variant in about a third of individuals.

Mutations in more than 250 genes are implicated in inherited retinal dystrophy; the encoded proteins are involved in a broad spectrum of pathways. The presence of unsolved families after highly parallel sequencing strategies suggests that further genes remain to be identified. Whole-exome and -genome sequencing studies employed here in large cohorts of affected individuals revealed biallelic mutations in ARHGEF18 in three such individuals. ARHGEF18 encodes ARHGEF18, a guanine nucleotide exchange factor that activates RHOA, a small GTPase protein that is a key component of tight junctions and adherens junctions. This biological pathway is known to be important for retinal development and function, as mutation of CRB1, encoding another component, causes retinal dystrophy. The retinal structure in individuals with ARHGEF18 mutations resembled that seen in subjects with CRB1 mutations. Five mutations were found on six alleles in the three individuals: c.808A>G (p.Thr270Ala), c.1617+5G>A (p.Asp540Glyfs∗63), c.1996C>T (p.Arg666∗), c.2632G>T (p.Glu878∗), and c.2738_2761del (p.Arg913_Glu920del). Functional tests suggest that each disease genotype might retain some ARHGEF18 activity, such that the phenotype described here is not the consequence of nullizygosity. In particular, the p.Thr270Ala missense variant affects a highly conserved residue in the DBL homology domain, which is required for the interaction and activation of RHOA. Previously, knock-out of Arhgef18 in the medaka fish has been shown to cause larval lethality which is preceded by retinal defects that resemble those seen in zebrafish Crumbs complex knock-outs. The findings described here emphasize the peculiar sensitivity of the retina to perturbations of this pathway, which is highlighted as a target for potential therapeutic strategies.

Inherited monogenic disease has an enormous impact on the well-being of children and their families. Over half of the children living with one of these conditions are without a molecular diagnosis because of the rarity of the disease, the marked clinical heterogeneity, and the reality that there are thousands of rare diseases for which causative mutations have yet to be identified. It is in this context that in 2010 a Canadian consortium was formed to rapidly identify mutations causing a wide spectrum of pediatric-onset rare diseases by using whole-exome sequencing. The FORGE (Finding of Rare Disease Genes) Canada Consortium brought together clinicians and scientists from 21 genetics centers and three science and technology innovation centers from across Canada. From nation-wide requests for proposals, 264 disorders were selected for study from the 371 submitted; disease-causing variants (including in 67 genes not previously associated with human disease; 41 of these have been genetically or functionally validated, and 26 are currently under study) were identified for 146 disorders over a 2-year period. Here, we present our experience with four strategies employed for gene discovery and discuss FORGE's impact in a number of realms, from clinical diagnostics to the broadening of the phenotypic spectrum of many diseases to the biological insight gained into both disease states and normal human development. Lastly, on the basis of this experience, we discuss the way forward for rare-disease genetic discovery both in Canada and internationally.

Ciliopathies are a group of hereditary disorders associated with defects in cilia structure and function. The distal appendages (DAPs) of centrioles are involved in the docking and anchoring of the mother centriole to the cellular membrane during ciliogenesis. The molecular composition of DAPs was recently elucidated and mutations in two genes encoding DAPs components (CEP164/NPHP15, SCLT1) have been associated with human ciliopathies, namely nephronophthisis and orofaciodigital syndrome. To identify additional DAP components defective in ciliopathies, we independently performed targeted exon sequencing of 1,221 genes associated with cilia and 5 known DAP protein-encoding genes in 1,255 individuals with a nephronophthisis-related ciliopathy. We thereby detected biallelic mutations in a key component of DAP-encoding gene, CEP83, in seven families. All affected individuals had early-onset nephronophthisis and four out of eight displayed learning disability and/or hydrocephalus. Fibroblasts and tubular renal cells from affected individuals showed an altered DAP composition and ciliary defects. In summary, we have identified mutations in CEP83, another DAP-component-encoding gene, as a cause of infantile nephronophthisis associated with central nervous system abnormalities in half of the individuals.

Using a combination of exome sequencing and linkage analysis, we investigated an English family with two affected siblings in their 40s with recessive Charcot-Marie Tooth disease type 2 (CMT2). Compound heterozygous mutations in the immunoglobulin-helicase-μ-binding protein 2 (IGHMBP2) gene were identified. Further sequencing revealed a total of 11 CMT2 families with recessively inherited IGHMBP2 gene mutations. IGHMBP2 mutations usually lead to spinal muscular atrophy with respiratory distress type 1 (SMARD1), where most infants die before 1 year of age. The individuals with CMT2 described here, have slowly progressive weakness, wasting and sensory loss, with an axonal neuropathy typical of CMT2, but no significant respiratory compromise. Segregating IGHMBP2 mutations in CMT2 were mainly loss-of-function nonsense in the 5' region of the gene in combination with a truncating frameshift, missense, or homozygous frameshift mutations in the last exon. Mutations in CMT2 were predicted to be less aggressive as compared to those in SMARD1, and fibroblast and lymphoblast studies indicate that the IGHMBP2 protein levels are significantly higher in CMT2 than SMARD1, but lower than controls, suggesting that the clinical phenotype differences are related to the IGHMBP2 protein levels.

Cytokines are essential regulatory components of the immune system, and their aberrant levels have been linked to many disease states. Despite increasing evidence that cytokines operate in concert, many of the physiological interactions between cytokines, and the shared genetic architecture that underlies them, remain unknown. Here, we aimed to identify and characterize genetic variants with pleiotropic effects on cytokines. Using three population-based cohorts (n = 9,263), we performed multivariate genome-wide association studies (GWAS) for a correlation network of 11 circulating cytokines, then combined our results in meta-analysis. We identified a total of eight loci significantly associated with the cytokine network, of which two (PDGFRB and ABO) had not been detected previously. In addition, conditional analyses revealed a further four secondary signals at three known cytokine loci. Integration, through the use of Bayesian colocalization analysis, of publicly available GWAS summary statistics with the cytokine network associations revealed shared causal variants between the eight cytokine loci and other traits; in particular, cytokine network variants at the ABO, SERPINE2, and ZFPM2 loci showed pleiotropic effects on the production of immune-related proteins, on metabolic traits such as lipoprotein and lipid levels, on blood-cell-related traits such as platelet count, and on disease traits such as coronary artery disease and type 2 diabetes.

Proteoglycans (PGs) are a major component of the extracellular matrix in many tissues and function as structural and regulatory molecules. PGs are composed of core proteins and glycosaminoglycan (GAG) side chains. The biosynthesis of GAGs starts with the linker region that consists of four sugar residues and is followed by repeating disaccharide units. By exome sequencing, we found that B3GALT6 encoding an enzyme involved in the biosynthesis of the GAG linker region is responsible for a severe skeletal dysplasia, spondyloepimetaphyseal dysplasia with joint laxity type 1 (SEMD-JL1). B3GALT6 loss-of-function mutations were found in individuals with SEMD-JL1 from seven families. In a subsequent candidate gene study based on the phenotypic similarity, we found that B3GALT6 is also responsible for a connective tissue disease, Ehlers-Danlos syndrome (progeroid form). Recessive loss-of-function mutations in B3GALT6 result in a spectrum of disorders affecting a broad range of skeletal and connective tissues characterized by lax skin, muscle hypotonia, joint dislocation, and spinal deformity. The pleiotropic phenotypes of the disorders indicate that B3GALT6 plays a critical role in a wide range of biological processes in various tissues, including skin, bone, cartilage, tendon, and ligament.