The primary structure of the COOH-terminal half of cholera toxin subunit A1 containing the ADP-ribosylation site

The sequence of 96 amino acid residues from the COOH-terminus of the active subunit of cholera toxin, A₁, has been determined as PheAsnValAsnAspVal LeuGlyAlaTyrAlaProHisProAsxGluGlu GluValSerAlaLeuGlyGly IleProTyrSerGluIleTyrGlyTrpTyrArg ValHisPheGlyValLeuAsp GluGluLeuHisArgGlyTyrArgAspArgTyr TyrSerAsnLeuAspIleAla ProAlaAlaAspGlyTyrGlyLeuAlaGlyPhe ProProGluHisArgAlaTrp ArgGluGluProTrpIleHisHisAlaPro ProGlyCysGlyAsnAlaProArg(OH). This is the largest fragment obtained by BrCN cleavage of the subunit A₁ (M_r 23,000), and has previously been indicated to contain the active site for the adenylate cyclase-stimulating activity. Unequivocal identification of the COOH-terminal structure was achieved by separation and analysis of the terminal peptide after the specific chemical cleavage at the only cysteine residue in A₁ polypeptide. The site of self ADP-ribosylation in the A₁ subunit [C. Y. Lai, Q.-C. Xia, and P. T. Salotra (1983) Biochem. Biophys. Res. Commun.116, 341–348] has now been identified as Arg-50 of this peptide, 46 residues removed from the COOH-terminus. The cysteine that forms disulfide bridge to A₂ subunit in the holotoxin is at position 91.     

Resistance of Acta2R149C/+ mice to aortic disease is associated with defective release of mutant smooth muscle α-actin from the chaperonin-containing TCP1 folding complex

Pathogenic variants of the gene for smooth muscle α-actin (ACTA2), which encodes smooth muscle (SM) α-actin, predispose to heritable thoracic aortic disease. The ACTA2 variant p.Arg149Cys (R149C) is the most common alteration; however, only 60% of carriers have a dissection or undergo repair of an aneurysm by 70 years of age. A mouse model of ACTA2 p.Arg149Cys was generated using CRISPR/Cas9 technology to determine the etiology of reduced penetrance. Acta2^R149C/+ mice had significantly decreased aortic contraction compared with WT mice but did not form aortic aneurysms or dissections when followed to 24 months, even when hypertension was induced. In vitro motility assays found decreased interaction of mutant SM α-actin filaments with SM myosin. Polymerization studies using total internal reflection fluorescence microscopy showed enhanced nucleation of mutant SM α-actin by formin, which correlated with disorganized and reduced SM α-actin filaments in Acta2^R149C/+ smooth muscle cells (SMCs). However, the most prominent molecular defect was the increased retention of mutant SM α-actin in the chaperonin-containing t-complex polypeptide folding complex, which was associated with reduced levels of mutant compared with WT SM α-actin in Acta2^R149C/+ SMCs. These data indicate that Acta2^R149C/+ mice do not develop thoracic aortic disease despite decreased contraction of aortic segments and disrupted SM α-actin filament formation and function in Acta2^R149C/+ SMCs. Enhanced binding of mutant SM α-actin to chaperonin-containing t-complex polypeptide decreases the mutant actin versus WT monomer levels in Acta2^R149C/+ SMCs, thus minimizing the effect of the mutation on SMC function and potentially preventing aortic disease in the Acta2^R149C/+ mice.     

Amino acid sequence of seminalplasmin, an antimicrobial protein from bull semen

Analytical ultracentrifugation of highly purified seminalplasmin revealed a molecular mass of 6300. Amino acid analysis of the protein preparation indicated the absence of sulfur-containing amino acids cysteine and methionine. The amino acid sequence of seminalplasmin was determined by manual Edman degradation of peptides obtained by proteolytic enzymes trypsin, chymotrypsin and thermolysin: NH₂-Ser Asp Glu Lys Ala Ser Pro Asp Lys His His Arg Phe Ser Leu Ser Arg Tyr Ala Lys Leu Ala Asn Arg Leu Ser Lys Trp Ile Gly Asn Arg Gly Asn Arg Leu Ala Asn Pro Lys Leu Leu Glu Thr Phe Lys Ser Val-COOH. The number of amino acids according to the sequence were 48, the molecular mass 6385. As predicted from the sequence, seminalplasmin very likely contains two α-helical domains in which residues 8-17 and 40-48 are involved. No evidence for the existence of β-sheet structures was obtained. Treatment of seminalplasmin with the above proteases as well as with amino peptidase M and carboxypeptidase Y completely eliminated biological activity.     

