Meier-Gorlin syndrome (MGS) is a genetically heterogeneous primordial dwarfism syndrome known to be caused by biallelic loss-of-function mutations in one of five genes encoding pre-replication complex proteins: ORC1, ORC4, ORC6, CDT1, and CDC6. Mutations in these genes cause disruption of the origin of DNA replication initiation. To date, only an autosomal-recessive inheritance pattern has been described in individuals with this disorder, with a molecular etiology established in about three-fourths of cases. Here, we report three subjects with MGS and de novo heterozygous mutations in the 5' end of GMNN, encoding the DNA replication inhibitor geminin. We identified two truncating mutations in exon 2 (the 1(st) coding exon), c.16A>T (p.Lys6(∗)) and c.35_38delTCAA (p.Ile12Lysfs(∗)4), and one missense mutation, c.50A>G (p.Lys17Arg), affecting the second-to-last nucleotide of exon 2 and possibly RNA splicing. Geminin is present during the S, G2, and M phases of the cell cycle and is degraded during the metaphase-anaphase transition by the anaphase-promoting complex (APC), which recognizes the destruction box sequence near the 5' end of the geminin protein. All three GMNN mutations identified alter sites 5' to residue Met28 of the protein, which is located within the destruction box. We present data supporting a gain-of-function mechanism, in which the GMNN mutations result in proteins lacking the destruction box and hence increased protein stability and prolonged inhibition of replication leading to autosomal-dominant MGS.

The overall understanding of the molecular etiologies of intellectual disability (ID) and developmental delay (DD) is increasing as next-generation sequencing technologies identify genetic variants in individuals with such disorders. However, detailed analyses conclusively confirming these variants, as well as the underlying molecular mechanisms explaining the diseases, are often lacking. Here, we report on an ID syndrome caused by de novo heterozygous loss-of-function (LoF) mutations in SON. The syndrome is characterized by ID and/or DD, malformations of the cerebral cortex, epilepsy, vision problems, musculoskeletal abnormalities, and congenital malformations. Knockdown of son in zebrafish resulted in severe malformation of the spine, brain, and eyes. Importantly, analyses of RNA from affected individuals revealed that genes critical for neuronal migration and cortex organization (TUBG1, FLNA, PNKP, WDR62, PSMD3, and HDAC6) and metabolism (PCK2, PFKL, IDH2, ACY1, and ADA) are significantly downregulated because of the accumulation of mis-spliced transcripts resulting from erroneous SON-mediated RNA splicing. Our data highlight SON as a master regulator governing neurodevelopment and demonstrate the importance of SON-mediated RNA splicing in human development.

Genetic differences that specify unique aspects of human evolution have typically been identified by comparative analyses between the genomes of humans and closely related primates, including more recently the genomes of archaic hominins. Not all regions of the genome, however, are equally amenable to such study. Recurrent copy number variation (CNV) at chromosome 16p11.2 accounts for approximately 1% of cases of autism and is mediated by a complex set of segmental duplications, many of which arose recently during human evolution. Here we reconstruct the evolutionary history of the locus and identify bolA family member 2 (BOLA2) as a gene duplicated exclusively in Homo sapiens. We estimate that a 95-kilobase-pair segment containing BOLA2 duplicated across the critical region approximately 282 thousand years ago (ka), one of the latest among a series of genomic changes that dramatically restructured the locus during hominid evolution. All humans examined carried one or more copies of the duplication, which nearly fixed early in the human lineage--a pattern unlikely to have arisen so rapidly in the absence of selection (P < 0.0097). We show that the duplication of BOLA2 led to a novel, human-specific in-frame fusion transcript and that BOLA2 copy number correlates with both RNA expression (r = 0.36) and protein level (r = 0.65), with the greatest expression difference between human and chimpanzee in experimentally derived stem cells. Analyses of 152 patients carrying a chromosome 16p11. rearrangement show that more than 96% of breakpoints occur within the H. sapiens-specific duplication. In summary, the duplicative transposition of BOLA2 at the root of the H. sapiens lineage about 282 ka simultaneously increased copy number of a gene associated with iron homeostasis and predisposed our species to recurrent rearrangements associated with disease.

