Front. Genet. Frontiers in Genetics Front. Genet. 1664-8021 Frontiers Media S.A. 735795 10.3389/fgene.2021.735795 Genetics Original Research Comparative Genome and Transcriptome Integration Studies Reveal the Mechanism of Pectoral Muscle Development and Function in Pigeons Hou et al. Pigeon Selection Signal and Transcriptome Hou Haobin 1 2 † Wang Xiaoliang 1 2 † Yang Changsuo 1 2 Cai Xia 1 2 Lv Wenwei 2 Tu Yingying 2 Bao Aodungerile 3 Wu Quanli 3 Zhao Weimin 3 Yao Junfeng 1 2 * Ding Weixing 1 2 * 1 Shanghai Academy of Agricultural Sciences, Shanghai, China 2 National Poultry Engineer Research Center, Shanghai, China 3 Shanghai Jinhuang Pigeon Company, Shanghai, China

Edited by: Jiuzhou Song, University of Maryland, United States

Reviewed by: Clint Rhode, Stellenbosch University, South Africa

Long Jin, Sichuan Agricultural University, China

 *Correspondence: Junfeng Yao, yaobison1983@163.com; Weixing Ding, weixingding@163.com

These authors have contributed equally to this work

This article was submitted to Livestock Genomics, a section of the journal Frontiers in Genetics

 21 12 2021 2021 12 735795 03 07 2021 26 11 2021 Copyright © 2021 Hou, Wang, Yang, Cai, Lv, Tu, Bao, Wu, Zhao, Yao and Ding. 2021 Hou, Wang, Yang, Cai, Lv, Tu, Bao, Wu, Zhao, Yao and Ding

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Pigeon breed resources provide a genetic model for the study of phenomics. The pectoral muscles play a key role for the meat production performance of the meat pigeon and the athletic ability of the High flyers. Euro-pigeons and Silver King pigeons are commercial varieties that exhibit good meat production performance. In contrast to the domestication direction of meat pigeons, the traditional Chinese ornamental pigeon breed, High flyers, has a small and light body. Here, we investigate the molecular mechanism of the pectoral muscle development and function of pigeons using whole-genome and RNA sequencing data. The selective sweep analysis (F ST and log2 (θπ ratio)) revealed 293 and 403 positive selection genes in Euro-pigeons and Silver King, respectively, of which 65 genes were shared. With the Silver King and Euro-pigeon as the control group, the High flyers were selected for 427 and 566 genes respectively. There were 673 differentially expressed genes in the breast muscle transcriptome between the commercial meat pigeons and ornamental pigeons. Pigeon genome selection signal combined with the breast muscle transcriptome revealed that six genes (SLC16A10, S100B, SYNE1, HECW2, CASQ2 and LOC110363470) from commercial varieties of pigeons and five genes (INSC, CALCB, ZBTB21, B2M and LOC110356506) from Chinese traditional ornamental pigeons were positively selected which were involved in pathways related to muscle development and function. This study provides new insights into the selection of different directions and the genetic mechanism related to muscle development in pigeons.

 pigeons pectoral muscle selective sweep transcriptome candidate gene Introduction

Archaeological evidence indicates that pigeons were historically consumed as food for thousands of years (Blasco et al., 2014) and represented an important protein source for humans. In Europe, North America, and Asia, squab is considered a delicacy and is very popular among consumers (Kokoszyński et al., 2020). The pigeon industry gradually emerged in China during the early 1970s. After nearly 40 years of development, the number of pigeons in stock, out of stock, and the total production of pigeons ranked in China first in the world. According to the 2018 China Poultry Industry Development Report, there were 256,000 pairs of grandparent breeding pigeons, 41.2 million pairs of parent breeding pigeons, an increase of 5% over 2017, and an annual production of 643 million squabs (China Animal Agriculture Association, 2018). Similar to other types of livestock and poultry, breeders tend to cultivate large sized, full-breasted pigeons with high fecundity. There are substantial differences in body shape among different pigeon breeds. The largest pigeon can reach 1,000 g and small pigeons can reach to 250 g, which is nearly a four-fold difference in body mass (Stringham et al., 2012). For example, Euro-pigeon and Silver King weigh over 600 g at 4 weeks of age (Pomianowski et al., 2009), whereas some local varieties only weigh 250 g (Zhang et al., 2020). In contrast to meat pigeons, some performance and carrier pigeons are selected for athletic capability, and these pigeons are generally highly adept at flying. Carrier pigeons have the capacity for long-distance exploration and are known as messengers (Biro et al., 2007), whereas flipping pigeons are recognized as air dancers. High flying pigeons are more famous for their high altitude flying and also have important ornamental traits.

Most studies investigating quality traits using pigeon genomics have focused on the feather crown (Shapiro et al., 2013), feather color (Vickrey et al., 2018), and foot feathers (Domyan et al., 2016; Boer et al., 2019). Moreover, analyses of complex traits have mainly focused on competition ability (Shao et al., 2020). Using genomics and transcriptomics, researchers have revealed that the CASK gene is the key gene involved in the homing ability of carrier pigeons (Gazda et al., 2018). The relationship between candidate gene polymorphisms and meat quality traits has also been studied (Ye et al., 2016; Mao et al., 2018; Dong et al., 2019), and the physical and chemical characteristics of pigeon meat were previously reported (Kokoszyński et al., 2020). Many studies have found that the MSTN gene is related to muscle development in cattle (Berry et al., 2002; Miretti et al., 2013), pigs (Stinckens et al., 2008; Li et al., 2020), chickens (Guernec et al., 2004; Dushyanth et al., 2016), and sheep (Grochowska et al., 2019). In addition, MSTN gene expression in the breast muscle of pigeons was found to be significantly higher than that in other tissues, and to increase with age (Liu et al., 2019).

