Deletion of a conserved genomic region associated with adolescent idiopathic scoliosis leads to vertebral rotation in mice.

Adolescent idiopathic scoliosis (AIS) is the most common form of scoliosis, in which spinal curvature develops in adolescence, and 90% of patients are female. Scoliosis is a debilitating disease that often requires bracing or surgery in severe cases. AIS affects 2%-5.2% of the population; however, the biological origin of the disease remains poorly understood. In this study, we aimed to determine the function of a highly conserved genomic region previously linked to AIS using a mouse model generated by CRISPR-CAS9 gene editing to knockout this area of the genome to understand better its contribution to AIS, which we named AIS_CRMΔ. We also investigated the upstream factors that regulate the activity of this enhancer in vivo, whether the spatial expression of the LBX1 protein would change with the loss of AIS-CRM function, and whether any phenotype would arise after deletion of this region. We found a significant increase in mRNA expression in the developing neural tube at E10.5, and E12.5, for not only Lbx1 but also other neighboring genes. Adult knockout mice showed vertebral rotation and proprioceptive deficits, also observed in human AIS patients. In conclusion, our study sheds light on the elusive biological origins of AIS, by targeting and investigating a highly conserved genomic region linked to AIS in humans. These findings provide valuable insights into the function of the investigated region and contribute to our understanding of the underlying causes of this debilitating disease.

Genome-wide analysis of early vascular tunic repair and regeneration for Botrylloides digenesis reveals striking similarities to human wound healing.

The early stages of regeneration after injury are similar to those of wound healing. The ascidian Botrylloides diegensis can regenerate an entire adult from a small fragment of vascular tunic following the removal of all zooids in an injury-induced regeneration model. We investigated the molecular and cellular changes following injury to determine the differences between the healing process and the initiation of whole-body regeneration (WBR). We conducted transcriptome analysis at specific time points during regeneration and wound healing to identify differentially expressed genes (DEGs) and the unique biological processes associated with each state. Our findings revealed 296 DEGs at 10 h post-injury (hpi), with 71 highly expressed in healed tissue and 225 expressed during the WBR process. These DEGs were predicted to play roles in tissue reorganization, integrin signaling, extracellular matrix organization, and the innate immune system. Pathway analysis of the upregulated genes in the healed tunic indicated functional enrichment related to tissue repair, as has been observed in other species. Additionally, we examined the cell types in the tunic and ampullae in both tissue states using histology and in situ hybridization for six genes identified by transcriptome analysis. We observed strong mRNA expression in cells within the WBR tunic, and in small RNA-positive granules near the tunic edge. We hypothesized that many of these genes function in the compaction of the ampullae tunic, which is a pivotal process for WBR and dormancy in B. diegensis, and in an immune response. These findings establish surprising similarities between ascidian regeneration and human wound healing, emphasizing the potential for future investigations into human regenerative and repair mechanisms. This study provides valuable insights into the gene sets specifically activated during regeneration compared to wound healing, shedding light on the divergent activities of these processes.

Identification of novel microRNAs in the embryonic mouse brain using deep sequencing.

Many advances in small RNA-seq technology and bioinformatics pipelines have been made recently, permitting the discovery of novel miRNAs in the embryonic day 15.5 (E15.5) mouse brain. We aimed to improve miRNA discovery in this tissue to expand our knowledge of the regulatory networks that underpin normal neurodevelopment, find new candidates for neurodevelopmental disorder aetiology, and deepen our understanding of non-coding RNA evolution. A high-quality small RNA-seq dataset of 458 M reads was generated. An unbiased miRNA discovery pipeline identified fifty putative novel miRNAs, six of which were selected for further validation. A combination of conservation analysis and target functional prediction was used to determine the authenticity of novel miRNA candidates. These findings demonstrate that miRNAs remain to be discovered, particularly if they have the features of other small RNA species.

RNA sequencing and expression analysis reveal a role for Lhx9 in the haploinsufficient adult mouse ovary.

