Early transcriptional similarities between two distinct neural lineages during ascidian embryogenesis.

In chordates, the central nervous system arises from precursors that have distinct developmental and transcriptional trajectories. Anterior nervous systems are ontogenically associated with ectodermal lineages while posterior nervous systems are associated with mesoderm. Taking advantage of the well-documented cell lineage of ascidian embryos, we asked to what extent the transcriptional states of the different neural lineages become similar during the course of progressive lineage restriction. We performed single-cell RNA sequencing (scRNA-seq) analyses on hand-dissected neural precursor cells of the two distinct lineages, together with those of their sister cell lineages, with a high temporal resolution covering five successive cell cycles from the 16-cell to neural plate stages. A transcription factor binding site enrichment analysis of neural specific genes at the neural plate stage revealed limited evidence for shared transcriptional control between the two neural lineages, consistent with their different ontogenies. Nevertheless, PCA analysis and hierarchical clustering showed that, by neural plate stages, the two neural lineages cluster together. Consistent with this, we identified a set of genes enriched in both neural lineages at the neural plate stage, including miR-124, Celf3.a, Zic.r-b, and Ets1/2. Altogether, the current study has revealed genome-wide transcriptional dynamics of neural progenitor cells of two distinct developmental origins. Our scRNA-seq dataset is unique and provides a valuable resource for future analyses, enabling a precise temporal resolution of cell types not previously described from dissociated embryos.

Molecular characterisation of a cellular conveyor belt in Clytia medusae.

The tentacular system of Clytia hemisphaerica medusa (Cnidaria, Hydrozoa) has recently emerged as a promising experimental model to tackle the developmental mechanisms that regulate cell lineage progression in an early-diverging animal phylum. From a population of proximal stem cells, the successive steps of tentacle stinging cell (nematocyte) elaboration, are spatially ordered along a "cellular conveyor belt". Furthermore, the C. hemisphaerica tentacular system exhibits bilateral organisation, with two perpendicular polarity axes (proximo-distal and oral-aboral). We aimed to improve our knowledge of this cellular system by combining RNAseq-based differential gene expression analyses and expression studies of Wnt signalling genes. RNAseq comparisons of gene expression levels were performed (i) between the tentacular system and a control medusa deprived of all tentacles, nematogenic sites and gonads, and (ii) between three samples staggered along the cellular conveyor belt. The behaviour in these differential expression analyses of two reference gene sets (stem cell genes; nematocyte genes), as well as the relative representations of selected gene ontology categories, support the validity of the cellular conveyor belt model. Expression patterns obtained by in situ hybridisation for selected highly differentially expressed genes and for Wnt signalling genes are largely consistent with the results from RNAseq. Wnt signalling genes exhibit complex spatial deployment along both polarity axes of the tentacular system, with the Wnt/ß-catenin pathway probably acting along the oral-aboral axis rather than the proximo-distal axis. These findings reinforce the idea that, despite overall radial symmetry, cnidarians have a full potential for elaboration of bilateral structures based on finely orchestrated deployment of an ancient developmental gene toolkit.

Acoelomorph flatworms are deuterostomes related to Xenoturbella.

Xenoturbellida and Acoelomorpha are marine worms with contentious ancestry. Both were originally associated with the flatworms (Platyhelminthes), but molecular data have revised their phylogenetic positions, generally linking Xenoturbellida to the deuterostomes and positioning the Acoelomorpha as the most basally branching bilaterian group(s). Recent phylogenomic data suggested that Xenoturbellida and Acoelomorpha are sister taxa and together constitute an early branch of Bilateria. Here we assemble three independent data sets-mitochondrial genes, a phylogenomic data set of 38,330 amino-acid positions and new microRNA (miRNA) complements-and show that the position of Acoelomorpha is strongly affected by a long-branch attraction (LBA) artefact. When we minimize LBA we find consistent support for a position of both acoelomorphs and Xenoturbella within the deuterostomes. The most likely phylogeny links Xenoturbella and Acoelomorpha in a clade we call Xenacoelomorpha. The Xenacoelomorpha is the sister group of the Ambulacraria (hemichordates and echinoderms). We show that analyses of miRNA complements have been affected by character loss in the acoels and that both groups possess one miRNA and the gene Rsb66 otherwise specific to deuterostomes. In addition, Xenoturbella shares one miRNA with the ambulacrarians, and two with the acoels. This phylogeny makes sense of the shared characteristics of Xenoturbellida and Acoelomorpha, such as ciliary ultrastructure and diffuse nervous system, and implies the loss of various deuterostome characters in the Xenacoelomorpha including coelomic cavities, through gut and gill slits.

