Search

Sotos syndrome (SS), the most common overgrowth with intellectual disability (OGID) disorder, is caused by inactivating germline mutations of NSD1, which encodes a histone H3 lysine 36 methyltransferase. To understand how NSD1 inactivation deregulates transcription and DNA methylation (DNAm), and to explore how these abnormalities affect human development, we profiled transcription and DNAm in SS patients and healthy control individuals. We identified a transcriptional signature that distinguishes individuals with SS from controls and was also deregulated in NSD1-mutated cancers. Most abnormally expressed genes displayed reduced expression in SS; these downregulated genes consisted mostly of bivalent genes and were enriched for regulators of development and neural synapse function. DNA hypomethylation was strongly enriched within promoters of transcriptionally deregulated genes: overexpressed genes displayed hypomethylation at their transcription start sites while underexpressed genes featured hypomethylation at polycomb binding sites within their promoter CpG island shores. SS patients featured accelerated molecular aging at the levels of both transcription and DNAm. Overall, these findings indicate that NSD1-deposited H3K36 methylation regulates transcription by directing promoter DNA methylation, partially by repressing polycomb repressive complex 2 (PRC2) activity. These findings could explain the phenotypic similarity of SS to OGID disorders that are caused by mutations in PRC2 complex-encoding genes.

Despite all the advances in the production of transgenic mice, the production efficiency of these animal models is still low. Given that the expression of developmental genes has a critical role in growth and development of embryo, we determined the expression pattern of pluripotency, trophectoderm and imprinting genes in the Rag1 (recombination‐activating gene 1) knocked‐out blastocysts resulting from microinjection of CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR‐associated protein 9) constructs into the zygote cytoplasm of C57bl6 mice. Following microinjection, the embryos were cultured and the gene expression of developed blastocysts and natural blastocysts (Sham and control groups) were evaluated using real‐time PCR. The agarose gel to confirm the deletion in the Rag1 gene in Rag1 knocked‐out blastocyst. Our results showed that the expression of trophectoderm genes (‐TEAD‐4 and Cdx2), pluripotency genes (Nanog and Oct‐4) and imprinting gene (H19) in the Rag1 knocked‐out group was significantly lower compared with the embryos obtained from Natural fertilization. According to these findings, manipulation, embryo culture and microinjection of CRISPR constructs into the zygote cytoplasm of mice led to reduced expression of imprinting, pluripotency and trophectoderm genes. Therefore, the Rag1 knocked‐out embryos produced by the CRISPR/Cas9 system are of low quality, which reduces the chances of live birth in these animals and may cause various abnormalities in fetuses., Our results showed that the expression of trophectoderm genes, pluripotency genes and imprinting gene in the Rag1 knocked‐out group were significantly lower compared with the embryos obtained from Natural fertilization. Therefore, the Rag1 knocked‐out embryos produced by the CRISPR/Cas9 system are of low quality, which reduces the chances of live birth in these animals and may cause various abnormalities in fetuses

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated non-uniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

Background Non-genetic disease inheritance and offspring phenotype are substantially influenced by germline epigenetic programming, including genomic imprinting. Loss of Polycomb Repressive Complex 2 (PRC2) function in oocytes causes non-genetically inherited effects on offspring, including embryonic growth restriction followed by post-natal offspring overgrowth. While PRC2-dependent non-canonical imprinting is likely to contribute, less is known about germline epigenetic programming of non-imprinted genes during oocyte growth. In addition, de novo germline mutations in genes encoding PRC2 lead to overgrowth syndromes in human patients, but the extent to which PRC2 activity is conserved in human oocytes is poorly understood. Results In this study, we identify a discrete period of early oocyte growth during which PRC2 is expressed in mouse growing oocytes. Deletion of Eed during this window led to the de-repression of 343 genes. A high proportion of these were developmental regulators, and the vast majority were not imprinted genes. Many of the de-repressed genes were also marked by the PRC2-dependent epigenetic modification histone 3 lysine 27 trimethylation (H3K27me3) in primary–secondary mouse oocytes, at a time concurrent with PRC2 expression. In addition, we found H3K27me3 was also enriched on many of these genes by the germinal vesicle (GV) stage in human oocytes, strongly indicating that this PRC2 function is conserved in the human germline. However, while the 343 genes were de-repressed in mouse oocytes lacking EED, they were not de-repressed in pre-implantation embryos and lost H3K27me3 during pre-implantation development. This implies that H3K27me3 is a transient feature that represses a wide range of genes in oocytes. Conclusions Together, these data indicate that EED has spatially and temporally distinct functions in the female germline to repress a wide range of developmentally important genes and that this activity is conserved in the mouse and human germlines.

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.