Short-read sequencing has enabled the de novo assembly of several individual human genomes, but with inherent limitations in characterizing repeat elements. Here we sequence a Chinese individual HX1 by single-molecule real-time (SMRT) long-read sequencing, construct a physical map by NanoChannel arrays and generate a de novo assembly of 2.93 Gb (contig N50: 8.3 Mb, scaffold N50: 22.0 Mb, including 39.3 Mb N-bases), together with 206 Mb of alternative haplotypes. The assembly fully or partially fills 274 (28.4%) N-gaps in the reference genome GRCh38. Comparison to GRCh38 reveals 12.8 Mb of HX1-specific sequences, including 4.1 Mb that are not present in previously reported Asian genomes. Furthermore, long-read sequencing of the transcriptome reveals novel spliced genes that are not annotated in GENCODE and are missed by short-read RNA-Seq. Our results imply that improved characterization of genome functional variation may require the use of a range of genomic technologies on diverse human populations.

Hematopoietic stem cells (HSCs) underlie the production of blood and immune cells for the lifetime of an organism. In vertebrate embryos, HSCs arise from the unique transdifferentiation of hemogenic endothelium comprising the floor of the dorsal aorta during a brief developmental window. To date, this process has not been replicated in vitro from pluripotent precursors, partly because the full complement of required signaling inputs remains to be determined. Here, we show that TNFR2 via TNF? activates the Notch and NF-?B signaling pathways to establish HSC fate, indicating a requirement for inflammatory signaling in HSC generation. We determine that primitive neutrophils are the major source of TNF?, assigning a role for transient innate immune cells in establishing the HSC program. These results demonstrate that proinflammatory signaling, in the absence of infection, is utilized by the developing embryo to generate the lineal precursors of the adult hematopoietic system.

Haematopoietic stem cells (HSCs) derive from haemogenic endothelial cells of the primitive dorsal aorta (DA) during vertebrate embryogenesis. The molecular mechanisms governing this unique endothelial to haematopoietic transition remain unclear. Here, we demonstrate a novel requirement for fibroblast growth factor (FGF) signalling in HSC emergence. This requirement is non-cell-autonomous, and acts within the somite to bridge the Wnt and Notch signalling pathways. We previously demonstrated that Wnt16 regulates the somitic expression of two Notch ligands, deltaC (dlc) and deltaD (dld), whose combined function is required for HSC fate. How Wnt16 connects to Notch function has remained an open question. Our current studies demonstrate that FGF signalling, via FGF receptor 4 (Fgfr4), mediates a signal-transduction pathway between Wnt16 and Dlc, but not Dld, to regulate HSC specification. Our findings demonstrate that FGF signalling acts as a key molecular relay within the developmental HSC niche to instruct HSC fate.

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

To assess the relative impact of inherited and de novo variants on autism risk, we generated a comprehensive set of exonic single-nucleotide variants (SNVs) and copy number variants (CNVs) from 2,377 families with autism. We find that private, inherited truncating SNVs in conserved genes are enriched in probands (odds ratio = 1.14, P = 0.0002) in comparison to unaffected siblings, an effect involving significant maternal transmission bias to sons. We also observe a bias for inherited CNVs, specifically for small (<100 kb), maternally inherited events (P = 0.01) that are enriched in CHD8 target genes (P = 7.4 × 10(-3)). Using a logistic regression model, we show that private truncating SNVs and rare, inherited CNVs are statistically independent risk factors for autism, with odds ratios of 1.11 (P = 0.0002) and 1.23 (P = 0.01), respectively. This analysis identifies a second class of candidate genes (for example, RIMS1, CUL7 and LZTR1) where transmitted mutations may create a sensitized background but are unlikely to be completely penetrant.