Understanding the molecular pathways that underpin ovarian development and function is vital for improving the research approaches to investigating fertility. Despite a significant improvement in our knowledge of molecular activity in the ovary, many questions remain unanswered in the quest to understand factors influencing fertility and ovarian pathologies such as cancer. Here, we present an investigation into the expression and function of the developmental transcription factor LIM Homeobox 9 (LHX9) in the adult mouse ovary. We have characterized Lhx9 expression in several cell types of the mature ovary across follicle stages. To evaluate possible LHX9 function in the adult ovary, we investigated ovarian anatomy and transcription in an Lhx9+/- knockout mouse model displaying subfertility. Despite a lack of gross anatomical differences between genotypes, RNA-sequencing found that 90 differentially expressed genes between Lhx9+/ - and Lhx9+/+ mice. Gene ontology analyses revealed a reduced expression of genes with major roles in ovarian steroidogenesis and an increased expression of genes associated with ovarian cancer. Analysis of the ovarian epithelium revealed Lhx9+/ - mice have a disorganized epithelial phenotype, corresponding to a significant increase in epithelial marker gene expression. These results provide an analysis of Lhx9 in the adult mouse ovary, suggesting a role in fertility and ovarian epithelial cancer.

Characterization of a novel Lbx1 mouse loss of function strain.

Adolescent Idiopathic Scoliosis (AIS) is the most common type of spine deformity affecting 2-3% of the population worldwide. The etiology of this disease is still poorly understood. Several GWAS studies have identified single nucleotide polymorphisms (SNPs) located near the gene LBX1 that is significantly correlated with AIS risk. LBX1 is a transcription factor with roles in myocyte precursor migration, cardiac neural crest specification, and neuronal fate determination in the neural tube. Here, we further investigated the role of LBX1 in the developing spinal cord of mouse embryos using a CRISPR-generated mouse model expressing a truncated version of LBX1 (Lbx1Δ). Homozygous mice died at birth, likely due to cardiac abnormalities. To further study the neural tube phenotype, we used RNA-sequencing to identify 410 genes differentially expressed between the neural tubes of E12.5 wildtype and Lbx1Δ/Δ embryos. Genes with increased expression in the deletion line were involved in neurogenesis and those with broad roles in embryonic development. Many of these genes have also been associated with scoliotic phenotypes. In comparison, genes with decreased expression were primarily involved in skeletal development. Subsequent skeletal and immunohistochemistry analysis further confirmed these results. This study aids in understanding the significance of links between LBX1 function and AIS susceptibility.

Histone deacetylase activity is required for Botrylloides leachii whole-body regeneration.

The colonial tunicate Botrylloides leachii is exceptional at regenerating from a piece of vascular tunic after loss of all adults from the colony. Previous transcriptome analyses indicate a brief period of healing before regeneration of a new adult (zooid) in as little as 8-10 days. However, there is little understanding of how the resulting changes to gene expression, required to drive regeneration, are initiated and how the overall process is regulated. Rapid changes to transcription often occur in response to chromatin changes, mediated by histone modifications such as histone acetylation. Here, we investigated a group of key epigenetic modifiers, histone deacetylases (HDAC), which are known to play an important role in many biological processes such as development, healing and regeneration. Through our transcriptome data, we identified and quantified the expression levels of HDAC and histone acetyltransferase enzymes during whole-body regeneration (WBR). To determine whether HDAC activity is required for WBR, we inhibited its action using valproic acid and trichostatin A. HDAC inhibition prevented the final morphological changes normally associated with WBR and resulted in aberrant gene expression. Botrylloides leachii genes including Slit2, TGF-ß, Piwi and Fzd4 all showed altered mRNA levels upon HDAC inhibition in comparison with the control samples. Additionally, atypical expression of Bl_Piwi was found in immunocytes upon HDAC inhibition. Together, these results show that HDAC function, specifically HDAC I/IIa class enzymes, are vital for B. leachii to undergo WBR successfully.

Mutations Preventing Regulated Exon Skipping in MET Cause Osteofibrous Dysplasia.

