A human embryonic limb cell atlas resolved in space and time.

Human limbs emerge during the fourth post-conception week as mesenchymal buds, which develop into fully formed limbs over the subsequent months1. This process is orchestrated by numerous temporally and spatially restricted gene expression programmes, making congenital alterations in phenotype common2. Decades of work with model organisms have defined the fundamental mechanisms underlying vertebrate limb development, but an in-depth characterization of this process in humans has yet to be performed. Here we detail human embryonic limb development across space and time using single-cell and spatial transcriptomics. We demonstrate extensive diversification of cells from a few multipotent progenitors to myriad differentiated cell states, including several novel cell populations. We uncover two waves of human muscle development, each characterized by different cell states regulated by separate gene expression programmes, and identify musculin (MSC) as a key transcriptional repressor maintaining muscle stem cell identity. Through assembly of multiple anatomically continuous spatial transcriptomic samples using VisiumStitcher, we map cells across a sagittal section of a whole fetal hindlimb. We reveal a clear anatomical segregation between genes linked to brachydactyly and polysyndactyly, and uncover transcriptionally and spatially distinct populations of the mesenchyme in the autopod. Finally, we perform single-cell RNA sequencing on mouse embryonic limbs to facilitate cross-species developmental comparison, finding substantial homology between the two species.

The Cousa objective: a long-working distance air objective for multiphoton imaging in vivo.

Multiphoton microscopy can resolve fluorescent structures and dynamics deep in scattering tissue and has transformed neural imaging, but applying this technique in vivo can be limited by the mechanical and optical constraints of conventional objectives. Short working distance objectives can collide with compact surgical windows or other instrumentation and preclude imaging. Here we present an ultra-long working distance (20 mm) air objective called the Cousa objective. It is optimized for performance across multiphoton imaging wavelengths, offers a more than 4 mm2 field of view with submicrometer lateral resolution and is compatible with commonly used multiphoton imaging systems. A novel mechanical design, wider than typical microscope objectives, enabled this combination of specifications. We share the full optical prescription, and report performance including in vivo two-photon and three-photon imaging in an array of species and preparations, including nonhuman primates. The Cousa objective can enable a range of experiments in neuroscience and beyond.

Cortical response selectivity derives from strength in numbers of synapses.

Single neocortical neurons are driven by populations of excitatory inputs, which form the basis of neuronal selectivity to features of sensory input. Excitatory connections are thought to mature during development through activity-dependent Hebbian plasticity1, whereby similarity between presynaptic and postsynaptic activity selectively strengthens some synapses and weakens others2. Evidence in support of this process includes measurements of synaptic ultrastructure and in vitro and in vivo physiology and imaging studies3-8. These corroborating lines of evidence lead to the prediction that a small number of strong synaptic inputs drive neuronal selectivity, whereas weak synaptic inputs are less correlated with the somatic output and modulate activity overall6,7. Supporting evidence from cortical circuits, however, has been limited to measurements of neighbouring, connected cell pairs, raising the question of whether this prediction holds for a broad range of synapses converging onto cortical neurons. Here we measure the strengths of functionally characterized excitatory inputs contacting single pyramidal neurons in ferret primary visual cortex (V1) by combining in vivo two-photon synaptic imaging and post hoc electron microscopy. Using electron microscopy reconstruction of individual synapses as a metric of strength, we find no evidence that strong synapses have a predominant role in the selectivity of cortical neuron responses to visual stimuli. Instead, selectivity appears to arise from the total number of synapses activated by different stimuli. Moreover, spatial clustering of co-active inputs appears to be reserved for weaker synapses, enhancing the contribution of weak synapses to somatic responses. Our results challenge the role of Hebbian mechanisms in shaping neuronal selectivity in cortical circuits, and suggest that selectivity reflects the co-activation of large populations of presynaptic neurons with similar properties and a mixture of strengths.

Genomic Diagnosis of Rare Pediatric Disease in the United Kingdom and Ireland.