Clinical whole-genome sequencing in severe early-onset epilepsy reveals new genes and improves molecular diagnosis.

In severe early-onset epilepsy, precise clinical and molecular genetic diagnosis is complex, as many metabolic and electro-physiological processes have been implicated in disease causation. The clinical phenotypes share many features such as complex seizure types and developmental delay. Molecular diagnosis has historically been confined to sequential testing of candidate genes known to be associated with specific sub-phenotypes, but the diagnostic yield of this approach can be low. We conducted whole-genome sequencing (WGS) on six patients with severe early-onset epilepsy who had previously been refractory to molecular diagnosis, and their parents. Four of these patients had a clinical diagnosis of Ohtahara Syndrome (OS) and two patients had severe non-syndromic early-onset epilepsy (NSEOE). In two OS cases, we found de novo non-synonymous mutations in the genes KCNQ2 and SCN2A. In a third OS case, WGS revealed paternal isodisomy for chromosome 9, leading to identification of the causal homozygous missense variant in KCNT1, which produced a substantial increase in potassium channel current. The fourth OS patient had a recessive mutation in PIGQ that led to exon skipping and defective glycophosphatidyl inositol biosynthesis. The two patients with NSEOE had likely pathogenic de novo mutations in CBL and CSNK1G1, respectively. Mutations in these genes were not found among 500 additional individuals with epilepsy. This work reveals two novel genes for OS, KCNT1 and PIGQ. It also uncovers unexpected genetic mechanisms and emphasizes the power of WGS as a clinical tool for making molecular diagnoses, particularly for highly heterogeneous disorders.

Gene panel sequencing improves the diagnostic work-up of patients with idiopathic erythrocytosis and identifies new mutations.

Erythrocytosis is a rare disorder characterized by increased red cell mass and elevated hemoglobin concentration and hematocrit. Several genetic variants have been identified as causes for erythrocytosis in genes belonging to different pathways including oxygen sensing, erythropoiesis and oxygen transport. However, despite clinical investigation and screening for these mutations, the cause of disease cannot be found in a considerable number of patients, who are classified as having idiopathic erythrocytosis. In this study, we developed a targeted next-generation sequencing panel encompassing the exonic regions of 21 genes from relevant pathways (~79 Kb) and sequenced 125 patients with idiopathic erythrocytosis. The panel effectively screened 97% of coding regions of these genes, with an average coverage of 450×. It identified 51 different rare variants, all leading to alterations of protein sequence, with 57 out of 125 cases (45.6%) having at least one of these variants. Ten of these were known erythrocytosis-causing variants, which had been missed following existing diagnostic algorithms. Twenty-two were novel variants in erythrocytosis-associated genes (EGLN1, EPAS1, VHL, BPGM, JAK2, SH2B3) and in novel genes included in the panel (e.g. EPO, EGLN2, HIF3A, OS9), some with a high likelihood of functionality, for which future segregation, functional and replication studies will be useful to provide further evidence for causality. The rest were classified as polymorphisms. Overall, these results demonstrate the benefits of using a gene panel rather than existing methods in which focused genetic screening is performed depending on biochemical measurements: the gene panel improves diagnostic accuracy and provides the opportunity for discovery of novel variants.

Causes and consequences of chromatin variation between inbred mice.