The periosteum contributes to bone repair and maintenance of cortical bone mass. In contrast to the understanding of bone development within the epiphyseal growth plate, factors that regulate periosteal osteogenesis have not been studied as intensively. Osteofibrous dysplasia (OFD) is a congenital disorder of osteogenesis and is typically sporadic and characterized by radiolucent lesions affecting the cortical bone immediately under the periosteum of the tibia and fibula. We identified germline mutations in MET, encoding a receptor tyrosine kinase, that segregate with an autosomal-dominant form of OFD in three families and a mutation in a fourth affected subject from a simplex family and with bilateral disease. Mutations identified in all families with dominant inheritance and in the one simplex subject with bilateral disease abolished the splice inclusion of exon 14 in MET transcripts, which resulted in a MET receptor (MET(Δ14)) lacking a cytoplasmic juxtamembrane domain. Splice exclusion of this domain occurs during normal embryonic development, and forced induction of this exon-exclusion event retarded osteoblastic differentiation in vitro and inhibited bone-matrix mineralization. In an additional subject with unilateral OFD, we identified a somatic MET mutation, also affecting exon 14, that substituted a tyrosine residue critical for MET receptor turnover and, as in the case of the MET(Δ14) mutations, had a stabilizing effect on the mature protein. Taken together, these data show that aberrant MET regulation via the juxtamembrane domain subverts core MET receptor functions that regulate osteogenesis within cortical diaphyseal bone.

Uncovering the pathways underlying whole body regeneration in a chordate model, Botrylloides leachi using de novo transcriptome analysis.

BACKGROUND: Regenerative capacity differs greatly between animals. In vertebrates regenerative abilities are highly limited and tissue or organ specific. However the closest related chordate to the vertebrate clade, Botrylloides leachi, can undergo whole body regeneration (WBR). Therefore, research on WBR in B. leachi has focused on pathways known to be important for regeneration in vertebrates. To obtain a comprehensive vision of this unique process we have carried out the first de novo transcriptome sequencing for multiple stages of WBR occurring in B. leachi. The identified changes in gene expression during B. leachi WBR offer novel insights into this remarkable ability to regenerate. RESULTS: The transcriptome of B. leachi tissue undergoing WBR were analysed using differential gene expression, gene ontology and pathway analyses. We observed up-regulation in the expression of genes involved in wound healing and known developmental pathways including WNT, TGF-ß and Notch, during the earliest stages of WBR. Later in WBR, the expression patterns in several pathways required for protein synthesis, biogenesis and the organisation of cellular components were up-regulated. CONCLUSIONS: While the genes expressed early on are characteristic of a necessary wound healing response to an otherwise lethal injury, the subsequent vast increase in protein synthesis conceivably sustains the reestablishment of the tissue complexity and body axis polarity within the regenerating zooid. We have, for the first time, provided a global overview of the genes and their corresponding pathways that are modulated during WBR in B. leachi.

RNA localization in the honeybee (Apis mellifera) oocyte reveals insights about the evolution of RNA localization mechanisms.

Subcellular localization of RNAs is a critical biological process for generation of cellular asymmetries for many cell types and a critical step in axis determination during the early development of animals. We have identified transcripts localized to the anterior and posterior of honeybee oocyte using laser capture microscopy and microarray analysis. Analysis of orthologous transcripts in Drosophila indicates that many do not show a conserved pattern of localization. By microinjecting fluorescently labeled honeybee transcripts into Drosophila egg chambers we show that these RNAs become localized in a similar manner to their localization in honeybee oocytes, indicating conservation of the localization machinery. Thus while the mechanisms for localizing RNA are conserved, the complement of localized RNAs are not. We propose that this complement of localized RNAs may change relatively rapidly through the loss or evolution of signal sequences detected by the conserved localization machinery, and show this has occurred in one transcript that is localized in a novel way in the honeybee. Our proposal, that the acquisition of novel RNA localization is relatively easy to evolve, has implications for the evolution of symmetry breaking mechanisms that trigger axis formation and development in animal embryos.

Maternal PRKAA1 and EDNRA genotypes are associated with birth weight, and PRKAA1 with uterine artery diameter and metabolic homeostasis at high altitude.