<Emphasis Type="Italic">CDKN2A</Emphasis> and <Emphasis Type="Italic">MC1R</Emphasis> variants found in Cypriot patients diagnosed with cutaneous melanoma

The prevalence of genetic variants associated to cutaneous melanoma (CM) has never been determined within Cypriot melanomas. This study evaluates the frequency of variants in cyclin-dependent kinase inhibitor 2A (CDKN2A) and melanocortin-1 receptor (MC1R) in 32 patients diagnosed with CM. Other characteristics and risk factors were also assessed. CDKN2A p.Ala148Thr was detected in three of 32 patients, while the control group revealed no variations within CDKN2A. MC1R screening in 32 patients revealed the following variations: p.Val60Leu in 11 patients, p.Arg142His in four patients, p.Thr314Thr in one patient, p.Arg160Trp in one patient, p.Val92Met/p.Thr314Thr in one patient and p.Val92Met/p.Arg142His/p.Thr314Thr in one patient. The control group revealed only p.Val60Leu (in 10 of 45 individuals), which is frequently found in general populations. Two unrelated patients carried CDKN2A p.Ala148Thr in combination with MC1R p.Arg142His, suggesting digenic inheritance that may provide evidence of different gene variants acting synergistically to contribute for CM development. This study confirms the presence of CDKN2A and MC1R variants among Cypriot melanomas and supports existing evidence of a role for these variants in susceptibility to melanoma.     

The alpha4 residues of human DNA topoisomerase IIalpha function in enzymatic activity and anticancer drug sensitivity

We introduced a series of Pro substitutions within and near the α4 helix, a part of the breakage/rejoining region, in human DNA topoisomerase IIα, and analyzed if this region is involved in determination of anti-cancer drug sensitivity in a temperature- sensitive yeast strain (top2-4 allele). Among the 19 mutants generated, H759P and N770P showed resistance to etoposide and doxorubicin at the non-permissive temperature, where cell growth depends on activity of the human enzyme. For these residues, mutants with an Ala substitution were further created, in which H759A also showed resistance to etoposide. H759P, H759A and N770P were expressed, purified and subjected to in vitro measurement of drug sensitivity. They generated lower amounts of the etoposide-induced cleavable complexes, and were also found to have lower decatenation activity than the wild-type. In the crystal structure, the yeast equivalent of His759 is found in the vicinity of the Arg713, a putative anchoring residue of the 3′-side of cleaved DNA strands. These results suggest that His759 and the other α4 helix residues are involved in the enzymatic activity and drug sensitivity of human DNA topoisomerase IIα, via interaction with cleaved DNA.     

Association between Rare Variants in AP4E1, a Component of Intracellular Trafficking,and Persistent Stuttering

Stuttering is a common, highly heritable neurodevelopmental disorder characterized by deficits in the volitional control of speech. Whole-exome sequencing identified two heterozygous AP4E1 coding variants, c.1549G>A (p.Val517Ile) and c.2401G>A (p.Glu801Lys), that co-segregate with persistent developmental stuttering in a large Cameroonian family, and we observed the same two variants in unrelated Cameroonians with persistent stuttering. We found 23 other rare variants, including predicted loss-of-function variants, in AP4E1 in unrelated stuttering individuals in Cameroon, Pakistan, and North America. The rate of rare variants in AP4E1 was significantly higher in unrelated Pakistani and Cameroonian stuttering individuals than in population-matched control individuals, and coding variants in this gene are exceptionally rare in the general sub-Saharan West African, South Asian, and North American populations. Clinical examination of the Cameroonian family members failed to identify any symptoms previously reported in rare individuals carrying homozygous loss-of-function mutations in this gene. AP4E1 encodes the ε subunit of the heterotetrameric (ε-β4-μ4-σ4) AP-4 complex, involved in protein sorting at the trans-Golgi network. We found that the μ4 subunit of AP-4 interacts with NAGPA, an enzyme involved in the synthesis of the mannose 6-phosphate signal that targets acid hydrolases to the lysosome and the product of a gene previously associated with stuttering. These findings implicate deficits in intracellular trafficking in persistent stuttering.     

Loss-of-Function Mutations in the WNT Co-receptor LRP6 Cause Autosomal-Dominant Oligodontia