BACKGROUND: Pediatric disorders include a range of highly penetrant, genetically heterogeneous conditions amenable to genomewide diagnostic approaches. Finding a molecular diagnosis is challenging but can have profound lifelong benefits. METHODS: We conducted a large-scale sequencing study involving more than 13,500 families with probands with severe, probably monogenic, difficult-to-diagnose developmental disorders from 24 regional genetics services in the United Kingdom and Ireland. Standardized phenotypic data were collected, and exome sequencing and microarray analyses were performed to investigate novel genetic causes. We developed an iterative variant analysis pipeline and reported candidate variants to clinical teams for validation and diagnostic interpretation to inform communication with families. Multiple regression analyses were performed to evaluate factors affecting the probability of diagnosis. RESULTS: A total of 13,449 probands were included in the analyses. On average, we reported 1.0 candidate variant per parent-offspring trio and 2.5 variants per singleton proband. Using clinical and computational approaches to variant classification, we made a diagnosis in approximately 41% of probands (5502 of 13,449). Of 3599 probands in trios who received a diagnosis by clinical assertion, approximately 76% had a pathogenic de novo variant. Another 22% of probands (2997 of 13,449) had variants of uncertain significance in genes that were strongly linked to monogenic developmental disorders. Recruitment in a parent-offspring trio had the largest effect on the probability of diagnosis (odds ratio, 4.70; 95% confidence interval [CI], 4.16 to 5.31). Probands were less likely to receive a diagnosis if they were born extremely prematurely (i.e., 22 to 27 weeks' gestation; odds ratio, 0.39; 95% CI, 0.22 to 0.68), had in utero exposure to antiepileptic medications (odds ratio, 0.44; 95% CI, 0.29 to 0.67), had mothers with diabetes (odds ratio, 0.52; 95% CI, 0.41 to 0.67), or were of African ancestry (odds ratio, 0.51; 95% CI, 0.31 to 0.78). CONCLUSIONS: Among probands with severe, probably monogenic, difficult-to-diagnose developmental disorders, multimodal analysis of genomewide data had good diagnostic power, even after previous attempts at diagnosis. (Funded by the Health Innovation Challenge Fund and Wellcome Sanger Institute.).

Characterization of an eye field-like state during optic vesicle organoid development.

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors. This crucial event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

Evidence for 28 genetic disorders discovered by combining healthcare and research data.

De novo mutations in protein-coding genes are a well-established cause of developmental disorders1. However, genes known to be associated with developmental disorders account for only a minority of the observed excess of such de novo mutations1,2. Here, to identify previously undescribed genes associated with developmental disorders, we integrate healthcare and research exome-sequence data from 31,058 parent-offspring trios of individuals with developmental disorders, and develop a simulation-based statistical test to identify gene-specific enrichment of de novo mutations. We identified 285 genes that were significantly associated with developmental disorders, including 28 that had not previously been robustly associated with developmental disorders. Although we detected more genes associated with developmental disorders, much of the excess of de novo mutations in protein-coding genes remains unaccounted for. Modelling suggests that more than 1,000 genes associated with developmental disorders have not yet been described, many of which are likely to be less penetrant than the currently known genes. Research access to clinical diagnostic datasets will be critical for completing the map of genes associated with developmental disorders.

Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants.

Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.

Short-read whole genome sequencing identifies causative variants in most individuals with previously unexplained aniridia.

BACKGROUND: Classic aniridia is a highly penetrant autosomal dominant disorder characterised by congenital absence of the iris, foveal hypoplasia, optic disc anomalies and progressive opacification of the cornea. >90% of cases of classic aniridia are caused by heterozygous, loss-of-function variants affecting the PAX6 locus. METHODS: Short-read whole genome sequencing was performed on 51 (39 affected) individuals from 37 different families who had screened negative for mutations in the PAX6 coding region. RESULTS: Likely causative mutations were identified in 22 out of 37 (59%) families. In 19 out of 22 families, the causative genomic changes have an interpretable deleterious impact on the PAX6 locus. Of these 19 families, 1 has a novel heterozygous PAX6 frameshift variant missed on previous screens, 4 have single nucleotide variants (SNVs) (one novel) affecting essential splice sites of PAX6 5' non-coding exons and 2 have deep intronic SNV (one novel) resulting in gain of a donor splice site. In 12 out of 19, the causative variants are large-scale structural variants; 5 have partial or whole gene deletions of PAX6, 3 have deletions encompassing critical PAX6 cis-regulatory elements, 2 have balanced inversions with disruptive breakpoints within the PAX6 locus and 2 have complex rearrangements disrupting PAX6. The remaining 3 of 22 families have deletions encompassing FOXC1 (a known cause of atypical aniridia). Seven of the causative variants occurred de novo and one cosegregated with familial aniridia. We were unable to establish inheritance status in the remaining probands. No plausibly causative SNVs were identified in PAX6 cis-regulatory elements. CONCLUSION: Whole genome sequencing proves to be an effective diagnostic test in most individuals with previously unexplained aniridia.

F-box receptor mediated control of substrate stability and subcellular location organizes cellular development of Aspergillus nidulans.

Fungal growth and development are coordinated with specific secondary metabolism. This coordination requires 8 of 74 F-box proteins of the filamentous fungus Aspergillus nidulans. F-box proteins recognize primed substrates for ubiquitination by Skp1-Cul1-Fbx (SCF) E3 ubiquitin RING ligases and degradation by the 26S proteasome. 24 F-box proteins are found in the nuclear fraction as part of SCFs during vegetative growth. 43 F-box proteins interact with SCF proteins during growth, development or stress. 45 F-box proteins are associated with more than 700 proteins that have mainly regulatory roles. This corroborates that accurate surveillance of protein stability is prerequisite for organizing multicellular fungal development. Fbx23 combines subcellular location and protein stability control, illustrating the complexity of F-box mediated regulation during fungal development. Fbx23 interacts with epigenetic methyltransferase VipC which interacts with fungal NF-κB-like velvet domain regulator VeA that coordinates fungal development with secondary metabolism. Fbx23 prevents nuclear accumulation of methyltransferase VipC during early development. These results suggest that in addition to their role in protein degradation, F-box proteins also control subcellular accumulations of key regulatory proteins for fungal development.

Detecting cryptic clinically relevant structural variation in exome-sequencing data increases diagnostic yield for developmental disorders.

Structural variation (SV) describes a broad class of genetic variation greater than 50 bp in size. SVs can cause a wide range of genetic diseases and are prevalent in rare developmental disorders (DDs). Individuals presenting with DDs are often referred for diagnostic testing with chromosomal microarrays (CMAs) to identify large copy-number variants (CNVs) and/or with single-gene, gene-panel, or exome sequencing (ES) to identify single-nucleotide variants, small insertions/deletions, and CNVs. However, individuals with pathogenic SVs undetectable by conventional analysis often remain undiagnosed. Consequently, we have developed the tool InDelible, which interrogates short-read sequencing data for split-read clusters characteristic of SV breakpoints. We applied InDelible to 13,438 probands with severe DDs recruited as part of the Deciphering Developmental Disorders (DDD) study and discovered 63 rare, damaging variants in genes previously associated with DDs missed by standard SNV, indel, or CNV discovery approaches. Clinical review of these 63 variants determined that about half (30/63) were plausibly pathogenic. InDelible was particularly effective at ascertaining variants between 21 and 500 bp in size and increased the total number of potentially pathogenic variants identified by DDD in this size range by 42.9%. Of particular interest were seven confirmed de novo variants in MECP2, which represent 35.0% of all de novo protein-truncating variants in MECP2 among DDD study participants. InDelible provides a framework for the discovery of pathogenic SVs that are most likely missed by standard analytical workflows and has the potential to improve the diagnostic yield of ES across a broad range of genetic diseases.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome.

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Paediatric genomics: diagnosing rare disease in children.

The majority of rare diseases affect children, most of whom have an underlying genetic cause for their condition. However, making a molecular diagnosis with current technologies and knowledge is often still a challenge. Paediatric genomics is an immature but rapidly evolving field that tackles this issue by incorporating next-generation sequencing technologies, especially whole-exome sequencing and whole-genome sequencing, into research and clinical workflows. This complex multidisciplinary approach, coupled with the increasing availability of population genetic variation data, has already resulted in an increased discovery rate of causative genes and in improved diagnosis of rare paediatric disease. Importantly, for affected families, a better understanding of the genetic basis of rare disease translates to more accurate prognosis, management, surveillance and genetic advice; stimulates research into new therapies; and enables provision of better support.

Paediatric genomics: diagnosing rare disease in children.

This corrects the article DOI: 10.1038/nrg.2017.116.

Diagnosis and management of Cornelia de Lange syndrome: first international consensus statement.

Cornelia de Lange syndrome (CdLS) is an archetypical genetic syndrome that is characterized by intellectual disability, well-defined facial features, upper limb anomalies and atypical growth, among numerous other signs and symptoms. It is caused by variants in any one of seven genes, all of which have a structural or regulatory function in the cohesin complex. Although recent advances in next-generation sequencing have improved molecular diagnostics, marked heterogeneity exists in clinical and molecular diagnostic approaches and care practices worldwide. Here, we outline a series of recommendations that document the consensus of a group of international experts on clinical diagnostic criteria, both for classic CdLS and non-classic CdLS phenotypes, molecular investigations, long-term management and care planning.

Differential tuning of excitation and inhibition shapes direction selectivity in ferret visual cortex.

To encode specific sensory inputs, cortical neurons must generate selective responses for distinct stimulus features. In principle, a variety of factors can contribute to the response selectivity of a cortical neuron: the tuning and strength of excitatory1-3 and inhibitory synaptic inputs4-6, dendritic nonlinearities7-9 and spike threshold10,11. Here we use a combination of techniques including in vivo whole-cell recording, synaptic- and cellular-resolution in vivo two-photon calcium imaging, and GABA (γ-aminobutyric acid) neuron-selective optogenetic manipulation to dissect the factors that contribute to the direction-selective responses of layer 2/3 neurons in ferret visual cortex (V1). Two-photon calcium imaging of dendritic spines12,13 revealed that each neuron receives a mixture of excitatory synaptic inputs selective for the somatic preferred or null direction of motion. The relative number of preferred- and null-tuned excitatory inputs predicted a neuron's somatic direction preference, but failed to account for the degree of direction selectivity. By contrast, in vivo whole-cell patch-clamp recordings revealed a notable degree of direction selectivity in subthreshold responses that was significantly correlated with spiking direction selectivity. Subthreshold direction selectivity was predicted by the magnitude and variance of the response to the null direction of motion, and several lines of evidence, including conductance measurements, demonstrate that differential tuning of excitation and inhibition suppresses responses to the null direction of motion. Consistent with this idea, optogenetic inactivation of GABAergic neurons in layer 2/3 reduced direction selectivity by enhancing responses to the null direction. Furthermore, by optogenetically mapping connections of inhibitory neurons in layer 2/3 in vivo, we find that layer 2/3 inhibitory neurons make long-range, intercolumnar projections to excitatory neurons that prefer the opposite direction of motion. We conclude that intracortical inhibition exerts a major influence on the degree of direction selectivity in layer 2/3 of ferret V1 by suppressing responses to the null direction of motion.

De novo mutations in regulatory elements in neurodevelopmental disorders.

We previously estimated that 42% of patients with severe developmental disorders carry pathogenic de novo mutations in coding sequences. The role of de novo mutations in regulatory elements affecting genes associated with developmental disorders, or other genes, has been essentially unexplored. We identified de novo mutations in three classes of putative regulatory elements in almost 8,000 patients with developmental disorders. Here we show that de novo mutations in highly evolutionarily conserved fetal brain-active elements are significantly and specifically enriched in neurodevelopmental disorders. We identified a significant twofold enrichment of recurrently mutated elements. We estimate that, genome-wide, 1-3% of patients without a diagnostic coding variant carry pathogenic de novo mutations in fetal brain-active regulatory elements and that only 0.15% of all possible mutations within highly conserved fetal brain-active elements cause neurodevelopmental disorders with a dominant mechanism. Our findings represent a robust estimate of the contribution of de novo mutations in regulatory elements to this genetically heterogeneous set of disorders, and emphasize the importance of combining functional and evolutionary evidence to identify regulatory causes of genetic disorders.

Common genetic variants contribute to risk of rare severe neurodevelopmental disorders.

There are thousands of rare human disorders that are caused by single deleterious, protein-coding genetic variants1. However, patients with the same genetic defect can have different clinical presentations2-4, and some individuals who carry known disease-causing variants can appear unaffected5. Here, to understand what explains these differences, we study a cohort of 6,987 children assessed by clinical geneticists to have severe neurodevelopmental disorders such as global developmental delay and autism, often in combination with abnormalities of other organ systems. Although the genetic causes of these neurodevelopmental disorders are expected to be almost entirely monogenic, we show that 7.7% of variance in risk is attributable to inherited common genetic variation. We replicated this genome-wide common variant burden by showing, in an independent sample of 728 trios (comprising a child plus both parents) from the same cohort, that this burden is over-transmitted from parents to children with neurodevelopmental disorders. Our common-variant signal is significantly positively correlated with genetic predisposition to lower educational attainment, decreased intelligence and risk of schizophrenia. We found that common-variant risk was not significantly different between individuals with and without a known protein-coding diagnostic variant, which suggests that common-variant risk affects patients both with and without a monogenic diagnosis. In addition, previously published common-variant scores for autism, height, birth weight and intracranial volume were all correlated with these traits within our cohort, which suggests that phenotypic expression in individuals with monogenic disorders is affected by the same variants as in the general population. Our results demonstrate that common genetic variation affects both overall risk and clinical presentation in neurodevelopmental disorders that are typically considered to be monogenic.

EyeG2P: an automated variant filtering approach improves efficiency of diagnostic genomic testing for inherited ophthalmic disorders.

BACKGROUND: Genomic variant prioritisation is one of the most significant bottlenecks to mainstream genomic testing in healthcare. Tools to improve precision while ensuring high recall are critical to successful mainstream clinical genomic testing, in particular for whole genome sequencing where millions of variants must be considered for each patient. METHODS: We developed EyeG2P, a publicly available database and web application using the Ensembl Variant Effect Predictor. EyeG2P is tailored for efficient variant prioritisation for individuals with inherited ophthalmic conditions. We assessed the sensitivity of EyeG2P in 1234 individuals with a broad range of eye conditions who had previously received a confirmed molecular diagnosis through routine genomic diagnostic approaches. For a prospective cohort of 83 individuals, we assessed the precision of EyeG2P in comparison with routine diagnostic approaches. For 10 additional individuals, we assessed the utility of EyeG2P for whole genome analysis. RESULTS: EyeG2P had 99.5% sensitivity for genomic variants previously identified as clinically relevant through routine diagnostic analysis (n=1234 individuals). Prospectively, EyeG2P enabled a significant increase in precision (35% on average) in comparison with routine testing strategies (p<0.001). We demonstrate that incorporation of EyeG2P into whole genome sequencing analysis strategies can reduce the number of variants for analysis to six variants, on average, while maintaining high diagnostic yield. CONCLUSION: Automated filtering of genomic variants through EyeG2P can increase the efficiency of diagnostic testing for individuals with a broad range of inherited ophthalmic disorders.

Genomically Aided Diagnosis of Severe Developmental Disorders.

Our ability to make accurate and specific genetic diagnoses in individuals with severe developmental disorders has been transformed by data derived from genomic sequencing technologies. These data reveal both the patterns and rates of different mutational mechanisms and identify regions of the human genome with fewer mutations than would be expected. In outbred populations, the most common identifiable cause of severe developmental disorders is de novo mutation affecting the coding region in one of approximately 500 different genes, almost universally showing constraint. Simply combining the location of a de novo genomic event with its predicted consequence on the gene product gives significant diagnostic power. Our knowledge of the diversity of phenotypic consequences associated with comparable diagnostic genotypes at each locus is improving. Computationally useful phenotype data will improve diagnostic interpretation of ultrarare genetic variants and, in the long run, indicate which specific embryonic processes have been perturbed.

Efficacy and toxicity of primary re-irradiation for malignant spinal cord compression based on radiobiological modelling: a phase II clinical trial.

BACKGROUND: The efficacy and safety of primary re-irradiation for MSCC are not known. Our aim was to establish the efficacy and safety of biologically effective dose-based re-irradiation. METHODS: Patients presenting with MSCC at a previously irradiated spine segment, and not proceeding with surgical decompression, were eligible. A 3 Gray per fraction experimental schedule (minimum 18 Gy/6 fractions, maximum 30 Gy/10 fractions) was used, delivering a maximum cumulative spinal dose of 100 Gy2 if the interval since the last radiotherapy was within 6 months, or 130 Gy2 if longer. The primary outcome was a change in mobility from week 1 to week 5 post-treatment, as assessed by the Tomita score. The RTOG SOMA score was used to screen for spinal toxicity, and an MRI performed to assess for radiation-induced myelopathy (RIM). RESULTS: Twenty-two patients were enroled, of whom eleven were evaluable for the primary outcome. Nine of eleven (81.8%) had stable or improved Tomita scores at 5 weeks. One of eight (12.5%) evaluable for late toxicity developed RIM. CONCLUSIONS: Re-irradiation is an efficacious treatment for MSCC. There is a risk of RIM with a cumulative dose of 120 Gy2. CLINICAL TRIAL REGISTRATION: Cancer Trials Ireland (ICORG 07-11); NCT00974168.