Domestication Explains Two-Thirds of Differential-Gene-Expression Variance between Domestic and Wild Animals; The Remaining One-Third Reflects Intraspecific and Interspecific Variation

Belyaev's concept of destabilizing selection during domestication was a major achievement in the XX century. Its practical value has been realized in commercial colors of the domesticated fox that never occur in the wild and has been confirmed in a wide variety of pet breeds. Many human disease models involving animals allow to test drugs before human testing. Perhaps this is why investigators doing transcriptomic profiling of domestic versus wild animals have searched for breed-specific patterns. Here we sequenced hypothalamic transcriptomes of tame and aggressive rats, identified their differentially expressed genes (DEGs), and, for the first time, applied principal component analysis to compare them with all the known DEGs of domestic versus wild animals that we could find. Two principal components, PC1 and PC2, respectively explained 67% and 33% of differential-gene-expression variance (hereinafter: log2 value) between domestic and wild animals. PC1 corresponded to multiple orthologous DEGs supported by homologs; these DEGs kept the log2 value sign from species to species and from tissue to tissue (i.e., a common domestication pattern). PC2 represented stand-alone homologous DEG pairs reversing the log2 value sign from one species to another and from tissue to tissue (i.e., representing intraspecific and interspecific variation).

Identification of a 42-bp heart-specific enhancer of the notch1b gene in zebrafish embryos

BACKGROUND

NOTCH1 plays a key role in the differentiation of ventricles, and mutations are strongly associated with both sporadic and familial bicuspid aortic valves. However, few heart-specific enhancers have been identified to date.

RESULTS

In this study, we investigated evolutionary conserved regions (ECRs) that might act as potential enhancers within the region approximately 150-kb upstream and downstream of the NOTCH1 gene. Functional validation revealed that one 127-bp ECR located ~85-kb downstream of the NOTCH1 gene drives green fluorescent protein expression in the zebrafish embryo heart. Transcription factor (TF) prediction and core TF distribution analyses were performed to identify the core region. Dissection of ECR3 was performed to identify the 42-bp sequence, which is sufficient for heart-specific expression. In situ hybridization experiments showed that notch1b is expressed in the heart. Overexpression experiments in cells indicated that NKX2-5 is critical for enhancer activity. Mutation of the NKX-5 binding site significantly decreased reporter gene expression. Next, compared with the commonly used myocardium-labeled zebrafish transgenic strain Tg(cmlc2: mCherry), this 42-bp enhancer-labeled stable line mediated a similar expression pattern but with a smaller core region.

CONCLUSION

This study identified a 42-bp heart-specific enhancer near the NOTCH1 gene and further verified its functional targeting by NKX2-5.

Hoxa13 regulates expression of common Hox target genes involved in cartilage development to coordinate the expansion of the autopodal anlage

To elucidate the role of Hox genes in limb cartilage development, we identified the target genes of HOXA11 and HOXA13 by ChIP‐Seq. The ChIP DNA fragment contained evolutionarily conserved sequences and multiple highly conserved HOX binding sites. A substantial portion of the HOXA11 ChIP fragment overlapped with the HOXA13 ChIP fragment indicating that both factors share common targets. Deletion of the target regions neighboring Bmp2 or Tshz2 reduced their expression in the autopod suggesting that they function as the limb bud‐specific enhancers. We identified the Hox downstream genes as exhibiting expression changes in the Hoxa13 knock out (KO) and Hoxd11‐13 deletion double mutant (Hox13 dKO) autopod by Genechip analysis. The Hox downstream genes neighboring the ChIP fragment were defined as the direct targets of Hox. We analyzed the spatial expression pattern of the Hox target genes that encode two different categories of transcription factors during autopod development and Hox13dKO limb bud. (a) Bcl11a, encoding a repressor of cartilage differentiation, was expressed in the E11.5 autopod and was substantially reduced in the Hox13dKO. (b) The transcription factors Aff3, Bnc2, Nfib and Runx1t1 were expressed in the zeugopodal cartilage but not in the autopod due to the repressive or relatively weak transcriptional activity of Hox13 at E11.5. Interestingly, the expression of these genes was later observed in the autopodal cartilage at E12.5. These results indicate that Hox13 transiently suspends the cartilage differentiation in the autopodal anlage via multiple pathways until establishing the paddle‐shaped structure required to generate five digits.

Identification of a limb enhancer that is removed by pathogenic deletions downstream of the SHOX gene

