Open Access

Identification of embryonic lethal genes in humans by autozygosity mapping and exome sequencing in consanguineous families

  • Hanan E. Shamseldin1,
  • Maha Tulbah2,
  • Wesam Kurdi2,
  • Maha Nemer2,
  • Nada Alsahan2,
  • Elham Al Mardawi3,
  • Ola Khalifa4, 8,
  • Amal Hashem5,
  • Ahmed Kurdi6,
  • Zainab Babay7,
  • Dalal K. Bubshait4, 9,
  • Niema Ibrahim1,
  • Firdous Abdulwahab1,
  • Zuhair Rahbeeni4,
  • Mais Hashem1 and
  • Fowzan S. Alkuraya1, 10Email author
Genome Biology201516:116

DOI: 10.1186/s13059-015-0681-6

Received: 4 February 2015

Accepted: 21 May 2015

Published: 3 June 2015

Abstract

Background

Identifying genetic variants that lead to discernible phenotypes is the core of Mendelian genetics. An approach that considers embryonic lethality as a bona fide Mendelian phenotype has the potential to reveal novel genetic causes, which will further our understanding of early human development at a molecular level. Consanguineous families in which embryonic lethality segregates as a recessive Mendelian phenotype offer a unique opportunity for high throughput novel gene discovery as has been established for other recessive postnatal phenotypes.

Results

We have studied 24 eligible families using autozygosity mapping and whole-exome sequencing. In addition to revealing mutations in genes previously linked to embryonic lethality in severe cases, our approach revealed seven novel candidate genes (THSD1, PIGC, UBN1, MYOM1, DNAH14, GALNT14, and FZD6). A founder mutation in one of these genes, THSD1, which has been linked to vascular permeability, accounted for embryonic lethality in three of the study families. Unlike the other six candidate genes, we were able to identify a second mutation in THSD1 in a family with a less severe phenotype consisting of hydrops fetalis and persistent postnatal edema, which provides further support for the proposed link between this gene and embryonic lethality.

Conclusions

Our study represents an important step towards the systematic analysis of “embryonic lethal genes” in humans.

Background

The in utero development of humans from a single-celled zygote to a newborn with trillions of remarkably diverse cells with three-dimensional organization as organs is a highly complex process that remains incompletely understood. Much of our understanding of the molecular control of human development comes from the discipline of developmental genetics. The use of animal models as surrogates to understand these molecular events by virtue of forward and reverse genetics approaches that deconvolute these complex events to the level of individual genes has been helpful but does not obviate the need to study human embryos.

In humans, the ability to trace certain birth defects to individual genetic lesions has served as a “shortcut” to accelerate the discovery of the genetic control of many aspects of normal development. The premise of this approach is that genes linked to birth defects must be indispensable to the normal development of the respective organ(s), especially in the Mendelian context where the link is directly causal. In fact, it can be argued that the study of Mendelian forms of birth defects has provided the most compelling list of genes that are required for normal human development [1].

Embryonic lethality is a well-studied phenomenon in many model organisms, and the proportion of genes that are indispensable to development is remarkable. For example, 19 % of yeast genes are lethal when knocked out [2]. A similar systematic survey has been more challenging to conduct in higher organisms but available data suggest a lethal knockout phenotype for 30 % of mouse genes in the homozygous state [3]. Since systematic knockout of each human gene to test its viability cannot be done experimentally, naturally occurring mutations that effectively knock out the respective genes provide a practical alternative. Building on the successful approach of using Mendelian genetics of birth defects to understand normal morphogenesis of individual organs, we reckoned that Mendelian forms of embryonic lethality offer a window into the essential genetic components of early organismal development in humans. In this study, we report our genomic analysis of multiplex consanguineous pedigrees in which embryonic lethality appears to follow a Mendelian recessive pattern. This approach revealed known disease genes that have been reported to cause embryonic lethality in severe cases, in addition to genes that have not been previously linked to human diseases, representing novel candidates. The candidacy of one of these novel genes is further supported by positional mapping data as well as the identification of two different mutations.

