- Research article
- Open Access
A de novo silencer causes elimination of MITF-M expression and profound hearing loss in pigs
- Lei Chen†1, 3,
- Weiwei Guo†2,
- Lili Ren†2,
- Mingyao Yang4,
- Yaofeng Zhao1,
- Zongyi Guo3,
- Haijin Yi2,
- Mingzhou Li4,
- Yiqing Hu1,
- Xi Long3,
- Boyuan Sun4,
- Jinxiu Li1,
- Suoqiang Zhai2,
- Tinghuan Zhang3,
- Shilin Tian4,
- Qingyong Meng1,
- Ning Yu2,
- Dan Zhu3,
- Guoqing Tang4,
- Qianzi Tang4,
- Liming Ren1,
- Ke Liu2,
- Shihua Zhang3,
- Tiandong Che4,
- Zhengquan Yu1,
- Nan Wu2,
- Lan Jing3,
- Ran Zhang1,
- Tao Cong2,
- Siqing Chen3,
- Yiqiang Zhao1,
- Yue Zhang2,
- Xiaoqing Bai3,
- Ying Guo1,
- Lidong Zhao2,
- Fengming Zhang3,
- Hui Zhao2,
- Liang Zhang3,
- Zhaohui Hou2,
- Jiugang Zhao3,
- Jianan Li2,
- Lijuan Zhang3,
- Wei Sun5,
- Xiangang Zou3,
- Tao Wang3,
- Liangpeng Ge3,
- Zuohua Liu3,
- Xiaoxiang Hu1,
- Jingyong Wang3Email author,
- Shiming Yang2Email author and
- Ning Li1Email author
© Chen et al. 2016
Received: 15 February 2016
Accepted: 10 June 2016
Published: 27 June 2016
Genesis of novel gene regulatory modules is largely responsible for morphological and functional evolution. De novo generation of novel cis-regulatory elements (CREs) is much rarer than genomic events that alter existing CREs such as transposition, promoter switching or co-option. Only one case of de novo generation has been reported to date, in fish and without involvement of phenotype alteration. Yet, this event likely occurs in other animals and helps drive genetic/phenotypic variation.
Using a porcine model of spontaneous hearing loss not previously characterized we performed gene mapping and mutation screening to determine the genetic foundation of the phenotype. We identified a mutation in the non-regulatory region of the melanocyte-specific promoter of microphthalmia-associated transcription factor (MITF) gene that generated a novel silencer. The consequent elimination of expression of the MITF-M isoform led to early degeneration of the intermediate cells of the cochlear stria vascularis and profound hearing loss, as well as depigmentation, all of which resemble the typical phenotype of Waardenburg syndrome in humans. The mutation exclusively affected MITF-M and no other isoforms. The essential function of Mitf-m in hearing development was further validated using a knock-out mouse model.
Elimination of the MITF-M isoform alone is sufficient to cause deafness and depigmentation. To our knowledge, this study provides the first evidence of a de novo CRE in mammals that produces a systemic functional effect.
Genome variation in non-coding regions may not only disrupt regulatory elements but can also create them . Due to the sequence-specific nature of transcription factors’ binding to DNA, introduction of a single-nucleotide polymorphism (SNP) or a small-sized insertion or deletion (indel) in the target regulatory sequence can incapacitate or inhibit binding [1, 2]. Moreover, these relatively small mutations can be sufficient to generate a novel cis-regulatory element (CRE). A growing body of evidence has shown that new regulatory behavior can accompany the modification of existing functional elements [3–5], with point mutations and deletions being sufficient to generate cis-regulatory divergences [6, 7].
One of the most well studied genetic events capable of creating new regulatory behaviors is the transposon, a transposable genetic element that has been shown to generate functional changes in existing enhancers; the extensive studies on transposons have revealed their activities in insects and plants as well as mammals [8, 9]. In contrast, de novo generation of a novel CRE is less well studied. Only one case has been reported to date, that of an enhancer generated by whole-genome duplication in sequences in fish that were demonstrated as formerly lacking cis-regulatory activity . No such cases have been reported for mammals, and there are no reports of minor changes driving either de novo genesis of regulatory elements or systemic functional alterations.
Waardenburg syndrome type 2 (WS2; OMIM #193510) is a hereditary sensorineural deafness syndrome caused by gene mutations; additional physical traits of the disease in humans include heterochromia iridis and white forelock [11, 12]. A common mutation found in WS2 sufferers involves the microphthalmia-associated transcription factor (MITF), which plays a critical role in melanocytes and melanoma [13, 14]. The Mitf multi-promoter gene encodes at least seven isoforms of MITF, each with a distinct N-termini. These seven isoforms have been identified in humans and mice, and are known to be translated from different transcriptional variants with differing first and/or second 5’-end exons . The various promoters associated with each isoform contribute to their tissue-specific expression and functions .
Profound hearing loss in a porcine model
Albino pigs (Fig. 1a) spontaneously arising from a native breed of swine in Southwest China (Rongchang pigs ) are well-studied for their observed phenotypes of deafness and depigmentation, similar to the phenotype of WS2 [11, 12]. Results from auditory brainstem response (ABR) tests show that the albino pigs produced no recognizable waveforms up to 100 dB sound pressure level (SPL) stimuli in the range from 4–32 kHz, whereas normal pigs produced ABR thresholds at 5–10 dB SPL (Fig. 1b). Loss of hair cells and stereocilia bundles were observed in the cochleae of the albino pigs by scanning electron microscopy (SEM; Fig. 1c). Because the hearing loss observed in human cases of WS2 is attributed to abnormal cochlear stria vascularis (SV) , the morphology of SV was examined in our study’s albino pigs using light microscopy and transmission electron microscopy (TEM). As shown in Fig. 1d and e, the albino pigs lacked intermediate cells and had thinner SVs, consisting of two layers of cells only. Because the major functions of the SV are secretion of potassium ions and production of endolymphatic potential (EP), we recorded EPs and measured the [K+] in the scala media of the cochlea. The EP and [K+] were significantly lower than those of normal pigs (P < 0.001, Student’s t-test; Fig. 1f, g). Since EP and high [K+] in the endolymph are reportedly the driving force for mechanotransduction in cochlear hair cells , a reduction in EP can lead to profound hearing loss. All these phenotypes were tested at postnatal day 13. Thus, these results confirmed the phenotype of profound hearing loss related to cochlear morphology defects in our study’s albino pigs. Considering that eye defects have also been observed in some WS2 patients, we also assessed the morphology of porcine eyes. The irises of the albino pigs presented with pale coloration due to lack of pigmentation (Additional file 1: Figure S1A). The paraffin-embedded sections of retinae showed hypopigmentation in the choroid, but the retinal pigment epithelium was normal in both the normal and albino pigs (Additional file 1: Figure S1B).
Gene mapping and mutation screening
To achieve more fine mapping of the causative mutation, we carried out gene screening for 12 mutant (Mitf r/r ) and 12 wild-type (Mitf R/R ) pigs. The porcine Mitf gene was not completely annotated in the reference genome (Sus scrofa genome 10.2) due to the poor assembly quality in this region. Therefore, we performed a homology annotation using the human Mitf mRNA from the Reference Sequence (RefSeq) database and the scaffold (TP_scaffold_24421) from our recently reported genome of Tibetan wild boars . We found that 15 exons, spanning 243 kb of consecutive sequence, had a high syntenic relationship with the Mitf genes in both human and mouse (Additional file 2: Table S2). We next amplified and sequenced all exons, exon-intron boundaries and proximal promoters of the full Mitf gene, out to 10 kb upstream of the transcription start site of the encoded MITF-M isoform. In total, 21 co-segregated variants were identified, including 16 SNPs, four insertions, and one deletion (Additional file 2: Table S3). A C > A non-synonymous variant, which induced a N106K transition, was detected in the protein-coding region (exon 3). However, because only the Mitf R/R pigs carried this variant, it seemed unlikely to be correlated with the hearing loss phenotype. The other 20 variants were located in regulatory regions.
As the associated region we discovered was too large for complete screening, the causative mutation could have been missed due to incomplete coverage by Sanger sequencing. To help rule out this possibility, we used whole-genome re-sequencing data for three of the Mitf r/r pigs and three of the Mitf R/R pigs (Additional file 2: Table S4). A total of 1711 SNPs were detected in the associated region (SSC13: 56,170,062 to 56,933,573; Additional file 3: Table S5) and 961 of these co-segregated with the hearing loss phenotype, including 362 ambiguous SNPs with missing data (Additional file 3: Table S6). Furthermore, 103 open-access porcine whole genome re-sequencing data sets (Additional file 2: Table S4) were included in the analysis to filter out common variants, which are not expected to be related to the hearing loss phenotype. Because this phenotype has not been reported in any other pig breed beyond the Rongchang breed, the 103 pigs from the open-access database were regarded as wild-type pigs (MITF R/R ). Next, 946 of the abovementioned 961 co-segregated SNPs were sought in the 103 MITF R/R pigs (Additional file 3: Table S6); only 15 of the co-segregated SNPs were identified as carried exclusively by the Mitf r/r pigs (Additional file 3: Table S7). Additionally, nine of those 21 co-segregated variants detected in the mutation screening noted above were excluded due to their existence in any Mitf R/R pigs (marked with bold text in Additional file 2: Table S3). Combining the data from our mutation screening and re-sequencing analysis provided a total of 26 co-segregated variants that were deemed as candidate mutations for further study (Additional file 4: Table S8).
Mitf expression analysis
Moreover, we found that the Mitf-m transcript was normally expressed in Mitf R/R cochlea, but not in Mitf r/r cochlea at any developmental stage (Fig. 3b). No obvious differences were found in the expression levels of Mitf-a or Mitf-h between the Mitf R/R and Mitf r/r cochlea (Fig. 3b). RNA-seq assay was used to detect differences between the transcriptome profiles of the Mitf R/r and Mitf r/r SVs at embryonic day 85; at this time in the development, melanocytes remained in the Mitf r/r SVs (Additional file 1: Figure S3). A total of 28 genes showed two-fold differential expression, including some genes with known functions related to pigmentation and melanogenesis, but the Mitf gene was not among them (Additional file 4: Table S9). Because the algorithm used for calculating the gene RPKM values (reads per kilobase per million mapped reads, which serve to estimate gene expression) cannot separate the expression of Mitf-m from other equally expressed transcript variants, we obtained the normalized read count for each exon individually, and found that expression of the M-exon in wild-type SVs was approximately 11.5-fold higher than in mutant SVs (Fig. 3d). In addition, RNA-seq data indicated that expression of some marker genes of melanocytes or the SV intermediate cells remained in the MITF r/r SVs; these genes included S100, KIT, MUM1, and Kir1.2 (Additional file 4: Table S10). Most of the genes expressed in MITF r/r SVs, however, did not show significantly lower expression than the genes in the MITF R/r SVs. When we considered these results along with those from TEM analysis of prenatal SVs (Additional file 1: Figure S3), we determined that the intermediate cells remained in the MITF r/r SVs at the embryo stage and disappeared around birth. Thus, the expression differences between MITF R/r and MITF r/r SVs that we observed were indeed caused by Mitf mutation, rather than a lack of intermediate cells.
Subsequent immunoblotting assay showed an undetectable level of polypeptides of 55–70 kDa (the reported size range of the MITF-M isoform in melanoma ) in the Mitf r/r cochleae (Fig. 3c). Consistent with the previous data, the levels of MITF-A and MITF-H detected by immunoblot (Fig. 3c) were similar between the Mitf R/R and Mitf r/r cochleae. These results indicated that the expression of MITF-M was eliminated in the Mitf r/r pigs at both the transcript and protein levels, and this differential expression pattern itself indicated the existence of regulatory mutations in M isoform-specific regions. Consistently, eight of the 26 co-segregated variants were located in the M isoform-specific promoter (M-promoter; Additional file 4: Table S8).
Transcriptional activity analysis of the MITF M-promoter
To address whether these candidate mutations were able to alter the protein binding landscape in the involved sequence region, we used electrophoretic mobility shift assay (EMSA) to compare the binding capacity of the wild-type and mutant sequences. Two sets of oligonucleotides, R1 and r1 (Fig. 4c, red and yellow highlight, respectively), which differed only in the 9-bp insertion and the GC>TT replacement, and R2 and r2, which differed only in the 14-bp insertion (Fig. 4c, Additional file 4: Table S11), were incubated with nuclear extracts from mouse B16 melanoma cells. Only one differential complex (C1 in Fig. 4d) was formed, namely that with the r2 probe and lacking the R2 probe. The specificity of the complex was confirmed by competition EMSA, wherein a 50- and 100-fold molar excess of unlabeled r2 probe showed effective binding competition but a 50- and 2100-fold excess of unlabeled R2 probe did not. Another complex (C2 in Fig. 4d) formed was a non-specific complex, because no competition was observed with even 100-fold molar excess r2 cold probe. No difference in protein binding was observed between the R1 and r1 probes (Additional file 1: Figure S4). Collectively, these results show that only the 14-bp insertion can induce specific transcription factor binding events, while the 9-bp insertion and the two continuous SNPs do not. Thus, the 14-bp insertion is the only variation detected in our study that was considered as potentially responsible for the observed down-regulation of Mitf-m transcription activity.
SOX proteins (SOX2, SOX3 and SOX9) can regulate inner ear development [30–32]. Our RNA-seq data showed that Sox2 and Sox9, but not Sox3, were expressed in porcine SVs (Additional file 4: Table S12), suggesting that SOX2 and/or SOX9 is capable of binding ectopically to a new negative CRE (i.e. a silencer), which may have been generated by the 14-bp short insertion in the Mitf regulatory region and which may have caused the abrogation of Mitf-m expression (Fig. 6). Thus, the phenotypes of hearing loss and depigmentation in the albino pigs appear to be caused, at least partially, by the 14-bp insertion that is located at –7532 bp relative to the transcription start site of Mitf-m. Additionally, we performed a further genotyping analysis of the causative mutant in 311 individual Rongchang pigs (Additional file 1: Figure S5). Both the 14- and 9-bp insertions were found to be carried by all albino pigs, and to co-segregate with the albino phenotype completely (Additional file 4: Table S13). No recombination event involving the two insertions was observed.
These results provide convincing evidence that even a small insertion in a region that lacks regulatory activity in the melanocyte lineage can create a transcription factor binding site (TFBS). Therefore, we next investigated whether this region was a non-regulatory sequence in other cell/species lineages. By exploiting the available data from the human ENCODE project, which had previously identified large numbers of regulatory elements, we were able to obtain the human ortholog of the mutant region identified in pigs and investigate whether this region overlapped with the CREs reported in the ENCODE data (Additional file 1: Figure S6A). The published RNA-seq data indicated that the wild-type M-exon was indeed transcribed in the human melanocyte lineage (Additional file 1: Figure S6C).
Using the UCSC browser to search the ENCODE data, we noted the following characteristics of these putative CREs. (1) The flanking regions of the causative mutant point has a relatively low level of conservation in mammals, suggesting a low probability of conserved functional CREs (Additional file 1: Figure S6B). (2) DNAse I hypersensitivity data indicate that chromatin accessibility around the mutant point is low in melanocytes and various cell lineages (Additional file 1: Figure S6D). (3) Chromatin immunoprecipitation and DNAse I footprinting data provide no evidence of protein binding sites near the mutant point (Additional file 1: Figure S6E). (4) The H3K4Me1 and H3K27Ac histone markers show no evidence of CREs in the flanking region of the mutant point (Additional file 1: Figure S6F). Because none of these data supported the existence of CREs near the mutant point, we concluded that the 14-bp insertion found in the albino pigs resulted in the de novo genesis of a silencer in the M-promoter (Fig. 6).
Mitf-m-specific mutations result in hearing loss in a mouse model
Genetic variations can create regulatory elements. Within the 10 million years of evolution involving modern organisms, hundreds of non-functional sequences have gained de novo regulatory functions . Induced de novo cis-regulatory behaviors have been reported in the literature as having emerged via the following four genetic mechanisms. (1) Duplicated coding genes may lose their coding functions and gain de novo regulatory functions through a whole-genome duplication event; these de novo enhancers are consequently known as “recycle regions” . (2) Point mutations can create new regulatory behaviors and are commonly observed, from flies to primates [6, 34], indicating their evolutionary benefit. (3) Deletions can also generate new enhancers by bringing together flanking sequences that create novel binding sites for activators . (4) Transposons can serve as a source for acquiring new regulatory behaviors [8, 9, 35]. However, to our knowledge, short insertions acting as the source of de novo CREs have not been reported previously. Thus, our results provide evidence that short insertions can also create systemic functional silencers by introducing novel binding sites.
Previous studies have shown that transposons containing existing TFBS can facilitate de novo regulatory behavior and increase the transcriptional activity of target genes [9, 35]. In the current study, we found a de novo silencer in a distal melanocyte-specific regulatory region that had been naturally created by short insertions. Moreover, this de novo silencer acted as a negative CRE, blocking the expression of Mitf-m. This novel enhancer feature, introduced by an insertion event, comprised an in silico predicted SOX protein-binding consensus sequence. Our evidence indicated that a transcription factor (predicted to be SOX2 or SOX9) interacts with this novel binding site to repress the expression of MITF-M. Both the 9-bp and the 14-bp insertions were predicted to create new TFBS, and although only the 14-bp insertion interacted directly with nucleoproteins in vitro, we cannot rule out the possibility that the 9-bp insertion is the causative mutation, because the distance between the two insertion regions is as close as 50 bp and the HMG binding domain of SOX family proteins can bind to a sequence as large as 100 bp . Thus, the 9-bp insertion region may act as a flanking sequence, affecting the robustness of the de novo element [37, 38].
Previous studies of the coding and regulatory regions of the Mitf gene have identified various functional mutations associated with congenital hearing loss in mammals [13, 20, 23, 24, 39]. Most of these mutations are located in the exons and introns that are conserved between species, and have been shown to cause equal genetic effects on multiple isoforms. Dysfunction in melanocytes in the SV is the major pathology of WS2. Apart from MITF-M, isoforms MITF-A and MITF-H are also expressed in melanocytes and melanoma cells [40–42]. The particular function(s) of each isoform in hearing loss remains to be fully elucidated. In a previous study of Mitf mi-bw mice, insertion of a retrotransposing L1 element into intron 3 of the MITF gene abolished MITF-M, but also affected the expression of MITF-A and MITF-H . Thus, the phenotypes on hearing loss and pigmentation cannot be exclusively attributable to the elimination of MITF-M, and there remains a possibility that MITF-A or MITF-H may play functional roles in the pathogenesis of WS2 and its manifested symptoms. Our study mapped a novel mutation located in the MITF-M-specific regulatory region that completely eliminated MITF-M expression, but which apparently had no effect on any other isoform (Fig. 3b, c). A search of the literature revealed that mutations in this region have also been shown to affect coat color in dogs  and horses . Furthermore, we generated the first M-exon specific knock-out mouse model, we found phenotypes that were consistent with Mitf r/r pigs, as well as similar Mitf expression profiles. Considering the collective results from the specific naturally-arising porcine mutation and the artificially-induced mouse mutation, it appears that MITF-M exerts a unique function in the inner ear and that dysfunction of the MITF-M isoform alone is sufficient to cause deafness.
Mitf is a well-described and frequent causative gene of WS2 [11, 12]. Most of the mutations found in Mitf are located in the coding region and the flanking splice donor/acceptor sites, accounting for the phenotype of hearing loss in approximately 15 % of WS2 cases; however, the mutation profile has not been determined in approximately 70 % of WS2 patients . Our results suggest that the upstream region of the M-exon, which serves to regulate the expression of Mitf-m, is a novel potential mutation region in WS2 patients. Identifying and characterizing such regulatory regions will have clinical implications for identifying causative mutations for genetic diseases such as WS2. New techniques, such as targeted, capture-based sequencing and whole-genome re-sequencing, provide the possibility of mutation screening for large-scale regulatory regions. Finally, genetic hearing loss is almost exclusively studied in mouse models. The lack of a large animal model for genetic hearing loss has impaired the development of gene therapy for clinical application. Our Mitf r/r pigs should serve as a valuable model to be tested for such therapeutic intervention both by traditional gene therapy and by new CRISPR/Cas9-mediated genome editing .
Study population of Rongchang pigs
Chinese Rongchang pigs were chosen from the animal breeding facility of the Chongqing Academy of Animal Science. Pigs from 1 month to 1.5 years old (male and female) were used. All Rongchang pigs underwent auditory electrophysiology diagnoses to determine the hearing phenotype. Pigs with normal hearing were assigned to the control group and albino pigs with no ABR response were used as the case group. Use of the pigs and mice was approved by the Institutional Animal Care and Use Committee of China Agricultural University in Beijing.
ABR, EP and [K+] recording were used to evaluate hearing function of pigs and mice. Details of these tests have been described in our previous publications . Briefly, ABRs were evoked with clicks and tone pips at 4, 8, 16, and 20 kHz. The ABR threshold was determined by visual inspection. EP and [K+] in the endolymph were recorded using double-barreled microelectrodes filled with 150 mM KCl and a potassium ion exchanger. The reference barrel was filled with 0.5 M NaCl for EP recording and the electrodes were calibrated in a 37 °C water bath chamber using a series of solutions of KCl and NaCl that had [K+] of 1, 5, 25, 50, 100, and 150 mM. At least four independent animals were tested to confirm the results.
The morphology of the stereocilia of cochlear hair cells was examined using SEM. Samples for SEM were prepared following procedures described in our previous publications . Briefly, the cochleae were fixed with 2.5 % glutaraldehyde. After dehydration, the samples were critical-point dried, mounted on aluminum stubs, sputter-coated with gold particles, and examined using a Hitachi S-3700N SEM (Japan). The cochlear SV morphology was observed on cochlear sections using light microscopy and TEM. For TEM, the cochleae after overnight fixation were decalcified, embedded in Epon resin, and sectioned on a Reichert Ultracut E ultramicrotome (Boeckeler Instruments, Tucson, NM). Ultrathin sections were mounted on formvar-coated slot grids, stained with lead citrate and uranyl acetate, and examined using a Philips CM 120 TEM. At least three independent samples were performed to confirm the results.
Genetic analysis and whole genome association study
Eleven pairs of putative heterozygous boars and sows were selected from the normal herd if they had albino offspring. In total, 74 piglets from 11 litters of heterozygous × heterozygous matings were phenotyped to investigate the Mendel segregation ratio. DNA from 28 pigs showing profound hearing loss and 73 pigs with normal hearing was isolated from tissue samples taken from pinnas. The DNA samples were genotyped using the Illumina PorcineSNP60 BeadChip. The PLINK software package was used for genome-wide mapping; all SNP markers with MAF > 0.05 and call rate > 75 % were used for a case-control association analysis. Genome-wide P values were ascertained through phenotype permutation testing (n = 100,000).
All exons, exon-intron boundaries and promoters (10-kb fragment upstream of TSS) were sequenced in 12 normal pigs and 12 albino pigs using the primers listed in Additional file 4: Table S11. Long-range PCR with LongAmpTaq DNA polymerase (NEB) was used to amplify these fragments. PCR fragments were gel-purified with an EZNA Gel Extraction kit (Omega Biotek) and then subjected to Sanger sequencing.
The mRNA expression levels of MITF transcript variants and β-actin were analyzed by a qPCR method using the primers in Additional file 4: Table S11. The cDNA samples (100 ng) and primers for the target genes were mixed with Power SYBR Green PCR Master Mix (Applied Biosystems) in 25-μL final volumes, and amplified using an ABI7900 instrument (Applied Biosystems). All samples were analyzed in triplicate. Protein expression levels of MITF isoforms and tubulin were examined by immunoblotting assays using lysed tissue from MITF R/R and MITF r/r cochleae or skin. MITF antibody C5 (#ab80651, Abcam) was used. The loading amounts were verified by determining levels of a housekeeping protein, β-tublin (Chemicon, Temecula, CA).
Transcriptome analysis of SVs
The lateral walls of the cochlea containing the SVs were dissected from pigs, and preserved with the RNAlater reagent (Ambion, Austin, TX). Total RNA was extracted using the RNeasy Micro Kit (#74004, Qiagen, GRE). RNA concentration and integrity was measured using a Qubit 2.0 Flurometer (Life Technologies, CA, USA) and a Bioanalyzer 2100 system (Agilent Technologies, CA, USA), respectively. The IlluminaTruSeq RNA Sample Preparation Kit (Illumina, San Diego, USA) was used to generate sequencing libraries following the manufacturer’s recommendations. The clustering of samples was performed with the cBot Cluster Generation System according to the manufacturer’s instructions. Clustered libraries were sequenced on a Hiseq 2000 platform (Illumina) and 100-bp paired-end reads were generated. After quality control, clean reads with high quality were obtained for downstream analysis. Genome assembly Sscrofa10.2 was used as the reference genome for read mapping. Clean reads were aligned to the reference genome using TopHat (ver. 2.0.7) [33, 47]. RPKM values (reads per kilobase of exon model per million mapped reads) were calculated to measure the expression level of genes.
Transcription activity analysis
The 7.8-kb fragments of the M-promoter of the R and r alleles were inserted into a luciferase reporter plasmid pGL3-basic (Promega, Madison, USA) to yield the pGL3-R-7.8 k and pGL3-r-7.8 k constructs, respectively. Similarly, truncated constructs including pGL3-R-6.4 k, pGL3-r-6.4 k, pGL3-R-1.2 k, and pGL3-r-1.2 k were generated. Murine B16 melanoma cells with high endogenous MITF-M expression were cultured for plasmid transfection. Luciferase reporter constructs (2.0 μg) and Renilla luciferase vector (0.5 μg) were co-transfected into B16 cells using Lipofectamine 2000 Transfection Reagent according to the manufacturer’s protocol (#11668-019, Invitrogen). At 48 h after transfection, luciferase activity was determined with a Dual-Luciferase Reporter Assay Kit (#E2920, Promega, Madison, USA). The Renilla luciferase activity was used to normalize the transfection efficiency. At least three independent experiments were performed.
Nuclear extracts from B16 cells were prepared using the Nuclear and Cytoplasmic Protein Extraction Kit (#P0027, Beyotime). Oligonucleotides representing the R and r allele fragments (Additional file 4: Table S11) were 3’ end-labeled with biotin and incubated with nuclear extract in the absence or presence of homologous unlabeled DNA (50-fold molar excess). The products were resolved by electrophoresis on an 8 % polyacrylamide gel with × 0.5 TBE at room temperature for 2 h at 150 V. At least three independent replicates were performed to confirm the results.
Mitf mi-ΔM/mi-ΔM -targeted mice were generated using the “recombineering” technology. Briefly, a targeting construct with flanking regions (10,868 bp) of the M-promoter/M-exon was used for standard targeting of mouse embryonic stem cells. Positive cell clones were microinjected into eight-cell embryos to obtain chimeric mice. Chimeric mice that could transmit the modified Mitf m-exon allele to their progeny were crossed with wild-type mice to generate Mitf mi-ΔM/+ mice. Mitf mi-ΔM/mi-ΔM -targeted mice were obtained by Mitf mi-ΔM/+ × Mitf mi-ΔM/+ mating. Genotyping of the mice was performed by PCR using primers flanking the m-exon of Mitf (as shown in Additional file 4: Table S11). Use of the mice was approved by the Institutional Animal Care and Use Committee of the General Hospital of PLA in Beijing.
ABR, auditory brainstem response; CRE, cis-regulatory element; EMSA, electrophoretic mobility shift assays; EP, endolymphatic potential; MITF, microphthalmia-associated transcription factor; SEM, scanning electron microscopy; SNP, single-nucleotide polymorphism; SPL, sound pressure level; SSC 13, Sus scrofa chromosome 13; SV, stria vascularis; TEM, transmission electron microscopy; TFBS, transcription factor binding site; WGD, whole-genome duplication; WS2, Waardenburg syndrome type 2
We thank Dr. David ZZ He at Creighton University and Dr. Zhengyi Chen at Harvard Medical School for critical review of the manuscript. This work was supported by grants the National Basic Research Program of China (973 Program) (#2012CB967900) to SY, the National High Technology Research and Development Program of China (863 Program) (#2013AA102502), the 973 Program (#2011CBA01000) to NL, the 863 Program (#2014AA020510 and #2015SKLAB6-16), the National Transgenic Research Project (#2011ZX08009-001-003), the National Science Foundation of China (#31101701, #81470700, #31372284, #81400472 and #81271082), the China Postdoctoral Science Foundation to LC and the Chongqing Fund of application and development (#cstc2013yykfC80003) to LC.
Availability of data and materials
Supporting data of quantification results are available in Additional file 5. The high-throughput data sets supporting the results of this article are available in the NCBI Sequence Read Archive and Gene Expression Omnibus repository. The Rongchang pig whole-genome re-sequencing data has been deposited at the NCBI Sequence Read Archive, (SRA:SRX397138, SRA:SRX397142, SRA:SRX397141, SRA:SRX397138, SRA:SRX397139, SRA:SRX397137). Raw whole-genome genotype data of the mapping population have been deposited in the Gene Expression Omnibus (GEO:GSE67618). Data of porcine SV transcriptome have been deposited in GEO: GSE80665.
LC, WG, LR, JW, SY, and NL led the experiments. ZG, SC, XB, LZ, XZ, LG, ZL, and LZ performed animal work and prepared biological samples. MY, YZ, YH, BS, JL, and XL performed transcriptional activity analysis. DZ, LR, SZ, ZY, TZ, and JL performed sanger sequencing. ML, ST, QM, QT, YZ, and TC performed the high-throughput sequencing. HY, SZ, KL, NW, TC, YZ, and LZ performed the cochlear morphology experiments. NY, RZ, YG, FZ, HZ, ZH, JL, and WS performed the auditory electrophysiology experiments. JZ, TW and XH mapped the causative gene. LC, WG and WS wrote the paper. JW, SY, MY, ZY, and NL revised the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Additional information is available in the online version of the paper.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Wray GA. The evolutionary significance of cis-regulatory mutations. Nat Rev Genet. 2007;8(3):206–16.View ArticlePubMedGoogle Scholar
- Herz HM, Hu D, Shilatifard A. Enhancer malfunction in cancer. Mol Cell. 2014;53(6):859–66.View ArticlePubMedPubMed CentralGoogle Scholar
- Rebeiz M, Jikomes N, Kassner VA, Carroll SB. Evolutionary origin of a novel gene expression pattern through co-option of the latent activities of existing regulatory sequences. Proc Natl Acad Sci U S A. 2011;108(25):10036–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Frankel N, Erezyilmaz DF, McGregor AP, Wang S, Payre F, Stern DL. Morphological evolution caused by many subtle-effect substitutions in regulatory DNA. Nature. 2011;474(7353):598–603.View ArticlePubMedPubMed CentralGoogle Scholar
- Wittkopp PJ, Kalay G. Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence. Nat Rev Genet. 2012;13(1):59–69.View ArticleGoogle Scholar
- Prabhakar S, Visel A, Akiyama JA, Shoukry M, Lewis KD, Holt A, et al. Human-specific gain of function in a developmental enhancer. Science. 2008;321(5894):1346–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Shirangi TR, Dufour HD, Williams TM, Carroll SB. Rapid evolution of sex pheromone-producing enzyme expression in Drosophila. PLoS Biol. 2009;7(8):e1000168.View ArticlePubMedPubMed CentralGoogle Scholar
- Schlenke TA, Begun DJ. Strong selective sweep associated with a transposon insertion in Drosophila simulans. Proc Natl Acad Sci U S A. 2004;101(6):1626–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Studer A, Zhao Q, Ross-Ibarra J, Doebley J. Identification of a functional transposon insertion in the maize domestication gene tb1. Nat Genet. 2011;43(11):1160–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Eichenlaub MP, Ettwiller L. De novo genesis of enhancers in vertebrates. PLoS Biol. 2011;9(11):e1001188.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang T, Hashimoto K, Bawle EV. Spontaneous contraction of leukodermic patches in Waardenburg syndrome. J Dermatol. 1993;20(11):707–11.View ArticlePubMedGoogle Scholar
- Read AP, Newton VE. Waardenburg syndrome. J Med Genet. 1997;34(8):656–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Pingault V, Ente D, Dastot-Le Moal F, Goossens M, Marlin S, Bondurand N. Review and update of mutations causing Waardenburg syndrome. Hum Mutat. 2010;31(4):391–406.View ArticlePubMedGoogle Scholar
- Levy C, Khaled M, Fisher DE. MITF: master regulator of melanocyte development and melanoma oncogene. Trends Mol Med. 2006;12(9):406–14.View ArticlePubMedGoogle Scholar
- Steingrimsson E, Copeland NG, Jenkins NA. Melanocytes and the microphthalmia transcription factor network. Annu Rev Genet. 2004;38:365–411.View ArticlePubMedGoogle Scholar
- Vachtenheim J, Borovansky J. "Transcription physiology" of pigment formation in melanocytes: central role of MITF. Exp Dermatol. 2010;19(7):617–27.View ArticlePubMedGoogle Scholar
- Li M, Wu H, Luo Z, Xia Y, Guan J, Wang T, et al. An atlas of DNA methylomes in porcine adipose and muscle tissues. Nat Commun. 2012;3:850.View ArticlePubMedPubMed CentralGoogle Scholar
- Steel KP, Barkway C. Another role for melanocytes: their importance for normal stria vascularis development in the mammalian inner ear. Development. 1989;107(3):453–63.PubMedGoogle Scholar
- Koshikawa S, Giorgianni MW, Vaccaro K, Kassner VA, Yoder JH, Werner T, et al. Gain of cis-regulatory activities underlies novel domains of wingless gene expression in Drosophila. Proc Natl Acad Sci U S A. 2015;112(24):7524–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Tassabehji M, Newton VE, Read AP. Waardenburg syndrome type 2 caused by mutations in the human microphthalmia (MITF) gene. Nat Genet. 1994;8(3):251–5.View ArticlePubMedGoogle Scholar
- Karlsson EK, Baranowska I, Wade CM, Salmon Hillbertz NH, Zody MC, Anderson N, et al. Efficient mapping of Mendelian traits in dogs through genome-wide association. Nat Genet. 2007;39(11):1321–8.View ArticlePubMedGoogle Scholar
- Hauswirth R, Haase B, Blatter M, Brooks SA, Burger D, Drogemuller C, et al. Mutations in MITF and PAX3 cause "splashed white" and other white spotting phenotypes in horses. PLoS Genet. 2012;8(4):e1002653.View ArticlePubMedPubMed CentralGoogle Scholar
- Philipp U, Lupp B, Momke S, Stein V, Tipold A, Eule JC, et al. A MITF mutation associated with a dominant white phenotype and bilateral deafness in German Fleckvieh cattle. PLoS One. 2011;6(12):e28857.View ArticlePubMedPubMed CentralGoogle Scholar
- Tachibana M, Kobayashi Y, Matsushima Y. Mouse models for four types of Waardenburg syndrome. Pigm Cell Res. 2003;16(5):448–54.View ArticleGoogle Scholar
- Li M, Tian S, Jin L, Zhou G, Li Y, Zhang Y, et al. Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars. Nat Genet. 2013;45(12):1431–8.View ArticlePubMedGoogle Scholar
- Okumura N, Hayashi T, Sekikawa H, Matsumoto T, Mikawa A, Hamasima N, et al. Sequencing, mapping and nucleotide variation of porcine coat colour genes EDNRB, MYO5A, KITLG, SLC45A2, RAB27A, SILV and MITF. Anim Genet. 2006;37(1):80–2.View ArticlePubMedGoogle Scholar
- Li KK, Goodall J, Goding CR, Liao SK, Wang CH, Lin YC, et al. The melanocyte inducing factor MITF is stably expressed in cell lines from human clear cell sarcoma. Brit J Cancer. 2003;89(6):1072–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Akiyama Y. TFSEARCH: searching transcription factor binding sites. Japan: Real World Computing Partnership; 1995.Google Scholar
- Chew LJ, Gallo V. The Yin and Yang of Sox proteins: activation and repression in development and disease. J Neurosci Res. 2009;87(15):3277–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Uwanogho D, Rex M, Cartwright EJ, Pearl G, Healy C, Scotting PJ, et al. Embryonic expression of the chicken Sox2, Sox3 and Sox11 genes suggests an interactive role in neuronal development. Mech Develop. 1995;49(1-2):23–36.View ArticleGoogle Scholar
- Neves J, Kamaid A, Alsina B, Giraldez F. Differential expression of Sox2 and Sox3 in neuronal and sensory progenitors of the developing inner ear of the chick. J Comp Neurol. 2007;503(4):487–500.View ArticlePubMedGoogle Scholar
- Trowe MO, Shah S, Petry M, Airik R, Schuster-Gossler K, Kist R, et al. Loss of Sox9 in the periotic mesenchyme affects mesenchymal expansion and differentiation, and epithelial morphogenesis during cochlea development in the mouse. Dev Biol. 2010;342(1):51–62.View ArticlePubMedGoogle Scholar
- Arnold CD, Gerlach D, Spies D, Matts JA, Sytnikova YA, Pagani M, et al. Quantitative genome-wide enhancer activity maps for five Drosophila species show functional enhancer conservation and turnover during cis-regulatory evolution. Nat Genet. 2014;46(7):685–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Gompel N, Prud'homme B, Wittkopp PJ, Kassner VA, Carroll SB. Chance caught on the wing: cis-regulatory evolution and the origin of pigment patterns in Drosophila. Nature. 2005;433(7025):481–7.View ArticlePubMedGoogle Scholar
- Chung H, Bogwitz MR, McCart C, Andrianopoulos A, Ffrench-Constant RH, Batterham P, et al. Cis-regulatory elements in the Accord retrotransposon result in tissue-specific expression of the Drosophila melanogaster insecticide resistance gene Cyp6g1. Genetics. 2007;175(3):1071–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Pil PM, Chow CS, Lippard SJ. High-mobility-group 1 protein mediates DNA bending as determined by ring closures. Proc Natl Acad Sci U S A. 1993;90(20):9465–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Goode DK, Callaway HA, Cerda GA, Lewis KE, Elgar G. Minor change, major difference: divergent functions of highly conserved cis-regulatory elements subsequent to whole genome duplication events. Development. 2011;138(5):879–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Erceg J, Saunders TE, Girardot C, Devos DP, Hufnagel L, Furlong EE. Subtle changes in motif positioning cause tissue-specific effects on robustness of an enhancer's activity. PLoS Genet. 2014;10(1):e1004060.View ArticlePubMedPubMed CentralGoogle Scholar
- Markakis MN, Soedring VE, Dantzer V, Christensen K, Anistoroaei R. Association of MITF gene with hearing and pigmentation phenotype in Hedlund white American mink (Neovison vison). J Genet. 2014;93(2):477–81.View ArticlePubMedGoogle Scholar
- Amae S, Fuse N, Yasumoto K, Sato S, Yajima I, Yamamoto H, et al. Identification of a novel isoform of microphthalmia-associated transcription factor that is enriched in retinal pigment epithelium. Biochem Bioph Res Commun. 1998;247(3):710–5.View ArticleGoogle Scholar
- Fuse N, Yasumoto K, Takeda K, Amae S, Yoshizawa M, Udono T, et al. Molecular cloning of cDNA encoding a novel microphthalmia-associated transcription factor isoform with a distinct amino-terminus. J Biochem. 1999;126(6):1043–51.View ArticlePubMedGoogle Scholar
- Vachtenheim J, Novotna H. Expression of genes for microphthalmia isoforms, Pax3 and MSG1, in human melanomas. Cell Mol Biol. 1999;45(7):1075–82.PubMedGoogle Scholar
- Yajima I, Sato S, Kimura T, Yasumoto K, Shibahara S, Goding CR, et al. An L1 element intronic insertion in the black-eyed white (Mitf[mi-bw]) gene: the loss of a single Mitf isoform responsible for the pigmentary defect and inner ear deafness. Hum Mol Genet. 1999;8(8):1431–41.View ArticlePubMedGoogle Scholar
- Zuris JA, Thompson DB, Shu Y, Guilinger JP, Bessen JL, Hu JH, et al. Cationic lipid-mediated delivery of proteins enables efficient protein-based genome editing in vitro and in vivo. Nat Biotechnol. 2015;33(1):73–80.View ArticlePubMedGoogle Scholar
- Guo W, Yi H, Ren L, Chen L, Zhao L, Sun W, et al. The morphology and electrophysiology of the cochlea of the miniature pig. Anat Res. 2015;298(3):494–500.Google Scholar
- Yang SM, Guo WW, Hu YY, Sun YX, Hou ZH, Sun JH, et al. Smad5 haploinsufficiency leads to hair cell and hearing loss. Dev Neurobiol. 2009;69(2-3):153–61.View ArticlePubMedGoogle Scholar
- Carroll SB. Evolution at two levels: on genes and form. PLoS Biol. 2005;3(7):e245.View ArticlePubMedPubMed CentralGoogle Scholar