Low birth weight and intrauterine growth restriction (IUGR) increase the risk of mortality and morbidity during the perinatal period as well as in adulthood. Environmental and genetic factors contribute to IUGR, but the influence of maternal genetic variation on birth weight is largely unknown. We implemented a gene-by-environment study wherein we utilized the growth restrictive effects of high altitude. Multigenerational high-altitude residents (Andeans) are protected from altitude-associated IUGR compared with recent migrants (Europeans). Using a combined cohort of low- and high-altitude European and Andean women, we tested 63 single nucleotide polymorphisms (SNPs) from 16 natural selection-nominated candidate gene regions for associations with infant birth weight. We identified significant SNP associations with birth weight near coding regions for two genes involved in oxygen sensing and vascular control, PRKAA1 and EDNRA, respectively. Next, we identified a significant association for the PRKAA1 SNP with an intermediate phenotype, uterine artery diameter, which has been shown to be related to Andean protection from altitude-associated reductions in fetal growth. To explore potential functional relationships for the effect of maternal SNP genotype on birth weight, we evaluated the relationship between maternal PRKAA1 SNP genotype and gene expression patterns in general and, in particular, of key pathways involved in metabolic homeostasis that have been proposed to play a role in the pathophysiology of IUGR. Our observations suggest that maternal genetic variation within genes that regulate oxygen sensing, metabolic homeostasis, and vascular control influence fetal growth and birth weight outcomes and hence Andean adaptation to high altitude.

Diversity in insect axis formation: two orthodenticle genes and hunchback act in anterior patterning and influence dorsoventral organization in the honeybee (Apis mellifera).

Axis formation is a key step in development, but studies indicate that genes involved in insect axis formation are relatively fast evolving. Orthodenticle genes have conserved roles, often with hunchback, in maternal anterior patterning in several insect species. We show that two orthodenticle genes, otd1 and otd2, and hunchback act as maternal anterior patterning genes in the honeybee (Apis mellifera) but, unlike other insects, act to pattern the majority of the anteroposterior axis. These genes regulate the expression domains of anterior, central and posterior gap genes and may directly regulate the anterior gap gene giant. We show otd1 and hunchback also influence dorsoventral patterning by regulating zerknült (zen) as they do in Tribolium, but that zen does not regulate the expression of honeybee gap genes. This suggests that interactions between anteroposterior and dorsal-ventral patterning are ancestral in holometabolous insects. Honeybee axis formation, and the function of the conserved anterior patterning gene orthodenticle, displays unique characters that indicate that, even when conserved genes pattern the axis, their regulatory interactions differ within orders of insects, consistent with relatively fast evolution in axis formation pathways.

Little effect of gestation at 3,100 m on fetal fat accretion or the fetal circulation.

OBJECTIVE: While chronic hypoxia has been recognized as the principal causative factor for decreasing birth weight at high altitude, unknown is whether fetal fat accretion and vascular function are affected. METHODS: Colorado women with normal singleton pregnancies (18 Denver residents, 1,600 m; 24 Leadville residents, 3,100 m) were studied longitudinally from 20 to 36 weeks gestation. Fetal biometry was used to obtain axial images for assessing mid-upper arm and mid-thigh subcutaneous tissue mass (MUA and MUL SQ) and Doppler waveform analysis conducted to measure indices of vascular function in the fetal umbilical arteries (UmbA), umbilical vein (UmbV), middle cerebral artery (MCA), and ductus venosus (DV). SAS PROC MIXED was used to compare altitudes with P < 0.05 considered significant and trends present when 0.05 < P < 0.10. RESULTS: The 3,100 m vs. 1,600 m babies weighed less at birth. Third trimester fetal biometry, MUA SQ and MUL SQ were somewhat lower, but neither the biometry nor the SQ altitudinal differences attained statistical significance. Greater prepregnant maternal BMI tended to decrease MUA SQ (P = 0.07) and increase MUL SQ (P = 0.07). UmbA S/D ratios decreased and UmbV flow increased with advancing gestation (both P < 0.001). Altitude did not affect the UmbA or MCA systolic/diastolic ratios (S/D), MCA peak-systolic velocity, UmbV flow, or the DV systolic/atrial flow ratio. CONCLUSION: The hypoxia of residence at high compared to moderate altitude lowered birth weight but did not significantly alter MUA or mid-thigh fetal subcutaneous tissue mass or Doppler indices of vascular function.

Andean and Tibetan patterns of adaptation to high altitude.

