Investigations of humans with disorders of sex development (DSDs) resulted in the discovery of many of the now-known mammalian sex-determining genes, including SRY, RSPO1, SOX9, NR5A1, WT1, NR0B1, and WNT4. Here, the locus for an autosomal sex-determining gene was mapped via linkage analysis in two families with 46,XY DSD to the long arm of chromosome 5 with a combined, multipoint parametric LOD score of 6.21. A splice-acceptor mutation (c.634-8T>A) in MAP3K1 segregated with the phenotype in the first family and disrupted RNA splicing. Mutations were demonstrated in the second family (p.Gly616Arg) and in two of 11 sporadic cases (p.Leu189Pro, p.Leu189Arg)-18% prevalence in this cohort of sporadic cases. In cultured primary lymphoblastoid cells from family 1 and the two sporadic cases, these mutations altered the phosphorylation of the downstream targets, p38 and ERK1/2, and enhanced binding of RHOA to the MAP3K1 complex. Map3k1 within the syntenic region was expressed in the embryonic mouse gonad prior to, and after, sex determination. Thus, mutations in MAP3K1 that result in 46,XY DSD with partial or complete gonadal dysgenesis implicate this pathway in normal human sex determination.

Orofaciodigital syndromes (OFDSs) consist of a group of heterogeneous disorders characterized by abnormalities in the oral cavity, face, and digits and associated phenotypic abnormalities that lead to the delineation of 13 OFDS subtypes. Here, by a combined approach of homozygozity mapping and exome ciliary sequencing, we identified truncating TCTN3 mutations as the cause of an extreme form of OFD associated with bone dysplasia, tibial defect, cystic kidneys, and brain anomalies (OFD IV, Mohr-Majewski syndrome). Analysis of 184 individuals with various ciliopathies (OFD, Meckel, Joubert, and short rib polydactyly syndromes) led us to identify four additional truncating TCTN3 mutations in unrelated fetal cases with overlapping Meckel and OFD IV syndromes and one homozygous missense mutation in a family with Joubert syndrome. By exploring roles of TCTN3 in human ciliary related functions, we found that TCTN3 is necessary for transduction of the sonic hedgehog (SHH) signaling pathway, as revealed by abnormal processing of GLI3 in patient cells. These results are consistent with the suggested role of its murine ortholog, which forms a complex at the ciliary transition zone with TCTN1 and TCTN2, both of which are also implicated in the transduction of SHH signaling. Overall, our data show the involvement of the transition zone protein TCTN3 in the regulation of the key SHH signaling pathway and that its disruption causes a severe form of ciliopathy, combining features of Meckel and OFD IV syndromes.

Variants in genes encoding ribosomal proteins have thus far been associated with Diamond-Blackfan anemia, a rare inherited bone marrow failure, and isolated congenital asplenia. Here, we report one de novo missense variant and three de novo splice variants in RPL13, which encodes ribosomal protein RPL13 (also called eL13), in four unrelated individuals with a rare bone dysplasia causing severe short stature. The three splice variants (c.477+1G>T, c.477+1G>A, and c.477+2 T>C) result in partial intron retention, which leads to an 18-amino acid insertion. In contrast to observations from Diamond-Blackfan anemia, we detected no evidence of significant pre-rRNA processing disturbance in cells derived from two affected individuals. Consistently, we showed that the insertion-containing protein is stably expressed and incorporated into 60S subunits similar to the wild-type protein. Erythroid proliferation in culture and ribosome profile on sucrose gradient are modified, suggesting a change in translation dynamics. We also provide evidence that RPL13 is present at high levels in chondrocytes and osteoblasts in mouse growth plates. Taken together, we show that the identified RPL13 variants cause a human ribosomopathy defined by a rare skeletal dysplasia, and we highlight the role of this ribosomal protein in bone development.

Mutations in the GLI3 zinc-finger transcription factor gene cause Greig cephalopolysyndactyly syndrome (GCPS) and Pallister-Hall syndrome (PHS), which are variable but distinct clinical entities. We hypothesized that GLI3 mutations that predict a truncated functional repressor protein cause PHS and that functional haploinsufficiency of GLI3 causes GCPS. To test these hypotheses, we screened patients with PHS and GCPS for GLI3 mutations. The patient group consisted of 135 individuals: 89 patients with GCPS and 46 patients with PHS. We detected 47 pathological mutations (among 60 probands); when these were combined with previously published mutations, two genotype-phenotype correlations were evident. First, GCPS was caused by many types of alterations, including translocations, large deletions, exonic deletions and duplications, small in-frame deletions, and missense, frameshift/nonsense, and splicing mutations. In contrast, PHS was caused only by frameshift/nonsense and splicing mutations. Second, among the frameshift/nonsense mutations, there was a clear genotype-phenotype correlation. Mutations in the first third of the gene (from open reading frame [ORF] nucleotides [nt] 1-1997) caused GCPS, and mutations in the second third of the gene (from ORF nt 1998-3481) caused primarily PHS. Surprisingly, there were 12 mutations in patients with GCPS in the 3' third of the gene (after ORF nt 3481), and no patients with PHS had mutations in this region. These results demonstrate a robust correlation of genotype and phenotype for GLI3 mutations and strongly support the hypothesis that these two allelic disorders have distinct modes of pathogenesis.

