Comparative analysis of embryo proper and suspensor transcriptomes in plant embryos with different morphologies.

An important question is what genes govern the differentiation of plant embryos into suspensor and embryo proper regions following fertilization and division of the zygote. We compared embryo proper and suspensor transcriptomes of four plants that vary in embryo morphology within the suspensor region. We determined that genes encoding enzymes in several metabolic pathways leading to the formation of hormones, such as gibberellic acid, and other metabolites are up-regulated in giant scarlet runner bean and common bean suspensors. Genes involved in transport and Golgi body organization are up-regulated within the suspensors of these plants as well, strengthening the view that giant specialized suspensors serve as a hormone factory and a conduit for transferring substances to the developing embryo proper. By contrast, genes controlling transcriptional regulation, development, and cell division are up-regulated primarily within the embryo proper. Transcriptomes from less specialized soybean and Arabidopsis suspensors demonstrated that fewer genes encoding metabolic enzymes and hormones are up-regulated. Genes active in the embryo proper, however, are functionally similar to those active in scarlet runner bean and common bean embryo proper regions. We uncovered a set of suspensor- and embryo proper-specific transcription factors (TFs) that are shared by all embryos irrespective of morphology, suggesting that they are involved in early differentiation processes common to all plants. Chromatin immunoprecipitation sequencing (ChIP-Seq) experiments with scarlet runner bean and soybean WOX9, an up-regulated suspensor TF, gained entry into a regulatory network important for suspensor development irrespective of morphology.

Genome-wide Hi-C analyses in wild-type and mutants reveal high-resolution chromatin interactions in Arabidopsis.

Chromosomes form 3D structures that are critical to the regulation of cellular and genetic processes. Here, we present a study of global chromatin interaction patterns in Arabidopsis thaliana. Our genome-wide approach confirmed interactions that were previously observed by other methods as well as uncovered long-range interactions such as those among small heterochromatic regions embedded in euchromatic arms. We also found that interactions are correlated with various epigenetic marks that are localized in active or silenced chromatin. Arabidopsis chromosomes do not contain large local interactive domains that resemble the topological domains described in animals but, instead, contain relatively small interactive regions scattered around the genome that contain H3K27me3 or H3K9me2. We generated interaction maps in mutants that are defective in specific epigenetic pathways and found altered interaction patterns that correlate with changes in the epigenome. These analyses provide further insights into molecular mechanisms of epigenetic regulation of the genome.

Large-scale heterochromatin remodeling linked to overreplication-associated DNA damage.

Previously, we have shown that loss of the histone 3 lysine 27 (H3K27) monomethyltransferases ARABIDOPSIS TRITHORAX-RELATED 5 (ATXR5) and ATXR6 (ATXR6) results in the overreplication of heterochromatin. Here we show that the overreplication results in DNA damage and extensive chromocenter remodeling into unique structures we have named "overreplication-associated centers" (RACs). RACs have a highly ordered structure with an outer layer of condensed heterochromatin, an inner layer enriched in the histone variant H2AX, and a low-density core containing foci of phosphorylated H2AX (a marker of double-strand breaks) and the DNA-repair enzyme RAD51. atxr5,6 mutants are strongly affected by mutations in DNA repair, such as ATM and ATR. Because of its dense packaging and repetitive DNA sequence, heterochromatin is a challenging environment in which to repair DNA damage. Previous work in animals has shown that heterochromatic breaks are translocated out of the heterochromatic domain for repair. Our results show that atxr5,6 mutants use a variation on this strategy for repairing heterochromatic DNA damage. Rather than being moved to adjacent euchromatic regions, as in animals, heterochromatin undergoes large-scale remodeling to create a compartment with low chromatin density.

Chromosome-level genome assembly and transcriptome of the green alga Chromochloris zofingiensis illuminates astaxanthin production.