Results

Identification of Mendelian forms of embryonic lethality

We have identified 24 consanguineous families in which two or more pregnancies were diagnosed with lethal non-immune hydrops fetalis (NIHF; Figure S1 in Additional file 1). High quality DNA was available from at least one affected fetus in 19 of the 24 families and these were subjected to whole-exome sequencing (WES). Additional file 2 summarizes the raw data characteristics from these 50× coverage exomes. This was followed by gene discovery using the algorithm outlined in Fig. 1 and detailed in Additional file 3. In seven of these families, WES revealed a likely pathogenic variant in a gene known to cause NIHF in severe cases (Table 1). Family 13DG1635 had a homozygous truncating mutation in CTSA, which predicts galactosialidosis, a known cause of NIHF [46]. Similarly, families 14DG1727, 14DG1819, 13DG0042, 13DG2155, and 14DG0947 had likely pathogenic mutations in GUSB as the only variants that remained after applying the various filters. These mutations predict GUSB-related mucopolysaccharidosis VII, another established cause of NIHF in severe cases [7, 8]. Family 13DG0975 was found to harbor a homozygous stopgain mutation in NEB (c.21076C > T:p.R7026X), which predicts a severe form of nemaline myopathy, which has rarely been reported to cause NIHF [9]. Family 10DG0827 was also found to harbor a homozygous mutation in the muscle-related gene CHRNA1 (c.762C > T:p.R254C) and we reported this family as the first application of WES to identify the cause of recurrent fetal loss in humans [10]. The candidate causal variants in the remaining families, on the other hand, involve genes that have not been previously linked to human phenotypes, thus highlighting them as potential candidates for embryonic lethality.
Fig. 1

Workflow of the study

Table 1

Summary of the variants identified by WES of families with recurrent fetal loss

Patient ID

Variant remaining after WES filtering

Effect of the mutation

Relevance to embryonic lethality

14DG1727

GUSB:NM_000181:exon2:c.307C > T:p.R103W

Replaces a highly conserved amino acid (PolyPhen score 1; SIFT score 0)

Known cause of NIHF in severe cases

14DG1819

GUSB:NM_000181:exon10:c.1586A > G:p.Y529C

Replaces a highly conserved amino acid (PolyPhen score 1; SIFT score 0)

Known cause of NIHF in severe cases

10DG0827

CHRNA1:NM_001039523:exon7:c.762C > T:p.R254C

Replaces a highly conserved amino acid (PolyPhen score 1; SIFT score 0)

Known cause of NIHF in severe cases

09DG01201

PIGC:NM_002642:exon2:c.659 T > C:p.L220P

Replaces a highly conserved amino acid (PolyPhen score 0.998; SIFT score 0.01)

Novel candidate

12DG2262

UBN1:NM_016936:exon14:c.2356 T > A:p.L786M

Replaces a highly conserved amino acid (PolyPhen score 0.915; SIFT score 0.02)

Novel candidate

13DG0042

GUSB:NM_000181:exon7:c.1144C > T:p.R382C

Replaces a highly conserved amino acid (PolyPhen score 0.995; SIFT score 0)

Known cause of NIHF in severe cases

13DG0259

None

-

-

13DG0447

DNAH14:NM_001373:exon23:c.3755 T > A:p.M1252K

Replaces a highly conserved amino acid (PolyPhen score 0.805; SIFT score 0)

Novel candidate

13DG0806

THSD1:NM_018676:exon3:c.617G > A:p.C206Y

Replaces a highly conserved amino acid (PolyPhen score 0.999; SIFT score 0)

Novel candidate

13DG0975

NEB:NM_001164507:exon140:c.21076C > T:p.R7026X

Truncation of >50 % of the protein

Known cause of NIHF in severe cases

