Multi-Omic Graph Diagnosis (MOGDx): A data integration tool to perform classification tasks for heterogeneous diseases.

MOTIVATION: Heterogeneity in human diseases presents challenges in diagnosis and treatments due to the broad range of manifestations and symptoms. With the rapid development of labelled multi-omic data, integrative machine learning methods have achieved breakthroughs in treatments by redefining these diseases at a more granular level. These approaches often have limitations in scalability, oversimplification, and handling of missing data. RESULTS: In this study, we introduce Multi-Omic Graph Diagnosis (MOGDx), a flexible command line tool for the integration of multi-omic data to perform classification tasks for heterogeneous diseases. MOGDx has a network taxonomy. It fuses patient similarity networks, augments this integrated network with a reduced vector representation of genomic data and performs classification using a graph convolutional network. MOGDx was evaluated on three datasets from the cancer genome atlas for breast invasive carcinoma, kidney cancer, and low grade glioma. MOGDx demonstrated state-of-the-art performance and an ability to identify relevant multi-omic markers in each task. It integrated more genomic measures with greater patient coverage compared to other network integrative methods. Overall, MOGDx is a promising tool for integrating multi-omic data, classifying heterogeneous diseases, and aiding interpretation of genomic marker data. AVAILABILITY: MOGDx source code is available from https://github.com/biomedicalinformaticsgroup/MOGDxhttps://github.com/biomedicalinformaticsgroup/MOGDx. SUPPLEMENTARY: Supplementary material is available in the accompanying file SupplementaryMaterial.pdf.

A phospho-dawn of protein modification anticipates light onset in the picoeukaryote Ostreococcus tauri.

Diel regulation of protein levels and protein modification had been less studied than transcript rhythms. Here, we compare transcriptome data under light-dark cycles with partial proteome and phosphoproteome data, assayed using shotgun MS, from the alga Ostreococcus tauri, the smallest free-living eukaryote. A total of 10% of quantified proteins but two-thirds of phosphoproteins were rhythmic. Mathematical modelling showed that light-stimulated protein synthesis can account for the observed clustering of protein peaks in the daytime. Prompted by night-peaking and apparently dark-stable proteins, we also tested cultures under prolonged darkness, where the proteome changed less than under the diel cycle. Among the dark-stable proteins were prasinophyte-specific sequences that were also reported to accumulate when O. tauri formed lipid droplets. In the phosphoproteome, 39% of rhythmic phospho-sites reached peak levels just before dawn. This anticipatory phosphorylation suggests that a clock-regulated phospho-dawn prepares green cells for daytime functions. Acid-directed and proline-directed protein phosphorylation sites were regulated in antiphase, implicating the clock-related casein kinases 1 and 2 in phase-specific regulation, alternating with the CMGC protein kinase family. Understanding the dynamic phosphoprotein network should be facilitated by the minimal kinome and proteome of O. tauri. The data are available from ProteomeXchange, with identifiers PXD001734, PXD001735, and PXD002909.

Presence of DNA from Chlamydia-like organisms in the nasal cavities of grey seal pups (Halichoerus grypus) and three different substrates present in a breeding colony.

BACKGROUND: Chlamydia-like organisms (CLO) have been found to be present in many environmental niches, including human sewage and agricultural run-off, as well as in a number of aquatic species worldwide. Therefore, monitoring their presence in sentinel wildlife species may be useful in assessing the wider health of marine food webs in response to habitat loss, pollution and disease. We used nasal swabs from live (n = 42) and dead (n = 50) pre-weaned grey seal pups and samples of differing natal substrates (n = 8) from an off-shore island devoid of livestock and permanent human habitation to determine if CLO DNA is present in these mammals and to identify possible sources. RESULTS: We recovered CLO DNA from 32/92 (34.7%) nasal swabs from both live (n = 17) and dead (n = 15) seal pups that clustered most closely with currently recognised species belonging to three chlamydial families: Parachlamydiaceae (n = 22), Rhabdochlamydiaceae (n = 6), and Simkaniaceae (n = 3). All DNA positive sediment samples (n = 7) clustered with the Rhabdochlamydiaceae. No difference was found in rates of recovery of CLO DNA in live versus dead pups suggesting the organisms are commensal but their potential as opportunistic secondary pathogens could not be determined. CONCLUSION: This is the first report of CLO DNA being found in marine mammals. This identification warrants further investigation in other seal populations around the coast of the UK and in other areas of the world to determine if this finding is unique or more common than shown by this data. Further investigation would also be warranted to determine if they are present as purely commensal organisms or whether they could also be opportunistic pathogens in seals, as well as to investigate possible sources of origin, including whether they originated as a result of anthropogenic impacts, including human waste and agricultural run-off.

