Cloning of HOXD1 from unfertilised human oocytes and expression analyses during murine oogenesis and embryogenesis

Diversity of human and mouse homeobox gene expression in development and adult tissues

Background

Homeobox genes encode a diverse set of transcription factors implicated in a vast range of biological processes including, but not limited to, embryonic cell fate specification and patterning. Although numerous studies report expression of particular sets of homeobox genes, a systematic analysis of the tissue specificity of homeobox genes is lacking.

Results

Here we analyse publicly-available transcriptome data from human and mouse developmental stages, and adult human tissues, to identify groups of homeobox genes with similar expression patterns. We calculate expression profiles for 242 human and 278 mouse homeobox loci across a combination of 59 human and 12 mouse adult tissues, early and late developmental stages. This revealed 20 human homeobox genes with widespread expression, primarily from the TALE, CERS and ZF classes. Most homeobox genes, however, have greater tissue-specificity, allowing us to compile homeobox gene expression lists for neural tissues, immune tissues, reproductive and developmental samples, and for numerous organ systems. In mouse development, we propose four distinct phases of homeobox gene expression: oocyte to zygote; 2-cell; 4-cell to blastocyst; early to mid post-implantation. The final phase change is marked by expression of ANTP class genes. We also use these data to compare expression specificity between evolutionarily-based gene classes, revealing that ANTP, PRD, LIM and POU homeobox gene classes have highest tissue specificity while HNF, TALE, CUT and CERS are most widely expressed.

Conclusions

The homeobox genes comprise a large superclass and their expression patterns are correspondingly diverse, although in a broad sense related to an evolutionarily-based classification. The ubiquitous expression of some genes suggests roles in general cellular processes; in contrast, most human homeobox genes have greater tissue specificity and we compile useful homeobox datasets for particular tissues, organs and developmental stages. The identification of a set of eutherian-specific homeobox genes peaking from human 8-cell to morula stages suggests co-option of new genes to new developmental roles in evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-016-0140-y) contains supplementary material, which is available to authorized users.

Transcript analysis reveals a specific HOX signature associated with positional identity of human endothelial cells

The endothelial cell has a remarkable ability for sub-specialisation, adapted to the needs of a variety of vascular beds. The role of developmental programming versus the tissue contextual environment for this specialization is not well understood. Here we describe a hierarchy of expression of HOX genes associated with endothelial cell origin and location. In initial microarray studies, differential gene expression was examined in two endothelial cell lines: blood derived outgrowth endothelial cells (BOECs) and pulmonary artery endothelial cells. This suggested shared and differential patterns of HOX gene expression between the two endothelial lines. For example, this included a cluster on chromosome 2 of HOXD1, HOXD3, HOXD4, HOXD8 and HOXD9 that was expressed at a higher level in BOECs. Quantative PCR confirmed the higher expression of these HOXs in BOECs, a pattern that was shared by a variety of microvascular endothelial cell lines. Subsequently, we analysed publically available microarrays from a variety of adult cell and tissue types using the whole “HOX transcriptome” of all 39 HOX genes. Using hierarchical clustering analysis the HOX transcriptome was able to discriminate endothelial cells from 61 diverse human cell lines of various origins. In a separate publically available microarray dataset of 53 human endothelial cell lines, the HOX transcriptome additionally organized endothelial cells related to their organ or tissue of origin. Human tissue staining for HOXD8 and HOXD9 confirmed endothelial expression and also supported increased microvascular expression of these HOXs. Together these observations suggest a significant involvement of HOX genes in endothelial cell positional identity.

Transcript analysis reveals a specific HOX signature associated with positional identity of human endothelial cells

