Domestication modulates the expression of genes involved in neurogenesis in high-quality eggs of Sander lucioperca.

Pikeperch, Sander lucioperca, is a species of high interest to the aquaculture. The expansion of its production can only be achieved by furthering domestication level. However, the mechanisms driving the domestication process in finfishes are poorly understood. Transcriptome profiling of eggs was found to be a useful tool allowing understanding of the domestication process in teleosts. In this study, using next-generation sequencing, the first pikeperch transcriptome has been generated as well as pikeperch-specific microarray comprising 35,343 unique probes. Next, we performed transcriptome profiling of eggs obtained from wild and domesticated populations. We found 710 differentially expressed genes that were linked mostly to nervous system development. These results provide new insights into processes that are directly involved in the domestication of finfishes. It can be suggested that all the identified processes were predetermined by the maternally derived set of genes contained in the unfertilized eggs. This allows us to suggest that fish behavior, along with many other processes, can be predetermined at the cellular level and may have significant implications on the adaptation of cultured fish to the natural environment. This also allows to suggest that fish behavior should be considered as a very important pikeperch aquaculture selection trait.

Evaluation and classification of severity for 176 genes on an expanded carrier screening panel.

BACKGROUND: Disease severity is important when considering genes for inclusion on reproductive expanded carrier screening (ECS) panels. We applied a validated and previously published algorithm that classifies diseases into four severity categories (mild, moderate, severe, and profound) to 176 genes screened by ECS. Disease traits defining severity categories in the algorithm were then mapped to four severity-related ECS panel design criteria cited by the American College of Obstetricians and Gynecologists (ACOG). METHODS: Eight genetic counselors (GCs) and four medical geneticists (MDs) applied the severity algorithm to subsets of 176 genes. MDs and GCs then determined by group consensus how each of these disease traits mapped to ACOG severity criteria, enabling determination of the number of ACOG severity criteria met by each gene. RESULTS: Upon consensus GC and MD application of the severity algorithm, 68 (39%) genes were classified as profound, 71 (40%) as severe, 36 (20%) as moderate, and one (1%) as mild. After mapping of disease traits to ACOG severity criteria, 170 out of 176 genes (96.6%) were found to meet at least one of the four criteria, 129 genes (73.3%) met at least two, 73 genes (41.5%) met at least three, and 17 genes (9.7%) met all four. CONCLUSION: This study classified the severity of a large set of Mendelian genes by collaborative clinical expert application of a trait-based algorithm. Further, it operationalized difficult to interpret ACOG severity criteria via mapping of disease traits, thereby promoting consistency of ACOG criteria interpretation.

Developmental gene expression as a phylogenetic data class: support for the monophyly of Arachnopulmonata.

Despite application of genome-scale datasets, the phylogenetic placement of scorpions within arachnids remains contentious between two different phylogenetic data classes. Paleontologists continue to recover scorpions in a basally branching position, partly owing to their morphological similarity to extinct marine orders like Eurypterida (sea scorpions). Phylogenomic datasets consistently recover scorpions in a derived position, as the sister group of Tetrapulmonata (a clade of arachnids that includes spiders). To adjudicate between these hypotheses using a rare genomic change (RGC), we leveraged the recent discovery of ancient paralogy in spiders and scorpions to assess phylogenetic placement. We identified homologs of four transcription factors required for appendage patterning (dachshund, homothorax, extradenticle, and optomotor blind) in arthropods that are known to be duplicated in spiders. Using genomic resources for a spider, a scorpion, and a harvestman, we conducted gene tree analyses and assayed expression patterns of scorpion gene duplicates. Here we show that scorpions, like spiders, retain two copies of all four transcription factors, whereas arachnid orders like mites and harvestmen bear a single copy. A survey of embryonic expression patterns of the scorpion paralogs closely matches those of their spider counterparts, with one paralog consistently retaining the putatively ancestral pattern found in the harvestman, as well as the mite, and/or other outgroups. These data comprise a rare genomic change in chelicerate phylogeny supporting the inference of a distal placement of scorpions. Beyond demonstrating the diagnostic power of developmental genetic data as a phylogenetic data class, a derived placement of scorpions within the arachnids, together with an array of stem-group Paleozoic scorpions that occupied marine habitats, effectively rules out a scenario of a single colonization of terrestrial habitat within Chelicerata, even in tree topologies contrived to recover the monophyly of Arachnida.