Tooth agenesis is one of the most common developmental anomalies in man. Oligodontia, a severe form of tooth agenesis, occurs both as an isolated anomaly and as a syndromal feature. We performed exome sequencing on 20 unrelated individuals with apparent non-syndromic oligodontia and failed to detect mutations in genes previously associated with oligodontia. In three of the probands, we detected heterozygous variants in LRP6, and sequencing of additional oligodontia-affected individuals yielded one additional mutation in LRP6. Three mutations (c.1144_1145dupAG [p.Ala383Glyfs^∗8], c.1779dupT [p.Glu594^∗], and c.2224_2225dupTT [p.Leu742Phefs^∗7]) are predicted to truncate the protein, whereas the fourth (c.56C>T [p.Ala19Val]) is a missense variant of a conserved residue located at the cleavage site of the protein’s signal peptide. All four affected individuals harboring a LRP6 mutation had a family history of tooth agenesis. LRP6 encodes a transmembrane cell-surface protein that functions as a co-receptor with members from the Frizzled protein family in the canonical Wnt/β-catenin signaling cascade. In this same pathway, WNT10A was recently identified as a major contributor in the etiology of non-syndromic oligodontia. We show that the LRP6 missense variant (c.56C>T) results in altered glycosylation and improper subcellular localization of the protein, resulting in abrogated activation of the Wnt pathway. Our results identify LRP6 variants as contributing to the etiology of non-syndromic autosomal-dominant oligodontia and suggest that this gene is a candidate for screening in DNA diagnostics.     

Identification by Mn2+ rescue of two residues essential for the proton transfer of tRNase Z catalysis

Thermotoga maritima tRNase Z cleaves pre-tRNAs containing the ₇₄CCA₇₆ sequence precisely after the A₇₆ residue to create the mature 3′ termini. Its crystal structure has revealed a four-layer αβ/βα sandwich fold that is typically found in the metallo-β-lactamase superfamily. The well-conserved six histidine and two aspartate residues together with metal ions are assumed to form the tRNase Z catalytic center. Here, we examined tRNase Z variants containing single amino acid substitutions in the catalytic center for pre-tRNA cleavage. Cleavage by each variant in the presence of Mg²⁺ was hardly detected, although it is bound to pre-tRNA. Surprisingly, however, Mn²⁺ ions restored the lost Mg²⁺-dependent activity with two exceptions of the Asp52Ala and His222Ala substitutions, which abolished the activity almost completely. These results provide a piece of evidence that Asp-52 and His-222 directly contribute the proton transfer for the catalysis.     

Mutations in DDX58, which Encodes RIG-I,Cause Atypical Singleton-Merten Syndrome

Singleton-Merten syndrome (SMS) is an autosomal-dominant multi-system disorder characterized by dental dysplasia, aortic calcification, skeletal abnormalities, glaucoma, psoriasis, and other conditions. Despite an apparent autosomal-dominant pattern of inheritance, the genetic background of SMS and information about its phenotypic heterogeneity remain unknown. Recently, we found a family affected by glaucoma, aortic calcification, and skeletal abnormalities. Unlike subjects with classic SMS, affected individuals showed normal dentition, suggesting atypical SMS. To identify genetic causes of the disease, we performed exome sequencing in this family and identified a variant (c.1118A>C [p.Glu373Ala]) of DDX58, whose protein product is also known as RIG-I. Further analysis of DDX58 in 100 individuals with congenital glaucoma identified another variant (c.803G>T [p.Cys268Phe]) in a family who harbored neither dental anomalies nor aortic calcification but who suffered from glaucoma and skeletal abnormalities. Cys268 and Glu373 residues of DDX58 belong to ATP-binding motifs I and II, respectively, and these residues are predicted to be located closer to the ADP and RNA molecules than other nonpathogenic missense variants by protein structure analysis. Functional assays revealed that DDX58 alterations confer constitutive activation and thus lead to increased interferon (IFN) activity and IFN-stimulated gene expression. In addition, when we transduced primary human trabecular meshwork cells with c.803G>T (p.Cys268Phe) and c.1118A>C (p.Glu373Ala) mutants, cytopathic effects and a significant decrease in cell number were observed. Taken together, our results demonstrate that DDX58 mutations cause atypical SMS manifesting with variable expression of glaucoma, aortic calcification, and skeletal abnormalities without dental anomalies.     

The Domain-Specific and Temperature-Dependent Protein Misfolding Phenotype of Variant Medium-Chain acyl-CoA Dehydrogenase