13DG1635

CTSA:NM_000308:exon6:c.649delC:p.L217fs

Truncation of >75 % of the protein sequence

Known cause of NIHF in severe cases

13DG1885

MYOM1:NM_003803:exon38:c.4987G > A:p.V1663M

Replaces a highly conserved amino acid (PolyPhen score 0.999; SIFT score 0)

Novel candidate

13DG2155

GUSB:NM_000181:exon3:c.398G > C:p.W133S,

Replaces a highly conserved amino acid (PolyPhen score 0.984; SIFT score 0)

Known cause of NIHF in severe cases

14DG0052

FZD6:NM_001164616:exon5:c.773A > G:p.Y258C

Replaces a highly conserved amino acid (PolyPhen score 1; SIFT score 0)

Novel candidate

14DG0946

THSD1:NM_018676:exon3:c.617G > A:p.C206Y

Replaces a highly conserved amino acid (PolyPhen score 0.999; SIFT score 0)

Novel candidate

14DG0947

GUSB:NM_000181:exon7:c.1069C > T:p.R357X

Truncation of >75 % of the protein sequence

Known cause of NIHF in severe cases

14DG1037

GALNT14:NM_024572:exon13:c.C1273T:p.R425X

Truncation of >75 % of the protein sequence

Novel candidate

14DG1695

THSD1:NM_018676:exon3:c.617G > A:p.C206Y

Replaces a highly conserved amino acid (PolyPhen score 0.999; SIFT score 0)

Novel candidate

14DG11738

THSD1:NM_018676:exon3:c.G670A:p.R224X

Truncation of >75 % of the protein sequence

Novel candidate

The single candidate variant that survived the various filters in two of the remaining families was identical — a missense mutation in THSD1 (c.617G > A:p.C206Y) — and linkage analysis combining the two families revealed a single linkage peak that overlaps with THSD1 (Fig. 2). Surprisingly, in family 13DG0806, the missense mutation in THSD1 was also found in the homozygous state in two siblings who are said to have had severe edema on prenatal ultrasound and on postnatal examination but that gradually resolved with age. Since this suggests that THSD1 mutation can be compatible with an attenuated phenotype, we decided to enroll two consanguineous families in which NIHF was not associated with a lethal outcome. In family 14DG1695, the mother presented during pregnancy and her fetus had typical findings of NIHF but gave history of two previous children with a similar presentation who later improved spontaneously. Sequencing of THSD1 in the fetus confirmed the presence of the same missense mutation. Unfortunately, this fetus died shortly after birth due to severe edema with respiratory compromise. Family 14DG1738 is another consanguineous family in which two siblings presented with persistent lymphedema and have history of NIHF during their pregnancies. Sequencing of THSD1 revealed a homozygous truncating mutation (c.G670A:p.R224X; Fig. 2; Figure S2 in Additional file 1). This variant has been reported to be heterozygous in two individuals in the ExAC browser (rs9536062) for an allele frequency of 0.00001648, which is compatible with this being a disease-causing allele. Reassuringly, WES in these two families independently revealed the same mutations by applying the filters we used in the other study families.
Fig. 2

a Pedigrees of the two families that generated linkage to THSD1 (index boxed in red). b Genome-wide linkage analysis showing a single significant peak (red arrow) that spans the THSD1 locus. c Cartoon of the THSD1 gene and protein with the location of the two identified mutations marked. The missense mutation replaces a highly conserved amino acid across species

In addition to THSD1, the candidate genes identified in the remaining families have not been previously linked to human diseases and thus represent novel candidates for embryonic lethality (PIGC, UBN1, DNAH14, MYOM1, FZD6, and GALNT14).

Discussion

