Hubrecht Institute for Developmental Biology and Stem Cell Research

A high-quality sequence assembly of the zebrafish genome reveals the largest gene set of any vertebrate and provides information on key genomic features, and comparison to the human reference genome shows that approximately 70% of human protein-coding genes have at least one clear zebrafish orthologue. The genome of the zebrafish — a key model organism for the study of development and human disease — has now been sequenced and published as a well-annotated reference genome. Zebrafish turns out to have the largest gene set of any vertebrate so far sequenced, and few pseudogenes. Importantly for disease studies, comparison between human and zebrafish sequences reveals that 70% of human genes have at least one obvious zebrafish orthologue. A second paper reports on an ongoing effort to identify and phenotype disruptive mutations in every zebrafish protein-coding gene. Using the reference genome sequence along with high-throughput sequencing and efficient chemical mutagenesis, the project's initial results — covering 38% of all known protein-coding genes — they describe phenotypic consequences of more than 1,000 alleles. The long-term goal is the creation of a knockout allele in every protein-coding gene in the zebrafish genome. All mutant alleles and data are freely available at go.nature.com/en6mos . Zebrafish have become a popular organism for the study of vertebrate gene function1,2. The virtually transparent embryos of this species, and the ability to accelerate genetic studies by gene knockdown or overexpression, have led to the widespread use of zebrafish in the detailed investigation of vertebrate gene function and increasingly, the study of human genetic disease3,4,5. However, for effective modelling of human genetic disease it is important to understand the extent to which zebrafish genes and gene structures are related to orthologous human genes. To examine this, we generated a high-quality sequence assembly of the zebrafish genome, made up of an overlapping set of completely sequenced large-insert clones that were ordered and oriented using a high-resolution high-density meiotic map. Detailed automatic and manual annotation provides evidence of more than 26,000 protein-coding genes6, the largest gene set of any vertebrate so far sequenced. Comparison to the human reference genome shows that approximately 70% of human genes have at least one obvious zebrafish orthologue. In addition, the high quality of this genome assembly provides a clearer understanding of key genomic features such as a unique repeat content, a scarcity of pseudogenes, an enrichment of zebrafish-specific genes on chromosome 4 and chromosomal regions that influence sex determination.

The rat, mouse and human estrogen receptor (ER) exists as two subtypes, ER alpha and ER beta, which differ in the C-terminal ligand-binding domain and in the N-terminal transactivation domain. In this study, we investigated the estrogenic activity of environmental chemicals and phytoestrogens in competition binding assays with ER alpha or ER beta protein, and in a transient gene expression assay using cells in which an acute estrogenic response is created by cotransfecting cultures with recombinant human ER alpha or ER beta complementary DNA (cDNA) in the presence of an estrogen-dependent reporter plasmid. Saturation ligand-binding analysis of human ER alpha and ER beta protein revealed a single binding component for [3H]-17beta-estradiol (E2) with high affinity [dissociation constant (Kd) = 0.05 - 0.1 nM]. All environmental estrogenic chemicals [polychlorinated hydroxybiphenyls, dichlorodiphenyltrichloroethane (DDT) and derivatives, alkylphenols, bisphenol A, methoxychlor and chlordecone] compete with E2 for binding to both ER subtypes with a similar preference and degree. In most instances the relative binding affinities (RBA) are at least 1000-fold lower than that of E2. Some phytoestrogens such as coumestrol, genistein, apigenin, naringenin, and kaempferol compete stronger with E2 for binding to ER beta than to ER alpha. Estrogenic chemicals, as for instance nonylphenol, bisphenol A, o, p'-DDT and 2',4',6'-trichloro-4-biphenylol stimulate the transcriptional activity of ER alpha and ER beta at concentrations of 100-1000 nM. Phytoestrogens, including genistein, coumestrol and zearalenone stimulate the transcriptional activity of both ER subtypes at concentrations of 1-10 nM. The ranking of the estrogenic potency of phytoestrogens for both ER subtypes in the transactivation assay is different; that is, E2 >> zearalenone = coumestrol > genistein > daidzein > apigenin = phloretin > biochanin A = kaempferol = naringenin > formononetin = ipriflavone = quercetin = chrysin for ER alpha and E2 >> genistein = coumestrol > zearalenone > daidzein > biochanin A = apigenin = kaempferol = naringenin > phloretin = quercetin = ipriflavone = formononetin = chrysin for ER beta. Antiestrogenic activity of the phytoestrogens could not be detected, except for zearalenone which is a full agonist for ER alpha and a mixed agonist-antagonist for ER beta. In summary, while the estrogenic potency of industrial-derived estrogenic chemicals is very limited, the estrogenic potency of phytoestrogens is significant, especially for ER beta, and they may trigger many of the biological responses that are evoked by the physiological estrogens.

