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.

Single-cell analyses reveal transient retinal progenitor cells in the ciliary margin of developing human retina.

The emergence of retinal progenitor cells and differentiation to various retinal cell types represent fundamental processes during retinal development. Herein, we provide a comprehensive single cell characterisation of transcriptional and chromatin accessibility changes that underline retinal progenitor cell specification and differentiation over the course of human retinal development up to midgestation. Our lineage trajectory data demonstrate the presence of early retinal progenitors, which transit to late, and further to transient neurogenic progenitors, that give rise to all the retinal neurons. Combining single cell RNA-Seq with spatial transcriptomics of early eye samples, we demonstrate the transient presence of early retinal progenitors in the ciliary margin zone with decreasing occurrence from 8 post-conception week of human development. In retinal progenitor cells, we identified a significant enrichment for transcriptional enhanced associate domain transcription factor binding motifs, which when inhibited led to loss of cycling progenitors and retinal identity in pluripotent stem cell derived organoids.

The effect of a minimum price per unit of alcohol in Scotland on alcohol-related ambulance call-outs: A controlled interrupted time-series analysis.

BACKGROUND AND AIMS: On 1 May 2018, Scotland introduced a minimum unit price (MUP) of £0.50 for alcohol, with one UK unit of alcohol being 10 ml of pure ethanol. This study measured the association between MUP and changes in the volume of alcohol-related ambulance call-outs in the overall population and in call-outs subsets (night-time call-outs and subpopulations with higher incidence of alcohol-related harm). DESIGN: An interrupted time-series (ITS) was used to measure variations in the daily volume of alcohol-related call-outs. We performed uncontrolled ITS on both the intervention and control group and a controlled ITS built on the difference between the two series. Data were from electronic patient clinical records from the Scottish Ambulance Service. SETTING AND CASES: Alcohol-related ambulance call-outs (intervention group) and total ambulance call-outs for people aged under 13 years (control group) in Scotland, from December 2017 to March 2020. MEASUREMENTS: Call-outs were deemed alcohol-related if ambulance clinicians indicated that alcohol was a 'contributing factor' in the call-out and/or a validated Scottish Ambulance Service algorithm determined that the call-out was alcohol-related. FINDINGS: No statistically significant association in the volume of call-outs was found in both the uncontrolled series [step change = 0.062, 95% confidence interval (CI) = -0.012, 0.0135 P = 0.091; slope change = -0.001, 95% CI = -0.001, 0.1 × 10-3 P = 0.139] and controlled series (step change = -0.01, 95% CI = -0.317, 0.298 P = 0.951; slope change = -0.003, 95% CI = -0.008, 0.002 P = 0.257). Similarly, no significant changes were found for the night-time series or for any population subgroups. CONCLUSIONS: There appears to be no statistically significant association between the introduction of minimum unit pricing for alcohol in Scotland and the volume of alcohol-related ambulance call-outs. This was observed overall, across subpopulations and at night-time.

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.

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.

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.

Gliotoxin-mediated bacterial growth inhibition is caused by specific metal ion depletion.

Overcoming antimicrobial resistance represents a formidable challenge and investigating bacterial growth inhibition by fungal metabolites may yield new strategies. Although the fungal non-ribosomal peptide gliotoxin (GT) is known to exhibit antibacterial activity, the mechanism(s) of action are unknown, although reduced gliotoxin (dithiol gliotoxin; DTG) is a zinc chelator. Furthermore, it has been demonstrated that GT synergises with vancomycin to inhibit growth of Staphylococcus aureus. Here we demonstrate, without precedent, that GT-mediated growth inhibition of both Gram positive and negative bacterial species is reversed by Zn2+ or Cu2+ addition. Both GT, and the known zinc chelator TPEN, mediate growth inhibition of Enterococcus faecalis which is reversed by zinc addition. Moreover, zinc also reverses the synergistic growth inhibition of E. faecalis observed in the presence of both GT and vancomycin (4 µg/ml). As well as zinc chelation, DTG also appears to chelate Cu2+, but not Mn2+ using a 4-(2-pyridylazo)resorcinol assay system and Zn2+ as a positive control. DTG also specifically reacts in Fe3+-containing Siderotec™ assays, most likely by Fe3+ chelation from test reagents. GSH or DTT show no activity in these assays. Confirmatory high resolution mass spectrometry, in negative ion mode, confirmed, for the first time, the presence of both Cu[DTG] and Fe[DTG]2 chelates. Label free quantitative proteomic analysis further revealed major intracellular proteomic remodelling within E. faecalis in response to GT exposure for 30-180 min. Globally, 4.2-7.2% of detectable proteins exhibited evidence of either unique presence/increased abundance or unique absence/decreased abundance (n = 994-1160 total proteins detected), which is the first demonstration that GT affects the bacterial proteome in general, and E. faecalis, specifically. Unique detection of components of the AdcABC and AdcA-II zinc uptake systems was observed, along with apparent ribosomal reprofiling to zinc-free paralogs in the presence of GT. Overall, we hypothesise that GT-mediated bacterial growth inhibition appears to involve intracellular zinc depletion or reduced bioavailability, and based on in vitro chelate formation, may also involve dysregulation of Cu2+ homeostasis.