Microalgae have potential to help meet energy and food demands without exacerbating environmental problems. There is interest in the unicellular green alga Chromochloris zofingiensis, because it produces lipids for biofuels and a highly valuable carotenoid nutraceutical, astaxanthin. To advance understanding of its biology and facilitate commercial development, we present a C. zofingiensis chromosome-level nuclear genome, organelle genomes, and transcriptome from diverse growth conditions. The assembly, derived from a combination of short- and long-read sequencing in conjunction with optical mapping, revealed a compact genome of ∼58 Mbp distributed over 19 chromosomes containing 15,274 predicted protein-coding genes. The genome has uniform gene density over chromosomes, low repetitive sequence content (∼6%), and a high fraction of protein-coding sequence (∼39%) with relatively long coding exons and few coding introns. Functional annotation of gene models identified orthologous families for the majority (∼73%) of genes. Synteny analysis uncovered localized but scrambled blocks of genes in putative orthologous relationships with other green algae. Two genes encoding beta-ketolase (BKT), the key enzyme synthesizing astaxanthin, were found in the genome, and both were up-regulated by high light. Isolation and molecular analysis of astaxanthin-deficient mutants showed that BKT1 is required for the production of astaxanthin. Moreover, the transcriptome under high light exposure revealed candidate genes that could be involved in critical yet missing steps of astaxanthin biosynthesis, including ABC transporters, cytochrome P450 enzymes, and an acyltransferase. The high-quality genome and transcriptome provide insight into the green algal lineage and carotenoid production.

A naïve Bayesian classifier for identifying plant microRNAs.

MicroRNAs (miRNAs) are important regulatory molecules in eukaryotic organisms. Existing methods for the identification of mature miRNA sequences in plants rely extensively on the search for stem-loop structures, leading to high false negative rates. Here, we describe a probabilistic method for ranking putative plant miRNAs using a naïve Bayes classifier and its publicly available implementation. We use a number of properties to construct the classifier, including sequence length, number of observations, existence of detectable predicted miRNA* sequences, the distribution of nearby reads and mapping multiplicity. We apply the method to small RNA sequence data from soybean, peach, Arabidopsis and rice and provide experimental validation of several predictions in soybean. The approach performs well overall and strongly enriches for known miRNAs over other types of sequences. By utilizing a Bayesian approach to rank putative miRNAs, our method is able to score miRNAs that would be eliminated by other methods, such as those that have low counts or lack detectable miRNA* sequences. As a result, we are able to detect several soybean miRNA candidates, including some that are 24 nucleotides long, a class that is almost universally eliminated by other methods.

Genome-wide erasure of DNA methylation in mouse primordial germ cells is affected by AID deficiency.

Epigenetic reprogramming including demethylation of DNA occurs in mammalian primordial germ cells (PGCs) and in early embryos, and is important for the erasure of imprints and epimutations, and the return to pluripotency. The extent of this reprogramming and its molecular mechanisms are poorly understood. We previously showed that the cytidine deaminases AID and APOBEC1 can deaminate 5-methylcytosine in vitro and in Escherichia coli, and in the mouse are expressed in tissues in which demethylation occurs. Here we profiled DNA methylation throughout the genome by unbiased bisulphite next generation sequencing in wild-type and AID-deficient mouse PGCs at embryonic day (E)13.5. Wild-type PGCs revealed marked genome-wide erasure of methylation to a level below that of methylation deficient (Np95(-/-), also called Uhrf1(-/-)) embryonic stem cells, with female PGCs being less methylated than male ones. By contrast, AID-deficient PGCs were up to three times more methylated than wild-type ones; this substantial difference occurred throughout the genome, with introns, intergenic regions and transposons being relatively more methylated than exons. Relative hypermethylation in AID-deficient PGCs was confirmed by analysis of individual loci in the genome. Our results reveal that erasure of DNA methylation in the germ line is a global process, hence limiting the potential for transgenerational epigenetic inheritance. AID deficiency interferes with genome-wide erasure of DNA methylation patterns, indicating that AID has a critical function in epigenetic reprogramming and potentially in restricting the inheritance of epimutations in mammals.

Relationship between nucleosome positioning and DNA methylation.

Nucleosomes compact and regulate access to DNA in the nucleus, and are composed of approximately 147 bases of DNA wrapped around a histone octamer. Here we report a genome-wide nucleosome positioning analysis of Arabidopsis thaliana using massively parallel sequencing of mononucleosomes. By combining this data with profiles of DNA methylation at single base resolution, we identified 10-base periodicities in the DNA methylation status of nucleosome-bound DNA and found that nucleosomal DNA was more highly methylated than flanking DNA. These results indicate that nucleosome positioning influences DNA methylation patterning throughout the genome and that DNA methyltransferases preferentially target nucleosome-bound DNA. We also observed similar trends in human nucleosomal DNA, indicating that the relationships between nucleosomes and DNA methyltransferases are conserved. Finally, as has been observed in animals, nucleosomes were highly enriched on exons, and preferentially positioned at intron-exon and exon-intron boundaries. RNA polymerase II (Pol II) was also enriched on exons relative to introns, consistent with the hypothesis that nucleosome positioning regulates Pol II processivity. DNA methylation is also enriched on exons, consistent with the targeting of DNA methylation to nucleosomes, and suggesting a role for DNA methylation in exon definition.