Network analyses reveal novel aspects of ALS pathogenesis.

Amyotrophic Lateral Sclerosis (ALS) is a fatal neurodegenerative disease characterized by selective loss of motor neurons, muscle atrophy and paralysis. Mutations in the human VAMP-associated protein B (hVAPB) cause a heterogeneous group of motor neuron diseases including ALS8. Despite extensive research, the molecular mechanisms underlying ALS pathogenesis remain largely unknown. Genetic screens for key interactors of hVAPB activity in the intact nervous system, however, represent a fundamental approach towards understanding the in vivo function of hVAPB and its role in ALS pathogenesis. Targeted expression of the disease-causing allele leads to neurodegeneration and progressive decline in motor performance when expressed in the adult Drosophila, eye or in its entire nervous system, respectively. By using these two phenotypic readouts, we carried out a systematic survey of the Drosophila genome to identify modifiers of hVAPB-induced neurotoxicity. Modifiers cluster in a diverse array of biological functions including processes and genes that have been previously linked to hVAPB function, such as proteolysis and vesicular trafficking. In addition to established mechanisms, the screen identified endocytic trafficking and genes controlling proliferation and apoptosis as potent modifiers of ALS8-mediated defects. Surprisingly, the list of modifiers was mostly enriched for proteins linked to lipid droplet biogenesis and dynamics. Computational analysis reveals that most modifiers can be linked into a complex network of interacting genes, and that the human genes homologous to the Drosophila modifiers can be assembled into an interacting network largely overlapping with that in flies. Identity markers of the endocytic process were also found to abnormally accumulate in ALS patients, further supporting the relevance of the fly data for human biology. Collectively, these results not only lead to a better understanding of hVAPB function but also point to potentially relevant targets for therapeutic intervention.

Analysis of compound heterozygotes reveals that the mouse floxed Pax6 (tm1Ued) allele produces abnormal eye phenotypes.

Analysis of abnormal phenotypes produced by different types of mutations has been crucial for our understanding of gene function. Some floxed alleles that retain a neomycin-resistance selection cassette (neo cassette) are not equivalent to wild-type alleles and provide useful experimental resources. Pax6 is an important developmental gene and the aim of this study was to determine whether the floxed Pax6 (tm1Ued) (Pax6 (fl) ) allele, which has a retained neo cassette, produced any abnormal eye phenotypes that would imply that it differs from the wild-type allele. Homozygous Pax6 (fl/fl) and heterozygous Pax6 (fl/+) mice had no overt qualitative eye abnormalities but morphometric analysis showed that Pax6 (fl/fl) corneas tended be thicker and smaller in diameter. To aid identification of weak effects, we produced compound heterozygotes with the Pax6 (Sey-Neu) (Pax6 (-)) null allele. Pax6 (fl/-) compound heterozygotes had more severe eye abnormalities than Pax6 (+/-) heterozygotes, implying that Pax6 (fl) differs from the wild-type Pax6 (+) allele. Immunohistochemistry showed that the Pax6 (fl/-) corneal epithelium was positive for keratin 19 and negative for keratin 12, indicating that it was abnormally differentiated. This Pax6 (fl) allele provides a useful addition to the existing Pax6 allelic series and this study demonstrates the utility of using compound heterozygotes with null alleles to unmask cryptic effects of floxed alleles.