An adaptive behavioral control motif mediated by cortical axo-axonic inhibition.

Genetically defined subgroups of inhibitory interneurons are thought to play distinct roles in learning, but heterogeneity within these subgroups has limited our understanding of the scope and nature of their specific contributions. Here we reveal that the chandelier cell (ChC), an interneuron type that specializes in inhibiting the axon-initial segment (AIS) of pyramidal neurons, establishes cortical microcircuits for organizing neural coding through selective axo-axonic synaptic plasticity. We found that organized motor control is mediated by enhanced population coding of direction-tuned premotor neurons, with tuning refined through suppression of irrelevant neuronal activity. ChCs contribute to learning-dependent refinements by providing selective inhibitory control over individual pyramidal neurons rather than global suppression. Quantitative analysis of structural plasticity across axo-axonic synapses revealed that ChCs redistributed inhibitory weights to individual pyramidal neurons during learning. These results demonstrate an adaptive logic of the inhibitory circuit motif responsible for organizing distributed neural representations. Thus, ChCs permit efficient cortical computation in a targeted cell-specific manner.

Developmental alignment of feedforward inputs and recurrent network activity drives increased response selectivity and reliability in primary visual cortex following the onset of visual experience.

Selective and reliable cortical sensory representations depend on synaptic interactions between feedforward inputs, conveying information from lower levels of the sensory pathway, and recurrent networks that reciprocally connect neurons functioning at the same hierarchical level. Here we explore the development of feedforward/recurrent interactions in primary visual cortex of the ferret that is responsible for the representation of orientation, focusing on the feedforward inputs from cortical layer 4 and its relation to the modular recurrent network in layer 2/3 before and after the onset of visual experience. Using simultaneous laminar electrophysiology and calcium imaging we found that in experienced animals, individual layer 4 and layer 2/3 neurons exhibit strongly correlated responses with the modular recurrent network structure in layer 2/3. Prior to experience, layer 2/3 neurons exhibit comparable modular correlation structure, but this correlation structure is missing for individual layer 4 neurons. Further analysis of the receptive field properties of layer 4 neurons in naïve animals revealed that they exhibit very poor orientation tuning compared to layer 2/3 neurons at this age, and this is accompanied by the lack of spatial segregation of ON and OFF subfields, the definitive property of layer 4 simple cells in experienced animals. Analysis of the response dynamics of layer 2/3 neurons with whole-cell patch recordings confirms that individual layer 2/3 neurons in naïve animals receive poorly-selective feedforward input that does not align with the orientation preference of the layer 2/3 responses. Further analysis reveals that the misaligned feedforward input is the underlying cause of reduced selectivity and increased response variability that is evident in the layer 2/3 responses of naïve animals. Altogether, our experiments indicate that the onset of visual experience is accompanied by a critical refinement in the responses of layer 4 neurons and the alignment of feedforward and recurrent networks that increases the selectivity and reliability of the representation of orientation in V1.

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.

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.).

Recombinant production, characterization and industrial application testing of a novel acidic exo/endo-chitinase from Rasamsonia emersonii.

An acid-active exo/endo-chitinase; comprising a GH18 catalytic domain and substrate insertion domain; originating from the thermophilic filamentous fungus Rasamsonia emersonii, was expressed in Pichia pastoris. In silico analysis including phylogenetic analysis, and recombinant production, purification, biochemical characterisation, and industrial application testing, was carried out. The expressed protein was identified by SDS-PAGE as a smear from 56.3 to 125.1 kDa, which sharpens into bands at 46.0 kDa, 48.4 kDa and a smear above 60 kDa when treated with PNGase F. The acid-active chitinase was primarily a chitobiosidase but displayed some endo-chitinase and acetyl-glucosamidase activity. The enzyme was optimally active at 50 °C, and markedly low pH of 2.8. As far as the authors are aware, this is the lowest pH optima reported for any fungal chitinase. The acid-active chitinase likely plays a role in chitin degradation for cell uptake in its native environment, perhaps in conjunction with a chitin deacetylase. Comparative studies with other R. emersonii chitinases indicate that they may play a synergistic role in this. The acid-active chitinase displayed some efficacy against non-treated substrates; fungal chitin and chitin from shrimp. Thus, it may be suited to industrial chitin hydrolysis reactions for extraction of glucosamine and chitobiose at low pH.