The prediction of the quaternary structure of biomolecular macromolecules is of paramount importance for fundamental understanding of cellular processes and drug design. In the era of integrative structural biology, one way of increasing the accuracy of modeling methods used to predict the structure of biomolecular complexes is to include as much experimental or predictive information as possible in the process. This has been at the core of our information-driven docking approach HADDOCK. We present here the updated version 2.2 of the HADDOCK portal, which offers new features such as support for mixed molecule types, additional experimental restraints and improved protocols, all of this in a user-friendly interface. With well over 6000 registered users and 108,000 jobs served, an increasing fraction of which on grid resources, we hope that this timely upgrade will help the community to solve important biological questions and further advance the field. The HADDOCK2.2 Web server is freely accessible to non-profit users at http://haddock.science.uu.nl/services/HADDOCK2.2.

The recent advent of methods for high-throughput single-cell molecular profiling has catalyzed a growing sense in the scientific community that the time is ripe to complete the 150-year-old effort to identify all cell types in the human body. The Human Cell Atlas Project is an international collaborative effort that aims to define all human cell types in terms of distinctive molecular profiles (such as gene expression profiles) and to connect this information with classical cellular descriptions (such as location and morphology). An open comprehensive reference map of the molecular state of cells in healthy human tissues would propel the systematic study of physiological states, developmental trajectories, regulatory circuitry and interactions of cells, and also provide a framework for understanding cellular dysregulation in human disease. Here we describe the idea, its potential utility, early proofs-of-concept, and some design considerations for the Human Cell Atlas, including a commitment to open data, code, and community.

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly ‘housekeeping’, whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research. A study from the FANTOM consortium using single-molecule cDNA sequencing of transcription start sites and their usage in human and mouse primary cells, cell lines and tissues reveals insights into the specificity and diversity of transcription patterns across different mammalian cell types. FANTOM5 (standing for functional annotation of the mammalian genome 5) is the fifth major stage of a major international collaboration that aims to dissect the transcriptional regulatory networks that define every human cell type. Two Articles in this issue of Nature present some of the project's latest results. The first paper uses the FANTOM5 panel of tissue and primary cell samples to define an atlas of active, in vivo bidirectionally transcribed enhancers across the human body. These authors show that bidirectional capped RNAs are a signature feature of active enhancers and identify more than 40,000 enhancer candidates from over 800 human cell and tissue samples. The enhancer atlas is used to compare regulatory programs between different cell types and identify disease-associated regulatory SNPs, and will be a resource for studies on cell-type-specific enhancers. In the second paper, single-molecule sequencing is used to map human and mouse transcription start sites and their usage in a panel of distinct human and mouse primary cells, cell lines and tissues to produce the most comprehensive mammalian gene expression atlas to date. The data provide a plethora of insights into open reading frames and promoters across different cell types in addition to valuable annotation of mammalian cell-type-specific transcriptomes.

MicroRNAs (miRNAs) are small noncoding RNAs, about 21 nucleotides in length, that can regulate gene expression by base-pairing to partially complementary mRNAs. Regulation by miRNAs can play essential roles in embryonic development. We determined the temporal and spatial expression patterns of 115 conserved vertebrate miRNAs in zebrafish embryos by microarrays and by in situ hybridizations, using locked-nucleic acid-modified oligonucleotide probes. Most miRNAs were expressed in a highly tissue-specific manner during segmentation and later stages, but not early in development, which suggests that their role is not in tissue fate establishment but in differentiation or maintenance of tissue identity.