Exploration of key regulators driving primary feather follicle induction in goose skin.

The primary feather follicles are universal skin appendages widely distributed in the skin of feathered birds. The morphogenesis and development of the primary feather follicles in goose skin remain largely unknown. Here, the induction of primary feather follicles in goose embryonic skin (pre-induction vs induction) was investigated by de novo transcriptome analyses to reveal 409 differentially expressed genes (DEGs). The DEGs were characterized to potentially regulate the de novo formation of feather follicle primordia consisting of placode (4 genes) and dermal condensate (12 genes), and the thickening of epidermis (5 genes) and dermal fibroblasts (17 genes), respectively. Further analyses enriched DEGs into GO terms represented as cell adhesion and KEGG pathways including Wnt and Hedgehog signaling pathways that are highly correlated with cell communication and molecular regulation. Six selected Wnt pathway genes were detected by qPCR with up-regulation in goose skin during the induction of primary feather follicles. The localization of WNT16, SFRP1 and FRZB by in situ hybridization showed weak expression in the primary feather primordia, whereas FZD1, LEF1 and DKK1 were expressed initially in the inter-follicular skin and feather follicle primordia, then mainly restricted in the feather primordia. The spatial-temporal expression patterns indicate that Wnt pathway genes DKK1, FZD1 and LEF1 are the important regulators functioned in the induction of primary feather follicle in goose skin. The dynamic molecular changes and specific gene expression patterns revealed in this report provide the general knowledge of primary feather follicle and skin development in waterfowl, and contribute to further understand the diversity of hair and feather development beyond the mouse and chicken models.

The Developmental Gene Hypothesis for Punctuated Equilibrium: Combined Roles of Developmental Regulatory Genes and Transposable Elements.

Theories of the genetics underlying punctuated equilibrium (PE) have been vague to date. Here the developmental gene hypothesis is proposed, which states that: 1) developmental regulatory (DevReg) genes are responsible for the orchestration of metazoan morphogenesis and their extreme conservation and mutation intolerance generates the equilibrium or stasis present throughout much of the fossil record and 2) the accumulation of regulatory elements and recombination within these same genes-often derived from transposable elements-drives punctuated bursts of morphological divergence and speciation across metazoa. This two-part hypothesis helps to explain the features that characterize PE, providing a theoretical genetic basis for the once-controversial theory. Also see the video abstract here https://youtu.be/C-fu-ks5yDs.

A dual role of dLsd1 in oogenesis: regulating developmental genes and repressing transposons.

The histone demethylase LSD1 is a key chromatin regulator that is often deregulated in cancer. Its ortholog, dLsd1 plays a crucial role in Drosophila oogenesis; however, our knowledge of dLsd1 function is insufficient to explain its role in the ovary. Here, we have performed genome-wide analysis of dLsd1 binding in the ovary, and we document that dLsd1 is preferentially associated to the transcription start site of developmental genes. We uncovered an unanticipated interplay between dLsd1 and the GATA transcription factor Serpent and we report an unexpected role for Serpent in oogenesis. Besides, our transcriptomic data show that reducing dLsd1 levels results in ectopic transposable elements (TE) expression correlated with changes in H3K4me2 and H3K9me2 at TE loci. In addition, our results suggest that dLsd1 is required for Piwi dependent TE silencing. Hence, we propose that dLsd1 plays crucial roles in establishing specific gene expression programs and in repressing transposons during oogenesis.

Programmed DNA elimination of germline development genes in songbirds.