The distinct topographic Hox expression patterns observed in vascular smooth muscle cells (VSMCs) of the adult cardiovascular system suggest that these transcriptional regulators are critical for maintaining region-specific physiological properties of blood vessels. To test this proposition, we expanded the vascular Hoxc11 expression domain normally restricted to the lower limbs by utilizing an innovative integrated tetracycline regulatory system and Transgelin promoter elements to induce Hoxc11 expression universally in VSMCs of transgenic mice. Ectopic Hoxc11 expression in carotid arteries, aortic arch and descending aorta resulted in drastic vessel wall remodeling involving elastic laminae fragmentation, medial smooth muscle cell loss, and intimal lesion formation. None of these alterations were observed upon induction of Hoxc11 transgene expression in the femoral artery, i.e. the natural Hoxc11 activity domain, although this vessel was greatly enlarged, comparable to the topographically restricted vascular changes seen in Hoxc11^−/− mice. To begin defining Hoxc11-controlled pathways of vascular remodeling, we performed immunolabeling studies in conjunction with co-transfection and chromatin immunoprecipitation (ChIP) assays using mouse vascular smooth muscle (MOVAS) cells. The results suggest direct transcriptional control of two members of the matrix metalloproteinase (Mmp) family, including Mmp2 and Mmp9 that are known as key players in the inception and progression of vascular remodeling events. In summary, the severe vascular abnormalities resulting from the induced dysregulated expression of a Hox gene with regional vascular patterning functions suggest that proper Hox function and regulation is critical for maintaining vascular functional integrity.

Fras1, a basement membrane-associated protein mutated in Fraser syndrome, mediates both the initiation of the mammalian kidney and the integrity of renal glomeruli

FRAS1 is mutated in some individuals with Fraser syndrome (FS) and the encoded protein is expressed in embryonic epidermal cells, localizing in their basement membrane (BM). Syndactyly and cryptophthalmos in FS are sequelae of skin fragility but the bases for associated kidney malformations are unclear. We demonstrate that Fras1 is expressed in the branching ureteric bud (UB), and that renal agenesis occurs in homozygous Fras1 null mutant blebbed (bl) mice on a C57BL6J background. In vivo, the bl/bl bud fails to invade metanephric mesenchyme which undergoes involution, events replicated in organ culture. The expression of glial cell line-derived neurotrophic factor and growth-differentiation factor 11 was defective in bl/bl renal primordia in vivo, whereas, in culture, the addition of either growth factor restored bud invasion into the mesenchyme. Mutant primordia also showed deficient expression of Hoxd11 and Six2 transcription factors, whereas the activity of bone morphogenetic protein 4, an anti-branching molecule, was upregulated. In wild types, Fras1 was also expressed by nascent nephrons. Foetal glomerular podocytes expressed Fras1 transcripts and Fras1 immunolocalized in a glomerular BM-like pattern. On a mixed background, bl mutants, and also compound mutants for bl and my, another bleb strain, sometimes survive into adulthood. These mice have two kidneys, which contain subsets of glomeruli with perturbed nephrin, podocin, integrin α3 and fibronectin expression. Thus, Fras1 protein coats branching UB epithelia and is strikingly upregulated in the nephron lineage after mesenchymal/epithelial transition. Fras1 deficiency causes defective interactions between the bud and mesenchyme, correlating with disturbed expression of key nephrogenic molecules. Furthermore, Fras1 may also be required for the formation of normal glomeruli.

Generation of amplified RNAs and cDNA libraries from single mammalian cells

With the near completion of the human genome sequencing effort, it is now possible to analyze the expression of the entire human gene complement. However, a major obstacle in performing such analysis is the ability to successfully generate enough cDNA or amplified RNA from a limited number of cells, such as biopsies, blood smears, cells obtained by laser capture microscopy, and preimplantation embryonic cells and germ cells. Because these samples yield extremely small amounts of RNA, reproducible methods are needed to amplify this RNA while maintaining the original message profile. A detailed description is given for generating pools of cDNA libraries containing a high proportion of cDNAs enriched with 5'-coding sequences from as little as 1 ng of total RNA using a modified switching mechanism at 5' end of RNA transcript protocol. In addition, the T7-promoter-linked double-stranded cDNAs can be in vitro transcribed linearly using T7-RNA polymerase to generate amplified RNA that is mRNA derived. The cDNA pools can be used directly for gene-specific reverse transcriptase polymerase chain reaction or processed for ligation into vectors of choice whereas the amplified RNA can be used for microarray-based expression profiling.