The reduced kinome of Ostreococcus tauri: core eukaryotic signalling components in a tractable model species.

BACKGROUND: The current knowledge of eukaryote signalling originates from phenotypically diverse organisms. There is a pressing need to identify conserved signalling components among eukaryotes, which will lead to the transfer of knowledge across kingdoms. Two useful properties of a eukaryote model for signalling are (1) reduced signalling complexity, and (2) conservation of signalling components. The alga Ostreococcus tauri is described as the smallest free-living eukaryote. With less than 8,000 genes, it represents a highly constrained genomic palette. RESULTS: Our survey revealed 133 protein kinases and 34 protein phosphatases (1.7% and 0.4% of the proteome). We conducted phosphoproteomic experiments and constructed domain structures and phylogenies for the catalytic protein-kinases. For each of the major kinases families we review the completeness and divergence of O. tauri representatives in comparison to the well-studied kinomes of the laboratory models Arabidopsis thaliana and Saccharomyces cerevisiae, and of Homo sapiens. Many kinase clades in O. tauri were reduced to a single member, in preference to the loss of family diversity, whereas TKL and ABC1 clades were expanded. We also identified kinases that have been lost in A. thaliana but retained in O. tauri. For three, contrasting eukaryotic pathways - TOR, MAPK, and the circadian clock - we established the subset of conserved components and demonstrate conserved sites of substrate phosphorylation and kinase motifs. CONCLUSIONS: We conclude that O. tauri satisfies our two central requirements. Several of its kinases are more closely related to H. sapiens orthologs than S. cerevisiae is to H. sapiens. The greatly reduced kinome of O. tauri is therefore a suitable model for signalling in free-living eukaryotes.

The gene regulatory cascade linking proneural specification with differentiation in Drosophila sensory neurons.

In neurogenesis, neural cell fate specification is generally triggered by proneural transcription factors. Whilst the role of proneural factors in fate specification is well studied, the link between neural specification and the cellular pathways that ultimately must be activated to construct specialised neurons is usually obscure. High-resolution temporal profiling of gene expression reveals the events downstream of atonal proneural gene function during the development of Drosophila chordotonal (mechanosensory) neurons. Among other findings, this reveals the onset of expression of genes required for construction of the ciliary dendrite, a key specialisation of mechanosensory neurons. We determine that atonal activates this cellular differentiation pathway in several ways. Firstly, atonal directly regulates Rfx, a well-known highly conserved ciliogenesis transcriptional regulator. Unexpectedly, differences in Rfx regulation by proneural factors may underlie variations in ciliary dendrite specialisation in different sensory neuronal lineages. In contrast, fd3F encodes a novel forkhead family transcription factor that is exclusively expressed in differentiating chordotonal neurons. fd3F regulates genes required for specialized aspects of chordotonal dendrite physiology. In addition to these intermediate transcriptional regulators, we show that atonal directly regulates a novel gene, dilatory, that is directly associated with ciliogenesis during neuronal differentiation. Our analysis demonstrates how early cell fate specification factors can regulate structural and physiological differentiation of neuronal cell types. It also suggests a model for how subtype differentiation in different neuronal lineages may be regulated by different proneural factors. In addition, it provides a paradigm for how transcriptional regulation may modulate the ciliogenesis pathway to give rise to structurally and functionally specialised ciliary dendrites.

Consolidation and translation regulation.

