Honey bees provide essential pollination services to the terrestrial ecosystem and produce important agricultural products. As a beneficial lactic acid bacterium, Enterococcus faecium is often supplied as a probiotic for honey bees and other animals. However, the underlying mechanisms of its actions and possible safety risks are not well understood. We present the first complete genome sequence of E. faecium isolated from the honey bee gut using nanopore sequencing, and investigate the effects and mechanisms of interactions between E. faecium and honey bees via transcriptome and miRNA analysis. E. faecium colonization increased honey bee gut weight. Transcriptome analysis showed that developmental genes were up-regulated. In accordance, the target genes of the down-regulated miRNAs were enriched in developmental pathways. We describe how E. faecium increases honey bee gut weight at the transcriptional and post-transcriptional levels, and add insights about how miRNAs mediate host and bacteria interactions.

Cestoda is a class of endoparasitic worms in the flatworm phylum (Platyhelminthes). During the course of their evolution cestodes have evolved some interesting aspects, such as their increased reproductive capacity. In this sense, they have serial repetition of their reproductive organs in the adult stage, which is often associated with external segmentation in a developmental process called strobilation. However, the molecular basis of strobilation is poorly understood. To assess this issue, an evolutionary comparative study among strobilated and non-strobilated flatworm species was conducted to identify genes and proteins related to the strobilation process.

The genes that are required for organismal survival are annotated as 'essential genes'. Identifying all the essential genes of an animal species can reveal critical functions that are needed during the development of the organism. To inform studies on mouse development, we developed a supervised machine learning classifier based on phenotype data from mouse knockout experiments. We used this classifier to predict the essentiality of mouse genes lacking experimental data. Validation of our predictions against a blind test set of recent mouse knockout experimental data indicated a high level of accuracy (>80%). We also validated our predictions for other mouse mutagenesis methodologies, demonstrating that the predictions are accurate for lethal phenotypes isolated in random chemical mutagenesis screens and embryonic stem cell screens. The biological functions that are enriched in essential and non-essential genes have been identified, showing that essential genes tend to encode intracellular proteins that interact with nucleic acids. The genome distribution of predicted essential and non-essential genes was analysed, demonstrating that the density of essential genes varies throughout the genome. A comparison with human essential and non-essential genes was performed, revealing conservation between human and mouse gene essentiality status. Our genome-wide predictions of mouse essential genes will be of value for the planning of mouse knockout experiments and phenotyping assays, for understanding the functional processes required during mouse development, and for the prioritisation of disease candidate genes identified in human genome and exome sequence datasets.

Autism is a complex disease involving both environmental and genetic factors. Recent efforts have implicated the correlation of genomic imprinting and brain development in autism, however the pathogenesis of autism is not completely clear. Here, we used bioinformatic tools to provide a comprehensive analysis of the autism-related genes, genomic imprinted genes and the spatially and temporally differentially expressed genes of human brain, aiming to explore the relationship between autism, brain development and genomic imprinting.

Fos-like antigens (Fosl) including Fosl1 and Fosl2 exclusively heterodimerize with Jun members to form AP-1 complex, thereby participating in various cellular progresses including cell cycle regulation. However, expression patterns of these two genes during embryonic development remains largely unknown. In the present study, both temporal and spatial expression patterns of fosl1 and fosl2 were examined during embryonic development of Xenopus tropicalis. Real-time quantitative PCR results showed that the expression of the two genes was increased from stage 2 to stage 42. However, expression level of fosl1 is much higher than that of fosl2 at stage 42. Whole-mount in situ hybridization showed that fosl1 was expressed in eyes, branchial arch, notochord, otic vesicle, and liver. However, fosl2 was expressed in lung primordium from stage 34 to stage 38, in addition to the moderate expression in eyes and branchial arch at stage 42. Thus, the developmental expression patterns of these two fosl genes is different in Xenopus embryos. These results provide a basis for further functional study of these two genes.

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.