Haploinsufficiency of the human SHOX gene causes Léri-Weill dyschondrosteosis (LWD), characterized by shortening of the middle segments of the limbs and Madelung deformity of the wrist. As many as 35% of LWD cases are caused by deletions of non-coding sequences downstream of SHOX that presumably remove an enhancer or enhancers necessary for SHOX expression in developing limbs. We searched for these active sequences using a transgenic mouse assay and identified a 563 basepair (bp) enhancer with specific activity in the limb regions where SHOX functions. This enhancer has previously escaped notice because of its poor evolutionary conservation, although it does contain 100 bp that are conserved in non-rodent mammals. A primary cell luciferase assay confirmed the enhancer activity of the conserved core sequence and demonstrated that putative HOX binding sites are required for its activity. This enhancer is removed in most non-coding deletions that cause LWD. However, we did not identify any likely pathogenic variants of the enhancer in a screen of 124 LWD individuals for whom no causative mutation had been found, suggesting that only larger deletions in the region commonly cause LWD. We hypothesize that loss of this enhancer contributes to the pathogenicity of deletions downstream of SHOX.

Enhancer redundancy provides phenotypic robustness in mammalian development

Distant-acting tissue-specific enhancers vastly outnumber protein-coding genes in mammalian genomes, but the functional significance of this regulatory complexity remains insufficiently understood^1,2. Here we show that the pervasive presence of multiple enhancers with similar activities near the same gene confers phenotypic robustness to loss-of-function mutations in individual enhancers. We used genome editing to create 23 mouse deletion lines and inter-crosses, including both single and combinatorial enhancer deletions at seven distinct loci required for limb development. Surprisingly, none of ten deletions of individual enhancers caused noticeable changes in limb morphology. In contrast, removal of pairs of limb enhancers near the same gene resulted in discernible phenotypes, indicating that enhancers function redundantly in establishing normal morphology. In a genetic background sensitized by reduced baseline expression of the target gene, even single enhancer deletions caused limb abnormalities, suggesting that functional redundancy is conferred by additive effects of enhancers on gene expression levels. A genome-wide analysis integrating epigenomic and transcriptomic data from 29 developmental mouse tissues revealed that mammalian genes are very commonly associated with multiple enhancers that have similar spatiotemporal activity. Systematic exploration of three representative developmental structures (limb, brain, heart) uncovered more than a thousand cases in which five or more enhancers with redundant activity patterns were found near the same gene. Taken together, our data indicate that enhancer redundancy is a remarkably widespread feature of mammalian genomes and provides an effective regulatory buffer preventing deleterious phenotypic consequences upon loss of individual enhancers.

Comparison of SHOX and associated elements duplications distribution between patients (Lėri-Weill dyschondrosteosis/idiopathic short stature) and population sample

The effect of heterozygous duplications of SHOX and associated elements on Lėri-Weill dyschondrosteosis (LWD) and idiopathic short stature (ISS) development is less distinct when compared to reciprocal deletions. The aim of our study was to compare frequency and distribution of duplications within SHOX and associated elements between population sample and LWD (ISS) patients. A preliminary analysis conducted on Czech population sample of 250 individuals compared to our previously reported sample of 352 ISS/LWD Czech patients indicated that rather than the difference in frequency of duplications it is the difference in their distribution. Particularly, there was an increased frequency of duplications residing to the CNE-9 enhancer in our LWD/ISS sample. To see whether the obtained data are consistent across published studies we made a literature survey to get published cases with SHOX or associated elements duplication and formed the merged LWD, the merged ISS, and the merged population samples. Relative frequency of particular region duplication in each of those merged samples were calculated. There was a significant difference in the relative frequency of CNE-9 enhancer duplications (11 vs. 3) and complete SHOX (exon1-6b) duplications (4 vs. 24) (p-value 0.0139 and p-value 0.000014, respectively) between the merged LWD sample and the merged population sample. We thus propose that partial SHOX duplications and small duplications encompassing CNE-9 enhancer could be highly penetrant alleles associated with ISS and LWD development.

A unique stylopod patterning mechanism by Shox2-controlled osteogenesis

Vertebrate appendage patterning is programmed by Hox-TALE factor-bound regulatory elements. However, it remains unclear which cell lineages are commissioned by Hox-TALE factors to generate regional specific patterns and whether other Hox-TALE co-factors exist. In this study, we investigated the transcriptional mechanisms controlled by the Shox2 transcriptional regulator in limb patterning. Harnessing an osteogenic lineage-specific Shox2 inactivation approach we show that despite widespread Shox2 expression in multiple cell lineages, lack of the stylopod observed upon Shox2 deficiency is a specific result of Shox2 loss of function in the osteogenic lineage. ChIP-Seq revealed robust interaction of Shox2 with cis-regulatory enhancers clustering around skeletogenic genes that are also bound by Hox-TALE factors, supporting a lineage autonomous function of Shox2 in osteogenic lineage fate determination and skeleton patterning. Pbx ChIP-Seq further allowed the genome-wide identification of cis-regulatory modules exhibiting co-occupancy of Pbx, Meis and Shox2 transcriptional regulators. Integrative analysis of ChIP-Seq and RNA-Seq data and transgenic enhancer assays indicate that Shox2 patterns the stylopod as a repressor via interaction with enhancers active in the proximal limb mesenchyme and antagonizes the repressive function of TALE factors in osteogenesis.