Evolutionary insights from de novo transcriptome assembly and SNP discovery in California white oaks.

BACKGROUND: Reference transcriptomes provide valuable resources for understanding evolution within and among species. We de novo assembled and annotated a reference transcriptome for Quercus lobata and Q. garryana and identified single-nucleotide polymorphisms (SNPs) to provide resources for forest genomicists studying this ecologically and economically important genus. We further performed preliminary analyses of genes important in interspecific divergent (positive) selection that might explain ecological differences among species, estimating rates of nonsynonymous to synonymous substitutions (d N/d S) and Fay and Wu's H. Functional classes of genes were tested for unusually high d N/d S or low H consistent with divergent positive selection. RESULTS: Our draft transcriptome is among the most complete for oaks, including 83,644 contigs (23,329 ≥ 1 kbp), 14,898 complete and 13,778 partial gene models, and functional annotations for 9,431 Arabidopsis orthologs and 19,365 contigs with Pfam hits. We identified 1.7 million possible sequence variants including 1.1 million high-quality diallelic SNPs - among the largest sets identified in any tree. 11 of 18 functional categories with significantly elevated d N/d S are involved in disease response, including 50+ genes with d N/d S > 1. Other high-d N/d S genes are involved in biotic response, flowering and growth, or regulatory processes. In contrast, median d N/d S was low (0.22), suggesting that purifying selection influences most genes. No functional categories have unusually low H. CONCLUSIONS: These results offer preliminary support for the hypothesis that divergent selection at pathogen resistance are important factors in species divergence in these hybridizing California oaks. Our transcriptome provides a solid foundation for future studies of gene expression, natural selection, and speciation in Quercus.

Systems biology approach in Chlamydomonas reveals connections between copper nutrition and multiple metabolic steps.

In this work, we query the Chlamydomonas reinhardtii copper regulon at a whole-genome level. Our RNA-Seq data simulation and analysis pipeline validated a 2-fold cutoff and 10 RPKM (reads per kilobase of mappable length per million mapped reads) (~1 mRNA per cell) to reveal 63 CRR1 targets plus another 86 copper-responsive genes. Proteomic and immunoblot analyses captured 25% of the corresponding proteins, whose abundance was also dependent on copper nutrition, validating transcriptional regulation as a major control mechanism for copper signaling in Chlamydomonas. The impact of copper deficiency on the expression of several O2-dependent enzymes included steps in lipid modification pathways. Quantitative lipid profiles indicated increased polyunsaturation of fatty acids on thylakoid membrane digalactosyldiglycerides, indicating a global impact of copper deficiency on the photosynthetic apparatus. Discovery of a putative plastid copper chaperone and a membrane protease in the thylakoid suggest a mechanism for blocking copper utilization in the chloroplast. We also found an example of copper sparing in the N assimilation pathway: the replacement of copper amine oxidase by a flavin-dependent backup enzyme. Forty percent of the targets are previously uncharacterized proteins, indicating considerable potential for new discovery in the biology of copper.

Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning.

Cytosine DNA methylation is important in regulating gene expression and in silencing transposons and other repetitive sequences. Recent genomic studies in Arabidopsis thaliana have revealed that many endogenous genes are methylated either within their promoters or within their transcribed regions, and that gene methylation is highly correlated with transcription levels. However, plants have different types of methylation controlled by different genetic pathways, and detailed information on the methylation status of each cytosine in any given genome is lacking. To this end, we generated a map at single-base-pair resolution of methylated cytosines for Arabidopsis, by combining bisulphite treatment of genomic DNA with ultra-high-throughput sequencing using the Illumina 1G Genome Analyser and Solexa sequencing technology. This approach, termed BS-Seq, unlike previous microarray-based methods, allows one to sensitively measure cytosine methylation on a genome-wide scale within specific sequence contexts. Here we describe methylation on previously inaccessible components of the genome and analyse the DNA methylation sequence composition and distribution. We also describe the effect of various DNA methylation mutants on genome-wide methylation patterns, and demonstrate that our newly developed library construction and computational methods can be applied to large genomes such as that of mouse.