BackgroundSeverity is an important factor for inclusion of diseases on expanded carrier screening (ECS) panels. Here, we applied a validated algorithm that objectively classifies diseases into severity categories to 176 genes on a clinically available ECS panel. We then mapped disease traits from the algorithm to severity-related criteria cited by the American College of Obstetricians and Gynecologists (ACOG).MethodsEight genetic counselors (GCs), followed by four medical geneticists (MDs), applied the algorithm to subsets of the 176 genes. MDs and GCs then determined which disease traits met ACOG severity criteria.ResultsUpon initial GC and MD review, 107/176 genes (61%) and 133/176 genes (76%), respectively, had concordant classifications, with consensus reached for all genes. Final severity classifications were 68 (39%) profound, 71 (40%) severe, 36 (20%) moderate, and one (1%) mild. The vast majority of genes (170 out of 176) met at least one of ACOG’s four severity criteria.ConclusionThis study classified the severity of a large set of Mendelian genes by collaborative clinical expert application of an algorithm. Further, it clarified and operationalized difficult to interpret ACOG severity criteria via mapping of disease traits, thereby promoting consistency of ACOG criteria interpretation across laboratories.What’s already known about this topic?Disease severity is an important consideration for disease inclusion on expanded carrier screening panels.An algorithm that objectively classifies diseases into severity categories has been published and validated.What does this study add?176 genes were classified into severity categories.The algorithm was used to bring clarity to American College of Obstetricians and Gynecologist’s (ACOG’s) severity criteria that are not easily interpretable.170 of 176 genes met at least one of ACOG’s severity criteria.Data Availability StatementThe data that support the findings of this study have been completely reported in this manuscript and shared in the Figures and Supplementary Material.

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.

Appropriate developmental gene regulation relies on the capacity of gene promoters to integrate inputs from distal regulatory elements, yet how this is achieved remains poorly understood. In embryonic stem cells (ESCs), a subset of silent developmental gene promoters are primed for activation by FBXL19, a CpG island binding protein, through its capacity to recruit CDK-Mediator. How mechanistically these proteins function together to prime genes for activation during differentiation is unknown. Here we discover that in mouse ESCs FBXL19 and CDK-Mediator support long-range interactions between silent gene promoters that rely on FBXL19 for their induction during differentiation and gene regulatory elements. During gene induction, these distal regulatory elements behave in an atypical manner, in that the majority do not acquire histone H3 lysine 27 acetylation and no longer interact with their target gene promoter following gene activation. Despite these atypical features, we demonstrate by targeted deletions that these distal elements are required for appropriate gene induction during differentiation. Together these discoveries demonstrate that CpG-island associated gene promoters can prime genes for activation by communicating with atypical distal gene regulatory elements to achieve appropriate gene expression.

Transposable elements (TEs) are the major sources of lineage-specific genomic innovation and comprise nearly half of the human genome, but most of their functions remain unclear. Here, we identify that a series of endogenous retroviruses (ERVs), a TE subclass, regulate the transcriptome at the definitive endoderm stage with in vitro differentiation model from human embryonic stem cell. Notably, these ERVs perform as enhancers containing binding sites for critical transcription factors for endoderm lineage specification. Genome-wide methylation analysis shows most of these ERVs are derepressed by TET1-mediated DNA demethylation. LTR6B, a representative definitive endoderm activating ERV, contains binding sites for FOXA2 and GATA4 and governs the primate-specific expression of its neighboring developmental genes such as ERBB4 in definitive endoderm. Together, our study proposes evidence that recently evolved ERVs represent potent de novo developmental regulatory elements, which, in turn, fine-tune species-specific transcriptomes during endoderm and embryonic development.

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.