The genetic mechanism of body weight and size traits in pigs (Gusev et al., 2016), horses (Srikanth et al., 2019), and ducks (Wang L. et al., 2017) was revealed using multi-omics; however, few genome-wide studies have been conducted on gene mapping of pigeon growth traits, particularly the development of pectoral muscles. Pigeon breast muscle accounted for approximately 30% of the slaughter weight (Kokoszyński et al., 2020). Moreover, the growth and development of breast muscle is extremely important for meat performance and the flying ability of pigeons. In this study, we selected three breeds and resequenced a total of 23 individuals, including eight Euro-pigeons (EU), eight Silver King pigeons (SK), and seven High flyers. European meat pigeons and SK pigeons are larger and have better meat performance. HF are good at flying at high altitudes, with the characteristics of a small size and light posture. The distinct purposes of the pigeons diverged these breeds via human selection. The selection signal revealed genes that were positively selected among each of the breeds. Next, 28-day-old pigeon breast muscle transcriptome differences combined with the selection signal revealed the key genes required for pigeon muscle growth. This study is of great significance for the discovery of candidate genes that affect pigeon growth and development traits, and provides new insight into muscle development and function of altricial birds using functional omics research.

 Materials and Methods Pigeons and Ethics Approval

The whole genome of 23 pigeons, including eight European meat pigeons, eight SK pigeons, and seven HF pigeons were resequenced. All pigeons are from Shanghai Jinhuang Pigeon Industry Co., Ltd. After blood collection from the pterygoid vein, DNA was extracted using a Tiangen kit, and the DNA concentration was detected with a Nanodrop Spectrophotometer 2000. An Illumina Hiseq PE150 sequencing platform was used for sequencing. Pigeons used in this study were approved by the Ethics and Animal Welfare Committee of Shanghai Academy of Agricultural Sciences (No. SAASPZ0521012).

 Comparison of Reference Genomes

The effective high quality sequencing data were compared to the reference genome (reference genome download link: https://www.ncbi.nlm.nih.gov/genome/10719?genome_ assembly_ id = 39619) using BWA (Li and Durbin 2009) software (parameter: mem-t 4-K 32-m), and comparison results were removed with SAMTOOLS (Li et al., 2009) (parameter: rmdup).

 SNP Detection and Annotation

We used SAMTOOLS (Li et al., 2009) to detect population SNPs, and ANNOVAR (Wang et al., 2010) software was used for SNP annotation. The Bayesian model was used to detect polymorphic loci in the population, and quality control SNPs were obtained through the following filtering and screening method: 1) Q20 quality control (SNPs with a quality value of Q20 (i.e., the sequencing error rate was greater than 1% were filtered out); 2) the SNPs were at least five bp apart from each other (Since the probability of two SNPs being so close is extremely low, it is considered to be due to errors in sequencing, experimental factors, or analysis, and the two SNPs are removed); 3) the support number (coverage depth) of the SNP was between [1/3, five] times of the average depth.

 Analysis of Pigeon Population Structure