Microduplications at the pseudoautosomal SHOX locus in autism spectrum disorders and related neurodevelopmental conditions

BACKGROUND

The pseudoautosomal short stature homeobox-containing (SHOX) gene encodes a homeodomain transcription factor involved in cell-cycle and growth regulation. SHOX/SHOX enhancers deletions cause short stature and skeletal abnormalities in a female-dominant fashion; duplications appear to be rare. Neurodevelopmental disorders (NDDs), such as autism spectrum disorders (ASDs), are complex disorders with high heritability and skewed sex ratio; several rare (<1% frequency) CNVs have been implicated in risk.

METHODS

We analysed data from a discovery series of 90 adult ASD cases, who underwent clinical genetic testing by array-comparative genomic hybridisation (CGH). Twenty-seven individuals harboured CNV abnormalities, including two unrelated females with microduplications affecting SHOX. To determine the prevalence of SHOX duplications and delineate their associated phenotypic spectrum, we subsequently examined array-CGH data from a follow-up sample of 26 574 patients, including 18 857 with NDD (3541 with ASD).

RESULTS

We found a significant enrichment of SHOX microduplications in the NDD cases (p=0.00036; OR 2.21) and, particularly, in those with ASD (p=9.18×10(-7); OR 3.63) compared with 12 594 population-based controls. SHOX duplications affecting the upstream or downstream enhancers were enriched only in females with NDD (p=0.0043; OR 2.69/p=0.00020; OR 7.20), but not in males (p=0.404; OR 1.38/p=0.096; OR 2.21).

CONCLUSIONS

Microduplications at the SHOX locus are a low penetrance risk factor for ASD/NDD, with increased risk in both sexes. However, a concomitant duplication of SHOX enhancers may be required to trigger a NDD in females. Since specific SHOX isoforms are exclusively expressed in the developing foetal brain, this may reflect the pathogenic effect of altered SHOX protein dosage on neurodevelopment.

Short Stature in Isodicentric Y Chromosome and Three Copies of the SHOX Gene: Clinical Report and Review of Literature

Short stature homeobox gene (SHOX) mutations and pseudoautosomal region 1 (PAR1) deletions encompassing SHOX are known causes of Léri-Weill dyschondrosteosis and isolated short stature, while 3 copies of SHOX in cases with triple sex chromosome constitution are responsible for tall stature. Duplications involving SHOX have been rarely reported, and they were found in individuals with short, normal and tall stature. An adopted boy with short stature, isodicentric Y chromosome and 3 copies of SHOX is described. Normal growth hormone (GH) secretion and insulin-like growth factor 1 (IGF1) increase during an IGF1 generation test were found, ruling out impaired GH-IGF1 axis. No other organic or psychiatric causes of impaired growth were found. GH treatment improved linear growth, as reported in children with SHOX haploinsufficiency. This new report and the review of literature support that SHOX duplication may cause short stature, especially in those children with duplications of the 5'SHOX regulatory elements. Chromosome analysis and detailed molecular characterization of the duplicated region should be warranted in individuals with SHOX duplications in order to investigate the presence of occult chromosome imbalance. Additional reports and follow-up till adult height are needed to give conclusions on long-term efficacy and safety of GH treatment in short children with SHOX duplication.

Duplications upstream and downstream of SHOX identified as novel causes of Leri-Weill dyschondrosteosis or idiopathic short stature

Leri-Weill dyschondrosteosis is a pseudoautosomal dominantly-inherited skeletal dysplasia ascribed to haploinsufficiency of the SHOX gene caused by deletions, point mutations, or partial duplications of the gene, or to heterozygous deletions upstream or downstream of the intact SHOX gene involving conserved non-coding cis-regulatory DNA elements that show enhancer activity. Recently, two SHOX conserved non-coding element duplications, one upstream and one downstream, were reported in patients referred with idiopathic short stature. To further evaluate the role of these duplications in SHOX-related disorders, we describe seven patients (five with Leri-Weill dyschondrosteosis and two with short stature) all of whom have duplications of part of the upstream or downstream conserved non-coding element regions, identified by multiplex ligation-dependent probe amplification. In addition, we show data from 32 patients with an apparently identical downstream duplication that includes a proposed putative regulatory element (identified by multiplex ligation-dependent probe amplification or array comparative genome hybridization), which results in a variable phenotype from normal to mild Leri-Weill dyschondrosteosis. These additional data provide further evidence that duplications of upstream and downstream long range cis-regulatory DNA elements can result in a SHOX-related phenotype.

Functionally conserved enhancers with divergent sequences in distant vertebrates

BACKGROUND