The spectrum of phenotypes associated with Mendelian disorders is very wide, and embryonic lethality can be viewed as the severe end of that spectrum [1]. It has long been noted that embryonic lethality can be the presenting feature of known Mendelian disorders in severe cases. However, until recently it has not been possible to study embryonic lethality as the “entry” phenotype because the phenotype is usually too nonspecific to guide testing. WES represents an exciting methodology that bypasses this limitation since it targets the entire coding genome regardless of the phenotype [10]. In this study, we extended our experience in combining the power of WES with autozygosity mapping for determining the genetic causes of many recessive Mendelian conditions in our highly consanguineous population [11] to the study of embryonic lethality as another recessive Mendelian condition in selected families. Specifically, we targeted multiplex consanguineous families in which embryonic lethality likely follows a recessive pattern. Our choice of lethal NIHF was based on the longstanding experience with this phenotype as the final common pathway of diverse pathophysiological perturbations leading to fetal loss, many of which are genetic in etiology, especially in recurrent cases [12].

The mutations we identified in GUSB, CTSA, NEB and CHRNA1 can be viewed as “positive controls” in that they lend credence to the methodology used in this study to identify the cause of embryonic lethality. These genes are known to be mutated in established Mendelian diseases (mucopolysaccharidosis, galactosialidosis, nemaline myopathy and fetal akinesia, respectively) that present as embryonic lethality in severe cases. On the other hand, the genes identified in the other families have not been previously linked to any Mendelian phenotype so they represent novel candidate genes. In the case of THSD1, it appears that mutations in this gene cause a range of Mendelian phenotypes ranging from a lethal form of NIHF to self-limited lymphedema that resolves with age. For the other novel candidate genes, their identification in one family each precludes speculation on the range of phenotypes they may cause besides the proposed link to embryonic lethality.

THSD1 is the most compelling candidate we identified in this study in view of the linkage data and the identification of two independent mutations. It encodes a thrombospondin type 1 domain-containing protein of poorly understood function. It was first identified in 2006 as a marker of primitive hematopoietic stem cells and endothelial cells [13]. Recently, it has been proposed that THSD1 has a potential role in angiogenesis and maintenance of vascular integrity [14]. Thus, it is tempting to speculate that the mutations we identified in this gene compromise vascular integrity resulting in a range of phenotypes from embryonic lethality to persistent or self-limiting edema as we observed in the families we report in this study.

Unlike THSD1, however, the remaining novel genes remain candidates since they were observed to be mutated only once, and additional families will be required to confirm their proposed link to embryonic lethality. Nevertheless, available information on the function of these genes and their involvement in basic cellular functions appear suggestive. For example, PIGC encodes an endoplasmic reticulum membrane protein required for the first step of glycosylphosphatidylinositol (GPI) biosynthesis [15]. GPI is a major posttranslational modification of many eukaryotic proteins to anchor them to membranes. Disruption of GPI2 (the yeast homolog of PIGC in human) is lethal in yeasts [16]. Similarly, GALNT14 encodes an enzyme that catalyzes the first step in the O-glycosylation of mammalian proteins by transferring N-acetyl-D-galactosamine (GalNAc) to peptide substrates. The homozygous truncating mutation we identified in this gene will likely cause a severe form of O-glycosylation disorder that may have resulted in embryonic lethality in the respective family. It is worth highlighting that disruption of the genes involved in this pathway are lethal in yeast and mouse [17, 18]. Available literature supports the candidacy of the other candidate genes as well. UBN1 encodes an essential component of a chromatin remodeling complex that also includes HIRA [19]. Intact UBN1 is required for HIRA stabilization and proper localization, and HIRA −/− is embryonic lethal in mouse [20, 21]. MYOM1 encodes myomesin, a major binding partner of another sarcomeric protein Titin, a major protein component of striated muscles [22]. Homozygous deletion of the Titin domain that interacts with myomesin results in embryonic lethality in mouse [23]. For the remaining genes, there is insufficient evidence to support their potential role in embryonic lethality.

Conclusions

