Transcript Profiles of Stria Vascularis in Models of Waardenburg Syndrome

Background Waardenburg syndrome is an uncommon genetic condition characterized by at least some degree of congenital hearing loss and pigmentation deficiencies. However, the genetic pathway affecting the development of stria vascularis is not fully illustrated. Methods The transcript profile of stria vascularis of Waardenburg syndrome was studied using Mitf-M mutant pig and mice models. Therefore, GO analysis was performed to identify the differential gene expression caused by Mitf-M mutation. Results There were 113 genes in tyrosine metabolism, melanin formation, and ion transportations showed significant changes in pig models and 191 genes in mice models. In addition, there were some spice's specific gene changes in the stria vascularis in the mouse and porcine models. The expression of tight junction-associated genes, including Cadm1, Cldn11, Pcdh1, Pcdh19, and Cdh24 genes, were significantly higher in porcine models compared to mouse models. Vascular-related and ion channel-related genes in the stria vascularis were also shown significantly difference between the two species. The expression of Col2a1, Col3a1, Col11a1, and Col11a2 genes were higher, and the expression of Col8a2, Cd34, and Ncam genes were lower in the porcine models compared to mouse models. Conclusions Our data suggests that there is a significant difference on the gene expression and function between these two models.


Background
Waardenburg syndrome (WS) is a rare genetic condition characterized by at least some degree of congenital hearing loss and pigmentation deficiencies [1]. WS has more than 20 mutations in the Mitf allele [2,3], including Mitf mi-vga9 , Mitf mi-bw , and Mitf mi-ce , which have been identified to cause hearing loss and changes in pigmentation. Although the mouse model is widely used in disease phenotypes and pathogenic mechanisms of deafness-related research [4][5][6][7][8], many shortcomings have also been found in studying human genetic diseases. As there is a tremendous revolutionary difference between mouse and human, it may cause a huge biological difference in anatomy, energy metabolism, and auditory perception [9,10]. For example, the developmental patterns of auditory organs are different in mice and humans: human's hearing developed before birth while mouse's hearing did not fully developed until two weeks after birth [11]. Some studies [4,12] found that human embryonic developmental diseases are difficult to be replicated in some of the mouse models. Therefore, different animal models, such as cattle [13], horses [14], dogs [15,16], and pigs [17] were also necessary to be used to study genetic diseases. Pigs are precocial species with fully developed auditory system at birth. Recent studies also found the cochlear anatomy is very similar to human [18][19][20][21]. As pigs are large-scale animals with high reproductive efficiency and economical convenience, it is a good model for study auditory genetic diseases [22][23][24].
The stria vascularis plays an important role in maintaining the cochlear endolymphatic potentials (EP) which is essential for the mechanical electrical conduction of the hair cells [25][26][27][28][29]. The stria vascularis is composed with the macrophage-like melanocytes, which also called the intermediate cells [28,30,31]. The potassium ions in the Scala media are produced by the intermediate cells and several potassium channels and transporters in the lateral wall, such as KCNQ1/KCNE1, KCNQ4, KCNN2, KCNJ10, and SLC12A2, are also involved in maintaining the endolymphatic potentials [26,[32][33][34]. When malfunction of stria vascularis will result in hearing loss [35]. For example, Marcus [36] reported that Kcnj10 knockout can decrease EP value from +80 mV to 1 mV, and the K + concentration decrease from 110 mM to 60 mM in their mouse models. In our previous studies, we found that Mitf-M knockout can decrease the EP to 18 mV in the mouse model [37]. In the Mitf knockout pig model, we found that the Mitf mutation caused the value of EP dropped from +78 mV to +3 mV, which was lower than the mouse model [18]. In the wild type pig, the potassium concentration in the endolymph was 142 mM higher than those in the perilymph. Our previous study found the potassium concentration dropped to 0 mM in the Mitf mutant pigs [18]. We expect that there may be different genes in maintaining the EP in mice and pigs, and the mutation of Mitf gene may cause a different change in potassium channels. To answer these questions, this study attempted to detect changes in the genetic profiles of these two species caused by Mitf-M gene mutation in RNA transcriptome level. As most of the current researches only use mouse models, this paper will further detect the RNA transcriptome difference in the stria vascularis between the large animals and mouse models.

Gene Expression Changes Caused by Mitf-M Mutation.
The activation of different genes in the stria vascularis of pigs and mice caused by the Mitf-M mutation compared to their W/T controls were screened using the DESeq package software. The conditions for screening differential genes were corrected p value <0.05. The intersections of the differential genes in the pigs and mice were obtained using the Venn diagrams. The results were shown in Table 1. There were 14 common differential genes between the Mitf mutant animals and the controls. There were 177 specific differential genes in mouse model and 99 specific differential genes in pig models ( Figure 1).
The GO analysis and the KEGG pathway analysis were performed on the David Database (adjusted p value <0.05). The results were shown in Figure 2. The main pathway caused by Mitf mutation was the KEGG pathway, enriched in the tyrosine metabolism (mmu00350) and the melanogen-esis pathway (Melanogenesis, mmu04916). The GO analysis was mainly enriched in the biological process of ion transport (ion transport, GO: 0006811) and the integral component of plasma membrane (GO: 0005887) ( Figure 2).

Stria Vascularis Specific Ion Transport-Related Gene
Analysis. Ion transport-related genes were extracted from RNA transcriptome data from the normal and Mitf-m mutant pigs and mice samples for cluster analysis. The results showed many ion transport-related genes were highly expressed in both species through MeV cluster analysis ( Figure 3). The Mitf mutation was coaffected with Trpm1, Kcnj13, and Slc45a2 genes in both species. There were significant differences in ion channel regulation between pigs and mice. The expression of Kcnn1, Clcn2, and Trpm4 genes was higher in pigs than those genes in the mice, whereas the expression of Trpm7, Kcnq1, and Kcnj8 genes was found higher in mice compared to the pigs.

The Specific Tight Junction-Associated Genes in the Stria
Vascularis. The tight junction-associated genes were extracted from the RNA transcriptome data from the Mitf mutant and   Neural Plasticity normal pigs/mice for cluster analysis. The expression of tight junctions in the stria vascularis of the two species was different ( Figure 4). The Cadm1, Cldn11, Pcdh1, Pcdh19, and Cdh24 genes expressed higher in pigs compared to those genes in mice, whereas Ncam, Cldn6, Cldn9, and Cldn14 genes expressed higher in mice compared to pigs. And it was found that both the structures of the stria vascularis in two groups were intact. As the marginal nuclei and the cell connections were intact. The three layers of cells were obvious, and the basal cells were closely connected ( Figure 5).

Stria Vascularis Specific Vascular Development-Related
Genes. Extracted vascular development-related genes from the RNA transcriptome data of the Mitf mutant and normal pigs/mice were used for cluster analysis. There was a    GJB6  ATP1A1  KCNJ16  CLCNKA  KCNJ8  TRPM3  GJB2  KCNQ4  KCNE1  KCNJ10  KCNN4  SLC45A2  KCNJ13  TRPM1  KCNJ2  KCNN1  SLC26A2  CLCN2  ATP1A2  TRPM4  TRPM2  3 Neural Plasticity significant difference in the vascular developmental genes in the stria vascularis between these two species ( Figure 6). The Col2a1, Col3a1, Col11a1, and Col11a2 genes expressed higher in the pigs than the mice, whereas the Col8a2, Cd34, and Ncam genes expressed higher in the mice compared to the pigs.

Discussion
This study reviewed the changes of Mitf-M mutation on gene expression in the cochlea. Mitf has many subtypes [38], in which type M is specifically expressed in melanocytes [39], by direct association with related pigmentases such as tyrosinase (Tyr), dopachrome tautomerase (Dct), endothelin receptor type B (Ednrb), and solute carrier family 45 member 2 (Slc45a2), regulating the survival, migration, and differentiation of melanocytes [40]. Among them, Mitf-M gene [38,39] plays a key role in regulating tyrosine metabolic pathway and melanin production, mainly regulating downstream pigment-related enzymes such as Tyr, Dct, and Tyrp1. Mitf-M also controls cytoskeleton and intercellular tight protein to regulate morphology and migration of melanocytes. In this study, we found the main gene pathway caused by the Mitf-M mutation is on ion transport pathway, including the tyrosine, acid metabolism, and melanin forma-tion pathways, in the cochlear stria vascularis of both mice and pigs. Our data are consistent with previous reports [39].
Our previous studies reported that the Mitf-M gene mutation in the Waardenburg 2A pigs and mice through a deletion of the Mitf-M genes, which caused melanocytes failed to migrate to the cochlear stria vascularis. It can cause drops of EP and damages of cochlear hair cells. In this study, we also identified a significant decrease of the K + channelassociated genes, i.e., Trpm1, Kcnj13, Slc45a2, and Kcnj10. From our RNA-seq sequencing analysis, Clcnka and Kcnj15 genes showed a significant difference in pig models, whereas the Kcnq4, Kcnn4, Kcne1, and Kcnj2 genes are affected mainly in the mouse models. KCNJ13 (KIR7.1) and KCNJ10 (KIR4.1) belong to the inward rectifier potassium channel category. KCNJ10 is known as the key channel of potassium transport. It has been deeply studied in deafness-related diseases, and its deletion can lead to the reduction of EP and potassium ion concentration [36,[41][42][43][44][45]. However, Kcnj13, Trpm1, and Slc45a2 were rarely reported in auditory researches. In addition, TRPM1 is a nonselective voltagegated cation channel in the transient receptor potential (TRP) family, and the Mitf mutation can lead to the deletion of Trpm1 [46]. SLC45A2 is a cross-mediated melanin synthesis in membrane transporter [47], which is regulated by Mitf via the cAMP pathway through Tyr and Dct genes, the major pigment-related genes [48,49].

Neural Plasticity
The stria vascularis transcriptome data of the two species indicated that Mitf-M played an important role in regulating the expression of the Trpm1, Kcnj13, Slc45a2, and Kcnj10 genes in the stria vascularis. In both species, Mitf-M may play an important role in the auditory development and maintain the EP in the cochlea. Although there were huge biological differences between pigs and mice, we found that common gene changes in both species caused by Mitf-M. Mitf-M mutation induced a significant change in Clcn2, Kcnn1, and Trpm4 genes in the both models. CLCN2 [34] is an important component of chloride channel, which coordinates potassium and chloride exchanges. The function of KCNN1 has not been reported in the inner ears. KCNN1 belongs to the calcium ion-mediated potassium channels and plays an important role in the regulation of neural inflammation and nerve aging by microglia [50]. In mice, Kcnq1, Trpm7, and Kcnj8 were significantly affected. KCNQ1 is a calcium ion-dependent potassium channel [45]. When Kcnq1 is deleted, it will cause degeneration of the outer hair cells, which is clinically characterized as Jervell and Lange-Nielsen syndrome, one condition that causes profound hearing loss from birth and a disruption of the heart's normal rhythm. KCNE1 and KCNQ1 are important potassiumsecreting channels in the stria vascularis marginal cells [26,45,51]. KCNE1 regulates KCNQ1 expression and increases ion transport [52].
The tight junctions and vascular endothelial cells are important components of the blood labyrinth barrier as well as ion channels [30,[53][54][55][56][57]. Our cluster analysis of the RNA transcriptome data from both pigs and mice showed that the tight junctions were significantly different in the stria vascularis of these two species. The expression of Cadm1, Cldn11, Pcdh1, Pcdh19, and Cdh24 was found higher in pigs compared to mice, whereas the expression of Ncam, Cldn6, Cldn9, and Cldn14 genes were higher in mice compared to pigs. Cluster analysis of vascular-related genes revealed that it was significantly different in the stria vascularis of the two species. The higher expression in pigs is Col2a1, Col3a1, Col11a1, and Col11a2, whereas the expression of Col8a2, Cd34, and Ncam genes was higher in mice compared to pigs. These results may reveal that the two animals may invoke different genes to regulate the tight junction, just as the ion channels. The differences between the two species' evolutionary relationship, living habits, and anatomy may result in significant differences in these gene expressions [56,58,59]. It is more suitable to choose animal model closer to humans to study auditory related diseases.
In summary, we show that by leveraging RNA-seq for the analysis of the stria vascularis of the WS models, it helps to understand the regulatory mechanisms related to the loss of EP and deafness. These data provide insight into ion channel-defining genes and illustrate the possible genes associated with the WS hearing loss. These results may make a fundamental effect on the gene therapy, which used to rescue the elapse of the endocochlear potential in the stria vascularis.

Conclusion
Our research reveals that there exists a huge difference on the gene expression and function between these two models. According to the different expression in the genetic profiles of these two species caused by Mitf-M gene mutation in RNA transcriptome level, there may be different genes transcript pathway caused by mitf mutation in regulating the potassium channels in mice and pigs. And this transcriptome data may provide a basis for the gene therapy in treating the Waardenburg syndrome.

Animals.
Both Mitf mutant and normal pigs and mice have been used in this experiment. The generation of the Mitf mutant pigs and mice have been described in our previous publications [17,37]. The experimental protocols were approved by the ethics committee of the Chinese PLA Medical School.

RNA Isolation from Stria Vascularis Tissue.
Tissues of the stria vascularis of pigs were obtained from four normal pigs and four Mitf mutant pigs at E85 of embryonic stage as previous study described [17]. The tissues of stria vascularis of mice were obtained from ten normal mice and ten Mitf mutant mice at postnatal 30 days as previous studies described [7,37]. The total RNA of these tissues was extracted separately using Trizol reagent (Invitrogen, CA, USA) following the manufacturer's protocol. The quantity, purity, and integrity of the collected total RNA were analyzed with NanoPhotometer® spectrophotometer (IMPLEN, CA,

Neural Plasticity
USA), a Bioanalyzer 2100, and RNA Nano 6000 Assay Kit (Agilent, CA, USA). Approximately, 4 μg of total RNA was used for the RNA sample preparations as previous studies described [60,61].

Library Construction and Sequencing. The NEBNext®
Ultra TM RNA Library Prep Kit for Illumina® (NEB, USA) was used for the sequencing library preparation, which was conducted with an Illumina HiSeq TM 2000 system following the manufacturer's recommended protocol (Illumina Company Ltd, San Diego, CA, USA) as previous studies described [61,62].

RNA-Seq Reads
Mapping. The reference genome and gene model annotation files were obtained from Genome Web (http://asia.ensembl.org/index.html). The index of the reference genome was built using Hisat2 software (v2.0.5), and the paired-end clean reads were aligned to the reference genome. A database of potential splice junctions was built and confirmed by comparing the previously unmapped reads against the database of putative junctions. The aligned read files were processed by Cufflinks software, which used the normalized RNA-seq fragment counts to measure the relative abundances of the transcriptome. The unit of measure-ment was fragmented per kilobase of exons per million fragments mapped (FPKM). Reads were mapped into the mouse NCBIM38 (ensemble release 68) and the Sus scrofa 11.1 (sus scrofa ensembl release 94) using default options. 5.5. Gene Ontology (GO) and Pathway Enrichment Analysis of DEGs. Differential expression analysis of Mitf-M mutant and normal pigs/mice were performed using the DEseq R package [63]. Using the adjusted p values 0.05 and setting the absolute fold change of 2 as the threshold for significantly differential expression. Using Gene Ontology (GO) and KEGG to analyze high-throughput genome and transcriptome data in the DAVID database [64][65][66][67][68], which is an important online tool for these analyses. The DEGs list was uploaded to the DAVID [64] analysis tool, and p < 0:05 was considered statistically significant. The DEGs was uploaded to the MeV software (https://sourceforge.net/ projects/mev-tm4/) to get the relevant heat map. 5.6. Selecting Deafness Gene of SV Transcriptomes from RNA-Seq Data. In obtaining our data for known deafness genes, we used a database of known deafness genes from these sources: (1) Hereditary Hearing Loss homepage [69] and (2) Hereditary hearing loss and deafness overview [27,45].  Figure 6: Heat analysis of clustering analysis of genes associated with specific vascular endothelial cell in normal and Mitf mutant pigs/mice (hom-mouse and het-mouse represent Mitf-M knockout mice and normal mice, hom-pig and het-pig represent Mitf-M mutant pigs and normal pigs). 6 Neural Plasticity

The Transmission Electron Microscopy (TEM).
To prepare samples for TEM examination, the stria vascularis were washed with 0.1 M PBS and then postfixed in 1% osmium tetroxide and then dehydrated by a series of ethanol before embedded in plastic Agar 100 resin. After polymerization, the stria vascularis was cut into ultrathin sections (3 μm), stained with toluidine blue, were mounted on 0.7% formvar coated copper grids, contrasted by 0.5% uranyl acetate and lead citrate, then examined under a transmission electron microscopy (Philips Tecnai10) [70].

Data Availability
Readers can access additional experimental data in optional supplementary materials.

Conflicts of Interest
The authors declare no competing financial interests.