FOXA2 has been known to play important roles in liver functions in rodents. However, its role in human hepatocytes is not fully understood. Recently, we generated FOXA2 mutant induced pluripotent stem cell (FOXA2−/−iPSC) lines and illustrated that loss of FOXA2 results in developmental defects in pancreatic islet cells. Here, we used FOXA2−/−iPSC lines to understand the role of FOXA2 on the development and function of human hepatocytes. Lack of FOXA2 resulted in significant alterations in the expression of key developmental and functional genes in hepatic progenitors (HP) and mature hepatocytes (MH) as well as an increase in the expression of ER stress markers. Functional assays demonstrated an increase in lipid accumulation, bile acid synthesis and glycerol production, while a decrease in glucose uptake, glycogen storage, and Albumin secretion. RNA-sequencing analysis further validated the findings by showing a significant increase in genes associated with lipid metabolism, bile acid secretion, and suggested the activation of hepatic stellate cells and hepatic fibrosis in MH lacking FOXA2. Overexpression of FOXA2 reversed the defective phenotypes and improved hepatocyte functionality in iPSC-derived hepatic cells lacking FOXA2. These results highlight a potential role of FOXA2 in regulating human hepatic development and function and provide a human hepatocyte model, which can be used to identify novel therapeutic targets for FOXA2-associated liver disorders.

Animal and plant microRNAs (miRNAs) are essential for the spatio-temporal regulation of development. Together with this role, plant miRNAs have been proposed to target transposable elements (TEs) and stimulate the production of epigenetically active small interfering RNAs. This activity is evident in the plant male gamete containing structure, the male gametophyte or pollen grain. How the dual role of plant miRNAs, regulating both genes and TEs, is integrated during pollen development and which mRNAs are regulated by miRNAs in this cell type at a genome-wide scale are unknown. Here, we provide a detailed analysis of miRNA dynamics and activity during pollen development in Arabidopsis thaliana using small RNA and degradome parallel analysis of RNA end high-throughput sequencing. Furthermore, we uncover miRNAs loaded into the two main active Argonaute (AGO) proteins in the uninuclear and mature pollen grain, AGO1 and AGO5. Our results indicate that the developmental progression from microspore to mature pollen grain is characterized by a transition from miRNAs targeting developmental genes to miRNAs regulating TE activity.

Spiders constitute more than 49,000 described species distributed all over the world, and all ecological environments. Their order, Araneae, is defined by a set of characteristics with no parallel among their arachnid counterparts (e.g., spinnerets, silk glands, chelicerae that inoculate venom, among others). Changes in developmental pathways often underlie the evolution of morphological synapomorphies, and as such spiders are a promising model to study the role of developmental genes in the origin of evolutionary novelties. With that in mind, we investigated changes in the evolutionary regime of a set of six developmental genes, using spiders as our model. The genes were mainly chosen for their roles in spinneret ontogeny, yet they are pleiotropic, and it is likely that the origins of other unique morphological phenotypes are also linked to changes in their sequences. Our results indicate no great differences in the selective pressures on those genes when comparing spiders to other arachnids, but a few site-specific positive selection evidence were found in the Araneae lineage. These findings lead us to new insights on spider evolution that are to be further tested.

Preterm neonates are at a high risk for nephron loss under adverse clinical conditions. Renal damage potentially collides with postnatal nephrogenesis. Recent animal studies suggest that nephron loss within this vulnerable phase leads to renal damage later in life. Nephrogenic pathways are commonly reactivated after kidney injury supporting renal regeneration. We hypothesized that nephron loss during nephrogenesis affects renal development, which, in turn, impairs tissue repair after secondary injury. Neonates prior to 36 wk of gestation show an active nephrogenesis. In rats, nephrogenesis is ongoing until day 10 after birth. Mimicking the situation of severe nephron loss during nephrogenesis, male pups were uninephrectomized at day 1 of life (UNXd1). A second group of males was uninephrectomized at postnatal day 14 (UNXd14), after terminated nephrogenesis. Age-matched controls were sham operated. Three days after uninephrectomy transcriptional changes in the right kidney were analyzed by RNA-sequencing, followed by functional pathway analysis. In UNXd1, 1,182 genes were differentially regulated, but only 143 genes showed a regulation both in UNXd1 and UNXd14. The functional groups “renal development” and “kidney injury” were among the most differentially regulated groups and revealed distinctive alterations. Reduced expression of candidate genes concerning renal development ( Bmp7, Gdnf, Pdgf-B, Wt1) and injury ( nephrin, podocin, Tgf-β1) were detected. The downregulation of Bmp7 and Gdnf persisted until day 28. In UNXd14, Six2 was upregulated and Pax2 was downregulated. We conclude that nephron loss during nephrogenesis affects renal development and induces a specific regulation of genes that might hinder tissue repair after secondary kidney injury.