The transcriptome of human oocytes

The identification of genes and deduced pathways from the mature human oocyte can help us better understand oogenesis, folliculogenesis, fertilization, and embryonic development. Human metaphase II oocytes were used within minutes after removal from the ovary, and its transcriptome was compared with a reference sample consisting of a mixture of total RNA from 10 different normal human tissues not including the ovary. RNA amplification was performed by using a unique protocol. Affymetrix Human Genome U133 Plus 2.0 GeneChip arrays were used for hybridizations. Compared with reference samples, there were 5,331 transcripts significantly up-regulated and 7,074 transcripts significantly down-regulated in the oocyte. Of the oocyte up-regulated probe sets, 1,430 have unknown function. A core group of 66 transcripts was identified by intersecting significantly up-regulated genes of the human oocyte with those from the mouse oocyte and from human and mouse embryonic stem cells. GeneChip array results were validated using RT-PCR in a selected set of oocyte-specific genes. Within the up-regulated probe sets, the top overrepresented categories were related to RNA and protein metabolism, followed by DNA metabolism and chromatin modification. This report provides a comprehensive expression baseline of genes expressed in in vivo matured human oocytes. Further understanding of the biological role of these genes may expand our knowledge on meiotic cell cycle, fertilization, chromatin remodeling, lineage commitment, pluripotency, tissue regeneration, and morphogenesis.

Whole-genome approaches for large-scale gene identification and expression analysis in mammalian preimplantation embryos

The elucidation, unravelling and understanding of the molecular basis of transcriptional control during preimplantion development is of utmost importance if we are to intervene and eliminate or reduce abnormalities associated with growth, disease and infertility by applying assisted reproduction. Importantly, these studies should enhance our knowledge of basic reproductive biology and its application to regenerative medicine and livestock production. A major obstacle impeding progress in these areas is the ability to successfully generate molecular portraits of preimplantation embryos from their minute amounts of RNA. The present review describes the various approaches whereby classical embryology fuses with molecular biology, high-throughput genomics and systems biology to address and solve questions related to early development in mammals.

Differences in expression pattern and function between zebrafish hoxc13 orthologs: recruitment of Hoxc13b into an early embryonic role

Vertebrate Hox genes are generally believed to initiate expression at the primitive streak or early neural plate stages. The timing and spatial restrictions of the Hox expression patterns during these stages correlate well with their demonstrated role in axial patterning. Here we demonstrate that one zebrafish hoxc13 ortholog, hoxc13a, has an expression pattern in the developing tail bud that is consistent with the gene playing a role in axial patterning. However, the second hoxc13 ortholog, hoxc13b, is maternally expressed and is detectable in every cell of early cleavage embryos through gastrulae. In addition, both transcript and protein are detectable at these stages. At 19 h post fertilization (hpf), hoxc13b expression is up-regulated in the tail bud, becoming restricted to the tail bud by 24 hpf. Importantly, by 24 hpf, hoxc13b morphants show a specific developmental delay, which can be rescued by co-injecting synthetic capped hoxc13a or hoxc13b message. These data suggest some functional divergence due to altered expression patterns of the two hoxc13 orthologs after duplication. Further characterization of the hoxc13b morphant delay reveals that it is biphasic in nature, with the first phase of the delay occurring before gastrulation, suggesting a new role for vertebrate Hox genes before their conserved role in axial patterning. The extent of the delay does not change through 20 hpf; however, an additional delay emerges at this time. Notably, this second phase of the delay correlates with hoxc13b expression pattern becoming restricted to the tail bud.