Variation at regulatory elements, identified through hypersensitivity to digestion by DNase I, is believed to contribute to variation in complex traits, but the extent and consequences of this variation are poorly characterized. Analysis of terminally differentiated erythroblasts in eight inbred strains of mice identified reproducible variation at approximately 6% of DNase I hypersensitive sites (DHS). Only 30% of such variable DHS contain a sequence variant predictive of site variation. Nevertheless, sequence variants within variable DHS are more likely to be associated with complex traits than those in non-variant DHS, and variants associated with complex traits preferentially occur in variable DHS. Changes at a small proportion (less than 10%) of variable DHS are associated with changes in nearby transcriptional activity. Our results show that whilst DNA sequence variation is not the major determinant of variation in open chromatin, where such variants exist they are likely to be causal for complex traits.

Improving animal phylogenies with genomic data.

Since the first animal genomes were completely sequenced ten years ago, evolutionary biologists have attempted to use the encoded information to reconstruct different aspects of the earliest stages of animal evolution. One of the most important uses of genome sequences is to understand relationships between animal phyla. Despite the wealth of data available, ranging from primary sequence data to gene and genome structures, our lack of understanding of the modes of evolution of genomic characters means that using these data is fraught with potential difficulties, leading to errors in phylogeny reconstruction. Improved understanding of how different character types evolve, the use of this knowledge to develop more accurate models of evolution, and denser taxonomic sampling, are now minimizing the sources of error. The wealth of genomic data now being produced promises that a well-resolved tree of the animal phyla will be available in the near future.

Phylogenomic analysis of echinoderm class relationships supports Asterozoa.

While some aspects of the phylogeny of the five living echinoderm classes are clear, the position of the ophiuroids (brittlestars) relative to asteroids (starfish), echinoids (sea urchins) and holothurians (sea cucumbers) is controversial. Ophiuroids have a pluteus-type larva in common with echinoids giving some support to an ophiuroid/echinoid/holothurian clade named Cryptosyringida. Most molecular phylogenetic studies, however, support an ophiuroid/asteroid clade (Asterozoa) implying either convergent evolution of the pluteus or reversals to an auricularia-type larva in asteroids and holothurians. A recent study of 10 genes from four of the five echinoderm classes used 'phylogenetic signal dissection' to separate alignment positions into subsets of (i) suboptimal, heterogeneously evolving sites (invariant plus rapidly changing) and (ii) the remaining optimal, homogeneously evolving sites. Along with most previous molecular phylogenetic studies, their set of heterogeneous sites, expected to be more prone to systematic error, support Asterozoa. The homogeneous sites, in contrast, support an ophiuroid/echinoid grouping, consistent with the cryptosyringid clade, leading them to posit homology of the ophiopluteus and echinopluteus. Our new dataset comprises 219 genes from all echinoderm classes; analyses using probabilistic Bayesian phylogenetic methods strongly support Asterozoa. The most reliable, slowly evolving quartile of genes also gives highest support for Asterozoa; this support diminishes in second and third quartiles and the fastest changing quartile places the ophiuroids close to the root. Using phylogenetic signal dissection, we find heterogenous sites support an unlikely grouping of Ophiuroidea + Holothuria while homogeneous sites again strongly support Asterozoa. Our large and taxonomically complete dataset finds no support for the cryptosyringid hypothesis; in showing strong support for the Asterozoa, our preferred topology leaves the question of homology of pluteus larvae open.

Scaffolding low quality genomes using orthologous protein sequences.

MOTIVATION: The ready availability of next-generation sequencing has led to a situation where it is easy to produce very fragmentary genome assemblies. We present a pipeline, SWiPS (Scaffolding With Protein Sequences), that uses orthologous proteins to improve low quality genome assemblies. The protein sequences are used as guides to scaffold existing contigs, while simultaneously allowing the gene structure to be predicted by homology. RESULTS: To perform, SWiPS does not depend on a high N50 or whole proteins being encoded on a single contig. We tested our algorithm on simulated next-generation data from Ciona intestinalis, real next-generation data from Drosophila melanogaster, a complex genome assembly of Homo sapiens and the low coverage Sanger sequence assembly of Callorhinchus milii. The improvements in N50 are of the order of ∼20% for the C.intestinalis and H.sapiens assemblies, which is significant, considering the large size of intergenic regions in these eukaryotes. Using the CEGMA pipeline to assess the gene space represented in the genome assemblies, the number of genes retrieved increased by >110% for C.milii and from 20 to 40% for C.intestinalis. The scaffold error rates are low: 85-90% of scaffolds are fully correct, and >95% of local contig joins are correct. AVAILABILITY: SWiPS is available freely for download at http://www.well.ox.ac.uk/∼yli142/swips.html. CONTACT: yang.li@well.ox.ac.uk or copley@well.ox.ac.uk