Conditionally activated, caged morpholino antisense agents (cMOs) are tools that enable the temporal and spatial investigation of gene expression, regulation, and function during embryonic development. Cyclic MOs are conformationally gated oligonucleotide analogs that do not block gene expression until they are linearized through the application of an external trigger, such as light or enzyme activity. Here, we describe the first examples of small molecule-responsive cMOs, which undergo rapid and efficient decaging via a Staudinger reduction. This is enabled by a highly flexible linker design that offers opportunities for the installation of chemically activated, self-immolative motifs. We synthesized cyclic cMOs against two distinct, developmentally relevant genes and demonstrated phosphine-triggered knockdown of gene expression in zebrafish embryos. This represents the first report of a small molecule-triggered antisense agent for gene knockdown, adding another bioorthogonal entry to the growing arsenal of gene knockdown tools.

Sarcomas are difficult to treat and the therapy, even when effective, is associated with long-term and life-threatening side effects. In addition, the treatment regimens for many sarcomas, including Ewing sarcoma, rhabdomyosarcoma, and osteosarcoma, are relatively unchanged over the past two decades, indicating a critical lack of progress. Although differentiation-based therapies are used for the treatment of some cancers, the application of this approach to sarcomas has proven challenging. Here, using a CRISPR-mediated gene knockout approach, we show that Inhibitor of DNA Binding 2 (ID2) is a critical regulator of developmental-related genes and tumor growth in vitro and in vivo in Ewing sarcoma tumors. We also identified that homoharringtonine, which is an inhibitor of protein translation and FDA-approved for the treatment of leukemia, decreases the level of the ID2 protein and significantly reduces tumor growth and prolongs mouse survival in an Ewing sarcoma xenograft model. Furthermore, in addition to targeting ID2, homoharringtonine also reduces the protein levels of ID1 and ID3, which are additional members of the ID family of proteins with well-described roles in tumorigenesis, in multiple types of cancer. Overall, these results provide insight into developmental regulation in Ewing sarcoma tumors and identify a novel, therapeutic approach to target the ID family of proteins using an FDA-approved drug.

The pituitary gland, as a nodal component of the endocrine system, is responsible for the regulation of growth, reproduction, metabolism, and homeostasis. Although pituitary formation though the hierarchical action of different transcription factors is well studied in mouse models, there is little evidence of the analogous developmental processes in humans. Herein, we present a female patient with a phenotype that includes blepharoptosis-ptosis-epicanthus syndrome and premature ovarian failure. Clinical exome sequencing revealed two heterozygous variants in two genes, LHX4 (pathogenic) and NR5A1 (VUS) genes and no mutation in FOXL2 gene. We propose a model of genetic interaction between LHX4 and NR5A1 during pituitary and ovarian development that may lead to a similar phenotype mediated by reduced FOXL2 expression.

The histone acetyltransferase Gcn5 is critical for gene expression and development. In Drosophila, Gcn5 is part of four complexes (SAGA, ATAC, CHAT and ADA) that are essential for fly viability and have key roles in regulating gene expression. Here, we show that although the SAGA, ADA and CHAT complexes play redundant roles in embryonic gene expression, the insect-specific CHAT complex uniquely regulates expression of a subset of developmental genes. We also identify a substantial decrease in histone acetylation in chiffon mutant embryos that exceeds that observed in Ada2b, suggesting broader roles for Chiffon in regulating histone acetylation outside of the Gcn5 complexes. The chiffon gene encodes two independent polypeptides that nucleate formation of either the CHAT or Dbf4-dependent kinase (DDK) complexes. DDK includes the cell cycle kinase Cdc7, which is necessary for maternally driven DNA replication in the embryo. We identify a temporal switch between the expression of these chiffon gene products during a short window during the early nuclear cycles in embryos that correlates with the onset of zygotic genome activation, suggesting a potential role for CHAT in this process. This article has an associated First Person interview with the first author of the paper.