Three acyltransferases and nitrogen-responsive regulator are implicated in nitrogen starvation-induced triacylglycerol accumulation in Chlamydomonas.

Algae have recently gained attention as a potential source for biodiesel; however, much is still unknown about the biological triggers that cause the production of triacylglycerols. We used RNA-Seq as a tool for discovering genes responsible for triacylglycerol (TAG) production in Chlamydomonas and for the regulatory components that activate the pathway. Three genes encoding acyltransferases, DGAT1, DGTT1, and PDAT1, are induced by nitrogen starvation and are likely to have a role in TAG accumulation based on their patterns of expression. DGAT1 and DGTT1 also show increased mRNA abundance in other TAG-accumulating conditions (minus sulfur, minus phosphorus, minus zinc, and minus iron). Insertional mutants, pdat1-1 and pdat1-2, accumulate 25% less TAG compared with the parent strain, CC-4425, which demonstrates the relevance of the trans-acylation pathway in Chlamydomonas. The biochemical functions of DGTT1 and PDAT1 were validated by rescue of oleic acid sensitivity and restoration of TAG accumulation in a yeast strain lacking all acyltransferase activity. Time course analyses suggest than a SQUAMOSA promoter-binding protein domain transcription factor, whose mRNA increases precede that of lipid biosynthesis genes like DGAT1, is a candidate regulator of the nitrogen deficiency responses. An insertional mutant, nrr1-1, accumulates only 50% of the TAG compared with the parental strain in nitrogen-starvation conditions and is unaffected by other nutrient stresses, suggesting the specificity of this regulator for nitrogen-deprivation conditions.

BS-Seeker2: a versatile aligning pipeline for bisulfite sequencing data.

BACKGROUND: DNA methylation is an important epigenetic modification involved in many biological processes. Bisulfite treatment coupled with high-throughput sequencing provides an effective approach for studying genome-wide DNA methylation at base resolution. Libraries such as whole genome bisulfite sequencing (WGBS) and reduced represented bisulfite sequencing (RRBS) are widely used for generating DNA methylomes, demanding efficient and versatile tools for aligning bisulfite sequencing data. RESULTS: We have developed BS-Seeker2, an updated version of BS Seeker, as a full pipeline for mapping bisulfite sequencing data and generating DNA methylomes. BS-Seeker2 improves mappability over existing aligners by using local alignment. It can also map reads from RRBS library by building special indexes with improved efficiency and accuracy. Moreover, BS-Seeker2 provides additional function for filtering out reads with incomplete bisulfite conversion, which is useful in minimizing the overestimation of DNA methylation levels. We also defined CGmap and ATCGmap file formats for full representations of DNA methylomes, as part of the outputs of BS-Seeker2 pipeline together with BAM and WIG files. CONCLUSIONS: Our evaluations on the performance show that BS-Seeker2 works efficiently and accurately for both WGBS data and RRBS data. BS-Seeker2 is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker2/ and the Galaxy server.

RNA-seq analysis of sulfur-deprived Chlamydomonas cells reveals aspects of acclimation critical for cell survival.

The Chlamydomonas reinhardtii transcriptome was characterized from nutrient-replete and sulfur-depleted wild-type and snrk2.1 mutant cells. This mutant is null for the regulatory Ser-Thr kinase SNRK2.1, which is required for acclimation of the alga to sulfur deprivation. The transcriptome analyses used microarray hybridization and RNA-seq technology. Quantitative RT-PCR evaluation of the results obtained by these techniques showed that RNA-seq reports a larger dynamic range of expression levels than do microarray hybridizations. Transcripts responsive to sulfur deprivation included those encoding proteins involved in sulfur acquisition and assimilation, synthesis of sulfur-containing metabolites, Cys degradation, and sulfur recycling. Furthermore, we noted potential modifications of cellular structures during sulfur deprivation, including the cell wall and complexes associated with the photosynthetic apparatus. Moreover, the data suggest that sulfur-deprived cells accumulate proteins with fewer sulfur-containing amino acids. Most of the sulfur deprivation responses are controlled by the SNRK2.1 protein kinase. The snrk2.1 mutant exhibits a set of unique responses during both sulfur-replete and sulfur-depleted conditions that are not observed in wild-type cells; the inability of this mutant to acclimate to S deprivation probably leads to elevated levels of singlet oxygen and severe oxidative stress, which ultimately causes cell death. The transcriptome results for wild-type and mutant cells strongly suggest the occurrence of massive changes in cellular physiology and metabolism as cells become depleted for sulfur and reveal aspects of acclimation that are likely critical for cell survival.