Genetic dissection of a behavioral quantitative trait locus shows that Rgs2 modulates anxiety in mice.

Here we present a strategy to determine the genetic basis of variance in complex phenotypes that arise from natural, as opposed to induced, genetic variation in mice. We show that a commercially available strain of outbred mice, MF1, can be treated as an ultrafine mosaic of standard inbred strains and accordingly used to dissect a known quantitative trait locus influencing anxiety. We also show that this locus can be subdivided into three regions, one of which contains Rgs2, which encodes a regulator of G protein signaling. We then use quantitative complementation to show that Rgs2 is a quantitative trait gene. This combined genetic and functional approach should be applicable to the analysis of any quantitative trait.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes.

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Next-generation sequencing in health-care delivery: lessons from the functional analysis of rhodopsin.

PURPOSE: The interpretation of genetic information has always been challenging, but next-generation sequencing produces data on such a vast scale that many more variants of uncertain pathogenicity will be found. We exemplify this issue with reference to human rhodopsin, in which pathogenic mutations can lead to autosomal dominant retinitis pigmentosa. METHODS: Rhodopsin variants, with unknown pathogenicity, were found in patients by next-generation and Sanger sequencing and a multidisciplinary approach was used to determine their functional significance. RESULTS: Four variants in rhodopsin were identified: F45L, P53R, R69H, and M39R, with the latter two substitutions being novel. We investigated the cellular transport and photopigment function of all four human substitutions and found that the F45L and R69H variants behave like wild-type and are highly unlikely to be pathogenic. By contrast, P53R (a de novo change) and M39R were retained in the endoplasmic reticulum with significantly reduced functionality and are clearly pathogenic. CONCLUSION: Potential pathogenicity of variants requires careful assessment using clinical, genetic, and functional data. We suggest that a multidisciplinary pathway of assessment, using several functional assays, will be required if next-generation sequencing is to be used effectively, reliably, and safely in the clinical environment.

Beyond the enhanceosome: cluster of novel κB sites downstream of the human IFN-ß gene is essential for lipopolysaccharide-induced gene activation.

The expression of interferon-ß (IFN-ß) in virus-infected HeLa cells established a paradigm of multifactorial gene regulation, in which cooperative assembly of transcription factors (TFs) at the composite DNA element (enhanceosome), is central for amplification of weak activating signals provided by individual TFs. However, whether the same TFs and the same DNA element are essential for IFN-ß induction in response to bacterial stimuli are less well understood. Here we report that rapid and transient transcription of IFN-ß in response to TLR4 stimulation with bacterial lipopolysaccharide (LPS) follows nuclear factor-κB (NF-κB) RelA activation and recruitment to the IFN-ß genomic locus at multiple spatially separated regulatory regions. We demonstrate that the IFN-ß enhanceosome region is not sufficient for maximal gene induction in response to LPS and identify an essential cluster of homotypic κB sites in the 3' downstream of the gene. The cluster is characterized by elevated levels of histone 3 lysine 4 mono-methylation, a chromatin signature of enhancers, and efficiently binds RelA-containing NF-κB complexes in vitro and in vivo. These findings demonstrate that IFN-ß gene activation via multifactorial enhanceosome assembly is potentiated in LPS-stimulated cells by NF-κB interactions with all functional κB sites in the locus.

Deuterostome phylogeny reveals monophyletic chordates and the new phylum Xenoturbellida.