In sole aquaculture production, consistency in the quality of produced eggs throughout the year is unpredictable. Hox genes have a crucial role in controlling embryonic development and their genetic variation could alter the phenotype dramatically. In teleosts genome duplication led paralog hox genes to become diverged. Direct association of polymorphism in hoxa1a, hoxa2a & hoxa2b of Solea solea with egg viability indicates hoxa2b as a potential genetic marker. High Resolution Melt (HRM) analysis was carried out in 52 viable and 61 non-viable eggs collected at 54±6 hours post fertilization (hpf). Allelic and genotypic frequencies of polymorphism were analyzed and results illustrated a significantly increased risk for non-viability for minor alleles and their homozygous genotypes. Haplotype analysis demonstrated a significant recessive effect on the risk of non-viability, by increasing the odds of disrupting embryonic development up to three-fold. Phylogenetic analysis showed that the paralog genes hoxa2a and hoxa2b, are separated distinctly in two clades and presented a significant ω variation, revealing their diverged evolutionary rate.

The study of mutants to elucidate gene functions has a long and successful history; however, to discover causative mutations in mutants that were generated by random mutagenesis often takes years of laboratory work and requires previously generated genetic and/or physical markers, or resources like DNA libraries for complementation. Here, we present an alternative method to identify defective genes in developmental mutants of the filamentous fungus Sordaria macrospora through Illumina/Solexa whole-genome sequencing. We sequenced pooled DNA from progeny of crosses of three mutants and the wild type and were able to pinpoint the causative mutations in the mutant strains through bioinformatics analysis. One mutant is a spore color mutant, and the mutated gene encodes a melanin biosynthesis enzyme. The causative mutation is a G to A change in the first base of an intron, leading to a splice defect. The second mutant carries an allelic mutation in the pro41 gene encoding a protein essential for sexual development. In the mutant, we detected a complex pattern of deletion/rearrangements at the pro41 locus. In the third mutant, a point mutation in the stop codon of a transcription factor-encoding gene leads to the production of immature fruiting bodies. For all mutants, transformation with a wild type-copy of the affected gene restored the wild-type phenotype. Our data demonstrate that whole-genome sequencing of mutant strains is a rapid method to identify developmental genes in an organism that can be genetically crossed and where a reference genome sequence is available, even without prior mapping information.

Among developmental control genes, transcription factor-target gene "linkages"--the direct connections between target genes and the factors that control their patterns of expression--can show remarkable evolutionary stability. However, the specific binding sites that mediate and define these regulatory connections are themselves often subject to rapid turnover. Here we describe several instances in which particular transcription factor binding motif combinations have evidently been conserved upstream of orthologous target genes for extraordinarily long evolutionary periods. This occurs against a backdrop in which other binding sites for the same factors are coming and going rapidly. Our examples include a particular Dpp Silencer Element upstream of insect brinker genes, in combination with a novel motif we refer to as the Downstream Element; combinations of a Suppressor of Hairless Paired Site (SPS) and a specific proneural protein binding site associated with arthropod Notch pathway target genes; and a three-motif combination, also including an SPS, upstream of deuterostome Hes repressor genes, which are also Notch targets. We propose that these stable motif architectures have been conserved intact from a deep ancestor, in part because they mediate a special mode of regulation that cannot be supplied by the other, unstable motif instances.

Genome sequencing projects have presented the opportunity for analysis of developmental genes in three vector mosquito species: Aedes aegypti, Culex quinquefasciatus, and Anopheles gambiae. A comparative genomic analysis of developmental genes in Drosophila melanogaster and these three important vectors of human disease was performed in this investigation. While the study was comprehensive, special emphasis centered on genes that 1) are components of developmental signaling pathways, 2) regulate fundamental developmental processes, 3) are critical for the development of tissues of vector importance, 4) function in developmental processes known to have diverged within insects, and 5) encode microRNAs (miRNAs) that regulate developmental transcripts in Drosophila. While most fruit fly developmental genes are conserved in the three vector mosquito species, several genes known to be critical for Drosophila development were not identified in one or more mosquito genomes. In other cases, mosquito lineage-specific gene gains with respect to D. melanogaster were noted. Sequence analyses also revealed that numerous repetitive sequences are a common structural feature of Drosophila and mosquito developmental genes. Finally, analysis of predicted miRNA binding sites in fruit fly and mosquito developmental genes suggests that the repertoire of developmental genes targeted by miRNAs is species-specific. The results of this study provide insight into the evolution of developmental genes and processes in dipterans and other arthropods, serve as a resource for those pursuing analysis of mosquito development, and will promote the design and refinement of functional analysis experiments.

