- Research Article
- Open Access
Copy number variation in the MSRB3 gene enlarges porcine ear size through a mechanism involving miR-584-5p
- Congying Chen†1Email authorView ORCID ID profile,
- Chenlong Liu†1, 2,
- Xinwei Xiong1,
- Shaoming Fang1,
- Hui Yang1,
- Zhiyan Zhang1,
- Jun Ren1,
- Yuanmei Guo1 and
- Lusheng Huang1Email author
© The Author(s) 2018
- Received: 20 April 2018
- Accepted: 18 December 2018
- Published: 27 December 2018
The size and type of ears are important conformation characteristics that distinguish pig breeds. A significant quantitative trait locus (QTL) for ear size has been identified on SSC5 (SSC for Sus scrofa chromosome) but the underlying causative gene and mutation remain unknown. Thus, our aim was to identify the gene responsible for enlarged ears in pig.
First, we narrowed down the QTL region on SSC5 to a 137.85-kb interval that harbors only the methionine sulfoxide reductase B3 (MSRB3) gene. Then, we identified a 38.7-kb copy number variation (CNV) that affects the last two exons of MSRB3 and could be the candidate causative mutation for this QTL. This CNV showed complete concordance with genotype at the QTL of the founder animals in a white Duroc × Erhualian F2 intercross and was found only in pigs from six Chinese indigenous breeds with large ears and from the Landrace breed with half-floppy ears. Moreover, it accounted for the significant association with ear size on SSC5 across the five pig populations tested. eQTL mapping revealed that this CNV was significantly associated with the expression of the microRNA (miRNA) miR-584-5p, which interacts with MSRB3, one of its target genes. In vivo and in vitro experiments confirmed that miR-584-5p inhibits the translation of MSRB3 mRNA. Taken together, these results led us to conclude that presence of the 38.7-kb CNV in the genome of some pig breeds affects ear size by altering the expression of miR-584-5p, which consequently hinders the expression of one of its target genes (e.g. MSRB3).
Our findings shed insight into the underlying mechanism of development of external ears in mammals and contribute to a better understanding of how the presence of CNV can regulate gene expression.
The size and type of ears are important conformation characteristics that distinguish pig breeds. Many Chinese indigenous pig breeds, such as Erhualian and Min, show unusually large and floppy ears. In contrast, Wild boars and European commercial breeds such as Large White and Duroc have small and erect ears. In humans, congenital underdevelopment of the external ear or microtia affects about 0.83 to 17.40 newborns per 10,000 in different ethnic populations . The pig is an important biomedical model  and naturally occurring mutations that affect ear size in pigs offer opportunities to investigate the mechanisms that underlie development of human external ears. To date, quantitative trait loci (QTL) for porcine ear size have been identified on several pig chromosomes (SSC for Sus scrofa) [3–5], of which two major QTL have consistently been detected on SSC5 and 7 in multiple populations [5–7]. Previously, we showed that a missense mutation in the peroxisome proliferator activated receptor delta (PPARD) gene was the causative mutation for the QTL on SSC7 . For the QTL on SSC5, Zhang et al.  analyzed the association of two SNPs within the methionine sulfoxide reductase B3 (MSRB3) gene with ear size in a Large white × Min F2 intercross. We refined the QTL on SSC5 to a 8.7-cM interval and proposed HMGA2 as a candidate gene . Interestingly, the genomic region between 31.74 and 33.78 Mb on SSC5 shows a high level of genetic differentiation between breeds with floppy ears and breeds with prick/partly floppy ears . In addition, this region is orthologous to a QTL for ear type in dogs [11–13] and a strong selection signature within MSRB3 has been reported in large-eared sheep . However, the causative mutation(s) for this QTL remain(s) unexplored in mammals.
In this paper, we report (1) the identification of a 38.7-kb copy number variation (CNV), which involves part of MSRB3, as the suggestive causative mutation for this QTL on SSC5, and (2) the underlying biological mechanism that is responsible for its effect on ear size.
Animals and measurement of the ear phenotype
Five pig populations that included 912 White Duroc × Erhualian F2 pigs (hereafter referred to as F2), 403 Sutai (F19 of a Duroc × Meishan cross), 312 purebred Laiwu, 343 Duroc × (Landrace × Yorkshire) (DLY), and 331 purebred Erhualian pigs were phenotyped and used for genome-wide association studies (GWAS) of ear size. The F2 intercross was constructed and managed as described previously [6, 15]. Briefly, two White Duroc boars and 17 Erhualian sows were mated to produce F1 animals, and then nine F1 sires and 59 F1 dams were randomly intercrossed to generate 1912 F2 individuals. A total of 912 F2 individuals were phenotyped after slaughter at 240 ± 3 days of age. Sutai is a Chinese synthetic breed that was originally produced by intercrossing Chinese Taihu and Western Duroc . In this study, 435 Sutai pigs were generated by crossing four Sutai boars with 55 Sutai sows and slaughtered at 240 ± 3 days of age. Except for the DLY population that was raised on a commercial farm (Xiushui, Jiangxi) until slaughter (90–100 kg), the four other populations were raised on an experimental farm at Jiangxi Agricultural University in Nanchang, Jiangxi.
Distribution of the 38.7-kb CNV across diverse pig breeds with different ear size
Small and erect ears
Large and floppy ears
Chinese wild boar
Small and erect ears
White Duroc (F0 boar)
Half floppy ears
Two methods were used to determine ear size: (1) the ear size of 912 F2, 331 Erhualian, 435 Sutai, 343 DLY and 314 Laiwu pigs, was measured as previously described : after dissection from the slaughtered pig, each ear profile was drawn on a sulphate paper, which was scanned, and ear areas were calculated by the Leica Qwin area measurement software (Leica, USA); and (2) for the 144 Landrace pigs, ear size was measured on live animals by taking photographs : the ears of each animal were fixed and photographed with a ruler placed on the surface of the ear as an internal size reference, and the surface of each ear was calculated using the Leica Qwin software (Leica, USA). In addition, we recorded the body weight of each phenotyped pig to use as covariate in the statistical analysis.
Collection of ear tissue samples for gene expression analysis
Ear tissue was sampled from 26 Sutai two-day old piglets, i.e. 10 CNV/CNV, 10 wild-type, and six CNV/wild-type piglets, and used for gene expression analysis. Of these 26 samples: six (three CNV/CNV and three wild-type piglets) were used for qRT-PCR, RNA sequencing, and miR-584-5p expression, western blot, and immunohistochemistry analyses; six (three CNV/CNV and three wild-type) piglets were used for Western blot and miR-584-5p expression analysis; four (two CNV/CNV and two wild-type) piglets were used for miR-584-5p expression and qRT-PCR analysis; six (CNV/wild-type piglets) were used for western blot analyses only; and the remaining four samples (two CNV/CNV and two wild-type piglets) were used for qRT-PCR analysis only. Furthermore, ear tissue samples from 59 Landrace piglets were used for eQTL mapping of miR-584-5p. Storage of ear tissue samples differed depending on their future use, i.e. for RNA-based expression analysis, samples were stored in RNAlater (Ambion, UK) at − 80 °C until use, for Western blot analysis, samples were immersed in liquid nitrogen and then stored at − 80 °C until use, and for immunohistochemistry analysis, samples were embedded and frozen in OCT freezing medium (Tissue-Tek, Japan).
Genomic DNA was extracted from ear tissue by the standard phenol/chloroform method. DNA quality was evaluated using a Nanodrop-1000 spectrophotometer (Thermo Fisher, USA). The F2, Sutai, Laiwu, Erhualian and DLY populations were genotyped for 62,163 single nucleotide polymorphisms (SNPs) on the porcine SNP 60 K Beadchip (v1.0) (Illumina, USA) according to the manufacturer’s protocol. Quality control of SNP genotype data was performed separately for each population by the check marker function of GenABEL . SNPs with a call rate higher than 95%, a minor allele frequency lower than 1% and a P value for deviation from Hardy–Weinberg equilibrium (HWE) less than 10−6 as well as X-linked SNPs that were likely autosomal (odds > 1000) were removed from further analyses. Animals with a SNP call rate lower than 95% were also excluded.
First, we performed a GWAS of ear size for each of the F2, Sutai, DLY, Laiwu and Erhualian populations. Associations of SNPs with ear size were evaluated using a mixed linear model that included a random polygenic effect with the variance–covariance matrix proportional to genome-wide identity-by-state derived using SNP genotypes  to account for population stratification. Sex and batch were included as fixed effects, and carcass weight was treated as a covariate. The single-marker GWAS was performed using the GenABEL package [17, 19]. Presence of residual population stratification was evaluated by examining the distribution of test statistics generated from thousands of association tests and assessing their deviation from the null distribution in a quantile–quantile (Q–Q) plot , which was created using the R software.
A meta-analysis of the GWAS for the five populations was conducted by summing, for each SNP, the statistical Chi square (χ2) value of the SNP for each experimental population to calculate a new χ2 value with four degrees of freedom. Bonferroni correction was used to adjust for multiple-testing by using a threshold P < 0.05 divided by the number of SNPs used for GWAS analysis as a genome-wide significance level (see Additional file 2: Table S1). Linkage disequilibrium (r2) between SNPs within the QTL region was estimated by PLINK (v1.07) . To determine each confidence interval, a region was defined with all SNPs (r2 > 0.8) in strong LD with the peak SNP. Then, within that region, the confidence interval was defined by the 2-LOD drop-off method, i.e. all SNPs with a lod-score higher than the peak lod score minus 2 were retained . The phenotypic variances explained by the top SNPs were calculated by (Vreduce − Vfull)/Vreduce, where Vfull and Vreduce are residual variances of the models for the association analyses with and without the effect of the top SNP genotype, respectively.
Haplotype sharing analysis
Haplotypes were reconstructed with the dualPHASE software  for all F0, F1 and F2 animals in the F2 population using the SNPs on SSC5. Phenotypic differences between F2 animals that carried different haplotypes were compared by applying the t-test in the lm function of R as described previously . The QTL genotypes of the 19 F0 animals were inferred by integrating the analysis of the segregation of phenotypic values and haplotypes in the F2 pigs. To further identify the shared haplotype within the QTL region in large-eared pig breeds, haplotype sharing analysis was also performed between Erhualian and Min pigs using whole-genome resequencing data, as described in our previous study . We reconstructed the haplotypes within a 0.99 Mb region, which covers the QTL intervals on SSC5 that were detected in the GWAS in four Erhualian and six Min pigs using 118 SNPs from this resequencing data. Wuzhishan (n = 3) and Tibetan (n = 3) pigs with small and erect ears were used as controls using the whole-genome resequencing dataset .
The genomes of the 19 F0 pigs of the F2 population were resequenced at ~ 25 to 30X coverage on an Illumina Hiseq 2500 platform, as described previously . Clean reads were aligned to the Sscrofa11.1 pig reference genome by using the Burrows-Wheeler alignment software . Local realignment around potential insertions/deletions (InDel) and base quality recalibration were conducted with the GATK software . Variants within the QTL interval on SSC5 were called using the HaplotypeCaller algorithm implemented in GATK . Raw SNPs were filtered with the GATK Variant Filtration tool using the following criteria: QD < 2.0, FS > 60.0, MQ < 40.0, a haplotype score > 13.0, mapping quality rank sum < -12.5 and ReadPosRankSum < − 8.0. Raw InDel were filtered by setting the following criteria QD < 2.0, FS > 200.0 and ReadPosRankSum < − 20.0. CNV were called using the depth of coverage method described by Nord et al. .. Briefly, first, sequence coverage was normalized and corrected for capture bias associated with GC-content; then a sliding window was used to identify regions for which the majority of bases had a coverage ratio ≥ 1.4 or ≤ 0.6, which indicates a copy number gain or loss, respectively. A 20-bp window with a minimum of 18 bp meeting the criteria for either gain or loss was defined as a variant window. CNV regions were generated by extending variant windows and merging neighboring variant regions. CNV were examined for presence of partially mapped reads within the predicted breakpoint region.
Determination of the breakpoints of the causative CNV
The breakpoints of the 38.7-kb CNV were determined by PCR and Sanger sequencing with the primers 5′-CGAGGGAGTAAGGCAGACAG-3′ and 5′-GGCTCGGATCATCAGTATCG-3′. The PCR reaction contained 1 × buffer, 2 mM Mg2+, 2 mM dNTP, 200 nM of each primer, 2.5U rTaq polymerase (Takara, Japan), and 50 ng genomic DNA, and thermocycling conditions were as follows: 94 °C for 4 min, 40 cycles (94 °C for 30 s, 68–55 °C touchdown for 30 s, 72 °C for 1 min) and a final step at 72 °C for 10 min. PCR products were run on 1.5% agarose gels and sequenced after purification with the PCR products purification kit (Qiagen, Germany).
The number of 38.7-kb CNV was estimated using both the comparative Ct (2−ΔΔCt) relative quantification (qPCR) and the absolute quantification method (ddPCR for droplet digital PCR). In total, 435 Sutai, 314 Laiwu, 343 DLY, and 144 Landrace pigs were genotyped for the causative CNV by qPCR. Briefly, one pair of primers and a probe (MSRB3-CNV1) (see Additional file 3: Table S2) were designed with the Primer Express 3.0 software (Applied Biosystems, USA). The 20-μL-reaction system contained 1 × Premix Ex Taq (TaKaRa, China), 1 × ROX Reference Dye II (TaKaRa, China), 20 nM FAM-labeled MGB probe (Applied Biosystems, USA), 4 pmol each of the forward and reverse primers and 21.5 ng genomic DNA. The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene was used as an internal control. Real-time PCR was performed in a 7500 fast real-time PCR system instrument (Applied Biosystems, USA). Each sample was analyzed in triplicate. In order to validate the accuracy of genotyping, some samples were randomly drawn and genotyped by another pair of primers and probe using the same reaction system and parameters (MSRB3-CNV2) (see Additional file 3: Table S2).
The number of CNV in 12 samples from six large-eared pig breeds was determined by ddPCR, as described previously [28, 29]. Briefly, genomic DNA was digested with the restriction enzyme BcuI (NEB, UK) that cuts the DNA outside of the amplicons. Primers and TaqMan probes used for ddPCR were the same as those used for the qPCR (MSRB3-CNV1) (see Additional file 3: Table S2). The estrogen receptor 1 (ESR1) gene was used as reference and was labeled with VIC at the 5′ end. The ddPCR reaction mixture was prepared in a 22-μL volume containing 40 ng of digested DNA template, 1 × ddPCR master mix (Bio-Rad, USA), 1.1 μL of each target primer (0.9 μM) and probe (0.25 μM) mix, and 1.1 μL of each reference primer (0.9 μM) and probe (0.25 μM) mix. The disposable cartridge was loaded with the reaction mixture and droplet generation oil and then placed into a QX100 Droplet Generator (Bio-Rad, USA) to acquire droplets. The droplets were run on a Thermal Cycler (Bio-Rad, USA) as follows: 95 °C for 10 min, 40 cycles (94 °C for 30 s, 60 °C for 1 min), 98 °C for 10 min and then kept at 12 °C. After amplification, we loaded the plate on the QX100 Droplet Reader (Bio-Rad, USA) to read the droplets. The Droplet reader processed each sample independently and interrogated both FAM and VIC fluorescence under the guidance set in the Quanta soft software (V126.96.36.199). The software measured the numbers of positive and negative droplets per fluorophore per sample, which were converted into digital signals. The numbers of CNV copies in the tested samples were read directly from the software (Bio-Rad, USA).
Association and interaction analyses
Ear tissue samples from 10 CNV/CNV and 10 wild-type Sutai piglets were used for qRT-PCR analysis (as described above). Total RNA was extracted using the Rneasy Fibrous Tissue mini kit (Qiagen, Germany). The first-strand cDNA was synthesized from 2 μg of total RNA using the Omniscript reverse transcriptase kit (Qiagen, Germany). Expression profiles of the two transcript isoforms of MSRB3 were determined using the SYBR Green method. qRT-PCR was performed in a 20-μL reaction volume containing 50 ng of template cDNA, 1 × Power SYBR Green Master mix (Applied Biosystems, USA) and 5 pM of each primer, with the following thermocycling conditions: 15 min at 95 °C and 48 cycles (95 °C for 15 s, 60 °C for 1 min and 72 °C for 1 min). Then, a dissociation curve analysis was performed. Expression levels of HMGA2 and LEMD3 were determined by the Taqman method. The 20-μL reaction volume included 50 ng of template cDNA, 1 × Premix Ex Taq (TaKaRa, China), 1 × ROX Reference Dye II (TaKaRa, China), 4 pM of each primer and 2 nM of probe, with the following thermocycling conditions:15 min at 95 °C and 40 cycles (95 °C for 15 s and 60 °C for 1 min). β-actin was used as an internal control. Gene-specific primers and probes were designed using the primer premier 5.0 software (see Additional file 4: Table S3). qRT-PCR was performed using the 7500 Fast real-time PCR System (Applied Biosystems, USA). Each sample was analyzed in triplicate.
RNA sequencing was done on ear tissue from six Sutai piglets, including three CNV/CNV and three wild-type animals. DNA contamination was removed from total RNA by incubation with RNase-free DNase I (NEB, UK) at 37 °C for 30 min. The quality of total RNA was assessed by a 2100 Bioanalyzer (Agilent, USA) and 1% agarose gel electrophoresis. RNA sequencing was performed as described previously . Briefly, mRNA was isolated from total RNA with oligo (dT) magnetic beads (Invitrogen, USA) and then fragmented using the RNA fragmentation kit (Ambion, USA). The first-strand cDNA was synthesized using random hexamer primers and reverse transcriptase (Invitrogen, USA), while the second-strand cDNA was synthesized using RNase H (Invitrogen, USA) and DNA polymerase I (NEB, UK). The cDNA libraries were loaded onto flow cell channels of a HiSeq 2000 platform (Illumina, USA) for paired-end 90 bp × 2 sequencing. The average insert size for the paired-end libraries was 200 bp. A paired-end cDNA library was constructed for each of the six samples.
The raw sequence dataset was processed to produce clean sequence data, which were then aligned to the porcine reference genome sequence (Sscrofa 11.1) using TopHat . The gene expression value was determined and normalized using fragments per kilobase of exon per million fragments mapped (FPKM) . The Cuffdiff program in Cufflinks  was used to identify differential expression genes (DEG) based on |log2 (FC)| > 1, where FC is fold change, and q-value (an adjusted P value) < 0.01. Potential alternative splicing transcripts of the MSRB3 gene were estimated in all six sequenced samples using the Cuffdiff program. We refined the gene structure by assembling transcripts with clean reads with the Cufflinks software . The assembled transcripts were compared to gene annotations from the reference genome to detect the extensions at the 5′ and 3′ ends of the corresponding gene annotation. In this study, we focused on the 3′ gene boundary of MSRB3.
Total protein was extracted from ear tissue samples of 18 Sutai piglets (described above) and porcine fetal fibroblast (PFF) cells by using the total protein extraction kit (Applygen, China) and quantified with the BCA Protein Quantification kit (Vazyme, China). Then, the extracted proteins were separated by SDS-PAGE. After transfer, the nitrocellulose membranes were incubated with rabbit anti-MSRB3 primary antibody (orb373913, Biorbyt, UK) at a 1:400 dilution and mouse anti-β-actin primary antibody (A1978, as endogenous control) (Sigma, USA) at a 1:1000 dilution. After incubation with secondary antibodies (goat anti-rabbit IgG peroxidase antibody (A0545) and rabbit anti-mouse IgG peroxidase antibody (A9044) for MSRB3 and β-actin, respectively) that were conjugated with horseradish peroxidase (Sigma, USA), the nitrocellulose membranes were visualized by using a BeyoECL plus kit (Beyotime Biotechnology, China). The images were acquired with the GeneSnap software (SynGene, UK) and quantified by the ImageJ software .
Ear tissue samples from six Sutai piglets were used for immunohistochemistry analysis (described above). All ear tissue samples that had been embedded and frozen in OCT freezing medium (Tissue-Tek, Japan) were cut in sections of 4 μm with a cryostat. The sections were fixed in 4% paraformaldehyde in phosphate buffered saline (PBS) for 10 min at room temperature, and then permeated in PBS containing 0.2% Triton X-100. The sections were blocked for 30 min in PBS containing 10% bovine serum albumin and 0.3% Triton X-100, and then incubated overnight at 4 °C in blocking solution containing the primary antibodies (orb373913, Biorbyt, UK) at a 1:100 dilution. The secondary antibody, goat anti-Rabbit IgG (whole molecule) (F9887, Sigma, USA) conjugated with fluorescein isothiocyanate at a dilution, was used to incubate the slides at room temperature for 45 min. Immunohistochemistry images were captured using a Nikon ECLIPSE Ti inverted microscope equipped with NIS-Elements D software (Nikon, Japan).
Expression profiling analysis of miR-584-5p in vivo and eQTL mapping
The mature miRNA-584-5p sequence was predicted by using the miRBase database (release 20) . Expression profiles of miR-584-5p in ear tissue from eight CNV/CNV and eight wild-type Sutai piglets, along with 59 Landrace piglets were analyzed. Extraction and quality control of total RNA were performed as described above. A stem-loop RT-PCR method  was used to design the primers and probes to quantify the expression level of miR-584-5p using the 2−ΔΔCt method (see Additional file 5: Table S4). The U6 gene was chosen as endogenous control. The qRT-PCR reaction and conditions were the same as those described above.
eQTL mapping was performed for miR-584-5p in ear tissue samples from 59 Landrace two-day old piglets. Thirty-five SNPs within the QTL interval on SSC5 were genotyped by PCR amplification and Sanger sequencing. CNV genotypes were determined by qPCR on an ABI7900HT instrument (Applied Biosystems, USA), as described above. The expression values of miR-584-5p were adjusted for the effects of sex, batch, and kinship using a robust linear regression model . eQTL mapping was performed using a mixed linear model implemented by the mmscore function of GenABEL in R package. Bonferroni correction was applied to adjust for multiple-testing.
In silico prediction of the target genes of miR-584-5p
The DIANA-microT web server (v5.0) was used to predict the target genes of miR-584-5p based on the rule that the miRNA regulate gene expression mainly by base pairing with miRNA-recognition elements (MRE) in the target mRNA 3′UTR or coding sequences (CDS) . The miTG (miRNA targeted genes) score was calculated by microT-CDS using the weighted sum of the scores of all conserved and non-conserved MRE in the 3′UTR and CDS of the target mRNA. The higher the score, the higher the probability of targeting MRE that are within the 3′UTR or in the CDS of the mRNA.
Construction of plasmids and luciferase reporter assay
The 3′UTR of MSRB3 was synthesized and cloned into the vector downstream of the firefly luciferase ORF by cloneEZ (Genscript, China). We obtained a pGL3 luciferase reporter with the mutant 3′UTR of MSRB3 that has a 7-bp deletion in the target site. Porcine fetal fibroblast (PFF) cells were co-transfected with 100 ng of the wild-type or mutant 3′UTR luciferase reporter and 100 nM of the miR-584-5p mimics or the negative control duplexes using lipofectamine 3000 (Invitrogen, USA) in 24-well plates. Each treatment was replicated three times. After transfection for 48 h, cells were harvested by adding 200 μL of passive lysis buffer. Renilla and firefly luciferase activities were assayed on an Infinite 200 PRO multimode reader (Tecan, Switzerland) using the luciferase assay system from Promega. The relative luciferase activity was normalized to the Relina luciferase activity in each sample. We compared the normalized firefly luciferase activity between the wild-type and mutant reporters using the Student’s t test (P < 0.05). This assay was repeated in triplicate.
Transfection of the miR-584-5p mimic to PFF cells
PFF cells (2 × 105 cells seeded in 6-well plates) were cultured in DMEM (with high glucose and pyruvate) (Gibco, USA) containing 10% fetal bovine serum, 50 U/mL penicillin and 50 μg/mL streptomycin (Gibco, USA) and incubated in a humidified incubator with 5% CO2 at 37 °C. When the cells were 70% confluent, the mirVana miRNA Mimic negative control (Ambion, USA) and the miR-584-5p mimic (Ambion, USA, 5′-UUAUGGUUUGCCUGGGACUGAG-3′) were transfected into the PFF cells with 6 μL of 100 nM lipofectamine 3000 (Invitrogen, USA). After 7 h, the transfection mix was replaced by complete medium. Cells were harvested 24 h after transfection. Total proteins were extracted from the cultured cells and used for western blot analysis using the protocol described above.
GWAS in diverse populations confirms the QTL for ear size on SSC5
Haplotype-sharing analysis refined the 509.28-kb QTL interval on SSC5 to 137.85 kb
Chinese Min pigs, such as Erhualian, also have unusually large and floppy ears, and we identified the same QTL for ear size on SSC5 in a Large White × Min F2 resource population . This led us to assume that Min pigs carry the same causative mutation that affects ear size. Thus, we expected the presence of shared haplotypes around the causative mutation between Min and Erhualian pigs. To test this, we performed a haplotype-sharing analysis on 10 large-eared pigs (four Erhualian and six Min pigs) and six small-eared pigs (three Wuzhishan and three Tibetian pigs). We reconstructed the haplotypes of these pigs using the 118 SNPs in the 0.99-Mb QTL interval that contains the 509.28-kb critical region of the SSC5 QTL. As predicted, all large-eared pigs shared two chromosomal segments of 295.03 kb (29.58–29.88 Mb) and 133.31 kb (29.45–29.58 Mb) (Fig. 3c). One of the two segments, i.e. the 295.03-kb (29.58–29.88 Mb) segment overlapped with the 509.28-kb critical region (29.74–30.25 Mb). Thus, the 137.85-kb overlapping region (29.74–29.88 Mb) was the most likely confidence interval for the QTL on SSC5 on which we focused all further analyses. Exons 2–7 of the MSRB3 gene are located within this 137.85-kb interval (Fig. 3d).
Whole-genome sequencing identified a CNV within the MSRB3 gene as a candidate causative mutation
First, we performed a comparative genomic sequence analysis (see “Methods”) of the 137.85-kb QTL interval and found that the reference genome (Sscrofa 11.1) was misassembled for this region (see Additional file 7: Figure S3). We constructed a corrected contig that was submitted to GenBank in NCBI with accession number MK028166. Then, we isolated the full-length mRNA sequence of MSRB3 and identified two transcript isoforms that encode 188 and 183 amino acids, respectively (accession numbers: KX557289 and KX557290). Comparison of these two transcripts showed that the MSRB3 exon 2 was absent in isoform 1 due to alternative splicing (see Additional file 8: Figure S4).
In theory, the causative allele (most likely derived) for the increased ear size should not be detected or occur at a very low frequency in breeds with small and erect ears. We analyzed the allele frequencies of these 64 polymorphisms in 30 individuals from two small and erect-eared breeds (Wuzhishan and Tibetan) and in 20 individuals from two large and floppy-eared breeds (Min and Erhualian) using the dataset described in our previous publication . For all small-eared pigs, the alleles that presumably increases ear size occurred at high frequencies for all 64 polymorphisms (see Additional file 9: Table S5). In contrast, the 38.7-kb CNV was absent from all small-eared pigs and present at high frequencies in large-eared pigs. This provided important supporting evidence that the 38.7-kb CNV is the most likely causative mutation for the QTL on SSC5. Next, we determined the copy number of this CNV in the 19 F0 founder animals of the F2 population (see “Methods”) and found that the genotypes of the founder animals at this CNV were completely concordant with their genotypes at the QTL, i.e. two copies in two White Duroc sires (qDqD), three copies in three Erhualian founder sows (QEqE), and four copies in 14 Erhualian sows (QEQE) (Fig. 4c).
Association of the causative CNV with ear size in four pig populations
Frequency of CNV (%)
Phenotypic variance explained by CNV (%)
Phenotypic variance explained by top SNP (%)
2.20 × 10−16
8.72 × 10−10
2.00 × 10−16
1.37 × 10−15
To obtain additional supporting evidence, we performed association studies between this CNV and ear size in four pig populations, i.e. Sutai, Laiwu, DLY and Landrace. As expected, the CNV was significantly associated with ear size across the Sutai (P = 2.2 × 10−16), Laiwu (P = 8.72 × 10−10), DLY (P = 2.0 × 10−16) and Landrace pigs (P = 1.37 × 10−15) (Table 2). Compared to the top SNPs identified in the GWAS, this CNV explained more phenotypic variance for ear size in Sutai (29.4 vs. 27.4%), Laiwu (13.8 vs. 12.4%) and DLY pigs (33.4 vs. 27.9%). To further evaluate the causality of this CNV, we conducted conditional GWAS in the Sutai, Laiwu and DLY populations by including the genotypes of the CNV as a fixed effect in the GWAS statistical model. This resulted in the disappearance of the previously detected significant association signals on SSC5 across these populations (see Additional file 11: Figure S5). This finding strengthened our hypothesis that the CNV is the causative mutation underlying the QTL for ear size on SSC5.
The levels of MSRB3 in ear tissue differed significantly between small and large-eared pigs
The causative CNV affects the expression level of miR-584-5p, which suppresses the expression of MSRB3
Given that the MSRB3 protein level was lower in CNV/CNV compared to wild-type piglets, we hypothesized that miR-584-5p suppresses the expression of MSRB3. To test this hypothesis, we constructed luciferase reporters that included either the wild-type or the mutant 3′UTR of MSRB3 (Fig. 7d) and co-transfected them with a miR-584-5p mimic or a negative control into porcine fetal fibroblast cells. We found that miR-584-5p significantly reduced the luciferase activity of the wild-type MSRB3 reporter compared with the negative control (P < 0.01). In contrast, we did not observe reduced activity with the mutant luciferase reporter (Fig. 7e), which confirmed that miR-584-5p targets the 3′UTR of MSRB3 directly. We examined the repression of miR-584-5p in the endogenous MSRB3 protein by transfecting miR-584-5p mimics into porcine fetal fibroblast cells. The results of qRT-PCR indicated that the amount of MSRB3 mRNA was not decreased by miR-584-5p compared to the control (Fig. 7f). Western blot analysis revealed that miR-584-5p significantly decreased the MSRB3 protein level (P < 0.01) (Fig. 7g) (see Additional file 14: Figure S8).
Other genes that may affect porcine ear size by combining the predicted target genes of miR-584-5p with the differentially expressed genes
Because miR-584-5p can affect the expression of multiple genes, other than MSRB3, and because other target genes may be involved in porcine ear size, we predicted the target genes of miR-584-5p in silico with the DIANA-microT web server (v5.0) . At a significance threshold for the miTG score ≥ 0.8, we identified 388 predicted target genes (see Additional file 16: Table S8). To identify other possible target genes associated with porcine ear size, we examined the literature for phenotypic consequences in knockout mice for each of these 388 predicted target genes . As a result, we identified five genes associated with ear morphology, including FREM2 which is related with abnormal ear shape , LRP6 with abnormal outer ear morphology , FBXO11  and SIX1  with abnormal middle ear morphology, and FGF10 with abnormal inner ear morphology .
It is well known that miRNA can regulate gene expression by inhibiting target mRNA translation at some stage or by promoting mRNA degradation . To examine whether some of the target genes identified for miR-584-5p were differentially expressed genes (DEG) in the ear tissue between CNV/CNV and wild-type piglets, we isolated a list of 280 DEG with a > 2-fold change at an FDR ≤ 0.001 using RNA sequencing data. These DEG included 82 up-regulated genes and 198 down-regulated genes in the large-eared (CNV/CNV) pigs (see Additional file 17: Table S9). We also examined the literature for phenotypic consequences in knockout mice for these 280 DEG and found abnormal ear morphology in knockout mice for the COL9A2, BMP5, ILDR1, TYRP1 and CTGF genes [47–49]. However, none of these five DEG were target genes of miR-584-5p, and none of the five miR-584-5p target genes described above were DEG. The DEG identified in this study should be located within a transcriptional network involving the miR-584-5p target genes, e.g. MSRB3.
In this study, we have narrowed down the interval containing the QTL affecting ear size on SSC5 to 137.85 kb and identified a 38.7-kb CNV as the causative mutation underlying the QTL. We also suggest that this CNV influences porcine ear size by increasing expression of miR-584-5p, which suppresses the expression of its target genes, e.g. MSRB3. Considering the universality of this QTL for ear type in pig, sheep and dogs, our findings provide important clues for identifying the causative gene or mutations for ear type in mammals, as well as insights about the development of human external ears as a biomedical model.
Previous studies proposed LEMD3, WIF1, and HMGA2 as positional candidate genes for the QTL on SSC5 [5, 7, 10]. Here, fine-mapping analysis allowed us to exclude these three genes as the responsible gene. The use of multiple pig populations allowed us to fine-map the QTL to a small critical interval. In particular, haplotype-sharing analysis in Chinese indigenous pigs with unusually large and floppy ears refined the QTL to a 137.85 kb region. It was previously reported that Chinese indigenous pig breeds show a low level of linkage disequilibrium [50, 51], which facilitates QTL fine-mapping and identification of causative genes for porcine complex traits. In this paper, we provide several solid lines of evidence that support the CNV on SSC5 as a causative mutation: (1) association of this CNV with ear size was strongest in Sutai, Laiwu, DLY and Landrace pigs, which are expected to have a very low level of linkage disequilibrium; (2) the CNV genotypes were completely concordant with the inferred QTL genotypes of all F0 individuals in the F2 population; and (3) this CNV occurred at very high frequencies in Chinese indigenous pig breeds that have large and floppy ears, but was absent in breeds with small and erect-ears. No other mutations within the 137.85-kb region of the QTL on SSC5 showed such concordance patterns in the tested breeds.
CNV are increasingly shown to have causative effects on phenotypes in livestock [52–56]. Here, we provide another example that a CNV has a major QTL effect on porcine ear size. CNV can account for variation in gene expression in a variety of ways, including through a gene dosage effect, disruption of gene coding regions, and deletion or duplication of regulatory elements . In this study, the CNV is located in the 3′ region of miR-584-5p and was shown to be associated with the expression of miR-584-5p expression by eQTL mapping. We propose two possibilities to explain the mechanism underlying the effect of this CNV on ear size: (1) the presence of this CNV could duplicate the regulatory elements (e.g. enhancer) that affect the expression of miR-584-5p but analysis of the CNV sequence did not identify any potential regulatory elements; and (2) because this CNV is quite large, it could modify the chromatin structure in this region and thereby cause a change in the expression pattern of miR-584-5p.
To investigate how the increased expression of miR-584-5p can lead to enlarged ear size in pigs, we predicted the target genes of miR-584-5p in silico. In addition to MSRB3, five of the 388 predicted target genes were associated with abnormal ear morphology in knockout mice. MSRB3 is within the 137.85-kb QTL interval and functional null mutations of MSRB3 have been associated with deafness in both humans  and mouse . Moreover, knockout of MSRB3 in mouse inhibits cell growth through the activation of the p53-p21 and p27 pathways . As a proof of principle, we performed experiments in vivo and in vitro to establish that miR-584-5p suppresses the translation of MSRB3, which subsequently results in enlarged porcine ear size (both immunohistochemistry and western blot analyses showed that expression of MSRB3 is associated with ear size). We did not observe a significant difference in MSRB3 mRNA level between CNV/CNV and wild-type animals. A similar mechanism was also observed in Texel sheep, where a G to A transition in the 3′UTR of the myostatin gene creates an illegitimate microRNA target that causes translational inhibition of the gene, which contributes to muscular hypertrophy . Our findings contribute to a better understanding of how the presence of CNV can regulate gene expression. The effects of other target genes of miR-584-5p on porcine ear size requires further investigation. RNA sequencing analysis identified several DEG related to ear morphology in knockout mice but none of these was a target gene of miR-584-5p. These DEG may be located within a transcriptional network involving the miR-584-5p target genes (e.g. MSRB3).
We show that the 38.7-kb CNV in the 3′ region of the MSRB3 gene is the suggestive causative mutation for a QTL for ear size on SSC5. The presence of this CNV results in an increase in the expression of miR-584-5p, which inhibits the expression of its target gene MSRB3, which affects porcine ear size. As a proof of principle, we confirmed that miR-584-5p hinders the translation of MSRB3 and we found that the MSRB3 protein level differed between pigs with large and small ears. These findings not only shed insight into the underlying mechanisms for development of external ears in mammals but also improve our understanding of the regulation mechanism involving CNV, miRNA, and functional genes.
LH conceived and designed the experiments, and revised the manuscript; CC designed the experiments, analyzed the data, wrote, and revised the manuscript; CL performed the experiments, analyzed the data and wrote the manuscript; XX, SF and HY performed the experiments; ZZ and YG performed the data analysis; JR provided advice on the manuscript. All authors read and approved the final manuscript.
We are grateful to colleagues in the State Key Laboratory of Pig Genetic Improvement and Production Technology, Jiangxi Agricultural University for sample collection and phenotype measurement.
All authors declare that they have no competing interests.
Availability of data and materials
The sequences of the two transcript isoforms of MSRB3 were submitted to NCBI GenBank with Accession Number: KX557289 and KX557290; the RNA-seq clean data are accessible at [GEO: GSE80613]; and the corrected contig in the QTL region was submitted to GenBank with accession number MK028166. The other datasets supporting the conclusions of this article are included within the article and its additional files.
Consent for publication
All authors consent to publication.
Ethics approval and consent to participate
All analyses involving animal samples were conducted according to the care and use guidelines of experimental animals established by the Ministry of Agriculture of China. This study was approved by Animal Care and Use Committee (ACUC) of Jiangxi Agricultural University.
This work was supported by Jiangxi Provincial Program for cultivating young scientist (2010BQ01500).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Suutarla S, Rautio J, Ritvanen A, Ala-Mello S, Jero J, Klockars T. Microtia in Finland: comparison of characteristics in different populations. Int J Pediatr Otorhinolaryngol. 2007;71:1211–7.View ArticleGoogle Scholar
- Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491:393–8.View ArticleGoogle Scholar
- Guo XL, Looft C, Reinsch N, Ernst K. QTL mapping for ear shape based on a commercial pig population. Yi Chuan Xue Bao. 2004;31:819–21.PubMedGoogle Scholar
- Wei WH, de Koning DJ, Penman JC, Finlayson HA, Archibald AL, Haley CS. QTL modulating ear size and erectness in pigs. Anim Genet. 2007;38:222–6.View ArticleGoogle Scholar
- Zhang L, Liang J, Luo W, Liu X, Yan H, Zhao K, et al. Genome-wide scan reveals LEMD3 and WIF1 on SSC5 as the candidates for porcine ear size. PLoS One. 2014;9:e102085.View ArticleGoogle Scholar
- Ma J, Qi W, Ren D, Duan Y, Qiao R, Guo Y, et al. A genome scan for quantitative trait loci affecting three ear traits in a White Duroc × Chinese Erhualian resource population. Anim Genet. 2009;40:463–7.View ArticleGoogle Scholar
- Li P, Xiao S, Wei N, Zhang Z, Huang R, Gu Y, et al. Fine mapping of a QTL for ear size on porcine chromosome 5 and identification of high mobility group AT-hook 2 (HMGA2) as a positional candidate gene. Genet Sel Evol. 2012;44:6.View ArticleGoogle Scholar
- Ren J, Duan Y, Qiao R, Yao F, Zhang Z, Yang B, et al. A missense mutation in PPARD causes a major QTL effect on ear size in pigs. PLoS Genet. 2011;7:e1002043.View ArticleGoogle Scholar
- Zhang Y, Liang J, Zhang L, Wang L, Liu X, Yan H, et al. Porcine methionine sulfoxide reductase B3: molecular cloning, tissue-specific expression profiles, and polymorphisms associated with ear size in Sus scrofa. J Anim Sci Biotechnol. 2015;6:60.View ArticleGoogle Scholar
- Wilkinson S, Lu ZH, Megens HJ, Archibald AL, Haley C, Jackson IJ, et al. Signatures of diversifying selection in European pig breeds. PLoS Genet. 2013;9:e1003453.View ArticleGoogle Scholar
- Jones P, Chase K, Martin A, Davern P, Ostrander EA, Lark KG. Single-nucleotide-polymorphism-based association mapping of dog stereotypes. Genetics. 2008;179:1033–44.View ArticleGoogle Scholar
- Vaysse A, Ratnakumar A, Derrien T, Axelsson E, Rosengren Pielberg G, Sigurdsson S, et al. Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping. PLoS Genet. 2011;7:e1002316.View ArticleGoogle Scholar
- Boyko AR, Quignon P, Li L, Schoenebeck JJ, Degenhardt JD, Lohmueller KE, et al. A simple genetic architecture underlies morphological variation in dogs. PLoS Biol. 2010;8:e1000451.View ArticleGoogle Scholar
- Wei C, Wang H, Liu G, Wu M, Cao J, Liu Z, et al. Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds. BMC Genomics. 2015;16:194.View ArticleGoogle Scholar
- Guo Y, Mao H, Ren J, Yan X, Duan Y, Yang G, et al. A linkage map of the porcine genome from a large-scale White Duroc × Erhualian resource population and evaluation of factors affecting recombination rates. Anim Genet. 2009;40:47–52.View ArticleGoogle Scholar
- Bao WB, Ye L, Pan ZY, Zhu J, Zhu GQ, Huang XG, et al. Beneficial genotype of swine FUT1 gene governing resistance to E. coli F18 is associated with important economic traits. J Genet. 2011;90:315–8.View ArticleGoogle Scholar
- Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007;23:1294–6.View ArticleGoogle Scholar
- Svishcheva GR, Axenovich TI, Belonogova NM, van Duijn CM, Aulchenko YS. Rapid variance components-based method for whole-genome association analysis. Nat Genet. 2012;44:1166–70.View ArticleGoogle Scholar
- Aulchenko YS, de Koning DJ, Haley C. Genomewide rapid association using mixed model and regression: a fast and simple method for genomewide pedigree-based quantitative trait loci association analysis. Genetics. 2007;177:577–85.View ArticleGoogle Scholar
- Pearson TA, Manolio TA. How to interpret a genome-wide association study. JAMA. 2008;299:1335–44.View ArticleGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.View ArticleGoogle Scholar
- Karim L, Takeda H, Lin L, Druet T, Arias JA, Baurain D, et al. Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature. Nat Genet. 2011;43:405–13.View ArticleGoogle Scholar
- Druet T, Georges M. A hidden markov model combining linkage and linkage disequilibrium information for haplotype reconstruction and quantitative trait locus fine mapping. Genetics. 2010;184:789–98.View ArticleGoogle Scholar
- Ai H, Fang X, Yang B, Huang Z, Chen H, Mao L, et al. Adaptation and possible ancient interspecies introgression in pigs identified by whole-genome sequencing. Nat Genet. 2015;47:217–25.View ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010;26:589–95.View ArticleGoogle Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.View ArticleGoogle Scholar
- Nord AS, Lee M, King MC, Walsh T. Accurate and exact CNV identification from targeted high-throughput sequence data. BMC Genomics. 2011;12:184.View ArticleGoogle Scholar
- Hindson BJ, Ness KD, Masquelier DA, Belgrader P, Heredia NJ, Makarewicz AJ, et al. High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Anal Chem. 2011;83:8604–10.View ArticleGoogle Scholar
- Pinheiro LB, Coleman VA, Hindson CM, Herrmann J, Hindson BJ, Bhat S, et al. Evaluation of a droplet digital polymerase chain reaction format for DNA copy number quantification. Anal Chem. 2012;84:1003–11.View ArticleGoogle Scholar
- Chen C, Ai H, Ren J, Li W, Li P, Qiao R, et al. A global view of porcine transcriptome in three tissues from a full-sib pair with extreme phenotypes in growth and fat deposition by paired-end RNA sequencing. BMC Genomics. 2011;12:448.View ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.View ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.View ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.View ArticleGoogle Scholar
- Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011;27:2325–9.View ArticleGoogle Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.View ArticleGoogle Scholar
- Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:D68–73.View ArticleGoogle Scholar
- Salone V, Rederstorff M. Stem-loop RT-PCR based quantification of small non-coding RNAs. Methods Mol Biol. 2015;1296:103–8.View ArticleGoogle Scholar
- Lourenco VM, Pires AM, Kirst M. Robust linear regression methods in association studies. Bioinformatics. 2011;27:815–21.View ArticleGoogle Scholar
- Paraskevopoulou MD, Georgakilas G, Kostoulas N, Vlachos IS, Vergoulis T, Reczko M, et al. DIANA-microT web server v5.0: service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013;41:W169–73.View ArticleGoogle Scholar
- Georgiopoulos G, Evangelou E. Power considerations for lambda inflation factor in meta-analyses of genome-wide association studies. Genet Res (Camb). 2016;98:e9.View ArticleGoogle Scholar
- Anthwal N, Thompson H. The development of the mammalian outer and middle ear. J Anat. 2016;228:217–32.View ArticleGoogle Scholar
- The mouse genome informatics. http://www.informatics.jax.org/. Accessed 5 June 2017.
- Hardisty RE, Erven A, Logan K, Morse S, Guionaud S, Sancho-Oliver S, et al. The deaf mouse mutant Jeff (Jf) is a single gene model of otitis media. J Assoc Res Otolaryngol. 2003;4:130–8.View ArticleGoogle Scholar
- Zheng W, Huang L, Wei ZB, Silvius D, Tang B, Xu PX. The role of Six1 in mammalian auditory system development. Development. 2003;130:3989–4000.View ArticleGoogle Scholar
- Pauley S, Wright TJ, Pirvola U, Ornitz D, Beisel K, Fritzsch B. Expression and function of FGF10 in mammalian inner ear development. Dev Dyn. 2003;227:203–15.View ArticleGoogle Scholar
- Jackson RJ, Standart N. How do microRNAs regulate gene expression? Sci STKE. 2007;2007:re1.View ArticleGoogle Scholar
- Vogelstein B, Kinzler KW. Cancer genes and the pathways they control. Nat Med. 2004;10:789–99.View ArticleGoogle Scholar
- Morozko EL, Nishio A, Ingham NJ, Chandra R, Fitzgerald T, Martelletti E, et al. ILDR1 null mice, a model of human deafness DFNB42, show structural aberrations of tricellular tight junctions and degeneration of auditory hair cells. Hum Mol Genet. 2015;24:609–24.View ArticleGoogle Scholar
- Doherty HE, Kim HS, Hiller S, Sulik KK, Maeda N. A mouse strain where basal connective tissue growth factor gene expression can be switched from low to high. PLoS One. 2010;5:e12909.View ArticleGoogle Scholar
- Amaral AJ, Megens HJ, Crooijmans RP, Heuven HC, Groenen MA. Linkage disequilibrium decay and haplotype block structure in the pig. Genetics. 2008;179:569–79.View ArticleGoogle Scholar
- Ai H, Huang L, Ren J. Genetic diversity, linkage disequilibrium and selection signatures in chinese and Western pigs revealed by genome-wide SNP markers. PLoS One. 2013;8:e56001.View ArticleGoogle Scholar
- Wright D, Boije H, Meadows JR, Bed’hom B, Gourichon D, Vieaud A, et al. Copy number variation in intron 1 of SOX5 causes the Pea-comb phenotype in chickens. PLoS Genet. 2009;5:e1000512.View ArticleGoogle Scholar
- Rosengren Pielberg G, Golovko A, Sundstrom E, Curik I, Lennartsson J, Seltenhammer MH, et al. A cis-acting regulatory mutation causes premature hair graying and susceptibility to melanoma in the horse. Nat Genet. 2008;40:1004–9.View ArticleGoogle Scholar
- Salmon Hillbertz NH, Isaksson M, Karlsson EK, Hellmen E, Pielberg GR, Savolainen P, et al. Duplication of FGF3, FGF4, FGF19 and ORAOV1 causes hair ridge and predisposition to dermoid sinus in Ridgeback dogs. Nat Genet. 2007;39:1318–20.View ArticleGoogle Scholar
- Fontanesi L, Beretti F, Riggio V, Gomez Gonzalez E, Dall’Olio S, Davoli R, et al. Copy number variation and missense mutations of the agouti signaling protein (ASIP) gene in goat breeds with different coat colors. Cytogenet Genome Res. 2009;126:333–47.View ArticleGoogle Scholar
- Johansson A, Pielberg G, Andersson L, Edfors-Lilja I. Polymorphism at the porcine dominant white/KIT locus influence coat colour and peripheral blood cell measures. Anim Genet. 2005;36:288–96.View ArticleGoogle Scholar
- Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007;315:848–53.View ArticleGoogle Scholar
- Ahmed ZM, Yousaf R, Lee BC, Khan SN, Lee S, Lee K, et al. Functional null mutations of MSRB3 encoding methionine sulfoxide reductase are associated with human deafness DFNB74. Am J Hum Genet. 2011;88:19–29.View ArticleGoogle Scholar
- Kwon TJ, Cho HJ, Kim UK, Lee E, Oh SK, Bok J, et al. Methionine sulfoxide reductase B3 deficiency causes hearing loss due to stereocilia degeneration and apoptotic cell death in cochlear hair cells. Hum Mol Genet. 2014;23:1591–601.View ArticleGoogle Scholar
- Lee E, Kwak GH, Kamble K, Kim HY. Methionine sulfoxide reductase B3 deficiency inhibits cell growth through the activation of p53-p21 and p27 pathways. Arch Biochem Biophys. 2014;547:1–5.View ArticleGoogle Scholar
- Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibe B, et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat Genet. 2006;38:813–8.View ArticleGoogle Scholar