In some eukaryotes, germline and somatic genomes differ dramatically in their composition. Here we characterise a major germline-soma dissimilarity caused by a germline-restricted chromosome (GRC) in songbirds. We show that the zebra finch GRC contains >115 genes paralogous to single-copy genes on 18 autosomes and the Z chromosome, and is enriched in genes involved in female gonad development. Many genes are likely functional, evidenced by expression in testes and ovaries at the RNA and protein level. Using comparative genomics, we show that genes have been added to the GRC over millions of years of evolution, with embryonic development genes bicc1 and trim71 dating to the ancestor of songbirds and dozens of other genes added very recently. The somatic elimination of this evolutionarily dynamic chromosome in songbirds implies a unique mechanism to minimise genetic conflict between germline and soma, relevant to antagonistic pleiotropy, an evolutionary process underlying ageing and sexual traits.

Compound Dynamics and Combinatorial Patterns of Amino Acid Repeats Encode a System of Evolutionary and Developmental Markers.

Homopolymeric amino acid repeats (AARs) like polyalanine (polyA) and polyglutamine (polyQ) in some developmental proteins (DPs) regulate certain aspects of organismal morphology and behavior, suggesting an evolutionary role for AARs as developmental "tuning knobs." It is still unclear, however, whether these are occasional protein-specific phenomena or hints at the existence of a whole AAR-based regulatory system in DPs. Using novel approaches to trace their functional and evolutionary history, we find quantitative evidence supporting a generalized, combinatorial role of AARs in developmental processes with evolutionary implications. We observe nonrandom AAR distributions and combinations in HOX and other DPs, as well as in their interactomes, defining elements of a proteome-wide combinatorial functional code whereby different AARs and their combinations appear preferentially in proteins involved in the development of specific organs/systems. Such functional associations can be either static or display detectable evolutionary dynamics. These findings suggest that progressive changes in AAR occurrence/combination, by altering embryonic development, may have contributed to taxonomic divergence, leaving detectable traces in the evolutionary history of proteomes. Consistent with this hypothesis, we find that the evolutionary trajectories of the 20 AARs in eukaryotic proteomes are highly interrelated and their individual or compound dynamics can sharply mark taxonomic boundaries, or display clock-like trends, carrying overall a strong phylogenetic signal. These findings provide quantitative evidence and an interpretive framework outlining a combinatorial system of AARs whose compound dynamics mark at the same time DP functions and evolutionary transitions.

Accurate genome-wide predictions of spatio-temporal gene expression during embryonic development.

Comprehensive information on the timing and location of gene expression is fundamental to our understanding of embryonic development and tissue formation. While high-throughput in situ hybridization projects provide invaluable information about developmental gene expression patterns for model organisms like Drosophila, the output of these experiments is primarily qualitative, and a high proportion of protein coding genes and most non-coding genes lack any annotation. Accurate data-centric predictions of spatio-temporal gene expression will therefore complement current in situ hybridization efforts. Here, we applied a machine learning approach by training models on all public gene expression and chromatin data, even from whole-organism experiments, to provide genome-wide, quantitative spatio-temporal predictions for all genes. We developed structured in silico nano-dissection, a computational approach that predicts gene expression in >200 tissue-developmental stages. The algorithm integrates expression signals from a compendium of 6,378 genome-wide expression and chromatin profiling experiments in a cell lineage-aware fashion. We systematically evaluated our performance via cross-validation and experimentally confirmed 22 new predictions for four different embryonic tissues. The model also predicts complex, multi-tissue expression and developmental regulation with high accuracy. We further show the potential of applying these genome-wide predictions to extract tissue specificity signals from non-tissue-dissected experiments, and to prioritize tissues and stages for disease modeling. This resource, together with the exploratory tools are freely available at our webserver http://find.princeton.edu, which provides a valuable tool for a range of applications, from predicting spatio-temporal expression patterns to recognizing tissue signatures from differential gene expression profiles.

Retinoic acid signaling pathways.

Retinoic acid (RA), a metabolite of retinol (vitamin A), functions as a ligand for nuclear RA receptors (RARs) that regulate development of chordate animals. RA-RARs can activate or repress transcription of key developmental genes. Genetic studies in mouse and zebrafish embryos that are deficient in RA-generating enzymes or RARs have been instrumental in identifying RA functions, revealing that RA signaling regulates development of many organs and tissues, including the body axis, spinal cord, forelimbs, heart, eye and reproductive tract. An understanding of the normal functions of RA signaling during development will guide efforts for use of RA as a therapeutic agent to improve human health. Here, we provide an overview of RA signaling and highlight its key functions during development.