The recent boom in microfluidics and combinatorial indexing strategies, combined with low sequencing costs, has empowered single-cell sequencing technology. Thousands-or even millions-of cells analyzed in a single experiment amount to a data revolution in single-cell biology and pose unique data science problems. Here, we outline eleven challenges that will be central to bringing this emerging field of single-cell data science forward. For each challenge, we highlight motivating research questions, review prior work, and formulate open problems. This compendium is for established researchers, newcomers, and students alike, highlighting interesting and rewarding problems for the coming years.

The PTH/PTHrP receptor binds to two ligands with distinct functions: the calcium-regulating hormone, parathyroid hormone (PTH), and the paracrine factor, PTH-related protein (PTHrP). Each ligand, in turn, is likely to activate more than one receptor. The functions of the PTH/PTHrP receptor were investigated by deletion of the murine gene by homologous recombination. Most PTH/PTHrP receptor (-/-) mutant mice died in mid-gestation, a phenotype not observed in PTHrP (-/-) mice, perhaps because of the effects of maternal PTHrP. Mice that survived exhibited accelerated differentiation of chondrocytes in bone, and their bones, grown in explant culture, were resistant to the effects of PTHrP and Sonic hedgehog. These results suggest that the PTH/PTHrP receptor mediates the effects of Indian Hedgehog and PTHrP on chondrocyte differentiation.

In many organisms the allocation of primordial germ cells (PGCs) is determined by the inheritance of maternal factors deposited in the egg. However, in mammals, inductive cell interactions are required around gastrulation to establish the germ line. Here, we show that Bmp4 homozygous null embryos contain no PGCs. They also lack an allantois, an extraembryonic mesodermal tissue derived, like the PGCs, from precursors in the proximal epiblast. Heterozygotes have fewer PGCs than normal, due to a reduction in the size of the founding population and not to an effect on its subsequent expansion. Analysis of beta-galactosidase activity in Bmp4(lacZneo) embryos reveals that prior to gastrulation, Bmp4 is expressed in the extraembryonic ectoderm. Later, Bmp4 is expressed in the extraembryonic mesoderm, but not in PGCs. Chimera analysis indicates that it is the Bmp4 expression in the extraembryonic ectoderm that regulates the formation of allantois and primordial germ cell precursors, and the size of the founding population of PGCs. The initiation of the germ line in the mouse therefore depends on a secreted signal from the previously segregated, extraembryonic, trophectoderm lineage.

Adult stem cells are crucial for physiological tissue renewal and regeneration after injury. Prevailing models assume the existence of a single quiescent population of stem cells residing in a specialized niche of a given tissue. Emerging evidence indicates that both quiescent (out of cell cycle and in a lower metabolic state) and active (in cell cycle and not able to retain DNA labels) stem cell subpopulations may coexist in several tissues, in separate yet adjoining locations. Here, we summarize these findings and propose that quiescent and active stem cell populations have separate but cooperative functional roles.

Organoids are self‐organizing 3D structures grown from stem cells that recapitulate essential aspects of organ structure and function. Here, we describe a method to establish long‐term‐expanding human airway organoids from broncho‐alveolar resections or lavage material. The pseudostratified airway organoids consist of basal cells, functional multi‐ciliated cells, mucus‐producing secretory cells, and CC10‐secreting club cells. Airway organoids derived from cystic fibrosis (CF) patients allow assessment of CFTR function in an organoid swelling assay. Organoids established from lung cancer resections and metastasis biopsies retain tumor histopathology as well as cancer gene mutations and are amenable to drug screening. Respiratory syncytial virus (RSV) infection recapitulates central disease features, dramatically increases organoid cell motility via the non‐structural viral NS2 protein, and preferentially recruits neutrophils upon co‐culturing. We conclude that human airway organoids represent versatile models for the in vitro study of hereditary, malignant, and infectious pulmonary disease. To date, persistent in vitro culture of adult human lung epithelium remains elusive. In this methods resource article, culture conditions to maintain three‐dimensional pulmonary tissue long‐term are reported and applied to recapitulate related diseases. Three‐dimensional human pulmonary tissue culture allows for investigation of hereditary diseases.