We used neighbor joining methods to construct the cluster tree. After SNP detection, the individual SNPs can be used to calculate the distance between populations. TreeBeST-1.9.2 software (https://mybiosoftware.com/treebest-1-9-2-softwares-phylogenetic-trees.html) was used to calculate the distance matrix. Using the distance matrix, a cluster tree was constructed using the neighbor joining method. The bootstrap values were calculated 1,000 times. We used EIGENSOFT (v5.0; https://www.hsph.harvard.edu/alkes-price/software/) for principal component analysis (PCA) on an individual scale for the 23 pigeons.

 Analysis of Selection Signal Based on FST and θπ

We calculated the genome-wide distribution of F ST values (Weir and Cockerham 1984) and θπ ratios among the threepigeon breeds using a sliding-window approach (40-kb windows with 20-kb increments). The θπ ratios were log2(θπ ratio) transformed. We considered the windows with the top 5% values for the F ST and log2(θπ ratio) simultaneously as candidate outliers under strong selective sweeps. Furthermore, the overlap information of a selected signal was obtained using the “vennDiagram” package in R (https://www.omicstudio.cn/tool/6).

 Comparative Transcriptome Analysis of the Pectoral Muscle

Twelve healthy 28-day-old female pigeons were used as research materials. The breeding method was caged and free-feeding. The breast muscles of four EU meat pigeons, four SK pigeons, and four HF pigeons were collected for transcriptome analysis, and the slaughter traits were recorded (Table 1). Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) in accordance with the manufacturer’s procedure. The amount and purity of the RNA in each sample was quantified using a NanoDrop ND-1000 (NanoDrop, Wilmington, DE, United States). RNA integrity was assessed using a Bioanalyzer 2100 (Agilent, CA, United States) with an RIN number >7.0, and confirmed by electrophoresis with denaturing agarose gel. Finally, we performed 2 × 150 bp paired-end sequencing (PE150) with an Illumina Novaseq™ 6,000 (LC-Bio Technology Co., Ltd., Hangzhou, China) following the vendor’s recommended protocol.

Slaughter traits of three pigeon breeds.

Breeds Live weight (g) Slaughter weight (g) Eviscerated weight (g) Pectorales weight (g)
Euro-pigeon 599.25 ± 64.20a 523.00 ± 65.56a 404.50 ± 50.49a 119.5 ± 19.21a
Silver King 482.75 ± 34.54b 423.50 ± 34.38b 331.50 ± 21.27b 86.00 ± 5.89b
High flyer 311.00 ± 32.76c 262.50 ± 32.80c 177.00 ± 25.53c 64.95 ± 9.16c

The data in the same column are superscripted with different lowercase letters (a, b, and c), indicating significant differences (p < 0.05).

 Sequence and Primary Analysis

Cutadapt software (https://cutadapt.readthedocs.io/en/stable/, version: cutadapt-1.9) was used to remove the reads that contained adaptor contamination (command line: ∼cutadapt -a ADAPT1 -A ADAPT2 -o out1. fastq -p out2. fastq in1. fastq in2. fastq -O 5 -m 100). After removing the low quality and undetermined bases, HISAT2 software (https://daehwankimlab.github.io/hisat2/, version: hisat2-2.0.4) was used to map reads to the Cliv_1.0 rock pigeon reference genome (command line: ∼hisat2 -1 R1. fastq.gz -2 R1. fastq.gz -S sample_mapped.sam). The mapped reads of each sample were assembled using StringTie (http://ccb.jhu.edu/software/stringtie/, version: stringtie-1.3.4d. Linux_x86_64) with default parameters (command line: ∼stringtie -p 4 -G genome. gtf -o output. gtf -l sample input. bam). The transcriptomes from all of the samples were merged to reconstruct a comprehensive transcriptome using GffCompare software (http://ccb.jhu.edu/software/stringtie/gffcompare.shtml, version: gffcompare-0.9.8. Linux_x86_64). After the final transcriptome was generated, StringTie and Ballgown (http://www.bioconductor.org/packages/release/bioc/html/ballgown.html) were used to estimate the levels of transcript expression and determine the mRNA expression level by calculating FPKM (FPKM = [total_exon_fragments/mapped_reads (millions) × exon_length (kB)]) (command line: ∼stringtie -e -B -p 4 -G merged. gtf -o samples. gtf samples. bam). The differentially expressed mRNAs with a fold-change > 2 or fold-change < 0.5, and p-value < 0.05 were selected using DESeq2(Pertea et al., 2015) (http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html).

 Gene Ontology Enrichment and Pathway Analysis

The function of the differentially expressed genes in different types of pigeon breast muscle was investigated by comparing the following groups: EU vs SK, EU vs HF, and SK vs HF. Goseq was used to perform the enrichment analysis of the GO term (Young et al., 2010). Based on the KEGG database, KOBAS 2.0 software was used to analyze the pathways associated with the differentially expressed genes (Kanehisa et al., 2008). Ggplot two was used to analyze the enrichment of the GO and KEGG databases, and the results were displayed as a scatter plot (bubble chart).

 Real-Time Quantitative PCR Analysis

Using 12 RNA samples isolated from pigeons (four EU pigeons, four SK pigeons, and four HF pigeons), the RNA OD value was detected with an ultra-micro nucleic acid protein analyzer (scandrop100). SYBR Green I was used to detect the expression of target genes in the samples using the A260/A280 ratio, β-actin as an internal reference gene, and four target genes. The fluorescence quantitative PCR program and system was performed as follows: Step 1: 95°C for 3 min; Step 2: 95°C for 10 s; Step 3: 60°C for 30 s + plate read; repeat step 2 for 39 cycles. Melt curve analysis: (60°C–95 °C, +1°C/cycle, holding time of 4 s).

 Results

A total of 23 samples were sequenced from different pigeons, with high quality data volume of 199.3 Gb, high sequencing quality (Q20 ≥ 96.0%, Q30 ≥ 90.2%), and normal GC distribution. None of the 23 samples were contaminated (Supplementary Table S1). The library was successfully constructed and sequenced. The average map rate of the population samples was 97.6%, the average sequencing depth of the genome (excluding gap regions) was 7.30 (only reads with a comparison quality >0 are considered), and the average coverage was 98.8% (at least one base is covered) (Supplementary Table S2). A total of 5,673,290 SNPs were detected (Supplementary Table S3).

Population Structure

Based on the degree of SNP differences among individuals of different pigeon breeds, PCA analysis revealed that EU and SK had a closer genetic distance, while Chinese traditional ornamental pigeon and commercial meat pigeon had a longer genetic distance (Figure 1A). The cluster tree analysis revealed that commercial meat pigeon breeds and HF pigeons formed two independent branches, indicating that there was a large genetic distance between them. Although there are obvious differences in feather color between EU pigeons and SK pigeons, the genetic distance between them was relatively close (Figure 1B).

(A) The 3D PCA plot of pigeon population. EU (Euro-pigeon), SK (Silver King), HF (High flyers). (B) Cluster tree of three pigeon breeds.

 Positive Selection of Meat Pigeon and High Flyers

By analyzing selection signals, 293 genes were positively selected in the EU population and 403 genes were positively selected in the SK population. There were 65 overlapping genes (Supplementary Table S4; Figures 2A,B). When EU was used as a control group, 566 genes were selected by HF pigeons. When SK was used as the control group, 427 genes were selected by HF pigeons. There were 262 overlapping genes (Supplementary Table S4; Figures 2C,D).

Distribution of log2 values (θ π •control/θ π •selected) and the top 5% highest F ST values calculated in 40-kb sliding windows with 20-kb increments. (A) Compared with HF, EU are subject to positive selection of genomic regions. (B) Compared with HF, SK are subject to positive selection of genomic regions. (C) Compared with EU, HF are subject to positive selection of genomic regions. (D) Compared with SK, HF are subject to positive selection of genomic regions.

 Comparative Transcriptome of the Pectoral Muscle

A total of 12 breast muscle samples (4 EUs, 4 SKs and 4 HFs) of three pigeon breeds were sequenced according to the standard operation, and a total of 94.73 G of raw data was obtained. After performing quality control, 92.48 G of clean data were obtained, 97.62% of which were effective reads. The proportion of bases with quality values ≥20 (sequencing error rate less than 0.01) was 99.97%, the proportion of bases with quality values ≥30 (sequencing error rate less than 0.001) was 98.11%, and the proportion of GC content was 48.42% (Supplementary Table S5). Hisat was used to compare the reference genome of valid data after preprocessing, revealing a comparison rate of 90.86% (Supplementary Table S6). The percentage of exon annotated sequences was 91.63%, whereas the percentage of intron and intergenic reads was 4.78 and 3.59%, respectively (Supplementary Figure S1).

 Analysis of Differentially Expressed Genes (DEGs)

The expression profiles of 26,640 genes were obtained using Ballgown package to provide file input for fpkm quantification. The genes with FC > 2 times or FC < 0.5 times and p value <0.05 were defined as differentially expressed genes. A total of 1,016 differentially expressed genes were obtained by comparing the transcriptome of breast muscle samples from EU and HF pigeons, including 408 up-regulated genes and 608 down-regulated genes (Supplementary Table S7; Figure 3A). There were 1,294 differentially expressed genes between SK and HF pigeons, including 465 up-regulated genes and 829 down-regulated genes (Supplementary Table S8; Figure 3B). A total of 322 genes were identified in EU vs SK group, of which 195 were up-regulated and 127 were down regulated (Supplementary Table S9).

(A) Volcano plot of EU. VS. HF breast muscle differential gene expression level. (B) Volcano plot of SK. VS. HF breast muscle differential gene expression level. (C) Venn diagram of the commercial pigeon positive selection signal combined with DEGs analysis of the breast muscle transcriptome. (D) Venn diagram of the HF positive selection signal combined with DEGs analysis of the breast muscle transcriptome.

The results showed that there were more differences in breast muscle gene expression between meat pigeons and HF; however, there were fewer differences between meat pigeons, which was consistent with the phenotypic analysis results. Further analysis with a Wayne diagram showed that there were 673 overlapping genes in the EU vs. HF and SK vs HF groups (Figure 3). Pigeon genome selection signal combined with the breast muscle transcriptome revealed that six genes from commercial varieties of pigeons (Figure 3C) and five genes from Chinese traditional ornamental pigeons (Figure 3D) were positively selected.

 GO and KEGG Enrichment Analysis

A total of 305 significant GO terms were identified in the EU vs. HF group, the top 20 of which were related to muscle function, including actin binding (GO:0003779, p = 0.0003) and myosin filament (GO:0032982, p = 0.0003) (Supplementary Table S10; Figure 4A). Moreover, the insulin-like growth factor binding (GO:0005520, p = 0.00003) pathway, which is related to growth and development, was also enriched. For SK and HF, the differentially expressed genes were enriched in 358 significant pathways (Supplementary Table S11; Figure 4B). Similarly, the first 20 pathways included actin binding (GO:0003779, p = 0.000) and myosin filament (GO:0032982, p = 0.00009). In addition, muscle contraction (GO:0006936) was also related to muscle function. Although 516 GO terms were enriched in the EU vs. SK, the top 20 GO terms were not related to muscle development or function (Supplementary Table S12).

Top 20 significantly enriched GO and KEGG pathways of the differentially expressed genes. (A) The top 20 significant GO terms in the EU vs. HF group. (B) The top 20 significant GO terms in the SK vs. HF group. (C) The top 20 significant KEGG pathway in the EU vs. HF group. (D) The top 20 significant KEGG pathway in the SK vs HF group.

A KEGG pathway enrichment analysis of the differentially expressed genes in the breast muscle tissue of commercial meat pigeons and HF was carried out. The differentially expressed genes in the breast muscle of EU and HF were significantly enriched in 217 pathways, 21 of which were extremely significantly enriched (p < 0.05) (Supplementary Table S13; Figure 4C). The differentially expressed genes in the breast muscle of SK and HF were significantly enriched in 178 pathways, of which 31 were extremely significantly enriched (p < 0.05) (Supplementary Table S14; Figure 4D). The significantly enriched pathways exhibited by the two meat breeds included participation in cell proliferation, differentiation, metabolism, and synthesis. These pathways included cell adhesion molecules (CAMs), ECM-receptor interaction, glycolysis/gluconeogenesis, focal adhesion, insulin signaling pathway, mucin type O-glycan biosynthesis, biotin metabolism, and the adipocytokine signaling pathway. Differential gene expression in the breast muscle of meat pigeon breeds was enriched in 111 pathways, 19 of which were significantly enriched (p < 0.05) (Supplementary Table S15).

 Verification of Key Gene Expression

Four genes, including SYNE1, INSC, IGFBP1 and MAFF exhibited significant different expression patterns between the meat breeds and the HF breed. The expression patterns of these genes in pigeon breast muscle were consistent with the RNA sequencing results, which verified the accuracy of the comparative transcriptome sequencing results (Figure 5; Supplementary Table S16).

(A) Gene expression of pigeon breast muscle transcriptome sequencing. FPKM (Fragments Per Kilobase of exon model per Million mapped reads). (B) Expression of five significant differentially expressed genes validated by qRT-PCR. The same lowercase letters in the histogram represent no significant difference (p > 0.05), and different letters represent significant difference (p < 0.05).

 Discussion

In this study, a genomics perspective was used to analyze the genes of two meat pigeon breeds (EU and SK) and an ornamental breed. To further study the genetic mechanism associated with the phenotype differences in the pectoral muscle between the two different breeds, a comparative transcriptome was used to reveal the key genes related to growth and development and more importantly, to discover genes that regulate muscle function. The selection signal combined with the comparative transcriptome revealed INSC, CALCB, ZBTB21, B2M, and LOC110356506 as the genes that affect breast muscle development in meat pigeons. INSC can regulate mesoderm differentiation of mouse embryonic stem (ES) cells (Ishibashi et al., 2016). In addition, the human INSC gene, which is closely associated with the CALCB gene by an interval of about 30 kb, was assigned to human chromosome 11p15.2-p15.1 (Katoh and Katoh 2003). While a short-term RNA interference-mediated CALCB knockdown had no effect on the proliferation and clonogenic growth of ES cells in vitro, its long-term knockdown decreased ES growth both in vitro and in vivo. The INSC gene also plays an important role in the growth and development of meat pigeon breast muscle. In the differential expression analysis, 673 genes were identified, among which the genes related to growth and development included IGFBP1, IGFBP4, FOXO3, HMGA1, and FAM184B. Studies have shown that the IGFBP1 gene has an important regulatory effect on fetal growth and development. Both hypoxia and leucine deprivation can increase the level of the IGFBP1 gene expression and phosphorylation, inhibit the effect of IGF, and lead to impaired embryonic development (Kajimura et al., 2005; Seferovic et al., 2009). IGF-1 and IGFBP1 have anabolic effects on skeletal muscle and are related to the preservation of lean meat (Velloso 2008). Consequently, IGFBP1 play a role in the growth and metabolism of pigeons and regulate the development of the pectoral muscle.

Compared with commercial meat pigeons, HFs are selected for light weight individuals. At the same time, the flying ability of HFs is inferior to that of carrier pigeons, characterized by a short hovering distance and slow flying speed. Genomic and transcriptomic analyses revealed that SLC16A10, S100B, SYNE1, HECW2, and CASQ2 genes were positively selected. Studies in mice have shown that the SLC16A10 gene is involved in promoting the cellular transport of thyroid hormone (Müller et al., 2014). Studies in pigs have demonstrated that a low protein diet induces higher expression of the SLC16A10 gene, resulting in limited protein synthesis and growth of the longissimus dorsi (Wang D. et al., 2017). Mutations in the HECW2 gene can cause neurodevelopmental delay, and the clinical features shared by patients include severe developmental delay and hypotonia (Berko et al., 2017). In addition, this gene mutation has also been linked to epilepsy-associated developmental delay (Ullman et al., 2018). CASQ2 plays an important role in regulating Ca2+ release in the sarcoplasmic reticulum, buffering of Ca2+ in the sarcoplasmic reticulum, and promoting the closure of cardiac ryanodine receptors during diastole (Györke and Terentyev 2008). CASQ2 is expressed in slow muscle throughout the lifespan of mice, but only in fast muscle during the newborn stage and early development. Indeed, CASQ2 −/− mice display ultrastructural changes only in the rapid twitch muscle (Valle et al., 2016).

A defect in nesprin-1 encoded by the SYNE1 gene can cause Emery-Dreifuss muscular dystrophy (EDMD), which is characterized by joint contracture, myasthenia, and cardiac abnormalities (Zhang et al., 2007; Pillers and Von Bergen 2016; Chen et al., 2017; Heller et al., 2020). The protein encoded by SYNE1 is widely expressed in a variety of tissues and connects the outer membrane of the nuclear membrane with the cytoskeleton by interacting with F-actin (Rajgor and Shanahan 2013) and is highly expressed in striated muscle (Zhang et al., 2010). At the same time, several studies throughout the world have reported that SYNE1 gene mutations can cause autosomal receptive cerebellar ataxia type 8 (SCAR8) (Izumi et al., 2013; Noreau et al., 2013; Hamza et al., 2015; Mademan et al., 2016; Synofzik et al., 2016). The SYNE1 gene was selected in HF pigeons and found to be significantly up-regulated (p < 0.05) compared with that in meat pigeons. The body weight and breast muscle weight of HF pigeons are significantly lower (p < 0.05) than those of commercial meat pigeons (Table 1), indicating that this gene has a negative regulatory effect on the growth and development of pigeons. With the enrichment analysis, SYNE1 was found to be involved in actin binding (GO: 0003779) and muscle cell differentiation (GO: 0042692). The HF breed has a lighter body that is suitable for gliding in the air, but has a slower flight speed and is not suitable for long-distance flight. It is speculated that mutations in the SYNE1 gene will affect the development and function of HF breast muscle.

 Conclusion

This study revealed the selection of genome regions between commercial meat pigeons and Chinese traditional ornamental pigeons using genomics analyses. The results of transcriptome sequencing and differential expression analysis revealed the molecular mechanism of breast muscle phenotypic differences between commercial meat pigeons and Chinese traditional ornamental pigeons. Multi-omics further revealed multiple genes related to cell differentiation, muscle development, and skeletal muscle function. Among these genes, INSC and CALCB were related to cell differentiation and were positively selected and up-regulated in both EU pigeons and SK pigeons. Thus, these genes may be involved in promoting the growth and development of the pigeon breast muscle. More importantly, we found that the SYNE1 gene was related to actin binding (GO: 0003779) and muscle cell differentiation (GO: 0042692) in the ornamental pigeon population. This study provides a few genes that require further study of their molecular mechanisms related to muscle development, which will help understand the genomic imprints of different types of pigeons due to artificial selection.

 Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

 Ethics Statement

The animal study was reviewed and approved by The process of collecting pigeon samples strictly follows the animal ethics procedures and guidelines of the People’s Republic of China. Pigeons used in this study were approved by the Ethics and Animal Welfare Committee of Shanghai Academy of Agricultural Sciences (No. SAASPZ0521012).

 Author Contributions

Experimental design: JY and WD; Project leader: CY; Article writer: HH; Data analysis and mining: HH and XW; Sample collection: XC, WL, YT, and AB; Research materials: QW and WZ.

 Funding

This research has been supported by Shanghai Agriculture Applied Technology Development Program, China (Grant No. 2018-02-08-00-12-F01546) and SAAS Program for Excellent Research Team No. 2017 (B-03).

 Conflict of Interest

Authors AB, QW, and WZ, are employed by Shanghai Jinhuang Pigeon Company.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

The authors are grateful to The Shanghai Jinhuang Pigeon Co., Ltd., provided support in the process of collecting research materials on High Flyers and commercial pigeon breeds. The authors would like to thank Mr. Liu for providing technical support in bioinformatics analysis. We thank International Science Editing (http://www.internationalscienceediting.com) for editing this manuscript.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.735795/full#supplementary-material

References Berko E. R. Cho M. T. Eng C. Shao Y. Sweetser D. A. Waxler J. (2017). De Novo missense Variants inHECW2are Associated with Neurodevelopmental Delay and Hypotonia. J. Med. Genet. 54, 84–86. 10.1136/jmedgenet-2016-103943 Berry C. Thomas M. Langley B. Sharma M. Kambadur R. (2002). Single Cysteine to Tyrosine Transition Inactivates the Growth Inhibitory Function of Piedmontese Myostatin. Am. J. Physiology-Cell Physiol. 283, C135–C141. 10.1152/ajpcell.00458.2001 Biro D. Freeman R. Meade J. Roberts S. Guilford T. (2007). Pigeons Combine Compass and Landmark Guidance in Familiar Route Navigation. Pnas 104, 7471–7476. 10.1073/pnas.0701575104 Blasco R. Finlayson C. Rosell J. Marco A. S. Finlayson S. Finlayson G. (2014). The Earliest pigeon Fanciers. Sci. Rep. 4, 5971. 10.1038/srep05971 Boer E. F. Van Hollebeke H. F. Park S. Infante C. Menke D. B. Shapiro M. D. (2019). Pigeon Foot Feathering Reveals Conserved Limb Identity Networks. Dev. Biol. 454 (2), 128–144. 10.1016/j.ydbio.2019.06.015 Chen Z. Ren Z. Mei W. Ma Q. Shi Y. Zhang Y. (2017). A Novel SYNE1 Gene Mutation in a Chinese Family of Emery-Dreifuss Muscular Dystrophy-like. BMC Med. Genet. 18, 63. 10.1186/s12881-017-0424-5 China Animal Agriculture Association (2018). China Poultry Industry Development Report in 2018. Beijing. Domyan E. T. Kronenberg Z. Infante C. R. Vickrey A. I. Stringham S. A. Bruders R. (2016). Molecular Shifts in Limb Identity Underlie Development of Feathered Feet in Two Domestic Avian Species. Elife 5, e12115. 10.7554/eLife.12115 Dong X. Cao H. Mao H. Hong Q. Yin Z. (2019). Association of MyoD1 Gene Polymorphisms with Meat Quality Traits in Domestic Pigeons (Columba livia). J. Poult. Sci. 56, 20–26. 10.2141/jpsa.0170182 Dushyanth K. Bhattacharya T. K. Shukla R. Chatterjee R. N. Sitaramamma T. Paswan C. (2016). Gene Expression and Polymorphism of Myostatin Gene and its Association with Growth Traits in Chicken. Anim. Biotechnol. 27, 269–277. 10.1080/10495398.2016.1182541 Gazda M. A. Andrade P. Afonso S. Dilytė J. Archer J. P. Lopes R. J. (2018). Signatures of Selection on Standing Genetic Variation Underlie Athletic and Navigational Performance in Racing Pigeons. Mol. Biol. Evol. 35, 1176–1189. 10.1093/molbev/msy030 Grochowska E. Borys B. Lisiak D. Mroczkowski S. (2019). Genotypic and Allelic Effects of the Myostatin Gene (MSTN) on Carcass, Meat Quality, and Biometric Traits in Colored Polish Merino Sheep. Meat Sci. 151, 4–17. 10.1016/j.meatsci.2018.12.010 Guernec A. Chevalier B. Duclos M. J. (2004). Nutrient Supply Enhances Both IGF-I and MSTN mRNA Levels in Chicken Skeletal Muscle. Domest. Anim. Endocrinol. 26, 143–154. 10.1016/j.domaniend.2003.10.001 Gusev A. Ko A. Shi H. Bhatia G. Chung W. Penninx B. W. J. H. (2016). Integrative Approaches for Large-Scale Transcriptome-wide Association Studies. Nat. Genet. 48, 245–252. 10.1038/ng.3506 Györke S. Terentyev D. (2008). Modulation of Ryanodine Receptor by Luminal Calcium and Accessory Proteins in Health and Cardiac Disease. Cardiovasc. Res. 77, 245–255. 10.1093/cvr/cvm038 Hamza W. Ali Pacha L. Hamadouche T. Muller J. Drouot N. Ferrat F. (2015). Molecular and Clinical Study of a Cohort of 110 Algerian Patients with Autosomal Recessive Ataxia. BMC Med. Genet. 16, 36. 10.1186/s12881-015-0180-3 Heller S. A. Shih R. Kalra R. Kang P. B. (2020). Emery‐Dreifuss Muscular Dystrophy. Muscle Nerve 61, 436–448. 10.1002/mus.26782 Ishibashi R. Kozuki S. Kamakura S. Sumimoto H. Toyoshima F. (2016). c-Rel Regulates Inscuteable Gene Expression during Mouse Embryonic Stem Cell Differentiation. J. Biol. Chem. 291, 3333–3345. 10.1074/jbc.m115.679563 Izumi Y. Miyamoto R. Morino H. Yoshizawa A. Nishinaka K. Udaka F. (2013). Cerebellar Ataxia with SYNE1 Mutation Accompanying Motor Neuron Disease. Neurology 80, 600–601. 10.1212/wnl.0b013e3182815529 Kajimura S. Aida K. Duan C. (2005). Insulin-like Growth Factor-Binding Protein-1 (IGFBP-1) Mediates Hypoxia-Induced Embryonic Growth and Developmental Retardation. Proc. Natl. Acad. Sci. 102, 1240–1245. 10.1073/pnas.0407443102 Kanehisa M. Araki M. Goto S. Hattori M. Hirakawa M. Itoh M. (2008). KEGG for Linking Genomes to Life and the Environment. Nucleic Acids Res. 36, D480–D484. 10.1093/nar/gkm882 Katoh M. Katoh M. (2003). Identification and Characterization of Human Inscuteable Gene In Silico. Int. J. Mol. Med. 11, 111–116. 10.3892/ijmm.11.1.111 Kokoszyński D. Stęczny K. Żochowska-Kujawska J. Sobczak M. Kotowicz M. Saleh M. (2020). Carcass Characteristics, Physicochemical Properties, and Texture and Microstructure of the Meat and Internal Organs of Carrier and King Pigeons. Animals (Basel) 10, 1315. 10.3390/ani10081315 Li H. Durbin R. (2009). Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform. Bioinformatics 25, 1754–1760. 10.1093/bioinformatics/btp324 Li H. Handsaker B. Wysoker A. Fennell T. Ruan J. Homer N. (2009). Genome Project Data Processing SThe Sequence Alignment/Map Format and SAMtools. Bioinformatics 25, 2078–2079. 10.1093/bioinformatics/btp352 Li R. Zeng W. Ma M. Wei Z. Liu H. Liu X. (2020). Precise Editing of Myostatin Signal Peptide by CRISPR/Cas9 Increases the Muscle Mass of Liang Guang Small Spotted Pigs. Transgenic Res. 29, 149–163. 10.1007/s11248-020-00188-w Liu H. H. Mao H. G. Dong X. Y. Cao H. Y. Liu K. Yin Z. Z. (2019). Expression of MSTN Gene and its Correlation with Pectoralis Muscle Fiber Traits in the Domestic Pigeons (Columba livia). Poult. Sci. 98, 5265–5271. 10.3382/ps/pez399 Mademan I. Harmuth F. Giordano I. Timmann D. Magri S. Deconinck T. (2016). MultisystemicSYNE1ataxia: Confirming the High Frequency and Extending the Mutational and Phenotypic Spectrum. Brain 139, e46. 10.1093/brain/aww115 Mao H. G. Cao H. Y. Liu H. H. Dong X. Y. Xu N. Y. Yin Z. Z. (2018). Association of ADSL Gene Polymorphisms with Meat Quality and Carcass Traits in Domestic Pigeons (Columba livia). Br. Poult. Sci. 59, 604–607. 10.1080/00071668.2018.1493188 Miretti S. Martignani E. Accornero P. Baratta M. (2013). Functional Effect of Mir-27b on Myostatin Expression: a Relationship in Piedmontese Cattle with Double-Muscled Phenotype. BMC Genomics 14, 194. 10.1186/1471-2164-14-194 Müller J. Mayerl S. Visser T. J. Darras V. M. Boelen A. Frappart L. (2014). Tissue-specific Alterations in Thyroid Hormone Homeostasis in Combined Mct10 and Mct8 Deficiency. Endocrinology 155, 315–325. 10.1210/en.2013-1800 Noreau A. Bourassa C. V. Szuto A. Levert A. Dobrzeniecka S. Gauthier J. (2013). SYNE1 Mutations in Autosomal Recessive Cerebellar Ataxia. JAMA Neurol. 70, 1296–1301. 10.1001/jamaneurol.2013.3268 Pertea M. Pertea G. M. Antonescu C. M. Chang T.-C. Mendell J. T. Salzberg S. L. (2015). StringTie Enables Improved Reconstruction of a Transcriptome from RNA-Seq Reads. Nat. Biotechnol. 33, 290–295. 10.1038/nbt.3122 Pillers D.-A. Von Bergen N. (2016). Emery–Dreifuss Muscular Dystrophy: a Test Case for Precision Medicine. Tacg 9, 27–32. 10.2147/tacg.s75028 Pomianowski J. F. Mikulski D. Pudyszak K. Cooper R. G. Angowski M. Jóźwik A. (2009). Chemical Composition, Cholesterol Content, and Fatty Acid Profile of pigeon Meat as Influenced by Meat-type Breeds. Poult. Sci. 88, 1306–1309. 10.3382/ps.2008-00217 Rajgor D. Shanahan C. M. (2013). Nesprins: from the Nuclear Envelope and beyond. Expert Rev. Mol. Med. 15, e5. 10.1017/erm.2013.6 Seferovic M. D. Ali R. Kamei H. Liu S. Khosravi J. M. Nazarian S. (2009). Hypoxia and Leucine Deprivation Induce Human Insulin-like Growth Factor Binding Protein-1 Hyperphosphorylation and Increase its Biological Activity. Endocrinology 150, 220–231. 10.1210/en.2008-0657 Shao Y. Tian H.-Y. Zhang J.-J. Kharrati-Koopaee H. Guo X. Zhuang X.-L. (2020). Genomic and Phenotypic Analyses Reveal Mechanisms Underlying Homing Ability in Pigeon. Mol. Biol. Evol. 37, 134–148. 10.1093/molbev/msz208 Shapiro M. D. Kronenberg Z. Li C. Domyan E. T. Pan H. Campbell M. (2013). Genomic Diversity and Evolution of the Head Crest in the Rock pigeon. Science 339, 1063–1067. 10.1126/science.1230422 Srikanth K. Kim N. Y. Park W. Kim J. M. Kim K. D. Lee K. T. (2019). Comprehensive Genome and Transcriptome Analyses Reveal Genetic Relationship, Selection Signature, and Transcriptome Landscape of Small-Sized Korean Native Jeju Horse. Sci. Rep. 9, 16672. 10.1038/s41598-019-53102-8 Stinckens A. Luyten T. Bijttebier J. Van den Maagdenberg K. Dieltiens D. Janssens S. (2008). Characterization of the Complete porcineMSTNgene and Expression Levels in Pig Breeds Differing in Muscularity. Anim. Genet. 39, 586–596. 10.1111/j.1365-2052.2008.01774.x Stringham S. A. Mulroy E. E. Xing J. Record D. Guernsey M. W. Aldenhoven J. T. (2012). Divergence, Convergence, and the Ancestry of Feral Populations in the Domestic Rock pigeon. Curr. Biol. 22, 302–308. 10.1016/j.cub.2011.12.045 Synofzik M. Smets K. Mallaret M. Di Bella D. Gallenmüller C. Baets J. (2016). SYNE1 Ataxia Is a Common Recessive Ataxia with Major Non-cerebellar Features: a Large Multi-centre Study. Brain 139, 1378–1393. 10.1093/brain/aww079 Ullman N. L. Smith-Hicks C. L. Desai S. Stafstrom C. E. (2018). De Novo HECW2 Mutation Associated with Epilepsy, Developmental Decline, and Intellectual Disability: Case Report and Review of Literature. Pediatr. Neurol. 85, 76–78. 10.1016/j.pediatrneurol.2018.03.005 Valle G. Vergani B. Sacchetto R. Reggiani C. De Rosa E. Maccatrozzo L. (2016). Characterization of Fast-Twitch and Slow-Twitch Skeletal Muscles of Calsequestrin 2 (CASQ2)-Knock Out Mice: Unexpected Adaptive Changes of Fast-Twitch Muscles Only. J. Muscle Res. Cell Motil 37, 225–233. 10.1007/s10974-016-9463-3 Velloso C. P. (2008). Regulation of Muscle Mass by Growth Hormone and IGF-I. Br. J. Pharmacol. 154, 557–568. 10.1038/bjp.2008.153 Vickrey A. I. Bruders R. Kronenberg Z. Mackey E. Bohlender R. J. Maclary E. T. (2018). Introgression of Regulatory Alleles and a Missense Coding Mutation Drive Plumage Pattern Diversity in the Rock pigeon. Elife 7, e34803. 10.7554/eLife.34803 Wang D. Wan X. Peng J. Xiong Q. Niu H. Li H. (2017a). The Effects of Reduced Dietary Protein Level on Amino Acid Transporters and mTOR Signaling Pathway in Pigs. Biochem. Biophysical Res. Commun. 485, 319–327. 10.1016/j.bbrc.2017.02.084 Wang K. Li M. Hakonarson H. (2010). ANNOVAR: Functional Annotation of Genetic Variants from High-Throughput Sequencing Data. Nucleic Acids Res. 38, e164. 10.1093/nar/gkq603 Wang L. Li X. Ma J. Zhang Y. Zhang H. (2017b). Integrating Genome and Transcriptome Profiling for Elucidating the Mechanism of Muscle Growth and Lipid Deposition in Pekin Ducks. Sci. Rep. 7, 3837. 10.1038/s41598-017-04178-7 Weir B. S. Cockerham C. C. (1984). ESTIMATING F-STATISTICS FOR THE ANALYSIS OF POPULATION STRUCTURE. Evolution 38, 1358–1370. 10.1111/j.1558-5646.1984.tb05657.x Ye M. Zhou B. Wei S. Ding M. Lu X. Shi X. (2016). Transcriptomic Analysis Identifies Candidate Genes Related to Intramuscular Fat Deposition and Fatty Acid Composition in the Breast Muscle of Squabs (Columba). G3 (Bethesda) 6, 2081–2090. 10.1534/g3.116.029793 Young M. D. Wakefield M. J. Smyth G. K. Oshlack A. (2010). Gene Ontology Analysis for RNA-Seq: Accounting for Selection Bias. Genome Biol. 11, R14. 10.1186/gb-2010-11-2-r14 Zhang J. Felder A. Liu Y. Guo L. T. Lange S. Dalton N. D. (2010). Nesprin 1 Is Critical for Nuclear Positioning and anchorage. Hum. Mol. Genet. 19, 329–341. 10.1093/hmg/ddp499 Zhang L. X. Xu Z. Maitusun A. K. Fu R. Wang Z. W. Li F. M. (2020). Study on the Growth and Development of 0-28 Day Old Tarim pigeon. J. Xinjiang Agric. Univ. 43, 49–53. (in Chinese). Zhang Q. Bethmann C. Worth N. F. Davies J. D. Wasner C. Feuer A. (2007). Nesprin-1 and -2 Are Involved in the Pathogenesis of Emery-Dreifuss Muscular Dystrophy and Are Critical for Nuclear Envelope Integrity. Hum. Mol. Genet. 16, 2816–2833. 10.1093/hmg/ddm238