A large French family including members affected by nonspecific X-linked mental retardation, with or without autism or pervasive developmental disorder in affected male patients, has been found to have a 2-base-pair deletion in the Neuroligin 4 gene (NLGN4) located at Xp22.33. This mutation leads to a premature stop codon in the middle of the sequence of the normal protein and is thought to suppress the transmembrane domain and sequences important for the dimerization of neuroligins that are required for proper cell-cell interaction through binding to beta-neurexins. As the neuroligins are mostly enriched at excitatory synapses, these results suggest that a defect in synaptogenesis may lead to deficits in cognitive development and communication processes. The fact that the deletion was present in both autistic and nonautistic mentally retarded males suggests that the NLGN4 gene is not only involved in autism, as previously described, but also in mental retardation, indicating that some types of autistic disorder and mental retardation may have common genetic origins.

Autosomal-recessive cerebellar ataxias comprise a clinically and genetically heterogeneous group of neurodegenerative disorders. In contrast to their dominant counterparts, unraveling the molecular background of these ataxias has proven to be more complicated and the currently known mutations provide incomplete coverage for genotyping of patients. By combining SNP array-based linkage analysis and targeted resequencing of relevant sequences in the linkage interval with the use of next-generation sequencing technology, we identified a mutation in a gene and have shown its association with autosomal-recessive cerebellar ataxia. In a Dutch consanguineous family with three affected siblings a homozygous 12.5 Mb region on chromosome 3 was targeted by array-based sequence capture. Prioritization of all detected sequence variants led to four candidate genes, one of which contained a variant with a high base pair conservation score (phyloP score: 5.26). This variant was a leucine-to-arginine substitution in the DUF 590 domain of a 16K transmembrane protein, a putative calcium-activated chloride channel encoded by anoctamin 10 (ANO10). The analysis of ANO10 by Sanger sequencing revealed three additional mutations: a homozygous mutation (c.1150_1151del [p.Leu384fs]) in a Serbian family and a compound-heterozygous splice-site mutation (c.1476+1G>T) and a frameshift mutation (c.1604del [p.Leu535X]) in a French family. This illustrates the power of using initial homozygosity mapping with next-generation sequencing technology to identify genes involved in autosomal-recessive diseases. Moreover, identifying a putative calcium-dependent chloride channel involved in cerebellar ataxia adds another pathway to the list of pathophysiological mechanisms that may cause cerebellar ataxia.

Autosomal-recessive immunodeficiency, centromeric instability, and facial anomalies (ICF) syndrome is mainly characterized by recurrent, often fatal, respiratory and gastrointestinal infections. About 50% of patients carry mutations in the DNA methyltransferase 3B gene (DNMT3B) (ICF1). The remaining patients carry unknown genetic defects (ICF2) but share with ICF1 patients the same immunological and epigenetic features, including hypomethylation of juxtacentromeric repeat sequences. We performed homozygosity mapping in five unrelated ICF2 patients with consanguineous parents and then performed whole-exome sequencing in one of these patients and Sanger sequencing in all to identify mutations in the zinc-finger- and BTB (bric-a-bric, tramtrack, broad complex)-domain-containing 24 (ZBTB24) gene in four consanguineously descended ICF2 patients. Additionally, we found ZBTB24 mutations in an affected sibling pair and in one patient for whom it was not known whether his parents were consanguineous. ZBTB24 belongs to a large family of transcriptional repressors that include members, such as BCL6 and PATZ1, with prominent regulatory roles in hematopoietic development and malignancy. These data thus indicate that ZBTB24 is involved in DNA methylation of juxtacentromeric DNA and in B cell development and/or B and T cell interactions. Because ZBTB24 is a putative DNA-binding protein highly expressed in the lymphoid lineage, we predict that by studying the molecular function of ZBTB24, we will improve our understanding of the molecular pathophysiology of ICF syndrome and of lymphocyte biology in general.