OBJECTIVES: High-altitude hypoxia, or decreased oxygen levels caused by low barometric pressure, challenges the ability of humans to live and reproduce. Despite these challenges, human populations have lived on the Andean Altiplano and the Tibetan Plateau for millennia and exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. We and others have identified natural selection candidate genes and gene regions for these adaptations using dense genome scan data. One gene previously known to be important in cellular oxygen sensing, egl nine homolog 1 (EGLN1), shows evidence of positive selection in both Tibetans and Andeans. Interestingly, the pattern of variation for this gene differs between the two populations. Continued research among Tibetan populations has identified statistical associations between hemoglobin concentration and single nucleotide polymorphism (SNP) genotype at EGLN1 and a second gene, endothelial PAS domain protein 1 (EPAS1). METHODS: To measure for the effects of EGLN1 and EPAS1 altitude genotypes on hemoglobin concentration among Andean highlanders, we performed a multiple linear regression analysis of 10 candidate SNPs in or near these two genes. RESULTS: Our analysis did not identify significant associations between EPAS1 or EGLN1 SNP genotypes and hemoglobin concentration in Andeans. CONCLUSIONS: These results contribute to our understanding of the unique set of adaptations developed in different highland groups to the hypoxia of high altitude. Overall, the results provide key insights into the patterns of genetic adaptation to high altitude in Andean and Tibetan populations.

Identifying signatures of natural selection in Tibetan and Andean populations using dense genome scan data.

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude.

Sex-biased gene and microRNA expression in the developing mouse brain is associated with neurodevelopmental functions and neurological phenotypes.

BACKGROUND: Sex differences pose a challenge and an opportunity in biomedical research. Understanding how sex chromosomes and hormones affect disease-causing mechanisms will shed light on the mechanisms underlying predominantly idiopathic sex-biased neurodevelopmental disorders such as ADHD, schizophrenia, and autism. Gene expression is a crucial conduit for the influence of sex on developmental processes; therefore, this study focused on sex differences in gene expression and the regulation of gene expression. The increasing interest in microRNAs (miRNAs), small, non-coding RNAs, for their contribution to normal and pathological neurodevelopment prompted us to test how miRNA expression differs between the sexes in the developing brain. METHODS: High-throughput sequencing approaches were used to identify transcripts, including miRNAs, that showed significantly different expression between male and female brains on day 15.5 of development (E15.5). RESULTS: Robust sex differences were identified for some genes and miRNAs, confirming the influence of biological sex on RNA. Many miRNAs that exhibit the greatest differences between males and females have established roles in neurodevelopment, implying that sex-biased expression may drive sex differences in developmental processes. In addition to highlighting sex differences for individual miRNAs, gene ontology analysis suggested several broad categories in which sex-biased RNAs might act to establish sex differences in the embryonic mouse brain. Finally, mining publicly available SNP data indicated that some sex-biased miRNAs reside near the genomic regions associated with neurodevelopmental disorders. CONCLUSIONS: Together, these findings reinforce the importance of cataloguing sex differences in molecular biology research and highlight genes, miRNAs, and pathways of interest that may be important for sexual differentiation in the mouse and possibly the human brain.

In biomedical research, understanding the differences between males and females is essential for understanding diseases that affect one sex more than the other. This study focused on gene expression and regulation differences between male and female mouse brains during development. We found that many microRNAs, small molecules that play a role in development were expressed differently between male and female brains. These differences could be important in understanding why males and females develop differently, particularly regarding neurodevelopmental disorders like ADHD, schizophrenia, and autism. We also found that some microRNAs that differed between males and females were located near genes associated with these disorders. Overall, the study highlights the importance of understanding sex differences in molecular biology research and provides insights into potential genes and pathways of interest for further study.

In vitro effects of wool-derived keratin on human dental pulp-derived stem cells for endodontic applications.