We have developed a computational method based on polyploid phasing of long sequence reads to resolve collapsed regions of segmental duplications within genome assemblies. Segmental Duplication Assembler (SDA; https://github.com/mvollger/SDA ) constructs graphs in which paralogous sequence variants define the nodes and long-read sequences provide attraction and repulsion edges, enabling the partition and assembly of long reads corresponding to distinct paralogs. We apply it to single-molecule, real-time sequence data from three human genomes and recover 33-79 megabase pairs (Mb) of duplications in which approximately half of the loci are diverged (<99.8%) compared to the reference genome. We show that the corresponding sequence is highly accurate (>99.9%) and that the diverged sequence corresponds to copy-number-variable paralogs that are absent from the human reference genome. Our method can be applied to other complex genomes to resolve the last gene-rich gaps, improve duplicate gene annotation, and better understand copy-number-variant genetic diversity at the base-pair level.

Mendelian disorders of cholesterol biosynthesis typically result in multi-system clinical phenotypes, underlining the importance of cholesterol in embryogenesis and development. FDFT1 encodes for an evolutionarily conserved enzyme, squalene synthase (SS, farnesyl-pyrophosphate farnesyl-transferase 1), which catalyzes the first committed step in cholesterol biosynthesis. We report three individuals with profound developmental delay, brain abnormalities, 2-3 syndactyly of the toes, and facial dysmorphisms, resembling Smith-Lemli-Opitz syndrome, the most common cholesterol biogenesis defect. The metabolite profile in plasma and urine suggested that their defect was at the level of squalene synthase. Whole-exome sequencing was used to identify recessive disease-causing variants in FDFT1. Functional characterization of one variant demonstrated a partial splicing defect and altered promoter and/or enhancer activity, reflecting essential mechanisms for regulating cholesterol biosynthesis/uptake in steady state.

Certain human traits such as neurodevelopmental disorders (NDs) and congenital anomalies (CAs) are believed to be primarily genetic in origin. However, even after whole-genome sequencing (WGS), a substantial fraction of such disorders remain unexplained. We hypothesize that some cases of ND-CA are caused by aberrant DNA methylation leading to dysregulated genome function. Comparing DNA methylation profiles from 489 individuals with ND-CAs against 1534 controls, we identify epivariations as a frequent occurrence in the human genome. De novo epivariations are significantly enriched in cases, while RNAseq analysis shows that epivariations often have an impact on gene expression comparable to loss-of-function mutations. Additionally, we detect and replicate an enrichment of rare sequence mutations overlapping CTCF binding sites close to epivariations, providing a rationale for interpreting non-coding variation. We propose that epivariations contribute to the pathogenesis of some patients with unexplained ND-CAs, and as such likely have diagnostic relevance.

The recent introductions of low-cost, long-read, and read-cloud sequencing technologies coupled with intense efforts to develop efficient algorithms have made affordable, high-quality de novo sequence assembly a realistic proposition. The result is an explosion of new, ultracontiguous genome assemblies. To compare these genomes, we need robust methods for genome annotation. We describe the fully open source Comparative Annotation Toolkit (CAT), which provides a flexible way to simultaneously annotate entire clades and identify orthology relationships. We show that CAT can be used to improve annotations on the rat genome, annotate the great apes, annotate a diverse set of mammals, and annotate personal, diploid human genomes. We demonstrate the resulting discovery of novel genes, isoforms, and structural variants-even in genomes as well studied as rat and the great apes-and how these annotations improve cross-species RNA expression experiments.

We report on four families affected by a clinical presentation of complex hereditary spastic paraplegia (HSP) due to recessive mutations in DDHD2, encoding one of the three mammalian intracellular phospholipases A(1) (iPLA(1)). The core phenotype of this HSP syndrome consists of very early-onset (<2 years) spastic paraplegia, intellectual disability, and a specific pattern of brain abnormalities on cerebral imaging. An essential role for DDHD2 in the human CNS, and perhaps more specifically in synaptic functioning, is supported by a reduced number of active zones at synaptic terminals in Ddhd-knockdown Drosophila models. All identified mutations affect the protein's DDHD domain, which is vital for its phospholipase activity. In line with the function of DDHD2 in lipid metabolism and its role in the CNS, an abnormal lipid peak indicating accumulation of lipids was detected with cerebral magnetic resonance spectroscopy, which provides an applicable diagnostic biomarker that can distinguish the DDHD2 phenotype from other complex HSP phenotypes. We show that mutations in DDHD2 cause a specific complex HSP subtype (SPG54), thereby linking a member of the PLA(1) family to human neurologic disease.

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

We sequenced the genomes of a ∼7,000-year-old farmer from Germany and eight ∼8,000-year-old hunter-gatherers from Luxembourg and Sweden. We analysed these and other ancient genomes with 2,345 contemporary humans to show that most present-day Europeans derive from at least three highly differentiated populations: west European hunter-gatherers, who contributed ancestry to all Europeans but not to Near Easterners; ancient north Eurasians related to Upper Palaeolithic Siberians, who contributed to both Europeans and Near Easterners; and early European farmers, who were mainly of Near Eastern origin but also harboured west European hunter-gatherer related ancestry. We model these populations' deep relationships and show that early European farmers had ∼44% ancestry from a 'basal Eurasian' population that split before the diversification of other non-African lineages.

Hematopoietic stem cells (HSCs) require multiple molecular inputs for proper specification, including activity of the Notch signaling pathway. A requirement for the Notch1 and dispensability of the Notch2 receptor has been demonstrated in mice, but the role of the remaining Notch receptors has not been investigated. Here, we demonstrate that three of the four Notch receptors are independently required for the specification of HSCs in the zebrafish. The orthologues of the murine Notch1 receptor, Notch1a and Notch1b, are each required intrinsically to fate HSCs, just prior to their emergence from aortic hemogenic endothelium. By contrast, the Notch3 receptor is required earlier within the developing somite to regulate HSC emergence in a non-cell-autonomous manner. Epistatic analyses demonstrate that Notch3 function lies downstream of Wnt16, which is required for HSC specification through its regulation of two Notch ligands, dlc and dld. Collectively, these findings demonstrate for the first time that multiple Notch signaling inputs are required to specify HSCs and that Notch3 performs a novel role within the somite to regulate the neighboring precursors of hemogenic endothelium.

Asthma is a complex genetic disease caused by a combination of genetic and environmental risk factors. We sought to test classes of genetic variants largely missed by genome-wide association studies (GWAS), including copy number variants (CNVs) and low-frequency variants, by performing whole-genome sequencing (WGS) on 16 individuals from asthma-enriched and asthma-depleted families. The samples were obtained from an extended 13-generation Hutterite pedigree with reduced genetic heterogeneity due to a small founding gene pool and reduced environmental heterogeneity as a result of a communal lifestyle. We sequenced each individual to an average depth of 13-fold, generated a comprehensive catalog of genetic variants, and tested the most severe mutations for association with asthma. We identified and validated 1960 CNVs, 19 nonsense or splice-site single nucleotide variants (SNVs), and 18 insertions or deletions that were out of frame. As follow-up, we performed targeted sequencing of 16 genes in 837 cases and 540 controls of Puerto Rican ancestry and found that controls carry a significantly higher burden of mutations in IL27RA (2.0% of controls; 0.23% of cases; nominal p = 0.004; Bonferroni p = 0.21). We also genotyped 593 CNVs in 1199 Hutterite individuals. We identified a nominally significant association (p = 0.03; Odds ratio (OR) = 3.13) between a 6 kbp deletion in an intron of NEDD4L and increased risk of asthma. We genotyped this deletion in an additional 4787 non-Hutterite individuals (nominal p = 0.056; OR = 1.69). NEDD4L is expressed in bronchial epithelial cells, and conditional knockout of this gene in the lung in mice leads to severe inflammation and mucus accumulation. Our study represents one of the early instances of applying WGS to complex disease with a large environmental component and demonstrates how WGS can identify risk variants, including CNVs and low-frequency variants, largely untested in GWAS.

Despite the high heritability of autism spectrum disorders (ASD), characterized by persistent deficits in social communication and interaction and restricted, repetitive patterns of behavior, interests or activities, a genetic diagnosis can be established in only a minority of patients. Known genetic causes include chromosomal aberrations, such as the duplication of the 15q11-13 region, and monogenic causes, as in Rett and fragile-X syndromes. The genetic heterogeneity within ASD is striking, with even the most frequent causes responsible for only 1% of cases at the most. Even with the recent developments in next-generation sequencing, for the large majority of cases no molecular diagnosis can be established. Here, we report ten patients with ASD and other shared clinical characteristics, including intellectual disability and facial dysmorphisms caused by a mutation in ADNP, a transcription factor involved in the SWI/SNF remodeling complex. We estimate this gene to be mutated in at least 0.17% of ASD cases, making it one of the most frequent ASD-associated genes known to date.

Recently, we identified in two individuals with intellectual disability (ID) different de novo mutations in DEAF1, which encodes a transcription factor with an important role in embryonic development. To ascertain whether these mutations in DEAF1 are causative for the ID phenotype, we performed targeted resequencing of DEAF1 in an additional cohort of over 2,300 individuals with unexplained ID and identified two additional individuals with de novo mutations in this gene. All four individuals had severe ID with severely affected speech development, and three showed severe behavioral problems. DEAF1 is highly expressed in the CNS, especially during early embryonic development. All four mutations were missense mutations affecting the SAND domain of DEAF1. Altered DEAF1 harboring any of the four amino acid changes showed impaired transcriptional regulation of the DEAF1 promoter. Moreover, behavioral studies in mice with a conditional knockout of Deaf1 in the brain showed memory deficits and increased anxiety-like behavior. Our results demonstrate that mutations in DEAF1 cause ID and behavioral problems, most likely as a result of impaired transcriptional regulation by DEAF1.

Complement factor H shows very strong association with Age-related Macular Degeneration (AMD), and recent data suggest that multiple causal variants are associated with disease. To refine the location of the disease associated variants, we characterized in detail the structural variation at CFH and its paralogs, including two copy number polymorphisms (CNP), CNP147 and CNP148, and several rare deletions and duplications. Examination of 34 AMD-enriched extended families (N = 293) and AMD cases (White N = 4210 Indian = 134; Malay = 140) and controls (White N = 3229; Indian = 117; Malay = 2390) demonstrated that deletion CNP148 was protective against AMD, independent of SNPs at CFH. Regression analysis of seven common haplotypes showed three haplotypes, H1, H6 and H7, as conferring risk for AMD development. Being the most common haplotype H1 confers the greatest risk by increasing the odds of AMD by 2.75-fold (95% CI = [2.51, 3.01]; p = 8.31×10(-109)); Caucasian (H6) and Indian-specific (H7) recombinant haplotypes increase the odds of AMD by 1.85-fold (p = 3.52×10(-9)) and by 15.57-fold (P = 0.007), respectively. We identified a 32-kb region downstream of Y402H (rs1061170), shared by all three risk haplotypes, suggesting that this region may be critical for AMD development. Further analysis showed that two SNPs within the 32 kb block, rs1329428 and rs203687, optimally explain disease association. rs1329428 resides in 20 kb unique sequence block, but rs203687 resides in a 12 kb block that is 89% similar to a noncoding region contained in ΔCNP148. We conclude that causal variation in this region potentially encompasses both regulatory effects at single markers and copy number.