Bone morphology is regulated modularly by global and regional genetic programs.

Bone protrusions provide stable anchoring sites for ligaments and tendons and define the unique morphology of each long bone. Despite their importance, the mechanism by which superstructures are patterned is unknown. Here, we identify components of the genetic program that control the patterning of Sox9+/Scx+ superstructure progenitors in mouse and show that this program includes both global and regional regulatory modules. Using light-sheet fluorescence microscopy combined with genetic lineage labeling, we mapped the broad contribution of the Sox9+/Scx+ progenitors to the formation of bone superstructures. Then, by combining literature-based evidence, comparative transcriptomic analysis and genetic mouse models, we identified Gli3 as a global regulator of superstructure patterning, whereas Pbx1, Pbx2, Hoxa11 and Hoxd11 act as proximal and distal regulators, respectively. Moreover, by demonstrating a dose-dependent pattern regulation in Gli3 and Pbx1 compound mutations, we show that the global and regional regulatory modules work in a coordinated manner. Collectively, our results provide strong evidence for genetic regulation of superstructure patterning, which further supports the notion that long bone development is a modular process.This article has an associated 'The people behind the papers' interview.

The recognition of development-related genes in the testis and MAGs of time-series Harmonia axyridis adults using a time-series analysis by RNA-seq.

As an important natural enemy of aphids and other pests in agriculture, Harmonia axyridis has been widely used in classic biological control. The testis and male accessory gland (MAG) of H. axyridis are typically associated with the ability of egg-laying by multiple mating and influence mass artificial breeding for biological control. Development-related genes might impact the growth of adult testis and MAGs, but this has not yet been reported. Here, we use an integrative time-series analysis with RNA-seq in the testis and MAGs of H. axyridis adults to detect development-related genes. >10.3 Gb of clean data were obtained. All differentially expressed genes (DEG) between samples were analyzed with five databases for building a DEG annotation database. By combining previous reports, the DEG annotation database, and gene expression level changes in four different stages, the 26 DEGs related to testis and MAGs development were identified. To validate the expression profile, 15 random DEGs were chosen to perform RT-qPCR, and they were all in accordance with the RNA-Seq results. Taken together, this study established a robust pipeline for the discovery of key genes using RNA-seq and allowed for the class identification of development-related genes in the testis and MAGs for comprehensive characterization.

Zebrafish embryos exposed to deltamethrin exhibit abnormalities despite induced expression of related genes ( you, you-too, momo and u-boot).

Evaluation of the toxic effects of a widely used synthetic pyrethroid, deltamethrin (DM), was carried out in this study. This pesticide is preferred for pest control because of its low environmental persistence and toxicity. We investigated the expression pattern of four genes, namely, you ( you), yot ( you-too), momo ( mom) and ubo ( u-boot) during early development of zebrafish, that is, from 12 hpf to 48 hpf stages. These stages are selected as most of the important developmental aspects take place during this period. All four genes are known to play a vital role in development of notochord and somites. To understand the effect of DM on development, embryos of 4 hpf stage were exposed to two concentrations (100 and 200 µg/L) of DM, and observations were made at 12, 24 and 48 hpf stages. Our earlier studies have shown phenotypic abnormalities such as notochord bending, tail deformation, yolk sac and pericardial edema, lightening of body and eye pigmentation and interfered in somite patterning, during these stages of development. Understanding the relationship of phenotypic abnormalities with these four genes has been our primary objective. These four genes were analyzed by Reverse transcription (RT)-polymerase chain reaction and intensity of the bands has shown induction in their expression after exposure to the toxicant. In spite of the expression of genes, it was noticed that DM caused abnormalities. It can be said from the results that translational pathway could have been affected.

Frequent monoallelic or skewed expression for developmental genes in CNS-derived cells and evidence for balancing selection.