In this study, we developed a correlative approach that combined DNA immunoprecipitation-seq and RNA-seq analyses to define the regulon of the Chlamydia trachomatis transcription factor Euo. We confirmed the proposed role of Euo as a transcriptional repressor of late chlamydial genes but also showed that Euo activates transcription of a subset of midcycle genes and autoregulates its own expression via negative feedback. This study validates and expands the role of Euo as an important developmental regulator in C. trachomatis. In addition, this genome-wide correlative approach can be applied to study transcription factors in other pathogenic bacteria.

We examined the allelic expression and positioning of two pluripotency-associated genes, OCT4 and SOX2, and two housekeeping genes, ACTB and TUBA, in 4- and 8-cell porcine embryos utilizing RNA and DNA fluorescence in situ hybridization (FISH) in single blastomeres. The proportion of blastomeres expressing SOX2 bi-allelically increased from 45% at the 4-cell stage to 60% at the 8-cell stage. Moreover, in 8-cell embryos, SOX2 was expressed bi-allelically in significantly more blastomeres than was the case for OCT4, and this was associated with a tendency for SOX2 alleles to move toward the nuclear interior during 4- to 8-cell transition. However, the radial location of OCT4 alleles did not change significantly during this transition. The locations of active and inactive alleles based on DNA and RNA FISH signals were also calculated. Inactive OCT4 alleles were located in very close proximity to the nuclear membrane, whereas active OCT4 alleles were more centrally disposed in the nucleus. Nevertheless, the nuclear location of active and inactive SOX2 alleles did not change in either 4- or 8-cell blastomeres. Our RNA and DNA FISH data provide novel information on the allelic expression patterns and positioning of pluripotency-associated genes, OCT4 and SOX2, during embryonic genome activation in pigs.

Osteoarthritis (OA), also known as degenerative arthritis, affects millions of people all over the world. OA occurs when the cartilage wears down over time, which is a worldwide complaint. The aim of this study was to screen and verify hub genes involved in developmental chondrogenesis as well as to explore potential molecular mechanisms.The expression profiles of GSE51812 were downloaded from the Gene Expression Omnibus (GEO) database, which contained 9 samples, including 6-week pre-chondrocytes (PC, 6 independent specimens) and 17-week fetal periarticular resting chondrocytes (RC, 3 independent specimens). The raw data were integrated to obtain differentially expressed genes (DEGs) and were further analyzed with bioinformatics analysis. The Gene Ontology (GO) and pathway enrichment of DEGs were conducted via Database for Annotation, Visualization, and Integrated Discovery (DAVID). The protein-protein interaction (PPI) networks of the DEGs were constructed based on data from the search tool for the retrieval of interacting genes (STRING) database. An intersection figure was provided to show the relationship between the DEGs identified in this study and genes from any existed related studies.A total of 9486 DEGs, including 4821 upregulated genes and 4665 downregulated genes were observed. The top 30 developmental chondrogenesis associated genes were identified, including matrix metalloproteinase (MMP)1, MMP3, MMP13, prostaglandin-endoperoxide synthase 2 (PTGS2), and so on. The majority of DEGs, including PTGS2, CCL20, CHI3L1, LIF, CXCL8, and CXCL12 were intensively enriched in immune-associated biological process terms, including inflammatory, and immune responses. Additionally, the majority of DEGs were mainly enriched in NF-kappa β (NF-kβ) signaling pathway and tumor necrosis factor (TNF) signaling pathway. The hub genes identified in STRING and Cytoscape databases included MMP1, MMP3, MMP13, PTGS2 and so on. Among the top 30 upregulated and downregulated DEGs, there were 15 genes have been reported to be associated with OA or developmental chondrogenesis.This large scale gene expression study observed genes associated with human developmental chondrogenesis and their relative GO function, which may offer opportunities for the research for cartilage tissue engineering and novel insights into the prevention of OA in the near future.