Conservation and divergence of methylation patterning in plants and animals.

Cytosine DNA methylation is a heritable epigenetic mark present in many eukaryotic organisms. Although DNA methylation likely has a conserved role in gene silencing, the levels and patterns of DNA methylation appear to vary drastically among different organisms. Here we used shotgun genomic bisulfite sequencing (BS-Seq) to compare DNA methylation in eight diverse plant and animal genomes. We found that patterns of methylation are very similar in flowering plants with methylated cytosines detected in all sequence contexts, whereas CG methylation predominates in animals. Vertebrates have methylation throughout the genome except for CpG islands. Gene body methylation is conserved with clear preference for exons in most organisms. Furthermore, genes appear to be the major target of methylation in Ciona and honey bee. Among the eight organisms, the green alga Chlamydomonas has the most unusual pattern of methylation, having non-CG methylation enriched in exons of genes rather than in repeats and transposons. In addition, the Dnmt1 cofactor Uhrf1 has a conserved function in maintaining CG methylation in both transposons and gene bodies in the mouse, Arabidopsis, and zebrafish genomes.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks.

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Algal Functional Annotation Tool: a web-based analysis suite to functionally interpret large gene lists using integrated annotation and expression data.

BACKGROUND: Progress in genome sequencing is proceeding at an exponential pace, and several new algal genomes are becoming available every year. One of the challenges facing the community is the association of protein sequences encoded in the genomes with biological function. While most genome assembly projects generate annotations for predicted protein sequences, they are usually limited and integrate functional terms from a limited number of databases. Another challenge is the use of annotations to interpret large lists of 'interesting' genes generated by genome-scale datasets. Previously, these gene lists had to be analyzed across several independent biological databases, often on a gene-by-gene basis. In contrast, several annotation databases, such as DAVID, integrate data from multiple functional databases and reveal underlying biological themes of large gene lists. While several such databases have been constructed for animals, none is currently available for the study of algae. Due to renewed interest in algae as potential sources of biofuels and the emergence of multiple algal genome sequences, a significant need has arisen for such a database to process the growing compendiums of algal genomic data. DESCRIPTION: The Algal Functional Annotation Tool is a web-based comprehensive analysis suite integrating annotation data from several pathway, ontology, and protein family databases. The current version provides annotation for the model alga Chlamydomonas reinhardtii, and in the future will include additional genomes. The site allows users to interpret large gene lists by identifying associated functional terms, and their enrichment. Additionally, expression data for several experimental conditions were compiled and analyzed to provide an expression-based enrichment search. A tool to search for functionally-related genes based on gene expression across these conditions is also provided. Other features include dynamic visualization of genes on KEGG pathway maps and batch gene identifier conversion. CONCLUSIONS: The Algal Functional Annotation Tool aims to provide an integrated data-mining environment for algal genomics by combining data from multiple annotation databases into a centralized tool. This site is designed to expedite the process of functional annotation and the interpretation of gene lists, such as those derived from high-throughput RNA-seq experiments. The tool is publicly available at http://pathways.mcdb.ucla.edu.

Copy number variation influences gene expression and metabolic traits in mice.

Copy number variants (CNVs) are genomic segments which are duplicated or deleted among different individuals. CNVs have been implicated in both Mendelian and complex traits, including immune and behavioral disorders, but the study of the mechanisms by which CNVs influence gene expression and clinical phenotypes in humans is complicated by the limited access to tissues and by population heterogeneity. We now report studies of the effect of 19 CNVs on gene expression and metabolic traits in a mouse intercross between strains C57BL/6J and C3H/HeJ. We found that 83% of genes predicted to occur within CNVs were differentially expressed. The expression of most CNV genes was correlated with copy number, but we also observed evidence that gene expression was altered in genes flanking CNVs, suggesting that CNVs may contain regulatory elements for these genes. Several CNVs mapped to hotspots, genomic regions influencing expression of tens or hundreds of genes. Several metabolic traits including cholesterol, triglycerides, glucose and body weight mapped to three CNVs in the genome, in mouse chromosomes 1, 4 and 17. Predicted CNV genes, such as Itlna, Defcr-1, Trim12 and Trim34 were highly correlated with these traits. Our results suggest that CNVs have a significant impact on gene expression and that CNVs may be playing a role in the mechanisms underlying metabolic traits in mice.