Higher plants progress through a juvenile and an adult phase of development before they enter the reproductive phase. The transition from the juvenile to the adult phase is referred to as vegetative phase change, and this is signified by the production of trichomes on the abaxial side of leaf blades in Arabidopsis; however, the molecular mechanism underlying this process is poorly understood. We identified a dominant mutation (gl1-D) in a forward genetic screen that accelerates abaxial trichome production during shoot development. This phenotype is the result of a G-to-A substitution in the 3' noncoding region in the GLABRA1(GL1) gene. We show that TOE1, an AP2-like transcription factor that acts downstream of the miR156-SPL pathway, represses GL1 expression by directly binding to this site, and that gl1-D prevents TOE1 binding. Our work reveals a molecular link between vegetative phase change and abaxial trichome production in Arabidopsis, and answers a long-standing question of how the vegetative phase change pathway regulates vegetative traits.

13 pages, 4 figures, supplementary data https://doi.org/10.1093/molbev/msz153, Domestication of wild animals induces a set of phenotypic characteristics collectively known as the domestication syndrome. However, how this syndrome emerges is still not clear. Recently, the neural crest cell deficit hypothesis proposed that it is generated by a mildly disrupted neural crest cell developmental program, but clear support is lacking due to the difficulties of distinguishing pure domestication effects from preexisting genetic differences between farmed and wild mammals and birds. Here, we use a farmed fish as model to investigate the role of persistent changes in DNA methylation (epimutations) in the process of domestication. We show that early domesticates of sea bass, with no genetic differences with wild counterparts, contain epimutations in tissues with different embryonic origins. About one fifth of epimutations that persist into adulthood are established by the time of gastrulation and affect genes involved in developmental processes that are expressed in embryonic structures, including the neural crest. Some of these genes are differentially expressed in sea bass with lower jaw malformations, a key feature of domestication syndrome. Interestingly, these epimutations significantly overlap with cytosine-to-thymine polymorphisms after 25 years of selective breeding. Furthermore, epimutated genes coincide with genes under positive selection in other domesticates. We argue that the initial stages of domestication include dynamic alterations in DNA methylation of developmental genes that affect the neural crest. Our results indicate a role for epimutations during the beginning of domestication that could be fixed as genetic variants and suggest a conserved molecular process to explain Darwin’s domestication syndrome across vertebrates, Research supported by the Spanish Ministry of Economy and Competitiveness (MINECO) grants “Epifarm” (ref. AGL2013-41047-R) and “Epimark” (ref. AGL2016-78710-R) to F.P. D.A. was supported by a PhD scholarship from the Spanish Government (BES-2011-044860) and an Epimark contract

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.

The colonization of an animal’s tissues by its microbial partners creates networks of communication across the host’s body. We used the natural binary light-organ symbiosis between the squid Euprymna scolopes and its luminous bacterial partner, Vibrio fischeri , to define the impact of colonization on transcriptomic networks in the host. A night-active predator, E. scolopes coordinates the bioluminescence of its symbiont with visual cues from the environment to camouflage against moon and starlight. Like mammals, this symbiosis has a complex developmental program and a strong day/night rhythm. We determined how symbiont colonization impacted gene expression in the light organ itself, as well as in two anatomically remote organs: the eye and gill. While the overall transcriptional signature of light organ and gill were more alike, the impact of symbiosis was most pronounced and similar in light organ and eye, both in juvenile and adult animals. Furthermore, the presence of a symbiosis drove daily rhythms of transcription within all three organs. Finally, a single mutation in V. fischeri —specifically, deletion of the lux operon, which abrogates symbiont luminescence—reduced the symbiosis-dependent transcriptome of the light organ by two-thirds. In addition, while the gills responded similarly to light-organ colonization by either the wild-type or mutant, luminescence was required for all of the colonization-associated transcriptional responses in the juvenile eye. This study defines not only the impact of symbiont colonization on the coordination of animal transcriptomes, but also provides insight into how such changes might impact the behavior and ecology of the host.

Maize often exhibits asynchronous pollination under abiotic and biotic stress conditions; however, the molecular basis of this developmental deficiency has not been elucidated. Tassel development is a key process affecting the anthesis-silking interval (ASI) in maize. In this study, we showed that pollen shedding was delayed and ASI was significantly increased in B73 and Chang7-2 inbred lines under water deficit conditions, which resulted in longer barren tip length and decreased yields under both controlled and field conditions. Comparative transcriptome analysis performed on immature tassels derived from plants grown under well-watered and water deficit conditions identified 1931 and 1713 differentially expressed genes (DEGs) in B73 and Chang7-2, respectively. Further, 28 differentially co-expressed transcription factors were identified across both lines. Collectively, we demonstrated that the molecular regulation of tassel development is associated with water deficit stress at early vegetative stage in maize. This finding extends our understanding of the molecular basis of maize tassel development during abiotic stress.