To assess the critical role of Wnt signals in intestinal crypts, we generated transgenic mice ectopically expressing Dickkopf1 (Dkk1), a secreted Wnt inhibitor. We find that epithelial proliferation is greatly reduced coincidentally with the loss of crypts. Although enterocyte differentiation appears unaffected, secretory cell lineages are largely absent. Disrupted intestinal homeostasis is reflected by an absence of nuclear beta-catenin, inhibition of c-myc expression, and subsequent up-regulation of p21CIP1/WAF1. Thus, our data are the first to establish a direct requirement for Wnt ligands in driving proliferation in the intestinal epithelium, and also define an unexpected role for Wnts in controlling secretory cell differentiation.

BACKGROUND: Cardiomyocytes derived from human embryonic stem (hES) cells could be useful in restoring heart function after myocardial infarction or in heart failure. Here, we induced cardiomyocyte differentiation of hES cells by a novel method and compared their electrophysiological properties and coupling with those of primary human fetal cardiomyocytes. METHODS AND RESULTS: hES cells were cocultured with visceral-endoderm (VE)-like cells from the mouse. This initiated differentiation to beating muscle. Sarcomeric marker proteins, chronotropic responses, and ion channel expression and function were typical of cardiomyocytes. Electrophysiology demonstrated that most cells resembled human fetal ventricular cells. Real-time intracellular calcium measurements, Lucifer yellow injection, and connexin 43 expression demonstrated that fetal and hES-derived cardiomyocytes are coupled by gap junctions in culture. Inhibition of electrical responses by verapamil demonstrated the presence of functional alpha1c-calcium ion channels. CONCLUSIONS: This is the first demonstration of induction of cardiomyocyte differentiation in hES cells that do not undergo spontaneous cardiogenesis. It provides a model for the study of human cardiomyocytes in culture and could be a step forward in the development of cardiomyocyte transplantation therapies.

The spatial organization of chromosomes influences many nuclear processes including gene expression. The cohesin complex shapes the 3D genome by looping together CTCF sites along chromosomes. We show here that chromatin loop size can be increased and that the duration with which cohesin embraces DNA determines the degree to which loops are enlarged. Cohesin's DNA release factor WAPL restricts this loop extension and also prevents looping between incorrectly oriented CTCF sites. We reveal that the SCC2/SCC4 complex promotes the extension of chromatin loops and the formation of topologically associated domains (TADs). Our data support the model that cohesin structures chromosomes through the processive enlargement of loops and that TADs reflect polyclonal collections of loops in the making. Finally, we find that whereas cohesin promotes chromosomal looping, it rather limits nuclear compartmentalization. We conclude that the balanced activity of SCC2/SCC4 and WAPL enables cohesin to correctly structure chromosomes.

The fate of cells in the epiblast at prestreak and early primitive streak stages has been studied by injecting horseradish peroxidase (HRP) into single cells in situ of 6.7-day mouse embryos and identifying the labelled descendants at midstreak to neural plate stages after one day of culture. Ectoderm was composed of descendants of epiblast progenitors that had been located in the embryonic axis anterior to the primitive streak. Embryonic mesoderm was derived from all areas of the epiblast except the distal tip and the adjacent region anterior to it: the most anterior mesoderm cells originated posteriorly, traversing the primitive streak early; labelled cells in the posterior part of the streak at the neural plate stage were derived from extreme anterior axial and paraxial epiblast progenitors; head process cells were derived from epiblast at or near the anterior end of the primitive streak. Endoderm descendants were most frequently derived from a region that included, but extended beyond, the region producing the head process: descendants of epiblast were present in endoderm by the midstreak stage, as well as at later stages. Yolk sac and amnion mesoderm developed from posterolateral and posterior epiblast. The resulting fate map is essentially the same as those of the chick and urodele and indicates that, despite geometrical differences, topological fate relationships are conserved among these vertebrates. Clonal descendants were not necessarily confined to a single germ layer or to extraembryonic mesoderm, indicating that these lineages are not separated at the beginning of gastrulation. The embryonic axis lengthened up to the neural plate stage by (1) elongation of the primitive streak through progressive incorporation of the expanding lateral and initially more anterior regions of epiblast and, (2) expansion of the region of epiblast immediately cranial to the anterior end of the primitive streak. The population doubling time of labelled cells was 7.5 h; a calculated 43% were in, or had completed, a 4th cell cycle, and no statistically significant regional differences in the number of descendants were found. This clonal analysis also showed that (1) growth in the epiblast was noncoherent and in most regions anisotropic and directed towards the primitive streak and (2) the midline did not act as a barrier to clonal spread, either in the epiblast in the anterior half of the axis or in the primitive streak. These results taken together with the fate map indicate that, while individual cells in the epiblast sheet behave independently with respect to their neighbours, morphogenetic movement during germ layer formation is coordinated in the population as a whole.