The immunoglobulin heavy-chain locus (IGH) encodes variable (IGHV), diversity (IGHD), joining (IGHJ), and constant (IGHC) genes and is responsible for antibody heavy-chain biosynthesis, which is vital to the adaptive immune response. Programmed V-(D)-J somatic rearrangement and the complex duplicated nature of the locus have impeded attempts to reconcile its genomic organization based on traditional B-lymphocyte derived genetic material. As a result, sequence descriptions of germline variation within IGHV are lacking, haplotype inference using traditional linkage disequilibrium methods has been difficult, and the human genome reference assembly is missing several expressed IGHV genes. By using a hydatidiform mole BAC clone resource, we present the most complete haplotype of IGHV, IGHD, and IGHJ gene regions derived from a single chromosome, representing an alternate assembly of ∼1 Mbp of high-quality finished sequence. From this we add 101 kbp of previously uncharacterized sequence, including functional IGHV genes, and characterize four large germline copy-number variants (CNVs). In addition to this germline reference, we identify and characterize eight CNV-containing haplotypes from a panel of nine diploid genomes of diverse ethnic origin, discovering previously unmapped IGHV genes and an additional 121 kbp of insertion sequence. We genotype four of these CNVs by using PCR in 425 individuals from nine human populations. We find that all four are highly polymorphic and show considerable evidence of stratification (Fst = 0.3-0.5), with the greatest differences observed between African and Asian populations. These CNVs exhibit weak linkage disequilibrium with SNPs from two commercial arrays in most of the populations tested.

Rare mutations, including copy-number variants (CNVs), contribute significantly to autism spectrum disorder (ASD) risk. Although their importance has been established in families with only one affected child (simplex families), the contribution of both de novo and inherited CNVs to ASD in families with multiple affected individuals (multiplex families) is less well understood. We analyzed 1,532 families from the Autism Genetic Resource Exchange (AGRE) to assess the impact of de novo and rare CNVs on ASD risk in multiplex families. We observed a higher burden of large, rare CNVs, including inherited events, in individuals with ASD than in their unaffected siblings (odds ratio [OR] = 1.7), but the rate of de novo events was significantly lower than in simplex families. In previously characterized ASD risk loci, we identified 49 CNVs, comprising 24 inherited events, 19 de novo events, and 6 events of unknown inheritance, a significant enrichment in affected versus control individuals (OR = 3.3). In 21 of the 30 families (71%) in whom at least one affected sibling harbored an established ASD major risk CNV, including five families harboring inherited CNVs, the CNV was not shared by all affected siblings, indicating that other risk factors are contributing. We also identified a rare risk locus for ASD and language delay at chromosomal region 2q24 (implicating NR4A2) and another lower-penetrance locus involving inherited deletions and duplications of WWOX. The genetic architecture in multiplex families differs from that in simplex families and is complex, warranting more complete genetic characterization of larger multiplex ASD cohorts.

Charcot-Marie-Tooth disease (CMT) is the most common cause of inherited peripheral neuropathy, with an estimated frequency of 1/2500. We studied a large family with 17 patients affected by the axonal form of CMT (CMT2). Analysis of the 15 genes or loci known to date was negative. Genome-wide genotyping identified a CMT2 locus in 16q21-q23 between D16S3050 and D16S3106. The maximum two-point LOD score was 4.77 at theta = 0 for marker D16S3050. Sequencing of candidate genes identified a unique mutation, c.986G>A (p.Arg329His), affecting a totally conserved amino acid in the helical domain of cytoplasmic alanyl-tRNA synthetase (AlaRS). A second family with the same mutation and a different founder was then identified in a cohort of 91 CMT2 families. Although mislocation of mutant Arg329His-AlaRS in axons remains to be evaluated, experimental data point mostly to a quantitative reduction in tRNA(Ala) aminoacylation. Aminoacylation and editing functions closely cooperate in AlaRS, and Arg329His mutation could also lead to qualitative errors participating in neurodegeneration. Our report documents in 18 patients the deleterious impact of a mutation in human cytoplasmic AlaRS and broadens the spectrum of defects found in tRNA synthetases. Patients present with sensory-motor distal degeneration secondary to predominant axonal neuropathy, slight demyelination, and no atypical or additional CNS features.

