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.

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.

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.

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.

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.