The implementation of expanded newborn screening programs reduced mortality and morbidity in medium-chain acyl-CoA dehydrogenase deficiency (MCADD) caused by mutations in the ACADM gene. However, the disease is still potentially fatal. Missense induced MCADD is a protein misfolding disease with a molecular loss-of-function phenotype. Here we established a comprehensive experimental setup to analyze the structural consequences of eight ACADM missense mutations (p.Ala52Val, p.Tyr67His, p.Tyr158His, p.Arg206Cys, p.Asp266Gly, p.Lys329Glu, p.Arg334Lys, p.Arg413Ser) identified after newborn screening and linked the corresponding protein misfolding phenotype to the site of side-chain replacement with respect to the domain. With fever being the crucial risk factor for metabolic decompensation of patients with MCADD, special emphasis was put on the analysis of structural and functional derangements related to thermal stress. Based on protein conformation, thermal stability and kinetic stability, the molecular phenotype in MCADD depends on the structural region that is affected by missense-induced conformational changes with the central β-domain being particularly prone to structural derangement and destabilization. Since systematic classification of conformational derangements induced by ACADM mutations may be a helpful tool in assessing the clinical risk of patients, we scored the misfolding phenotype of the variants in comparison to p.Lys329Glu (K304E), the classical severe mutation, and p.Tyr67His (Y42H), discussed to be mild. Experiments assessing the impact of thermal stress revealed that mutations in the ACADM gene lower the temperature threshold at which MCAD loss-of-function occurs. Consequently, increased temperature as it occurs during intercurrent infections, significantly increases the risk of further conformational derangement and loss of function of the MCAD enzyme explaining the life-threatening clinical courses observed during fever episodes. Early and aggressive antipyretic treatment thus may be life-saving in patients suffering from MCADD.     

X-Linked Adrenoleukodystrophy: Molecular and Functional Analysis of the ABCD1 Gene in Argentinean Patients

X-linked adrenoleukodystrophy (X-ALD) is an inherited metabolic disease associated with mutations in the ABCD1 gene that encodes an ATP-binding cassette transporter protein, ALDP. The disease is characterized by increased concentrations of very long-chain fatty acids (VLCFAs) in plasma and in adrenal, testicular and nervous tissues, due to a defect in peroxisomal VLCFA β-oxidation. In the present study, we analyzed 10 male patients and 17 female carriers from 10 unrelated pedigrees with X-ALD from Argentina. By sequencing the ABCD1 we detected 9 different mutations, 8 of which were novel. These new mutations were verified by a combination of methods that included both functional (western blot and peroxisomal VLCFA β-oxidation) and bioinformatics analysis. The spectrum of novel mutations consists of 3 frameshift (p.Ser284fs*16, p.Glu380Argfs*21 and p.Thr254Argfs*82); a deletion (p.Ser572_Asp575del); a splicing mutation (c.1081+5G>C) and 3 missense mutations (p.Ala341Asp, p.His420Pro and p.Tyr547Cys). In one patient 2 changes were found: a known missense (p.His669Arg) and an unpublished amino acid substitution (p.Ala19Ser). In vitro studies suggest that p.Ala19Ser is a polymorphism. Moreover, we identified two novel intronic polymorphisms and two amino acid polymorphisms. In conclusion, this study extends the spectrum of mutation in X-ALD and facilitates the identification of heterozygous females.     

Variants in the ASB10 Gene Are Associated with Primary Open Angle Glaucoma

Background

Recently nonsynonymous coding variants in the ankyrin repeats and suppressor of cytokine signaling box-containing protein 10 (ASB10) gene were found to be associated with primary open angle glaucoma (POAG) in cohorts from Oregon and Germany, but this finding was not confirmed in an independent cohort from Iowa. The aim of the current study was to assess the role of ASB10 gene variants in Pakistani glaucoma patients.

Methods

Sanger sequencing of the coding exons and splice junctions of the ASB10 gene was performed in 30 probands of multiplex POAG families, 208 sporadic POAG patients and 151 healthy controls from Pakistan. Genotypic associations of individual variants with POAG were analyzed with the Fisher’s exact or Chi-square test.

Results