An adaptive behavioral control motif mediated by cortical axo-axonic inhibition.

Neural circuits are reorganized with specificity during learning. Genetically-defined subgroups of inhibitory interneurons are thought to play distinct roles in learning, but heterogeneity within these subgroups has limited our understanding of the scope and nature of their specific contributions to learning. Here we reveal that the chandelier cell (ChC), an interneuron type that specializes in inhibiting the axon-initial segment (AIS) of pyramidal neurons, establishes cortical microcircuits for organizing neural coding through selective axo-axonic synaptic plasticity. We find that organized motor control is mediated by enhanced population coding of direction-tuned premotor neurons, whose tuning is refined through suppression of irrelevant neuronal activity. ChCs are required for learning-dependent refinements via providing selective inhibitory control over pyramidal neurons rather than global suppression. Quantitative analysis on structural plasticity of axo-axonic synapses revealed that ChCs redistributed inhibitory weights to individual pyramidal neurons during learning. These results demonstrate an adaptive logic of the inhibitory circuit motif responsible for organizing distributed neural representations. Thus, ChCs permit efficient cortical computation in a target cell specific manner, which highlights the significance of interneuron diversity.

Comparative genomic study of the Penicillium genus elucidates a diverse pangenome and 15 lateral gene transfer events.

The Penicillia are known to produce a wide range natural products-some with devastating outcome for the agricultural industry and others with unexploited potential in different applications. However, a large-scale overview of the biosynthetic potential of different species has been lacking. In this study, we sequenced 93 Penicillium isolates and, together with eleven published genomes that hold similar assembly characteristics, we established a species phylogeny as well as defining a Penicillium pangenome. A total of 5612 genes were shared between ≥ 98 isolates corresponding to approximately half of the average number of genes a Penicillium genome holds. We further identified 15 lateral gene transfer events that have occurred in this collection of Penicillium isolates, which might have played an important role, such as niche adaption, in the evolution of these fungi. The comprehensive characterization of the genomic diversity in the Penicillium genus supersedes single-reference genomes, which do not necessarily capture the entire genetic variation.

Expanding the neurodevelopmental phenotype associated with HK1 de novo heterozygous missense variants.

Neurodevelopmental disorder with visual defects and brain anomalies (NEDVIBA) is a recently described genetic condition caused by de novo missense HK1 variants. Phenotypic data is currently limited; only seven patients have been published to date. This descriptive case series of a further four patients with de novo missense HK1 variants, alongside integration of phenotypic data with the reported cases, aims to improve our understanding of the associated phenotype. We provide further evidence that de novo HK1 variants located within the regulatory-terminal domain and alpha helix are associated with neurological problems and visual problems. We highlight for the first time an association with a raised cerebrospinal fluid lactate and specific abnormalities to the basal ganglia on brain magnetic resonance imaging, as well as associated respiratory issues and swallowing/feeding difficulties. We propose that this distinctive neurodevelopmental phenotype could arise through disruption of the regulatory glucose-6-phosphate binding site and subsequent gain of function of HK1 within the brain.

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.

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.

IMPROVE-DD: Integrating multiple phenotype resources optimizes variant evaluation in genetically determined developmental disorders.

Diagnosing rare developmental disorders using genome-wide sequencing data commonly necessitates review of multiple plausible candidate variants, often using ontologies of categorical clinical terms. We show that Integrating Multiple Phenotype Resources Optimizes Variant Evaluation in Developmental Disorders (IMPROVE-DD) by incorporating additional classes of data commonly available to clinicians and recorded in health records. In doing so, we quantify the distinct contributions of sex, growth, and development in addition to Human Phenotype Ontology (HPO) terms and demonstrate added value from these readily available information sources. We use likelihood ratios for nominal and quantitative data and propose a classifier for HPO terms in this framework. This Bayesian framework results in more robust diagnoses. Using data systematically collected in the Deciphering Developmental Disorders study, we considered 77 genes with pathogenic/likely pathogenic variants in ≥10 individuals. All genes showed at least a satisfactory prediction by receiver operating characteristic when testing on training data (AUC ≥ 0.6), and HPO terms were the best predictor for the majority of genes, though a minority (13/77) of genes were better predicted by other phenotypic data types. Overall, classifiers based upon multiple integrated phenotypic data sources performed better than those based upon any individual source, and importantly, integrated models produced notably fewer false positives. Finally, we show that IMPROVE-DD models with good predictive performance on cross-validation can be constructed from relatively few individuals. This suggests new strategies for candidate gene prioritization and highlights the value of systematic clinical data collection to support diagnostic programs.