Our study, to our knowledge, represents the first attempt to systematically catalogue monogenic causes of embryonic lethality in humans by combining WES and autozygome analysis in selected families. The apparently high yield of this approach in the study cohort justifies future expansion to hundreds of similarly selected families to accelerate the discovery of embryonic lethal genes in humans.

Materials and methods

Human subjects

Families were enrolled in this study if parents were consanguineous and had a history of at least two pregnancies diagnosed with lethal forms of NIHF using standard criteria used in maternal-fetal medicine. After signing written informed consent, as part of an institutional review board-approved research protocol in compliance with the Declaration of Helsinki (KFSHRC RAC# 2080006 and 2121053), a blood sample from at least one affected fetus was obtained from the umbilical cord and stored in EDTA tubes for further analysis. Detailed morphological examination by ultrasonography and, when applicable, postnatal clinical examination were recorded.

Autozygosity mapping

Determination of the autozygous genomic intervals (autozygome) was achieved by using runs of homozygosity (ROH) that are >2 Mb in size extracted from a genome-wide genotyping array as described before [24]. Briefly, DNA extracted from whole blood was run on Axiom SNP Chip platform following the manufacturer’s instructions (Affymetrix) followed by ROH determination using AutoSNPa software [25].

Whole-exome sequencing

Exome capture was performed using TruSeq Exome Enrichment kit (Illumina) following the manufacturer’s protocol. Samples were prepared as an Illumina sequencing library, and in the second step, the sequencing libraries were enriched for the desired target using the Illumina Exome Enrichment protocol. The captured libraries were sequenced using an Illumina HiSeq 2000 Sequencer. The reads were mapped against UCSC hg19 [26] by BWA [27]. The SNPs and Indels were detected by SAMTOOLS [28].

Mapping of likely causal variants

Variants from WES were filtered such that only novel (or very low frequency, 0.1 %), coding/splicing, homozygous variants that are within the autozygome of the affected fetus and are predicted to be pathogenic were considered as likely causal variants [29]. Frequency of variants was determined using publically available variant databases (1000 Genomes, Exome Variant Server and ExAC) as well as a database of 560 in-house ethnically matched exomes. Pathogenicity was likely if the mutation is loss-of-function (splicing/truncating) or, in the case of missense/in-frame indels, removes a highly conserved amino acid and is predicted to be pathogenic by the two in silico prediction modules PolyPhen and SIFT. All exomes described in this study can be accessed at [30].

Availability of supporting data

All novel variants in this study have been deposited in ClinVar with the following accession numbers: SCV000222665–SCV000222684. Access to the whole-exome data used in this study is provided through the TraBioS Data Warehouse [30].

Abbreviations

GPI: 

glycosylphosphatidylinositol

NIHF: 

non-immune hydrops fetalis

ROH: 

runs of homozygosity

WES: 

whole-exome sequencing

Declarations

Acknowledgement

We are very grateful to the families for their enthusiastic participation in this study. We also thank the Genotyping and Sequencing Core Facilities at KFSHRC for their technical help. This work was supported in part by KACST 13-BIO1113-20 (FSA).

Authors’ Affiliations

(1)
Department of Genetics, King Faisal Specialist Hospital and Research Center
(2)
Department of Obstetrics and Gynecology, King Faisal Specialist Hospital and Research Center
(3)
Department of Obstetrics and Gynecology, Security Forces Hospital
(4)
Department of Medical Genetics, King Faisal Specialist Hospital and Research Center
(5)
Department of Pediatrics, Prince Sultan Military Medical City
(6)
Department of Obstetrics and Gynecology, Prince Sultan Military Medical City
(7)
Department of Obstetrics and Gynecology, College of Medicine, King Saud University
(8)
Department of Pediatrics, Faculty of Medicine, Ain Shams University
(9)
Department of Pediatrics, King Fahd Hospital of the University, University of Dammam
(10)
Department of Anatomy and Cell Biology, College of Medicine, Alfaisal University

References

  1. Alkuraya FS. Human knockout research: new horizons and opportunities. Trends Genet. 2015;31:108–15.PubMedView ArticleGoogle Scholar
  2. Giaever G, Chu AM, Ni L, Connelly C, Riles L, Véronneau S, et al. Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002;418:387–91.PubMedView ArticleGoogle Scholar
  3. Ayadi A, Birling M-C, Bottomley J, Bussell J, Fuchs H, Fray M, et al. Mouse large-scale phenotyping initiatives: Overview of the European Mouse Disease Clinic (EUMODIC) and of the Wellcome Trust Sanger Institute Mouse Genetics Project. Mamm Genome. 2012;23:600–10.PubMed CentralPubMedView ArticleGoogle Scholar
  4. Kattner E, Schäfer A, Harzer K. Hydrops fetalis: manifestation in lysosomal storage diseases including Farber disease. Eur J Pediatr. 1997;156:292–5.PubMedView ArticleGoogle Scholar
  5. Carvalho S, Martins M, Fortuna A, Ramos U, Ramos C, Rodrigues MC. Galactosialidosis presenting as nonimmune fetal hydrops: a case report. Prenat Diagn. 2009;29:895–6.PubMedView ArticleGoogle Scholar
  6. Kleijer W, Geilen G, Janse H, Van Diggelen O, Zhou XY, Galjart N, et al. Cathepsin A deficiency in galactosialidosis: studies of patients and carriers in 16 families. Pediatr Res. 1996;39:1067–71.PubMedView ArticleGoogle Scholar
  7. Stangenberg M, Lingman G, Roberts G, Ozand P. Mucopolysaccharidosis VII as cause of fetal hydrops in early pregnancy. Am J Med Genet. 1992;44:142–4.PubMedView ArticleGoogle Scholar
  8. Schroer A, Stefanova I, Weichert J. P07. 07: Prenatal diagnosis of rare mucopolysaccharidosis type VII (MPS VII, Morbus Sly) associated with hydrops fetalis. Ultrasound Obstet Gynecol. 2011;38:192.View ArticleGoogle Scholar
  9. Vardon D, Chau C, Sigodi S, Figarella-Branger D, Boubli L. Congenital rapidly fatal form of nemaline myopathy with fetal hydrops and arthrogryposis. Fetal Diagn Ther. 1998;13:244–9.PubMedView ArticleGoogle Scholar
  10. Shamseldin HE, Swaid A, Alkuraya FS. Lifting the lid on unborn lethal Mendelian phenotypes through exome sequencing. Genet Med. 2012;15:307–9.PubMed CentralPubMedView ArticleGoogle Scholar
  11. Alkuraya FS. Genetics and genomic medicine in Saudi Arabia. Mol Genet Genomic Med. 2014;2:369–78.PubMed CentralPubMedView ArticleGoogle Scholar
  12. Bellini C, Hennekam R, Fulcheri E, Rutigliani M, Morcaldi G, Boccardo F, et al. Etiology of nonimmune hydrops fetalis: a systematic review. Am J Med Genet A. 2009;149:844–51.View ArticleGoogle Scholar
  13. S-i T, Hiroyama T, Yamazaki S, Nakajima T, Morita Y, Usui J, et al. Genetic marking of hematopoietic stem and endothelial cells: identification of the Tmtsp gene encoding a novel cell surface protein with the thrombospondin-1 domain. Blood. 2006;107:4317–25.View ArticleGoogle Scholar
  14. den Dekker WK, Haasdijk R, Tempel D, Houtgraaf JH, de Jong R, Cheng C, et al. Thsd-1 determines plaque phenotype by regulating vascular permeability and intraplaque haemorrhaging. Circulation. http://circ.ahajournals.org.ezp-prod1.hul.harvard.edu/cgi/content/meeting_abstract/124/21_MeetingAbstracts/A14844.
  15. Inoue N, Watanabe R, Takeda J, Kinoshita T. PIG-C, one of the three human genes involved in the first step of glycosylphosphatidylinositol biosynthesis is a homologue of Saccharomyces cerevisiae GPI2. Biochem Biophys Res Commun. 1996;226:193–9.PubMedView ArticleGoogle Scholar
  16. Leidich SD, Kostova Z, Latek RR, Costello LC, Drapp DA, Gray W, et al. Temperature-sensitive yeast GPI anchoring mutants gpi2 and gpi3 are defective in the synthesis of N-acetylglucosaminyl phosphatidylinositol. Cloning of the GPI2 gene. J Biol Chem. 1995;270:13029–35.PubMedView ArticleGoogle Scholar
  17. Gentzsch M, Tanner W. The PMT gene family: protein O-glycosylation in Saccharomyces cerevisiae is vital. EMBO J. 1996;15:5752.PubMed CentralPubMedGoogle Scholar
  18. Wells L, Whelan SA, Hart GW. O-GlcNAc: a regulatory post-translational modification. Biochem Biophys Res Commun. 2003;302:435–41.PubMedView ArticleGoogle Scholar
  19. Banumathy G, Somaiah N, Zhang R, Tang Y, Hoffmann J, Andrake M, et al. Human UBN1 is an ortholog of yeast Hpc2p and has an essential role in the HIRA/ASF1a chromatin-remodeling pathway in senescent cells. Mol Cell Biol. 2009;29:758–70.PubMed CentralPubMedView ArticleGoogle Scholar
  20. Tang Y, Puri A, Ricketts MD, Rai TS, Hoffmann J, Hoi E, et al. Identification of an ubinuclein 1 region required for stability and function of the human HIRA/UBN1/CABIN1/ASF1a histone H3. 3 chaperone complex. Biochemistry. 2012;51:2366–77.PubMed CentralPubMedView ArticleGoogle Scholar
  21. Roberts C, Sutherland HF, Farmer H, Kimber W, Halford S, Carey A, et al. Targeted mutagenesis of the Hira gene results in gastrulation defects and patterning abnormalities of mesoendodermal derivatives prior to early embryonic lethality. Mol Cell Biol. 2002;22:2318–28.PubMed CentralPubMedView ArticleGoogle Scholar
  22. Hein S, Schaper J. Weakness of a giant: mutations of the sarcomeric protein titin. Trends Mol Med. 2002;8:311–3.PubMedView ArticleGoogle Scholar
  23. Gotthardt M, Hammer RE, Hübner N, Monti J, Witt CC, McNabb M, et al. Conditional expression of mutant M-line titins results in cardiomyopathy with altered sarcomere structure. J Biol Chem. 2003;278:6059–65.PubMedView ArticleGoogle Scholar
  24. Alkuraya FS. Discovery of rare homozygous mutations from studies of consanguineous pedigrees. Curr Protoc Hum Genet. 2012:61211–16.12. 13.
  25. Carr IM, Flintoff KJ, Taylor GR, Markham AF, Bonthron DT. Interactive visual analysis of SNP data for rapid autozygosity mapping in consanguineous families. Hum Mutat. 2006;27:1041–6.PubMedView ArticleGoogle Scholar
  26. UCSC Genome Browser. http://genome.ucsc.edu/.
  27. Burrows-Wheeler Aligner. http://bio-bwa.sourceforge.net/.
  28. SAM (Sequence Alignment/Map) Tools. http://samtools.sourceforge.net/.
  29. Alkuraya FS. The application of next-generation sequencing in the autozygosity mapping of human recessive diseases. Hum Genet. 2013;132:1197–211.PubMedView ArticleGoogle Scholar
  30. TraBioS Data Warehouse. http://shgp.kfshrc.edu.sa/bioinf/db/variants/dg/details.html.

Copyright

© Shamseldin et al. 2015

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.