A dimerized HMX1 inhibits EPHA6/epha4b in mouse and zebrafish retinas

Significance of Premature Vertebral Mineralization in Zebrafish Models in Mechanistic and Pharmaceutical Research on Hereditary Multisystem Diseases

Zebrafish are increasingly becoming an important model organism for studying the pathophysiological mechanisms of human diseases and investigating how these mechanisms can be effectively targeted using compounds that may open avenues to novel treatments for patients. The zebrafish skeleton has been particularly instrumental in modeling bone diseases as-contrary to other model organisms-the lower load on the skeleton of an aquatic animal enables mutants to survive to early adulthood. In this respect, the axial skeletons of zebrafish have been a good read-out for congenital spinal deformities such as scoliosis and degenerative disorders such as osteoporosis and osteoarthritis, in which aberrant mineralization in humans is reflected in the respective zebrafish models. Interestingly, there have been several reports of hereditary multisystemic diseases that do not affect the vertebral column in human patients, while the corresponding zebrafish models systematically show anomalies in mineralization and morphology of the spine as their leading or, in some cases, only phenotype. In this review, we describe such examples, highlighting the underlying mechanisms, the already-used or potential power of these models to help us understand and amend the mineralization process, and the outstanding questions on how and why this specific axial type of aberrant mineralization occurs in these disease models.

Premature Vertebral Mineralization in hmx1-Mutant Zebrafish