In total 24 variants were identified in POAG probands and sporadic patients, including 11 novel variants and 13 known variants. 13 of the variants were nonsynonymous, 6 were synonymous, and 5 were intronic. Three nonsynonymous variants (p.Arg49Cys, p.Arg237Gly, p.Arg453Cys) identified in the probands were not segregating in the respective families. This is not surprising since glaucoma is a multifactorial disease, and multiple factors are likely to be involved in the disease manifestation in these families. However a nonsynonymous variant, p.Arg453Cys (rs3800791), was found in 6 sporadic POAG patients but not in controls, suggesting that it infers increased risk for the disease. In addition, one synonymous variant was found to be associated with sporadic POAG: p.Ala290Ala and the association of the variant with POAG remained significant after correction for multiple testing (uncorrected p-value 0.002, corrected p-value 0.047). The cumulative burden of rare, nonsynonymous variants was significantly higher in sporadic POAG patients compared to control individuals (p-value 0.000006).

Conclusions

Variants in ASB10 were found to be significantly associated with sporadic POAG in the Pakistani population. This supports previous findings that sequence variants in the ASB10 gene may act as a risk factor for glaucoma.     

Structural analysis of new compound heterozygous variants in PEPD gene identified in a patient with Prolidase Deficiency diagnosed by exome sequencing

Prolidase Deficiency (PD) is an autosomal recessive rare disorder caused by loss or reduction of prolidase enzymatic activity due to variants in the PEPD gene. PD clinical features vary among affected individuals: skin ulcerations, recurrent infections, and developmental delay are common. In this study, we describe a 16-year-old boy with a mild PD phenotype comprising chronic eczema, recurrent infections and elevated IgE. Whole exome sequencing analysis revealed three PEPD variants: c.575T>C p.(Leu192Pro) inherited from the mother, and c.692_694del p.(Tyr231del) and c.1409G>A p.(Arg470His), both inherited from the father. The variant p.(Tyr231del) has been previously characterized by high-resolution X-ray structure analysis as altering protein dynamics/flexibility. In order to study the effects of the other two prolidase variants, we performed site directed mutagenesis purification and crystallization studies. A high-resolution X-ray structure could only be obtained for the p.(Arg470His) variant, which showed no significant structural differences in comparison to WT prolidase. On the other hand, the p.(Leu192Pro) variant led to significant protein destabilization. Hence, we conclude that the maternal p.(Leu192Pro) variant was likely causally associated with the proband´s disease, together with the known pathogenic paternal variant p.(Tyr231del). Our results demonstrated the utility of exome sequencing to perform diagnosis in PD cases with mild phenotype.     

A Single Lys Residue on the First Intracellular Loop Modulates the Endoplasmic Reticulum Export and Cell-Surface Expression of α2A-Adrenergic Receptor

Export from the endoplasmic reticulum (ER) represents an initial step in intracellular trafficking of G protein-coupled receptors (GPCRs). However, the underlying molecular mechanisms remain poorly understood. We have previously demonstrated that a highly conserved Leu residue on the first intracellular loop (ICL1) is required for exit of several GPCRs from the ER. Here we found that, in addition to Leu64 residue in the ICL1, the neighboring positively charged residue Lys65also modulates the cell-surface transport of α_2A-adrenergic receptor (α_2A-AR). Mutation of Lys65 to Ala, Glu and Gln significantly attenuated, whereas mutation of Lys65 to Arg strongly augmented α_2A-AR expression at the cell surface. Consistent with the effects on the cell-surface expression of α_2A-AR, mutation of Lys65 to Ala and Arg produced opposing effects on α_2A-AR-mediated ERK1/2 activation. Furthermore, confocal microscopy revealed that the α_2A-AR mutant K65A displayed a strong intracellular expression pattern and was extensively co-localized with the ER marker DsRed2-ER, suggestive of ER accumulation. These data provide the first evidence indicating an important function for a single Lys residue on the ICL1 in the ER export and cell-surface expression of α_2A-AR. These data also suggest that the ICL1 may possess multiple signals that control the cell-surface targeting of GPCRs via distinct mechanisms.     

Defective Initiation of Glycosaminoglycan Synthesis due to B3GALT6 Mutations Causes a Pleiotropic Ehlers-Danlos-Syndrome-like Connective Tissue Disorder