Phenotypic variation is ubiquitous in biology and is often traceable to underlying genetic and environmental variation. However, even genetically identical cells in identical environments display variable phenotypes. Stochastic gene expression, or gene expression "noise," has been suggested as a major source of this variability, and its physiological consequences have been topics of intense research for the last decade. Several recent studies have measured variability in protein and messenger RNA levels, and they have discovered strong connections between noise and gene regulation mechanisms. When integrated with discrete stochastic models, measurements of cell-to-cell variability provide a sensitive "fingerprint" with which to explore fundamental questions of gene regulation. In this review, we highlight several studies that used gene expression variability to develop a quantitative understanding of the mechanisms and dynamics of gene regulation.

BACKGROUND & AIMS: We previously established long-term, 3-dimensional culture of organoids from mouse tissues (intestine, stomach, pancreas, and liver) and human intestine and pancreas. Here we describe conditions required for long-term 3-dimensional culture of human gastric stem cells. The technology can be applied to study the epithelial response to infection with Helicobacter pylori. METHODS: We generated organoids from surgical samples of human gastric corpus. Culture conditions were developed based on those for the mouse gastric and human intestinal systems. We used microinjection to infect the organoids with H pylori. Epithelial responses were measured using microarray and quantitative polymerase chain reaction analyses. RESULTS: Human gastric cells were expanded indefinitely in 3-dimensional cultures. We cultured cells from healthy gastric tissues, single-sorted stem cells, or tumor tissues. Organoids maintained many characteristics of their respective tissues based on their histology, expression of markers, and euploidy. Organoids from healthy tissue expressed markers of 4 lineages of the stomach and self-organized into gland and pit domains. They could be directed to specifically express either lineages of the gastric gland, or the gastric pit, by addition of nicotinamide and withdrawal of WNT. Although gastric pit lineages had only marginal reactions to bacterial infection, gastric gland lineages mounted a strong inflammatory response. CONCLUSIONS: We developed a system to culture human gastric organoids. This system can be used to study H pylori infection and other gastric pathologies.

MicroRNAs (miRNAs) are small non-coding RNA molecules that post-transcriptionally regulate gene expression by base-pairing to mRNAs. Hundreds of miRNAs have been identified in various multicellular organisms and many miRNAs are evolutionarily conserved. Although the biological functions of most miRNAs are unknown, miRNAs are predicted to regulate up to 30% of the genes within the human genome. Gradually, we are beginning to understand the functions of individual miRNAs and the general function of miRNA action. Here, we review the recent advances in miRNA biology in animals. Particularly, we focus on the roles of miRNAs in vertebrate development and disease.

Although Apc is well characterized as a tumor-suppressor gene in the intestine, the precise mechanism of this suppression remains to be defined. Using a novel inducible Ahcre transgenic line in conjunction with a loxP-flanked Apc allele we, show that loss of Apc acutely activates Wnt signaling through the nuclear accumulation of beta-catenin. Coincidentally, it perturbs differentiation, migration, proliferation, and apoptosis, such that Apc-deficient cells maintain a "crypt progenitor-like" phenotype. Critically, for the first time we confirm a series of Wnt target molecules in an in vivo setting and also identify a series of new candidate targets within the same setting.

Little is known about the signaling mechanisms that determine the highly regular patterning of the intestinal epithelium into crypts and villi. With the use of mouse models, we show that bone morphogenetic protein (BMP)-4 expression occurs exclusively in the intravillus mesenchyme. Villus epithelial cells respond to the BMP signal. Inhibition of BMP signaling by transgenic expression of noggin results in the formation of numerous ectopic crypt units perpendicular to the crypt-villus axis. These changes phenocopy the intestinal histopathology of patients with the cancer predisposition syndrome juvenile polyposis (JP), including the frequent occurrence of intraepithelial neoplasia. Many JP cases are known to harbor mutations in BMP pathway genes. These data indicate that intestinal BMP signaling represses de novo crypt formation and polyp growth.