H6 family homeobox 1 (HMX1) regulates multiple aspects of craniofacial development, and mutations in HMX1 are linked to an ocular defect termed oculoauricular syndrome of Schorderet–Munier–Franceschetti (OAS) (MIM #612109). Recently, additional altered orofacial features have been reported, including short mandibular rami, asymmetry of the jaws, and altered premaxilla. We found that in two mutant zebrafish lines termed hmx1mut10 and hmx1mut150, precocious mineralization of the proximal vertebrae occurred. Zebrafish hmx1mut10 and hmx1mut150 report mutations in the SD1 and HD domains, which are essential for dimerization and activity of hmx1. In hmx1mut10, the bone morphogenetic protein (BMP) antagonists chordin and noggin1 were downregulated, while bmp2b and bmp4 were highly expressed and specifically localized to the dorsal region prior to the initiation of the osteogenic process. The osteogenic promoters runx2b and spp1 were also upregulated. Supplementation with DMH1—an inhibitor of the BMP signaling pathway—at the specific stage in which bmp2b and bmp4 are highly expressed resulted in reduced vertebral mineralization, resembling the wildtype mineralization progress of the axial skeleton. These results point to a possible role of hmx1 as part of a complex gene network that inhibits bmp2b and bmp4 in the dorsal region, thus regulating early axial skeleton development.

Developmental gene regulatory network connections predicted by machine learning from gene expression data alone

Gene regulatory network (GRN) inference can now take advantage of powerful machine learning algorithms to complement traditional experimental methods in building gene networks. However, the dynamical nature of embryonic development-representing the time-dependent interactions between thousands of transcription factors, signaling molecules, and effector genes-is one of the most challenging arenas for GRN prediction. In this work, we show that successful GRN predictions for a developmental network from gene expression data alone can be obtained with the Priors Enriched Absent Knowledge (PEAK) network inference algorithm. PEAK is a noise-robust method that models gene expression dynamics via ordinary differential equations and selects the best network based on information-theoretic criteria coupled with the machine learning algorithm Elastic Net. We test our GRN prediction methodology using two gene expression datasets for the purple sea urchin, Stronglyocentrotus purpuratus, and cross-check our results against existing GRN models that have been constructed and validated by over 30 years of experimental results. Our results find a remarkably high degree of sensitivity in identifying known gene interactions in the network (maximum 81.58%). We also generate novel predictions for interactions that have not yet been described, which provide a resource for researchers to use to further complete the sea urchin GRN. Published ChIPseq data and spatial co-expression analysis further support a subset of the top novel predictions. We conclude that GRN predictions that match known gene interactions can be produced using gene expression data alone from developmental time series experiments.

Single cell transcriptomics of the developing zebrafish lens and identification of putative controllers of lens development

The vertebrate lens is a valuable model system for investigating the gene expression changes that coordinate tissue differentiation due to its inclusion of two spatially separated cell types, the outer epithelial cells and the deeper denucleated fiber cells that they support. Zebrafish are a useful model system for studying lens development given the organ's rapid development in the first several days of life in an accessible, transparent embryo. While we have strong foundational knowledge of the diverse lens crystallin proteins and the basic gene regulatory networks controlling lens development, no study has detailed gene expression in a vertebrate lens at single cell resolution. Here we report an atlas of lens gene expression in zebrafish embryos and larvae at single cell resolution through five days of development, identifying a number of novel putative regulators of lens development. Our data address open questions about the temperospatial expression of α-crystallins during lens development that will support future studies of their function and provide the first detailed view of β- and γ-crystallin expression in and outside the lens. We describe divergent expression in transcription factor genes that occur as paralog pairs in the zebrafish. Finally, we examine the expression dynamics of cytoskeletal, membrane associated, RNA-binding, and transcription factor genes, identifying a number of novel patterns. Overall these data provide a foundation for identifying and characterizing lens developmental regulatory mechanisms and revealing targets for future functional studies with potential therapeutic impact.

Hmx1 regulates urfh1 expression in the craniofacial region in zebrafish

H6 family homeobox 1 (HMX1) regulates multiple aspects of craniofacial development as it is widely expressed in the eye, peripheral ganglia and branchial arches. Mutations in HMX1 are linked to an ocular defect termed Oculo-auricular syndrome of Schorderet-Munier-Franceschetti (MIM #612109). We identified UHRF1 as a target of HMX1 during development. UHRF1 and its partner proteins actively regulate chromatin modifications and cellular proliferation. Luciferase assays and in situ hybridization analyses showed that HMX1 exerts a transcriptional inhibitory effect on UHRF1 and a modification of its expression pattern. Overexpression of hmx1 in hsp70-hmx1 zebrafish increased uhrf1 expression in the cranial region, while mutations in the hmx1 dimerization domains reduced uhrf1 expression. Moreover, the expression level of uhrf1 and its partner dnmt1 was increased in the eye field in response to hmx1 overexpression. These results indicate that hmx1 regulates uhrf1 expression and, potentially through regulating the expression of factors involved in DNA methylation, contribute to the development of the craniofacial region of zebrafish.

Hmx3a Has Essential Functions in Zebrafish Spinal Cord, Ear and Lateral Line Development

Transcription factors that contain a homeodomain DNA-binding domain have crucial functions in most aspects of cellular function and embryonic development in both animals and plants. Hmx proteins are a subfamily of NK homeodomain-containing proteins that have fundamental roles in development of sensory structures such as the eye and the ear. However, Hmx functions in spinal cord development have not been analyzed. Here, we show that zebrafish (Danio rerio) hmx2 and hmx3a are coexpressed in spinal dI2 and V1 interneurons, whereas hmx3b, hmx1, and hmx4 are not expressed in spinal cord. Using mutational analyses, we demonstrate that, in addition to its previously reported role in ear development, hmx3a is required for correct specification of a subset of spinal interneuron neurotransmitter phenotypes, as well as correct lateral line progression and survival to adulthood. Surprisingly, despite similar expression patterns of hmx2 and hmx3a during embryonic development, zebrafish hmx2 mutants are viable and have no obviously abnormal phenotypes in sensory structures or neurons that require hmx3a In addition, embryos homozygous for deletions of both hmx2 and hmx3a have identical phenotypes to severe hmx3a single mutants. However, mutating hmx2 in hypomorphic hmx3a mutants that usually develop normally, results in abnormal ear and lateral line phenotypes. This suggests that while hmx2 cannot compensate for loss of hmx3a, it does function in these developmental processes, although to a much lesser extent than hmx3a More surprisingly, our mutational analyses suggest that Hmx3a may not require its homeodomain DNA-binding domain for its roles in viability or embryonic development.

Breakpoint Mapping of Symptomatic Balanced Translocations Links the EPHA6, KLF13 and UBR3 Genes to Novel Disease Phenotype

De novo balanced chromosomal aberrations (BCAs), such as reciprocal translocations and inversions, are genomic aberrations that, in approximately 25% of cases, affect the human phenotype. Delineation of the exact structure of BCAs may provide a precise diagnosis and/or point to new disease loci. We report on six patients with de novo balanced chromosomal translocations (BCTs) and one patient with a de novo inversion, in whom we mapped breakpoints to a resolution of 1 bp, using shallow whole-genome mate pair sequencing. In all seven cases, a disruption of at least one gene was found. In two patients, the phenotypic impact of the disrupted genes is well known (NFIA, ATP7A). In five patients, the aberration damaged genes: PARD3, EPHA6, KLF13, STK24, UBR3, MLLT10 and TLE3, whose influence on the human phenotype is poorly understood. In particular, our results suggest novel candidate genes for retinal degeneration with anophthalmia (EPHA6), developmental delay with speech impairment (KLF13), and developmental delay with brain dysembryoplastic neuroepithelial tumor (UBR3). In conclusion, identification of the exact structure of symptomatic BCTs using next generation sequencing is a viable method for both diagnosis and finding novel disease candidate genes in humans.

Identification of Candidate Genes for Reactivity in Guzerat (Bos indicus) Cattle: A Genome-Wide Association Study

Temperament is fundamental to animal production due to its direct influence on the animal-herdsman relationship. When compared to calm animals, the aggressive, anxious or fearful ones exhibit less weight gain, lower reproductive efficiency, decreased milk production and higher herd maintenance costs, all of which contribute to reduced profits. However, temperament is a trait that is complex and difficult to assess. Recently, a new quantitative system, REATEST^®, for assessing reactivity, a phenotype of temperament, was developed. Herein, we describe the results of a Genome-wide association study for reactivity, assessed using REATEST^® with a sample of 754 females from five dual-purpose (milk and meat production) Guzerat (Bos indicus) herds. Genotyping was performed using a 50k SNP chip and a two-step mixed model approach (Grammar-Gamma) with a one-by-one marker regression was used to identify QTLs. QTLs for reactivity were identified on chromosomes BTA1, BTA5, BTA14, and BTA25. Five intronic and two intergenic markers were significantly associated with reactivity. POU1F1, DRD3, VWA3A, ZBTB20, EPHA6, SNRPF and NTN4 were identified as candidate genes. Previous QTL reports for temperament traits, covering areas surrounding the SNPs/genes identified here, further corroborate these associations. The seven genes identified in the present study explain 20.5% of reactivity variance and give a better understanding of temperament biology.

Retinoic Acid Signaling Regulates Differential Expression of the Tandemly-Duplicated Long Wavelength-Sensitive Cone Opsin Genes in Zebrafish

The signaling molecule retinoic acid (RA) regulates rod and cone photoreceptor fate, differentiation, and survival. Here we elucidate the role of RA in differential regulation of the tandemly-duplicated long wavelength-sensitive (LWS) cone opsin genes. Zebrafish embryos were treated with RA from 48 hours post-fertilization (hpf) to 75 hpf, and RNA was isolated from eyes for microarray analysis. ~170 genes showed significantly altered expression, including several transcription factors and components of cellular signaling pathways. Of interest, the LWS1 opsin gene was strongly upregulated by RA. LWS1 is the upstream member of the tandemly duplicated LWS opsin array and is normally not expressed embryonically. Embryos treated with RA 48 hpf to 100 hpf or beyond showed significant reductions in LWS2-expressing cones in favor of LWS1-expressing cones. The LWS reporter line, LWS-PAC(H) provided evidence that individual LWS cones switched from LWS2 to LWS1 expression in response to RA. The RA signaling reporter line, RARE:YFP indicated that increased RA signaling in cones was associated with this opsin switch, and experimental reduction of RA signaling in larvae at the normal time of onset of LWS1 expression significantly inhibited LWS1 expression. A role for endogenous RA signaling in regulating differential expression of the LWS genes in postmitotic cones was further supported by the presence of an RA signaling domain in ventral retina of juvenile zebrafish that coincided with a ventral zone of LWS1 expression. This is the first evidence that an extracellular signal may regulate differential expression of opsin genes in a tandemly duplicated array.