Cellular mosaicism due to monoallelic autosomal expression (MAE), with cell selection during development, is becoming increasingly recognized as prevalent in mammals, leading to interest in understanding its extent and mechanism(s). We report here use of clonal cell lines derived from the CNS of adult female [Formula: see text] hybrid (C57BL/6 X JF1) mice to characterize MAE as neural stem cells (nscs) differentiate to astrocyte-like cells (asls). We found that different subsets of genes show MAE in the two populations of cells; in each case, there is strong enrichment for genes specific to the respective developmental state. Genes that exhibit MAE are 22% of nsc-specific genes and 26% of asl-specific genes. Moreover, the promoters of genes with MAE have reduced CpG dinucleotides but increased CpG differences between the two parental mouse strains. Extending the study of variability to wild populations of mice, we found evidence for balancing selection as a contributing force in evolution of those genes showing developmental specificity (i.e., expressed in either nsc or asl), not just for genes showing MAE. Furthermore, we found that genes showing skewed allelic expression (SKE) were similarly enriched among cell type-specific genes and also showed a heightened probability of balancing selection. Thus, developmental stage-specific genes and genes with MAE or SKE seem to make up overlapping classes subject to selection for increased diversity. The implications of these results for development and evolution are discussed in the context of a model with stochastic epigenetic modifications taking place only during a relatively brief developmental window.

Identification of reference genes suitable for RT-qPCR studies of murine gastrulation and patterning.

Quantitative reverse transcriptase PCR (RT-qPCR), a powerful and efficient means of rapidly comparing gene expression between experimental conditions, is routinely used as a phenotyping tool in developmental biology. The accurate comparison of gene expression across multiple embryonic stages requires normalisation to reference genes that have stable expression across the time points to be examined. As the embryo and its constituent tissues undergo rapid growth and differentiation during development, reference genes known to be stable across some time points cannot be assumed to be stable across all developmental stages. The immediate post-implantation events of gastrulation and patterning are characterised by a rapid expansion in cell number and increasing specialisation of cells. The optimal reference genes for comparative gene expression studies at these specific stages have not been experimentally identified. In this study, the expression of five commonly used reference genes (H2afz, Ubc, Actb, Tbp and Gapdh) was measured across murine gastrulation and patterning (6.5-9.5 dpc) and analysed with the normalisation tools geNorm, Bestkeeper and Normfinder. The results, validated by RT-qPCR analysis of two genes with well-documented expression patterns across these stages, indicated the best strategy for RT-qPCR studies spanning murine gastrulation and patterning utilises the concurrent reference genes H2afz and Ubc.

Developmental origin and morphogenesis of the diaphragm, an essential mammalian muscle.

The diaphragm is a mammalian skeletal muscle essential for respiration and for separating the thoracic and abdominal cavities. Development of the diaphragm requires the coordinated development of muscle, muscle connective tissue, tendon, nerves, and vasculature that derive from different embryonic sources. However, defects in diaphragm development are common and the cause of an often deadly birth defect, Congenital Diaphragmatic Hernia (CDH). Here we comprehensively describe the normal developmental origin and complex spatial-temporal relationship between the different developing tissues to form a functional diaphragm using a developmental series of mouse embryos genetically and immunofluorescently labeled and analyzed in whole mount. We find that the earliest developmental events are the emigration of muscle progenitors from cervical somites followed by the projection of phrenic nerve axons from the cervical neural tube. Muscle progenitors and phrenic nerve target the pleuroperitoneal folds (PPFs), transient pyramidal-shaped structures that form between the thoracic and abdominal cavities. Subsequently, the PPFs expand across the surface of the liver to give rise to the muscle connective tissue and central tendon, and the leading edge of their expansion precedes muscle morphogenesis, formation of the vascular network, and outgrowth and branching of the phrenic nerve. Thus development and morphogenesis of the PPFs is critical for diaphragm formation. In addition, our data indicate that the earliest events in diaphragm development are critical for the etiology of CDH and instrumental to the evolution of the diaphragm. CDH initiates prior to E12.5 in mouse and suggests that defects in the early PPF formation or their ability to recruit muscle are an important source of CDH. Also, the recruitment of muscle progenitors from cervical somites to the nascent PPFs is uniquely mammalian and a key developmental innovation essential for the evolution of the muscularized diaphragm.