mRNA translation, or protein synthesis, is a major component of the transformation of the genetic code into any cellular activity. This complicated, multistep process is divided into three phases: initiation, elongation, and termination. Initiation is the step at which the ribosome is recruited to the mRNA, and is regarded as the major rate-limiting step in translation, while elongation consists of the elongation of the polypeptide chain; both steps are frequent targets for regulation, which is defined as a change in the rate of translation of an mRNA per unit time. In the normal brain, control of translation is a key mechanism for regulation of memory and synaptic plasticity consolidation, i.e., the off-line processing of acquired information. These regulation processes may differ between different brain structures or neuronal populations. Moreover, dysregulation of translation leads to pathological brain function such as memory impairment. Both normal and abnormal function of the translation machinery is believed to lead to translational up-regulation or down-regulation of a subset of mRNAs. However, the identification of these newly synthesized proteins and determination of the rates of protein synthesis or degradation taking place in different neuronal types and compartments at different time points in the brain demand new proteomic methods and system biology approaches. Here, we discuss in detail the relationship between translation regulation and memory or synaptic plasticity consolidation while focusing on a model of cortical-dependent taste learning task and hippocampal-dependent plasticity. In addition, we describe a novel systems biology perspective to better describe consolidation.

Improved prediction of hiking speeds using a data driven approach.

Hikers and hillwalkers typically use the gradient in the direction of travel (walking slope) as the main variable in established methods for predicting walking time (via the walking speed) along a route. Research into fell-running has suggested further variables which may improve speed algorithms in this context; the gradient of the terrain (hill slope) and the level of terrain obstruction. Recent improvements in data availability, as well as widespread use of GPS tracking now make it possible to explore these variables in a walking speed model at a sufficient scale to test statistical significance. We tested various established models used to predict walking speed against public GPS data from almost 88,000 km of UK walking / hiking tracks. Tracks were filtered to remove breaks and non-walking sections. A new generalised linear model (GLM) was then used to predict walking speeds. Key differences between the GLM and established rules were that the GLM considered the gradient of the terrain (hill slope) irrespective of walking slope, as well as the terrain type and level of terrain obstruction in off-road travel. All of these factors were shown to be highly significant, and this is supported by a lower root-mean-square-error compared to existing functions. We also observed an increase in RMSE between the GLM and established methods as hill slope increases, further supporting the importance of this variable.

SFARI genes and where to find them; modelling Autism Spectrum Disorder specific gene expression dysregulation with RNA-seq data.

Autism Spectrum Disorders (ASD) have a strong, yet heterogeneous, genetic component. Among the various methods that are being developed to help reveal the underlying molecular aetiology of the disease one approach that is gaining popularity is the combination of gene expression and clinical genetic data, often using the SFARI-gene database, which comprises lists of curated genes considered to have causative roles in ASD when mutated in patients. We build a gene co-expression network to study the relationship between ASD-specific transcriptomic data and SFARI genes and then analyse it at different levels of granularity. No significant evidence is found of association between SFARI genes and differential gene expression patterns when comparing ASD samples to a control group, nor statistical enrichment of SFARI genes in gene co-expression network modules that have a strong correlation with ASD diagnosis. However, classification models that incorporate topological information from the whole ASD-specific gene co-expression network can predict novel SFARI candidate genes that share features of existing SFARI genes and have support for roles in ASD in the literature. A statistically significant association is also found between the absolute level of gene expression and SFARI's genes and Scores, which can confound the analysis if uncorrected. We propose a novel approach to correct for this that is general enough to be applied to other problems affected by continuous sources of bias. It was found that only co-expression network analyses that integrate information from the whole network are able to reveal signatures linked to ASD diagnosis and novel candidate genes for the study of ASD, which individual gene or module analyses fail to do. It was also found that the influence of SFARI genes permeates not only other ASD scoring systems, but also lists of genes believed to be involved in other neurodevelopmental disorders.

Comparison of rule- and ordinary differential equation-based dynamic model of DARPP-32 signalling network.