Progression of development has to be insulated from the damaging impacts of environmental and genetic perturbations to produce highly predictable phenotypes. Molecular chaperones, such as the heat shock proteins (HSPs), are known to buffer various environmental stresses, and are deeply involved in protein homeostasis. These characteristics of HSPs imply that they might affect developmental buffering and canalization.

Mollusks display a striking morphological disparity, including, among others, worm-like animals (the aplacophorans), snails and slugs, bivalves, and cephalopods. This phenotypic diversity renders them ideal for studies into animal evolution. Despite being one of the most species-rich phyla, molecular and in silico studies concerning specific key developmental gene families are still scarce, thus hampering deeper insights into the molecular machinery that governs the development and evolution of the various molluscan class-level taxa.

The draft nuclear genome sequence of the snail-transmitted, dimorphic, parasitic, platyhelminth Schistosoma mansoni revealed eight genes encoding proteins that contain the Universal Stress Protein (USP) domain. Schistosoma mansoni is a causative agent of human schistosomiasis, a severe and debilitating Neglected Tropical Disease (NTD) of poverty, which is endemic in at least 76 countries. The availability of the genome sequences of Schistosoma species presents opportunities for bioinformatics and genomics analyses of associated gene families that could be targets for understanding schistosomiasis ecology, intervention, prevention and control. Proteins with the USP domain are known to provide bacteria, archaea, fungi, protists and plants with the ability to respond to diverse environmental stresses. In this research investigation, the functional annotations of the USP genes and predicted nucleotide and protein sequences were initially verified. Subsequently, sequence clusters and distinctive features of the sequences were determined. A total of twelve ligand binding sites were predicted based on alignment to the ATP-binding universal stress protein from Methanocaldococcus jannaschii. In addition, six USP sequences showed the presence of ATP-binding motif residues indicating that they may be regulated by ATP. Public domain gene expression data and RT-PCR assays confirmed that all the S. mansoni USP genes were transcribed in at least one of the developmental life cycle stages of the helminth. Six of these genes were up-regulated in the miracidium, a free-swimming stage that is critical for transmission to the snail intermediate host. It is possible that during the intra-snail stages, S. mansoni gene transcripts for universal stress proteins are low abundant and are induced to perform specialized functions triggered by environmental stressors such as oxidative stress due to hydrogen peroxide that is present in the snail hemocytes. This report serves to catalyze the formation of a network of researchers to understand the function and regulation of the universal stress proteins encoded in genomes of schistosomes and their snail intermediate hosts.

The RNA interference technique is a powerful tool to understand gene function. Intriguingly, RNA interference cannot only be used for cells in vitro, but also in living organisms. Here, we have adapted the method for use in the chick embryo. However, this technique is limited by the uncertainty in predicting the RNAi transfection efficiency and site in the embryo. Hence, we elaborated a modified vector system, pEGFP-shRNA, which can coexpress enhanced green fluorescent protein (EGFP) and short hairpin RNA (shRNA) simultaneously to facilitate analysis of gene silencing in chicken embryos. We tested the silencing of two highly conserved genes (cAxin2, cParaxis), which play crucial roles in chicken embryonic developmental processes. For each target gene, four to five small DNA inserts, each of them encoding one shRNA, were selected and cloned individually to the vector downstream of the Pol III promoter (either human H1 or U6 promoter), which shared with highly conserved motifs in human and chicken. The pEGFP-shRNA constructs were electroporated into the neural tube or somites. After subsequent re-incubation of 24 h, the EGFP expression, with green fluorescent signal, indicated the transfected regions in the neural tube or somites. The EGFP expressing embryos were further submitted into the process of in situ hybridization for examination of the silencing effects. The results show that the EGFP signal in transfected areas correlated with the silencing of the target genes (cAxin2, cParaxis). The cAxin2 expression was inhibited by shRNAs of either targeting the RGS domain or the DAX domain coding region. The cParaxis mRNA level in transgenic somites and the related migratory myogenic population was also reduced. The results suggest that our novel dual expression EGFP-shRNA system opens a new possibility to study gene function in a convenient and efficient way.