Arabidopsis HEAT SHOCK TRANSCRIPTION FACTORA1b regulates multiple developmental genes under benign and stress conditions.

In Arabidopsis thaliana, HEAT SHOCK TRANSCRIPTION FACTORA1b (HSFA1b) controls resistance to environmental stress and is a determinant of reproductive fitness by influencing seed yield. To understand how HSFA1b achieves this, we surveyed its genome-wide targets (ChIP-seq) and its impact on the transcriptome (RNA-seq) under non-stress (NS), heat stress (HS) in the wild type, and in HSFA1b-overexpressing plants under NS. A total of 952 differentially expressed HSFA1b-targeted genes were identified, of which at least 85 are development associated and were bound predominantly under NS. A further 1780 genes were differentially expressed but not bound by HSFA1b, of which 281 were classified as having development-associated functions. These genes are indirectly regulated through a hierarchical network of 27 transcription factors (TFs). Furthermore, we identified 480 natural antisense non-coding RNA (cisNAT) genes bound by HSFA1b, defining a further mode of indirect regulation. Finally, HSFA1b-targeted genomic features not only harboured heat shock elements, but also MADS box, LEAFY, and G-Box promoter motifs. This revealed that HSFA1b is one of eight TFs that target a common group of stress defence and developmental genes. We propose that HSFA1b transduces environmental cues to many stress tolerance and developmental genes to allow plants to adjust their growth and development continually in a varying environment.

Conserved non-coding elements: developmental gene regulation meets genome organization.

Comparative genomics has revealed a class of non-protein-coding genomic sequences that display an extraordinary degree of conservation between two or more organisms, regularly exceeding that found within protein-coding exons. These elements, collectively referred to as conserved non-coding elements (CNEs), are non-randomly distributed across chromosomes and tend to cluster in the vicinity of genes with regulatory roles in multicellular development and differentiation. CNEs are organized into functional ensembles called genomic regulatory blocks-dense clusters of elements that collectively coordinate the expression of shared target genes, and whose span in many cases coincides with topologically associated domains. CNEs display sequence properties that set them apart from other sequences under constraint, and have recently been proposed as useful markers for the reconstruction of the evolutionary history of organisms. Disruption of several of these elements is known to contribute to diseases linked with development, and cancer. The emergence, evolutionary dynamics and functions of CNEs still remain poorly understood, and new approaches are required to enable comprehensive CNE identification and characterization. Here, we review current knowledge and identify challenges that need to be tackled to resolve the impasse in understanding extreme non-coding conservation.

Bivalent complexes of PRC1 with orthologs of BRD4 and MOZ/MORF target developmental genes in Drosophila.

Regulatory decisions in Drosophila require Polycomb group (PcG) proteins to maintain the silent state and Trithorax group (TrxG) proteins to oppose silencing. Since PcG and TrxG are ubiquitous and lack apparent sequence specificity, a long-standing model is that targeting occurs via protein interactions; for instance, between repressors and PcG proteins. Instead, we found that Pc-repressive complex 1 (PRC1) purifies with coactivators Fs(1)h [female sterile (1) homeotic] and Enok/Br140 during embryogenesis. Fs(1)h is a TrxG member and the ortholog of BRD4, a bromodomain protein that binds to acetylated histones and is a key transcriptional coactivator in mammals. Enok and Br140, another bromodomain protein, are orthologous to subunits of a mammalian MOZ/MORF acetyltransferase complex. Here we confirm PRC1-Br140 and PRC1-Fs(1)h interactions and identify their genomic binding sites. PRC1-Br140 bind developmental genes in fly embryos, with analogous co-occupancy of PRC1 and a Br140 ortholog, BRD1, at bivalent loci in human embryonic stem (ES) cells. We propose that identification of PRC1-Br140 "bivalent complexes" in fly embryos supports and extends the bivalency model posited in mammalian cells, in which the coexistence of H3K4me3 and H3K27me3 at developmental promoters represents a poised transcriptional state. We further speculate that local competition between acetylation and deacetylation may play a critical role in the resolution of bivalent protein complexes during development.