Dynamic modelling has considerably improved our understanding of complex molecular mechanisms. Ordinary differential equations (ODEs) are the most detailed and popular approach to modelling the dynamics of molecular systems. However, their application in signalling networks, characterised by multi-state molecular complexes, can be prohibitive. Contemporary modelling methods, such as rule- based (RB) modelling, have addressed these issues. The advantages of RB modelling over ODEs have been presented and discussed in numerous reviews. In this study, we conduct a direct comparison of the time courses of a molecular system founded on the same reaction network but encoded in the two frameworks. To make such a comparison, a set of reactions that underlie an ODE model was manually encoded in the Kappa language, one of the RB implementations. A comparison of the models was performed at the level of model specification and dynamics, acquired through model simulations. In line with previous reports, we confirm that the Kappa model recapitulates the general dynamics of its ODE counterpart with minor differences. These occur when molecules have multiple sites binding the same interactor. Furthermore, activation of these molecules in the RB model is slower than in the ODE one. As reported for other molecular systems, we find that, also for the DARPP-32 reaction network, the RB representation offers a more expressive and flexible syntax that facilitates access to fine details of the model, easing model reuse. In parallel with these analyses, we report a refactored model of the DARPP-32 interaction network that can serve as a canvas for the development of more complex dynamic models to study this important molecular system.

Bio::Homology::InterologWalk--a Perl module to build putative protein-protein interaction networks through interolog mapping.

BACKGROUND: Protein-protein interaction (PPI) data are widely used to generate network models that aim to describe the relationships between proteins in biological systems. The fidelity and completeness of such networks is primarily limited by the paucity of protein interaction information and by the restriction of most of these data to just a few widely studied experimental organisms. In order to extend the utility of existing PPIs, computational methods can be used that exploit functional conservation between orthologous proteins across taxa to predict putative PPIs or 'interologs'. To date most interolog prediction efforts have been restricted to specific biological domains with fixed underlying data sources and there are no software tools available that provide a generalised framework for 'on-the-fly' interolog prediction. RESULTS: We introduce Bio::Homology::InterologWalk, a Perl module to retrieve, prioritise and visualise putative protein-protein interactions through an orthology-walk method. The module uses orthology and experimental interaction data to generate putative PPIs and optionally collates meta-data into an Interaction Prioritisation Index that can be used to help prioritise interologs for further analysis. We show the application of our interolog prediction method to the genomic interactome of the fruit fly, Drosophila melanogaster. We analyse the resulting interaction networks and show that the method proposes new interactome members and interactions that are candidates for future experimental investigation. CONCLUSIONS: Our interolog prediction tool employs the Ensembl Perl API and PSICQUIC enabled protein interaction data sources to generate up to date interologs 'on-the-fly'. This represents a significant advance on previous methods for interolog prediction as it allows the use of the latest orthology and protein interaction data for all of the genomes in Ensembl. The module outputs simple text files, making it easy to customise the results by post-processing, allowing the putative PPI datasets to be easily integrated into existing analysis workflows. The Bio::Homology::InterologWalk module, sample scripts and full documentation are freely available from the Comprehensive Perl Archive Network (CPAN) under the GNU Public license.

A unified resource and configurable model of the synapse proteome and its role in disease.

Genes encoding synaptic proteins are highly associated with neuronal disorders many of which show clinical co-morbidity. We integrated 58 published synaptic proteomic datasets that describe over 8000 proteins and combined them with direct protein-protein interactions and functional metadata to build a network resource that reveals the shared and unique protein components that underpin multiple disorders. All the data are provided in a flexible and accessible format to encourage custom use.

Merged consensus clustering to assess and improve class discovery with microarray data.

BACKGROUND: One of the most commonly performed tasks when analysing high throughput gene expression data is to use clustering methods to classify the data into groups. There are a large number of methods available to perform clustering, but it is often unclear which method is best suited to the data and how to quantify the quality of the classifications produced. RESULTS: Here we describe an R package containing methods to analyse the consistency of clustering results from any number of different clustering methods using resampling statistics. These methods allow the identification of the the best supported clusters and additionally rank cluster members by their fidelity within the cluster. These metrics allow us to compare the performance of different clustering algorithms under different experimental conditions and to select those that produce the most reliable clustering structures. We show the application of this method to simulated data, canonical gene expression experiments and our own novel analysis of genes involved in the specification of the peripheral nervous system in the fruitfly, Drosophila melanogaster. CONCLUSIONS: Our package enables users to apply the merged consensus clustering methodology conveniently within the R programming environment, providing both analysis and graphical display functions for exploring clustering approaches. It extends the basic principle of consensus clustering by allowing the merging of results between different methods to provide an averaged clustering robustness. We show that this extension is useful in correcting for the tendency of clustering algorithms to treat outliers differently within datasets. The R package, clusterCons, is freely available at CRAN and sourceforge under the GNU public licence.