In this study we examine the influence of wool-derived keratin intermediate filament proteins (kIFPs) on human dental pulp-derived stem cells (hDPSCs). kIFPs were diluted (10 mg/mL to 0.001 mg/mL) in cell culture media. Effects on hDPSCs proliferation were measured using Alamar blue assay. Keratin concentrations of 1 mg/mL and 0.1 mg/mL were tested for odontogenic differentiation and mineralisation. Alkaline phosphatase (ALP) quantification (7th, 14th, and 21st days), alizarin red S (AR-S) staining and calcium quantification (21st day), reverse transcription polymerase chain reaction (RT-PCR, collagen expression), and immunocytochemical staining for dentin matrix protein (DMP) were performed. hDPSCs showed higher proliferation with kIFPs of 0.1 mg/mL or less (p < 0.0001). The 0.1 mg/mL keratin concentration promoted odontogenic differentiation, confirmed by increased ALP activity, significant calcium deposits (AR-S staining, p < 0.05), up-regulated collagen expression (RT-PCR, p < 0.05), and positive DMP staining. These results suggest that kIFPs could be a potential biomaterial for pulp-dentin regeneration.

The CNS-penetrating taxane drug TPI 287 potentiates antiglioma activity of the AURKA inhibitor alisertib in vivo.

INTRODUCTION: Glioblastoma (GBM) has a very poor prognosis despite current treatment. We previously found cytotoxic synergy between the AURKA inhibitor alisertib and the CNS-penetrating taxane TPI 287 against GBM tumor cells in vitro. METHODS: We used an orthotopic human GBM xenograft mouse model to test if TPI 287 potentiates alisertib in vivo. Western blotting, immunohistochemistry, siRNA knockdown, annexin V binding, and 3-dimensional Matrigel invasion assays were used to investigate potential mechanisms of alisertib and TPI 287 treatment interactions. RESULTS: Alisertib + TPI 287 combination therapy significantly prolonged animal survival compared to vehicle (p = 0.011), but only marginally compared to alisertib alone. Alisertib, TPI 287, and combined alisertib + TPI 287 reduced animal tumor volume compared to vehicle-treated controls. This was statistically significant for the combination therapy at 4 weeks (p < 0.0001). Alisertib + TPI 287 treatment decreased anti-apoptotic Bcl-2 protein levels in vivo and in vitro. Expression of the pro-apoptotic protein Bak was significantly increased by combination treatment (p < 0.0001). Pro-apoptotic Bim and Bak knockdown by siRNA decreased apoptosis by alisertib + TPI 287 in GB9, GB30, and U87 cells (p = 0.0005 to 0.0381). Although alisertib and TPI 287 significantly reduced GBM cell invasion (p < 0.0001), their combination was no more effective than TPI 287 alone. CONCLUSIONS: Results suggest that apoptosis is the dominant mechanism of potentiation of GBM growth inhibition by alisertib + TPI 287, in part through effects on Bcl-2 family proteins, providing a rationale for further laboratory testing of an AURKA inhibitor plus TPI 287 as a potential therapy against GBM.

Deep sequencing and expression of microRNAs from early honeybee (Apis mellifera) embryos reveals a role in regulating early embryonic patterning.

BACKGROUND: Recent evidence supports the proposal that the observed diversity of animal body plans has been produced through alterations to the complexity of the regulatory genome rather than increases in the protein-coding content of a genome. One significant form of gene regulation is the contribution made by the non-coding content of the genome. Non-coding RNAs play roles in embryonic development of animals and these functions might be expected to evolve rapidly. Using next-generation sequencing and in situ hybridization, we have examined the miRNA content of early honeybee embryos. RESULTS: Through small RNA sequencing we found that 28% of known miRNAs are expressed in the early embryo. We also identified developmentally expressed microRNAs that are unique to the Apoidea clade. Examination of expression patterns implied these miRNAs have roles in patterning the anterior-posterior and dorso-ventral axes as well as the extraembryonic membranes. Knockdown of Dicer, a key component of miRNA processing, confirmed that miRNAs are likely to have a role in patterning these tissues. CONCLUSIONS: Examination of the expression patterns of novel miRNAs, some unique to the Apis group, indicated that they are likely to play a role in early honeybee development. Known miRNAs that are deeply conserved in animal phyla display differences in expression pattern between honeybee and Drosophila, particularly at early stages of development. This may indicate miRNAs play a rapidly evolving role in regulating developmental pathways, most likely through changes to the way their expression is regulated.