BS Seeker: precise mapping for bisulfite sequencing.

BACKGROUND: Bisulfite sequencing using next generation sequencers yields genome-wide measurements of DNA methylation at single nucleotide resolution. Traditional aligners are not designed for mapping bisulfite-treated reads, where the unmethylated Cs are converted to Ts. We have developed BS Seeker, an approach that converts the genome to a three-letter alphabet and uses Bowtie to align bisulfite-treated reads to a reference genome. It uses sequence tags to reduce mapping ambiguity. Post-processing of the alignments removes non-unique and low-quality mappings. RESULTS: We tested our aligner on synthetic data, a bisulfite-converted Arabidopsis library, and human libraries generated from two different experimental protocols. We evaluated the performance of our approach and compared it to other bisulfite aligners. The results demonstrate that among the aligners tested, BS Seeker is more versatile and faster. When mapping to the human genome, BS Seeker generates alignments significantly faster than RMAP and BSMAP. Furthermore, BS Seeker is the only alignment tool that can explicitly account for tags which are generated by certain library construction protocols. CONCLUSIONS: BS Seeker provides fast and accurate mapping of bisulfite-converted reads. It can work with BS reads generated from the two different experimental protocols, and is able to efficiently map reads to large mammalian genomes. The Python program is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker/BS_Seeker.html.

Whole-genome analysis of histone H3 lysine 27 trimethylation in Arabidopsis.

Trimethylation of histone H3 lysine 27 (H3K27me3) plays critical roles in regulating animal development, and in several cases, H3K27me3 is also required for the proper expression of developmentally important genes in plants. However, the extent to which H3K27me3 regulates plant genes on a genome-wide scale remains unknown. In addition, it is not clear whether the establishment and spreading of H3K27me3 occur through the same mechanisms in plants and animals. We identified regions containing H3K27me3 in the genome of the flowering plant Arabidopsis thaliana using a high-density whole-genome tiling microarray. The results suggest that H3K27me3 is a major silencing mechanism in plants that regulates an unexpectedly large number of genes in Arabidopsis (~4,400), and that the maintenance of H3K27me3 is largely independent of other epigenetic pathways, such as DNA methylation or RNA interference. Unlike in animals, where H3K27m3 occupies large genomic regions, in Arabidopsis, we found that H3K27m3 domains were largely restricted to the transcribed regions of single genes. Furthermore, unlike in animals systems, H3K27m3 domains were not preferentially associated with low-nucleosome density regions. The results suggest that different mechanisms may underlie the establishment and spreading of H3K27me3 in plants and animals.

Tissue-Specific Transcriptomes Reveal Gene Expression Trajectories in Two Maturing Skin Epithelial Layers in Zebrafish Embryos.

Epithelial cells are the building blocks of many organs, including skin. The vertebrate skin initially consists of two epithelial layers, the outer periderm and inner basal cell layers, which have distinct properties, functions, and fates. The embryonic periderm ultimately disappears during development, whereas basal cells proliferate to form the mature, stratified epidermis. Although much is known about mechanisms of homeostasis in mature skin, relatively little is known about the two cell types in pre-stratification skin. To define the similarities and distinctions between periderm and basal skin epithelial cells, we purified them from zebrafish at early development stages and deeply profiled their gene expression. These analyses identified groups of genes whose tissue enrichment changed at each stage, defining gene flow dynamics of maturing vertebrate epithelia. At each of 52 and 72 hr post-fertilization (hpf), more than 60% of genes enriched in skin cells were similarly expressed in both layers, indicating that they were common epithelial genes, but many others were enriched in one layer or the other. Both expected and novel genes were enriched in periderm and basal cell layers. Genes encoding extracellular matrix, junctional, cytoskeletal, and signaling proteins were prominent among those distinguishing the two epithelial cell types. In situ hybridization and BAC transgenes confirmed our expression data and provided new tools to study zebrafish skin. Collectively, these data provide a resource for studying common and distinguishing features of maturing epithelia.