Deuterostomes comprise vertebrates, the related invertebrate chordates (tunicates and cephalochordates) and three other invertebrate taxa: hemichordates, echinoderms and Xenoturbella. The relationships between invertebrate and vertebrate deuterostomes are clearly important for understanding our own distant origins. Recent phylogenetic studies of chordate classes and a sea urchin have indicated that urochordates might be the closest invertebrate sister group of vertebrates, rather than cephalochordates, as traditionally believed. More remarkable is the suggestion that cephalochordates are closer to echinoderms than to vertebrates and urochordates, meaning that chordates are paraphyletic. To study the relationships among all deuterostome groups, we have assembled an alignment of more than 35,000 homologous amino acids, including new data from a hemichordate, starfish and Xenoturbella. We have also sequenced the mitochondrial genome of Xenoturbella. We support the clades Olfactores (urochordates and vertebrates) and Ambulacraria (hemichordates and echinoderms). Analyses using our new data, however, do not support a cephalochordate and echinoderm grouping and we conclude that chordates are monophyletic. Finally, nuclear and mitochondrial data place Xenoturbella as the sister group of the two ambulacrarian phyla. As such, Xenoturbella is shown to be an independent phylum, Xenoturbellida, bringing the number of living deuterostome phyla to four.

Past, present and future of Clytia hemisphaerica as a laboratory jellyfish.

The hydrozoan species Clytia hemisphaerica was selected in the mid-2000s to address the cellular and molecular basis of body axis specification in a cnidarian, providing a reliable daily source of gametes and building on a rich foundation of experimental embryology. The many practical advantages of this species include genetic uniformity of laboratory jellyfish, derived clonally from easily-propagated polyp colonies. Phylogenetic distance from other laboratory models adds value in providing an evolutionary perspective on many biological questions. Here we outline the current state of the art regarding available experimental approaches and in silico resources, and illustrate the contributions of Clytia to understanding embryo patterning mechanisms, oogenesis and regeneration. Looking forward, the recent establishment of transgenesis methods is now allowing gene function and imaging studies at adult stages, making Clytia particularly attractive for whole organism biology studies across fields and extending its scientific impact far beyond the original question of interest.

Genomic patterns of divergence in the early and late steps of speciation of the deep-sea vent thermophilic worms of the genus Alvinella.

BACKGROUND: The transient and fragmented nature of the deep-sea hydrothermal environment made of ridge subduction, plate collision and the emergence of new rifts is currently acting to separate of vent populations, promoting local adaptation and contributing to bursts of speciation and species specialization. The tube-dwelling worms Alvinella pompejana called the Pompeii worm and its sister species A. caudata live syntopically on the hottest part of deep-sea hydrothermal chimneys along the East Pacific Rise. They are exposed to extreme thermal and chemical gradients, which vary greatly in space and time, and thus represent ideal candidates for understanding the evolutionary mechanisms at play in the vent fauna evolution. RESULTS: We explored genomic patterns of divergence in the early and late stages of speciation of these emblematic worms using transcriptome assemblies and the first draft genome to better understand the relative role of geographic isolation and habitat preference in their genome evolution. Analyses were conducted on allopatric populations of Alvinella pompejana (early stage of separation) and between A. pompejana and its syntopic species Alvinella caudata (late stage of speciation). We first identified divergent genomic regions and targets of selection as well as their position in the genome over collections of orthologous genes and, then, described the speciation dynamics by documenting the annotation of the most divergent and/or positively selected genes involved in the isolation process. Gene mapping clearly indicated that divergent genes associated with the early stage of speciation, although accounting for nearly 30% of genes, are highly scattered in the genome without any island of divergence and not involved in gamete recognition or mito-nuclear incompatibilities. By contrast, genomes of A. pompejana and A. caudata are clearly separated with nearly all genes (96%) exhibiting high divergence. This congealing effect however seems to be linked to habitat specialization and still allows positive selection on genes involved in gamete recognition, as a possible long-duration process of species reinforcement. CONCLUSION: Our analyses highlight the non-negligible role of natural selection on both the early and late stages of speciation in the iconic thermophilic worms living on the walls of deep-sea hydrothermal chimneys. They shed light on the evolution of gene divergence during the process of speciation and species specialization over a very long period of time.

A RING-type ubiquitin ligase family member required to repress follicular helper T cells and autoimmunity.