Heterozygous mutations in HNF1B cause the complex syndrome renal cysts and diabetes (RCAD), characterized by developmental abnormalities of the kidneys, genital tracts and pancreas, and a variety of renal, pancreas and liver dysfunctions. The pathogenesis underlying this syndrome remains unclear as mice with heterozygous null mutations have no phenotype, while constitutive/conditional Hnf1b ablation leads to more severe phenotypes. We generated a novel mouse model carrying an identified human mutation at the intron-2 splice donor site. Unlike heterozygous mice previously characterized, mice heterozygous for the splicing mutation exhibited decreased HNF1B protein levels and bilateral renal cysts from embryonic day 15, originated from glomeruli, early proximal tubules (PTs) and intermediate nephron segments, concurrently with delayed PT differentiation, hydronephrosis and rare genital tract anomalies. Consistently, mRNA sequencing showed that most downregulated genes in embryonic kidneys were primarily expressed in early PTs and the loop of Henle and involved in ion/drug transport, organic acid and lipid metabolic processes, while the expression of previously identified targets upon Hnf1b ablation, including cystic disease genes, was weakly or not affected. Postnatal analyses revealed renal abnormalities, ranging from glomerular cysts to hydronephrosis and, rarely, multicystic dysplasia. Urinary proteomics uncovered a particular profile predictive of progressive decline in kidney function and fibrosis, and displayed common features with a recently reported urine proteome in an RCAD pediatric cohort. Altogether, our results show that reduced HNF1B levels lead to developmental disease phenotypes associated with the deregulation of a subset of HNF1B targets. They further suggest that this model represents a unique clinical/pathological viable model of the RCAD disease., Summary: A novel established mouse model carrying a heterozygous splicing human mutation in the Hnf1b gene exhibits phenotypes similar to those of patients with renal cysts and diabetes disease.

BackgroundIt has been proposed that recently duplicated genes are more likely to be redundant with one another compared to ancient paralogues. The evolutionary logic underpinning this idea is simple, as the assumption is that recently derived paralogous genes are more similar in sequence compared to members of ancient gene families. We set out to test this idea by using molecular phylogenetics and exploiting the genetic tractability of the model nematode,Caenorhabditis elegans,in studying the nematode-specific family of Hedgehog-related genes, the Warthogs. Hedgehog is one of a handful of signal transduction pathways that underpins the development of bilaterian animals. While having lost abona fideHedgehog gene, most nematodes have evolved an expanded repertoire of Hedgehog-related genes, ten of which reside within the Warthog family.ResultsWe have characterised their evolutionary origin and their roles inC. elegansand found that these genes have adopted new functions in aspects of post-embryonic development, including left–right asymmetry and cell fate determination, akin to the functions of their vertebrate counterparts. Analysis of various double and triple mutants of the Warthog family reveals that more recently derived paralogues are not redundant with one another, while a pair of divergent Warthogs do display redundancy with respect to their function in cuticle biosynthesis.ConclusionsWe have shown that newer members of taxon-restricted gene families are not always functionally redundant despite their recent inception, whereas much older paralogues can be, which is considered paradoxical according to the current framework in gene evolution.

Vertebrate embryos achieve developmental competency during zygotic genome activation (ZGA) by establishing chromatin states that silence yet poise developmental genes for subsequent lineage-specific activation. Here, we reveal the order of chromatin states in establishing developmental gene poising in preZGA zebrafish embryos. Poising is established at promoters and enhancers that initially contain open/permissive chromatin with ‘Placeholder’ nucleosomes (bearing H2A.Z, H3K4me1, and H3K27ac), and DNA hypomethylation. Silencing is initiated by the recruitment of polycomb repressive complex 1 (PRC1), and H2Aub1 deposition by catalytic Rnf2 during preZGA and ZGA stages. During postZGA, H2Aub1 enables Aebp2-containing PRC2 recruitment and H3K27me3 deposition. Notably, preventing H2Aub1 (via Rnf2 inhibition) eliminates recruitment of Aebp2-PRC2 and H3K27me3, and elicits transcriptional upregulation of certain developmental genes during ZGA. However, upregulation is independent of H3K27me3 – establishing H2Aub1 as the critical silencing modification at ZGA. Taken together, we reveal the logic and mechanism for establishing poised/silent developmental genes in early vertebrate embryos.