MicroRNA-34a Acutely Regulates Synaptic Efficacy in the Adult Dentate Gyrus In Vivo.

Activity-dependent synaptic plasticity involves rapid regulation of neuronal protein synthesis on a time-scale of minutes. miRNA function in synaptic plasticity and memory formation has been elucidated by stable experimental manipulation of miRNA expression and activity using transgenic approaches and viral vectors. However, the impact of rapid miRNA modulation on synaptic efficacy is unknown. Here, we examined the effect of acute (12 min), intrahippocampal infusion of a miR-34a antagonist (antimiR) on medial perforant path-evoked synaptic transmission in the dentate gyrus of adult anesthetised rats. AntimiR-34a infusion acutely depressed medial perforant path-evoked field excitatory post-synaptic potentials (fEPSPs). The fEPSP decrease was detected within 9 min of infusion, lasted for hours, and was associated with knockdown of antimiR-34a levels. AntimiR-34a-induced synaptic depression was sequence-specific; no changes were elicited by infusion of scrambled or mismatch control. The rapid modulation suggests that a target, or set of targets, is regulated by miR-34a. Western blot analysis of dentate gyrus lysates revealed enhanced expression of Arc, a known miR-34a target, and four novel predicted targets (Ctip2, PKI-1α, TCF4 and Ube2g1). Remarkably, antimiR-34a had no effect when infused during the maintenance phase of long-term potentiation. We conclude that miR-34a regulates basal synaptic efficacy in the adult dentate gyrus in vivo. To our knowledge, these in vivo findings are the first to demonstrate acute (< 9 min) regulation of synaptic efficacy in the adult brain by a miRNA.

Mixed-species RNA-seq for elucidation of non-cell-autonomous control of gene transcription.

Transcriptomic changes induced in one cell type by another mediate many biological processes in the brain and elsewhere; however, achieving artifact-free physical separation of cell types to study them is challenging and generally allows for analysis of only a single cell type. We describe an approach using a co-culture of distinct cell types from different species that enables physical cell sorting to be replaced by in silico RNA sequencing (RNA-seq) read sorting, which is possible because of evolutionary divergence of messenger RNA (mRNA) sequences. As an exemplary experiment, we describe the co-culture of purified neurons, astrocytes, and microglia from different species (12-14 d). We describe how to use our Python tool, Sargasso, to separate the reads from conventional RNA-seq according to species and to eliminate any artifacts borne of imperfect genome annotation (10 h). We show how this procedure, which requires no special skills beyond those that might normally be expected of wet lab and bioinformatics researchers, enables the simultaneous transcriptomic profiling of different cell types, revealing the distinct influence of microglia on astrocytic and neuronal transcriptomes under inflammatory conditions.

Author Correction: Neurons and neuronal activity control gene expression in astrocytes to regulate their development and metabolism.

This corrects the article DOI: 10.1038/ncomms15132.

Neurons and neuronal activity control gene expression in astrocytes to regulate their development and metabolism.