Despite the sequencing of the human and mouse genomes, few genetic mechanisms for protecting against autoimmune disease are currently known. Here we systematically screen the mouse genome for autoimmune regulators to isolate a mouse strain, sanroque, with severe autoimmune disease resulting from a single recessive defect in a previously unknown mechanism for repressing antibody responses to self. The sanroque mutation acts within mature T cells to cause formation of excessive numbers of follicular helper T cells and germinal centres. The mutation disrupts a repressor of ICOS, an essential co-stimulatory receptor for follicular T cells, and results in excessive production of the cytokine interleukin-21. sanroque mice fail to repress diabetes-causing T cells, and develop high titres of autoantibodies and a pattern of pathology consistent with lupus. The causative mutation is in a gene of previously unknown function, roquin (Rc3h1), which encodes a highly conserved member of the RING-type ubiquitin ligase protein family. The Roquin protein is distinguished by the presence of a CCCH zinc-finger found in RNA-binding proteins, and localization to cytosolic RNA granules implicated in regulating messenger RNA translation and stability.

Lack of support for Deuterostomia prompts reinterpretation of the first Bilateria.

The bilaterally symmetric animals (Bilateria) are considered to comprise two monophyletic groups, Protostomia (Ecdysozoa and the Lophotrochozoa) and Deuterostomia (Chordata and the Xenambulacraria). Recent molecular phylogenetic studies have not consistently supported deuterostome monophyly. Here, we compare support for Protostomia and Deuterostomia using multiple, independent phylogenomic datasets. As expected, Protostomia is always strongly supported, especially by longer and higher-quality genes. Support for Deuterostomia, however, is always equivocal and barely higher than support for paraphyletic alternatives. Conditions that cause tree reconstruction errors-inadequate models, short internal branches, faster evolving genes, and unequal branch lengths-coincide with support for monophyletic deuterostomes. Simulation experiments show that support for Deuterostomia could be explained by systematic error. The branch between bilaterian and deuterostome common ancestors is, at best, very short, supporting the idea that the bilaterian ancestor may have been deuterostome-like. Our findings have important implications for the understanding of early animal evolution.

Whole-animal multiplexed single-cell RNA-seq reveals transcriptional shifts across Clytia medusa cell types.

We present an organism-wide, transcriptomic cell atlas of the hydrozoan medusa Clytia hemisphaerica and describe how its component cell types respond to perturbation. Using multiplexed single-cell RNA sequencing, in which individual animals were indexed and pooled from control and perturbation conditions into a single sequencing run, we avoid artifacts from batch effects and are able to discern shifts in cell state in response to organismal perturbations. This work serves as a foundation for future studies of development, function, and regeneration in a genetically tractable jellyfish species. Moreover, we introduce a powerful workflow for high-resolution, whole-animal, multiplexed single-cell genomics that is readily adaptable to other traditional or nontraditional model organisms.

Genome-wide association of hypoxia-inducible factor (HIF)-1alpha and HIF-2alpha DNA binding with expression profiling of hypoxia-inducible transcripts.

Hypoxia-inducible factor (HIF) controls an extensive range of adaptive responses to hypoxia. To better understand this transcriptional cascade we performed genome-wide chromatin immunoprecipitation using antibodies to two major HIF-alpha subunits, and correlated the results with genome-wide transcript profiling. Within a tiled promoter array we identified 546 and 143 sequences that bound, respectively, to HIF-1alpha or HIF-2alpha at high stringency. Analysis of these sequences confirmed an identical core binding motif for HIF-1alpha and HIF-2alpha (RCGTG) but demonstrated that binding to this motif was highly selective, with binding enriched at distinct regions both upstream and downstream of the transcriptional start. Comparison of HIF-promoter binding data with bidirectional HIF-dependent changes in transcript expression indicated that whereas a substantial proportion of positive responses (>20% across all significantly regulated genes) are direct, HIF-dependent gene suppression is almost entirely indirect. Comparison of HIF-1alpha- versus HIF-2alpha-binding sites revealed that whereas some loci bound HIF-1alpha in isolation, many bound both isoforms with similar affinity. Despite high-affinity binding to multiple promoters, HIF-2alpha contributed to few, if any, of the transcriptional responses to acute hypoxia at these loci. Given emerging evidence for biologically distinct functions of HIF-1alpha versus HIF-2alpha understanding the mechanisms restricting HIF-2alpha activity will be of interest.