Proteoglycans are important components of cell plasma membranes and extracellular matrices of connective tissues. They consist of glycosaminoglycan chains attached to a core protein via a tetrasaccharide linkage, whereby the addition of the third residue is catalyzed by galactosyltransferase II (β3GalT6), encoded by B3GALT6. Homozygosity mapping and candidate gene sequence analysis in three independent families, presenting a severe autosomal-recessive connective tissue disorder characterized by skin fragility, delayed wound healing, joint hyperlaxity and contractures, muscle hypotonia, intellectual disability, and a spondyloepimetaphyseal dysplasia with bone fragility and severe kyphoscoliosis, identified biallelic B3GALT6 mutations, including homozygous missense mutations in family 1 (c.619G>C [p.Asp207His]) and family 3 (c.649G>A [p.Gly217Ser]) and compound heterozygous mutations in family 2 (c.323_344del [p.Ala108Glyfs^∗163], c.619G>C [p.Asp207His]). The phenotype overlaps with several recessive Ehlers-Danlos variants and spondyloepimetaphyseal dysplasia with joint hyperlaxity. Affected individuals’ fibroblasts exhibited a large decrease in ability to prime glycosaminoglycan synthesis together with impaired glycanation of the small chondroitin/dermatan sulfate proteoglycan decorin, confirming β3GalT6 loss of function. Dermal electron microcopy disclosed abnormalities in collagen fibril organization, in line with the important regulatory role of decorin in this process. A strong reduction in heparan sulfate level was also observed, indicating that β3GalT6 deficiency alters synthesis of both main types of glycosaminoglycans. In vitro wound healing assay revealed a significant delay in fibroblasts from two index individuals, pointing to a role for glycosaminoglycan defect in impaired wound repair in vivo. Our study emphasizes a crucial role for β3GalT6 in multiple major developmental and pathophysiological processes.     

Somatic Activating PIK3CA Mutations Cause Venous Malformation

Somatic mutations in TEK, the gene encoding endothelial cell tyrosine kinase receptor TIE2, cause more than half of sporadically occurring unifocal venous malformations (VMs). Here, we report that somatic mutations in PIK3CA, the gene encoding the catalytic p110α subunit of PI3K, cause 54% (27 out of 50) of VMs with no detected TEK mutation. The hotspot mutations c.1624G>A, c.1633G>A, and c.3140A>G (p.Glu542Lys, p.Glu545Lys, and p.His1047Arg), frequent in PIK3CA-associated cancers, overgrowth syndromes, and lymphatic malformation (LM), account for >92% of individuals who carry mutations. Like VM-causative mutations in TEK, the PIK3CA mutations cause chronic activation of AKT, dysregulation of certain important angiogenic factors, and abnormal endothelial cell morphology when expressed in human umbilical vein endothelial cells (HUVECs). The p110α-specific inhibitor BYL719 restores all abnormal phenotypes tested, in PIK3CA- as well as TEK-mutant HUVECs, demonstrating that they operate via the same pathogenic pathways. Nevertheless, significant genotype-phenotype correlations in lesion localization and histology are observed between individuals with mutations in PIK3CA versus TEK, pointing to gene-specific effects.     

Biallelic SUN5 Mutations Cause Autosomal-Recessive Acephalic Spermatozoa Syndrome

Acephalic spermatozoa syndrome is a rare and severe form of teratozoospermia characterized by a predominance of headless spermatozoa in the ejaculate. Family clustering and consanguinity suggest a genetic origin; however, causative mutations have yet to be identified. We performed whole-exome sequencing in two unrelated infertile men and subsequent variant filtering identified one homozygous (c.824C>T [p.Thr275Met]) and one compound heterozygous (c.1006C>T [p.Arg356Cys] and c.485T>A [p.Met162Lys]) SUN5 (also named TSARG4) variants. Sanger sequencing of SUN5 in 15 additional unrelated infertile men revealed four compound heterozygous (c.381delA [p.Val128Serfs^∗7] and c.824C>T [p.Thr275Met]; c.381delA [p.Val128Serfs^∗7] and c.781G>A [p.Val261Met]; c.216G>A [p.Trp72^∗] and c.1043A>T [p.Asn348Ile]; c.425+1G>A/c.1043A>T [p.Asn348Ile]) and two homozygous (c.851C>G [p.Ser284^∗]; c.350G>A [p.Gly114Arg]) variants in six individuals. These 10 SUN5 variants were found in 8 of 17 unrelated men, explaining the genetic defect in 47.06% of the affected individuals in our cohort. These variants were absent in 100 fertile population-matched control individuals. SUN5 variants lead to absent, significantly reduced, or truncated SUN5, and certain variants altered SUN5 distribution in the head-tail junction of the sperm. In summary, these results demonstrate that biallelic SUN5 mutations cause male infertility due to autosomal-recessive acephalic spermatozoa syndrome.