Cell proliferation and differentiation require signalling pathways that enforce appropriate and timely gene expression. We find that Tor2, the catalytic subunit of the TORC1 complex in fission yeast, targets a conserved nuclear RNA elimination network, particularly the serine and proline-rich protein Pir1, to control gene expression through RNA decay and facultative heterochro-matin assembly. Phosphorylation by Tor2 protects Pir1 from degradation by the ubiquitin-proteasome system involving the polyubiquitin Ubi4 stress-response protein and the Cul4–Ddb1 E3 ligase. This pathway suppresses widespread and untimely gene expression and is critical for sustaining cell proliferation. Moreover, we find that the dynamic nature of Tor2-mediated control of RNA elimination machinery defines gene expression patterns that coordinate fundamental chromosomal events dur-ing gametogenesis, such as meiotic double-strand-break formation and chromosome segregation. These findings have impor-tant implications for understanding how the TOR signalling pathway reprogrammes gene expression patterns and contributes to diseases such as cancer.

Whole-genome duplications (WGDs) have occurred multiple times during animal evolution, including in lineages leading to vertebrates, teleosts, horseshoe crabs, and arachnopulmonates. These dramatic events initially produce a wealth of new genetic material, generally followed by extensive gene loss. It appears, however, that developmental genes such as homeobox genes, signaling pathway components and microRNAs are frequently retained as duplicates (so-called ohnologs) following WGD. These not only provide the best evidence for WGD, but an opportunity to study its evolutionary consequences. Although these genes are well studied in the context of vertebrate WGD, similar comparisons across the extant arachnopulmonate orders are patchy. We sequenced embryonic transcriptomes from two spider species and two amblypygid species and surveyed three important gene families, Hox, Wnt, and frizzled, across these and 12 existing transcriptomic and genomic resources for chelicerates. We report extensive retention of putative ohnologs, further supporting the ancestral arachnopulmonate WGD. We also found evidence of consistent evolutionary trajectories in Hox and Wnt gene repertoires across three of the six arachnopulmonate orders, with interorder variation in the retention of specific paralogs. We identified variation between major clades in spiders and are better able to reconstruct the chronology of gene duplications and losses in spiders, amblypygids, and scorpions. These insights shed light on the evolution of the developmental toolkit in arachnopulmonates, highlight the importance of the comparative approach within lineages, and provide substantial new transcriptomic data for future study.

Steps of mRNA maturation are important gene regulatory events that occur in distinct cellular locations. However, transcriptomic analyses often lose information on the subcellular distribution of processed and unprocessed transcripts. We generated extensive RNA-seq data sets to track mRNA maturation across subcellular locations in mouse embryonic stem cells, neuronal progenitor cells, and postmitotic neurons. We find disparate patterns of RNA enrichment between the cytoplasmic, nucleoplasmic, and chromatin fractions, with some genes maintaining more polyadenylated RNA in chromatin than in the cytoplasm. We bioinformatically defined four regulatory groups for intron retention, including complete cotranscriptional splicing, complete intron retention in the cytoplasmic RNA, and two intron groups present in nuclear and chromatin transcripts but fully excised in cytoplasm. We found that introns switch their regulatory group between cell types, including neuronally excised introns repressed by polypyrimidine track binding protein 1 (PTBP1). Transcripts for the neuronal gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1) are highly expressed in mESCs but are absent from the cytoplasm. Instead, incompletely spliced Gabbr1 RNA remains sequestered on chromatin, where it is bound by PTBP1, similar to certain long noncoding RNAs. Upon neuronal differentiation, Gabbr1 RNA becomes fully processed and exported for translation. Thus, splicing repression and chromatin anchoring of RNA combine to allow posttranscriptional regulation of Gabbr1 over development. For this and other genes, polyadenylated RNA abundance does not indicate functional gene expression. Our data sets provide a rich resource for analyzing many other aspects of mRNA maturation in subcellular locations and across development.