To examine the contributions of sequence and function conservation in the evolution of enhancers, we systematically identified enhancers whose sequences are not conserved among distant groups of vertebrate species, but have homologous function and are likely to be derived from a common ancestral sequence. Our approach combined comparative genomics and epigenomics to identify potential enhancer sequences in the genomes of three groups of distantly related vertebrate species.

RESULTS

We searched for sequences that were conserved within groups of closely related species but not between groups of more distant species, and were associated with an epigenetic mark of enhancer activity. To facilitate inferring orthology between non-conserved sequences, we limited our search to introns whose orthology could be unambiguously established by mapping the bracketing exons. We show that a subset of these non-conserved but syntenic sequences from the mouse and zebrafish genomes have homologous functions in a zebrafish transgenic enhancer assay. The conserved expression patterns driven by these enhancers are probably associated with short transcription factor-binding motifs present in the divergent sequences.

CONCLUSIONS

We have identified numerous potential enhancers with divergent sequences but a conserved function. These results indicate that selection on function, rather than sequence, may be a common mode of enhancer evolution; evidence for selection at the sequence level is not a necessary criterion to define a gene regulatory element.

Shox2 is required for the proper development of the facial motor nucleus and the establishment of the facial nerves

Background

Axons from the visceral motor neurons (vMNs) project from nuclei in the hindbrain to innervate autonomic ganglia and branchial arch-derived muscles. Although much is known about the events that govern specification of somatic motor neurons, the genetic pathways responsible for the development of vMNs are less well characterized. We know that vMNs, like all motor neurons, depend on sonic hedgehog signaling for their generation. Similarly, the paired-like homeobox 2b (Phox2b) gene, which is expressed in both proliferating progenitors and post-mitotic motor neurons, is essential for the development of vMNs. Given that our previous study identified a novel role for the short stature homeobox 2 (Shox2) gene in the hindbrain, and since SHOX2 has been shown to regulate transcription of islet 1 (Isl1), an important regulator of vMN development, we sought to determine whether Shox2 is required for the proper development of the facial motor nucleus.

Results

Using a Nestin-Cre driver, we show that elimination of Shox2 throughout the brain results in elevated cell death in the facial motor nucleus at embryonic day 12.5 (E12.5) and E14.5, which correlates with impaired axonal projection properties of vMNs. We also observed changes in the spatial expression of the vMN cell fate factors Isl1 and Phox2b, and concomitant defects in Shh and Ptch1 expression in Shox2 mutants. Furthermore, we demonstrate that elimination of Shox2 results in the loss of dorsomedial and ventromedial subnuclei by postnatal day 0 (P0), which may explain the changes in physical activity and impaired feeding/nursing behavior in Shox2 mutants.

Conclusions

Combined, our data show that Shox2 is required for development of the facial motor nucleus and its associated facial (VII) nerves, and serves as a new molecular tool to probe the genetic programs of this understudied hindbrain region.

Electronic supplementary material

The online version of this article (doi:10.1186/s12868-015-0176-0) contains supplementary material, which is available to authorized users.

A common Shox2-Nkx2-5 antagonistic mechanism primes the pacemaker cell fate in the pulmonary vein myocardium and sinoatrial node

In humans, atrial fibrillation is often triggered by ectopic pacemaking activity in the myocardium sleeves of the pulmonary vein (PV) and systemic venous return. The genetic programs that abnormally reinforce pacemaker properties at these sites and how this relates to normal sinoatrial node (SAN) development remain uncharacterized. It was noted previously that Nkx2-5, which is expressed in the PV myocardium and reinforces a chamber-like myocardial identity in the PV, is lacking in the SAN. Here we present evidence that in mice Shox2 antagonizes the transcriptional output of Nkx2-5 in the PV myocardium and in a functional Nkx2-5(+) domain within the SAN to determine cell fate. Shox2 deletion in the Nkx2-5(+) domain of the SAN caused sick sinus syndrome, associated with the loss of the pacemaker program. Explanted Shox2(+) cells from the embryonic PV myocardium exhibited pacemaker characteristics including node-like electrophysiological properties and the capability to pace surrounding Shox2(-) cells. Shox2 deletion led to Hcn4 ablation in the developing PV myocardium. Nkx2-5 hypomorphism rescued the requirement for Shox2 for the expression of genes essential for SAN development in Shox2 mutants. Similarly, the pacemaker-like phenotype induced in the PV myocardium in Nkx2-5 hypomorphs reverted back to a working myocardial phenotype when Shox2 was simultaneously deleted. A similar mechanism is also adopted in differentiated embryoid bodies. We found that Shox2 interacts with Nkx2-5 directly, and discovered a substantial genome-wide co-occupancy of Shox2, Nkx2-5 and Tbx5, further supporting a pivotal role for Shox2 in the core myogenic program orchestrating venous pole and pacemaker development.