Mesomelia-synostoses syndrome (MSS) or mesomelic dysplasia with acral synostoses Verloes-David-Pfeiffer type is a rare autosomal-dominant disorder characterized by mesomelic limb shortening, acral synostoses, and multiple congenital malformations. So far, five patients in four unrelated families have been reported worldwide with MMS. By using whole-genome oligonucleotide array CGH, we have identified an interstitial deletion at 8q13 in all patients. The deletions vary from 582 Kb to 738 Kb in size, but invariably encompass only two genes: SULF1, encoding the heparan sulfate 6-O-endosulfatase 1, and SLCO5A1, encoding the solute carrier organic anion transporter family member 5A1. SULF1 acts as a regulator of numerous growth factors in skeletal embryonic development whereas the function of SLCO5A1 is yet unknown. Breakpoint sequence analyses performed in two families showed nonrecurrent deletions. Real-time quantitative RT-PCR analysis showed the highest levels of SULF1 transcripts in human osteoblasts and cartilage whereas SLCO5A1 was highly expressed in human fetal and adult brain and heart. Our results strongly suggest that haploinsufficiency of SULF1 contributes to this mesomelic chondrodysplasia, highlighting the critical role of endosulfatase in human skeletal development. Codeletion of SULF1 and SLCO5A1--which does not result from a low-copy repeats (LCRs)-mediated recombination event in at least two families--was found in all patients, so we suggest that haploinsufficiency of SULF1 combined with haploinsufficiency of SLCO5A1 (or the altered expression of a neighboring gene through position effect) could be necessary in the pathogenesis of MSS.

Primary defects in lung branching morphogenesis, resulting in neonatal lethal pulmonary hypoplasias, are incompletely understood. To elucidate the pathogenetics of human lung development, we studied a unique collection of samples obtained from deceased individuals with clinically and histopathologically diagnosed interstitial neonatal lung disorders: acinar dysplasia (n = 14), congenital alveolar dysplasia (n = 2), and other lethal lung hypoplasias (n = 10). We identified rare heterozygous copy-number variant deletions or single-nucleotide variants (SNVs) involving TBX4 (n = 8 and n = 2, respectively) or FGF10 (n = 2 and n = 2, respectively) in 16/26 (61%) individuals. In addition to TBX4, the overlapping ∼2 Mb recurrent and nonrecurrent deletions at 17q23.1q23.2 identified in seven individuals with lung hypoplasia also remove a lung-specific enhancer region. Individuals with coding variants involving either TBX4 or FGF10 also harbored at least one non-coding SNV in the predicted lung-specific enhancer region, which was absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies. The occurrence of rare coding variants involving TBX4 or FGF10 with the putative hypomorphic non-coding SNVs implies a complex compound inheritance of these pulmonary hypoplasias. Moreover, they support the importance of TBX4-FGF10-FGFR2 epithelial-mesenchymal signaling in human lung organogenesis and help to explain the histopathological continuum observed in these rare lethal developmental disorders of the lung.

Congenital poikiloderma is characterized by a combination of mottled pigmentation, telangiectasia, and epidermal atrophy in the first few months of life. We have previously described a South African European-descent family affected by a rare autosomal-dominant form of hereditary fibrosing poikiloderma accompanied by tendon contracture, myopathy, and pulmonary fibrosis. Here, we report the identification of causative mutations in FAM111B by whole-exome sequencing. In total, three FAM111B missense mutations were identified in five kindreds of different ethnic backgrounds. The mutation segregated with the disease in one large pedigree, and mutations were de novo in two other pedigrees. All three mutations were absent from public databases and were not observed on Sanger sequencing of 388 ethnically matched control subjects. The three single-nucleotide mutations code for amino acid changes that are clustered within a putative trypsin-like cysteine/serine peptidase domain of FAM111B. These findings provide evidence of the involvement of FAM111B in congenital poikiloderma and multisystem fibrosis.

Genomic imprinting is a mechanism in which gene expression varies depending on parental origin. Imprinting occurs through differential epigenetic marks on the two parental alleles, with most imprinted loci marked by the presence of differentially methylated regions (DMRs). To identify sites of parental epigenetic bias, here we have profiled DNA methylation patterns in a cohort of 57 individuals with uniparental disomy (UPD) for 19 different chromosomes, defining imprinted DMRs as sites where the maternal and paternal methylation levels diverge significantly from the biparental mean. Using this approach we identified 77 DMRs, including nearly all those described in previous studies, in addition to 34 DMRs not previously reported. These include a DMR at TUBGCP5 within the recurrent 15q11.2 microdeletion region, suggesting potential parent-of-origin effects associated with this genomic disorder. We also observed a modest parental bias in DNA methylation levels at every CpG analyzed across ∼1.9 Mb of the 15q11-q13 Prader-Willi/Angelman syndrome region, demonstrating that the influence of imprinting is not limited to individual regulatory elements such as CpG islands, but can extend across entire chromosomal domains. Using RNA-seq data, we detected signatures consistent with imprinted expression associated with nine novel DMRs. Finally, using a population sample of 4,004 blood methylomes, we define patterns of epigenetic variation at DMRs, identifying rare individuals with global gain or loss of methylation across multiple imprinted loci. Our data provide a detailed map of parental epigenetic bias in the human genome, providing insights into potential parent-of-origin effects.