A challenge in clinical genomics is to predict whether copy number variation (CNV) affecting a gene or multiple genes will manifest as disease. Increasing recognition of gene dosage effects in neurodevelopmental disorders prompted us to develop a computational approach based on critical-exon (highly expressed in brain, highly conserved) examination for potential etiologic effects. Using a large CNV dataset, our updated analyses revealed significant (P < 1.64 × 10(-15)) enrichment of critical-exons within rare CNVs in cases compared to controls. Separately, we used a weighted gene co-expression network analysis (WGCNA) to construct an unbiased protein module from prenatal and adult tissues and found it significantly enriched for critical exons in prenatal (P < 1.15 × 10(-50), OR = 2.11) and adult (P < 6.03 × 10(-18), OR = 1.55) tissues. WGCNA yielded 1,206 proteins for which we prioritized the corresponding genes as likely to have a role in neurodevelopmental disorders. We compared the gene lists obtained from critical-exon and WGCNA analysis and found 438 candidate genes associated with CNVs annotated as pathogenic, or as variants of uncertain significance (VOUS), from among 10,619 developmental delay cases. We identified genes containing CNVs previously considered to be VOUS to be new candidate genes for neurodevelopmental disorders (GIT1, MVB12B and PPP1R9A) demonstrating the utility of this strategy to index the clinical effects of CNVs.

Polycomb repression in mouse embryonic stem cells (ESCs) is tightly associated with promoter co-occupancy of RNA polymerase II (RNAPII) which is thought to prime genes for activation during early development. However, it is unknown whether RNAPII poising is a general feature of Polycomb repression, or is lost during differentiation. Here, we map the genome-wide occupancy of RNAPII and Polycomb from pluripotent ESCs to non-dividing functional dopaminergic neurons. We find that poised RNAPII complexes are ubiquitously present at Polycomb-repressed genes at all stages of neuronal differentiation. We observe both loss and acquisition of RNAPII and Polycomb at specific groups of genes reflecting their silencing or activation. Strikingly, RNAPII remains poised at transcription factor genes which are silenced in neurons through Polycomb repression, and have major roles in specifying other, non-neuronal lineages. We conclude that RNAPII poising is intrinsically associated with Polycomb repression throughout differentiation. Our work suggests that the tight interplay between RNAPII poising and Polycomb repression not only instructs promoter state transitions, but also may enable promoter plasticity in differentiated cells.

Almond is one of the most featured nut crops owing to its high nutritional value. However, due to three different waves of flower and fruitlet drop, fruit drop is a major concern for growers. In this study, we carried out a time-course transcriptome analysis to investigate gene expression differences between normal and abnormal fruitlet development. By de novo assembly analysis, we identified 33,577 unigenes and provided their functional annotations. In total, we identified 7,469 differentially expressed genes and observed the most apparent difference between normal and abnormal fruits at 12 and 17 days after flowering. Their biological functions were enriched in carbon metabolism, carbon fixation in photosynthetic organisms and plant hormone signal transduction. RT-qPCR validated the expression pattern of 14 representative genes, including glycosyltransferase like family 2, MYB39, IAA13, gibberellin-regulated protein 11-like and POD44, which confirmed the reliability of our transcriptome data. This study provides an insight into the association between abnormal fruit development and carbohydrate signaling from the early developmental stages and could be served as useful information for understanding the regulatory mechanisms related to almond fruit drop.