Androgens/androgen receptor (AR)-mediated signaling pathways are essential for prostate development, morphogenesis and regeneration. Specifically, stromal AR signaling has been shown to be essential for prostatic initiation. However, the molecular mechanisms underlying AR-initiated mesenchymal-epithelial interactions in prostate development remain unclear. Here, using a newly generated mouse model, we have directly addressed the fate and role of genetically marked AR-expressing cells during embryonic prostate development. Androgen signaling-initiated signaling pathways were identified in mesenchymal niche populations at single-cell transcriptomic resolution. The dynamic cell-signaling networks regulated by stromal AR were additionally characterized in relation to prostatic epithelial bud formation. Pseudotime analyses further revealed the differentiation trajectory and fate of AR-expressing cells in both prostatic mesenchymal and epithelial cell populations. Specifically, the cellular properties of Zeb1-expressing progenitors were assessed. Selective deletion of AR signaling in a subpopulation of mesenchymal rather than epithelial cells dysregulated the expression of the master regulators and significantly impaired prostatic bud formation. These data provide novel, high-resolution evidence demonstrating the important role of mesenchymal androgen signaling in the cellular niche controlling prostate early development by initiating dynamic mesenchyme-epithelia cell interactions.

Forest fires generate a large amount of carbon that remains resident on the site as dead and partially 'pyrolysed' (i.e. burnt) material that has long residency times and constitutes a significant pool in fire-prone ecosystems. In addition, fire-induced hydrophobic soil layers, caused by condensation of pyrolysed waxes and lipids, increase post-fire erosion and can lead to long-term productivity losses. A small set of pyrophilous fungi dominate post-fire soils and are likely to be involved with the degradation of all these compounds, yet almost nothing is currently known about what these fungi do or the metabolic processes they employ. In this study, we sequenced and analysed genomes from fungi isolated after Rim fire near Yosemite National Park in 2013 and showed the enrichment/expansion of CAZymes and families known to be involved in fruiting body initiation when compared to other basidiomycete fungi. We found gene families potentially involved in the degradation of the hydrophobic layer and pyrolysed organic matter, such as hydrophobic surface binding proteins, laccases (AA1_1), xylanases (GH10, GH11), fatty acid desaturases and tannases. Thus, pyrophilous fungi are important actors to restate the soil's functional capabilities.

The indiscriminate use of nitrogenous fertilizers continues unabated for commercial crop production, resulting in air and water pollution. The development of rice varieties with enhanced nitrogen use efficiency (NUE) will require a thorough understanding of the molecular basis of a plant&rsquo, s response to low nitrogen (N) availability. The global expression profiles of root tissues collected from low and high N treatments at different time points in two rice genotypes, Pokkali and Bengal, with contrasting responses to N stress and contrasting root architectures were examined. Overall, the number of differentially expressed genes (DEGs) in Pokkali (indica) was higher than in Bengal (japonica) during low N and early N recovery treatments. Most low N DEGs in both genotypes were downregulated whereas early N recovery DEGs were upregulated. Of these, 148 Pokkali-specific DEGs might contribute to Pokkali&rsquo, s advantage under N stress. These DEGs included transcription factors and transporters and were involved in stress responses, growth and development, regulation, and metabolism. Many DEGs are co-localized with quantitative trait loci (QTL) related to root growth and development, chlorate-resistance, and NUE. Our findings suggest that the superior growth performance of Pokkali under low N conditions could be due to the genetic differences in a diverse set of genes influencing N uptake through the regulation of root architecture.

Appropriate flowering and maturity time are important for soybean production. Four maturity genes E1, E2, E3 and E4 have been molecularly identified and found to play major roles in the control of flowering and maturity of soybean. Here, to further investigate the effect of different allele combinations of E1-E4, we performed Kompetitive Allele Specific PCR (KASP) assays based on single nucleotide polymorphisms (SNPs) at these four E loci, and genotyped E1-E4 genes across 308 Chinese cultivars with a wide range of maturity groups. In total, twenty-one allele combinations for E1-E4 genes were identified across these Chinese cultivars. Various combinations of mutations at four E loci gave rise to the diversity of flowering and maturity time, which were associated with the adaptation of soybean cultivars to diverse geographic regions and farming systems. In particular, the cultivars with mutations at all four E loci reached flowering and maturity very early, and adapted to high-latitude cold regions. The allele combinations e1-as/e2-ns/e3-tr/E4, E1/e2-ns/E3/E4 and E1/E2/E3/E4 played important roles in the Northeast China, Huang-Huai-Hai (HHH) Rivers Valley and South China regions, respectively. Notably, E1 and E2, especially E2, affected flowering and maturity time of soybean significantly. Our study will be beneficial for germplasm evaluation, cultivar improvement and regionalization of cultivation in soybean production.