The influence that neurons exert on astrocytic function is poorly understood. To investigate this, we first developed a system combining cortical neurons and astrocytes from closely related species, followed by RNA-seq and in silico species separation. This approach uncovers a wide programme of neuron-induced astrocytic gene expression, involving Notch signalling, which drives and maintains astrocytic maturity and neurotransmitter uptake function, is conserved in human development, and is disrupted by neurodegeneration. Separately, hundreds of astrocytic genes are acutely regulated by synaptic activity via mechanisms involving cAMP/PKA-dependent CREB activation. This includes the coordinated activity-dependent upregulation of major astrocytic components of the astrocyte-neuron lactate shuttle, leading to a CREB-dependent increase in astrocytic glucose metabolism and elevated lactate export. Moreover, the groups of astrocytic genes induced by neurons or neuronal activity both show age-dependent decline in humans. Thus, neurons and neuronal activity regulate the astrocytic transcriptome with the potential to shape astrocyte-neuron metabolic cooperation.

Positive autoregulation of the transcription factor Pax6 in response to increased levels of either of its major isoforms, Pax6 or Pax6(5a), in cultured cells.

BACKGROUND: Pax6 is a transcription factor essential for normal development of the eyes and nervous system. It has two major isoforms, Pax6 and Pax6(5a), and the ratios between their expression levels vary within narrow limits. We tested the effects of overexpressing either one or other isoform on endogenous Pax6 expression levels in Neuro2A and NIH3T3 cells. RESULTS: We found that both isoforms caused an up-regulation of endogenous Pax6 expression in cells with (Neuro2A) or without (NIH3T3) constitutive Pax6 expression. Western blots showed that cells stably transfected with constructs expressing either Pax6 or Pax6(5a) contained raised levels of both Pax6 and Pax6(5a). Quantitative RT-PCR confirmed an increase in levels of Pax6(5a) mRNA in cells containing Pax6-expressing constructs and an increase in levels of Pax6 mRNA in cells containing Pax6(5a)-expressing constructs. The fact that the introduction of constructs expressing only one isoform increased the cellular levels of not only that isoform but also the other indicates that activation of the endogenous Pax6 locus occurred. The ratio between the levels of the two isoforms was maintained close to physiological values. The overexpression of either isoform in neuroblastoma (Neuro2A) cell lines also promoted morphological change and an increase in beta-III-tubulin expression, indicating an increase in neurogenesis. CONCLUSION: Our results demonstrate that Pax6 can up-regulate production of Pax6 protein from an entire intact endogenous Pax6 locus in its genomic environment. This adds to previous studies showing that Pax6 can up-regulate reporter expression driven by isolated Pax6 regulatory elements. Furthermore, our results suggest that an important function of positive feedback might be to stabilise the relative levels of Pax6 and Pax6(5a).

Penetrance of eye defects in mice heterozygous for mutation of Gli3 is enhanced by heterozygous mutation of Pax6.

BACKGROUND: Knowledge of the consequences of heterozygous mutations of developmentally important genes is important for understanding human genetic disorders. The Gli3 gene encodes a zinc finger transcription factor and homozygous loss-of-function mutations of Gli3 are lethal. Humans heterozygous for mutations in this gene suffer Greig cephalopolysyndactyly or Pallister-Hall syndromes, in which limb defects are prominent, and mice heterozygous for similar mutations have extra digits. Here we examined whether eye development, which is abnormal in mice lacking functional Gli3, is defective in Gli3+/- mice. RESULTS: We showed that Gli3 is expressed in the developing eye but that Gli3+/- mice have only very subtle eye defects. We then generated mice compound heterozygous for mutations in both Gli3 and Pax6, which encodes another developmentally important transcription factor known to be crucial for eye development. Pax6+/-; Gli3+/- eyes were compared to the eyes of wild-type, Pax6+/- or Gli3+/- siblings. They exhibited a range of abnormalities of the retina, iris, lens and cornea that was more extensive than in single Gli3+/- or Pax6+/- mutants or than would be predicted by addition of their phenotypes. CONCLUSION: These findings indicate that heterozygous mutations of Gli3 can impact on eye development. The importance of a normal Gli3 gene dosage becomes greater in the absence of a normal Pax6 gene dosage, suggesting that the two genes co-operate during eye morphogenesis.