Next Article in Journal
Phenolic Compounds Recovery from Pomegranate (Punica granatum L.) By-Products of Pressurized Liquid Extraction
Next Article in Special Issue
Interlaboratory Validation of a DNA Metabarcoding Assay for Mammalian and Poultry Species to Detect Food Adulteration
Previous Article in Journal
Applications of Plant Protein in the Dairy Industry
Previous Article in Special Issue
New Insights on the Impact of Cattle Handling on Post-Mortem Myofibrillar Muscle Proteome and Meat Tenderization
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Foods
Volume 11
Issue 8
10.3390/foods11081068
foods-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data

by
Zehu Yuan
1,
Ling Ge
2,
Weibo Zhang
2,
Xiaoyang Lv
1,
Shanhe Wang
2,
Xiukai Cao
1 and
Wei Sun
1,2,*
1
Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education, Yangzhou University, Yangzhou 225000, China
2
College of Animal Science and Technology, Yangzhou University, Yangzhou 225000, China
*
Author to whom correspondence should be addressed.
Foods 2022, 11(8), 1068; https://doi.org/10.3390/foods11081068
Submission received: 19 February 2022 / Revised: 29 March 2022 / Accepted: 4 April 2022 / Published: 7 April 2022
(This article belongs to the Special Issue Recent Advances in OMICs Technologies and Application for Ensuring Meat Quality, Safety and Authenticity)
Browse Figures
Versions Notes

Abstract

:
Tenderness is an important indicator of meat quality. Novel isoforms associated with meat tenderness and the role of the CCCTC-binding factor (CTCF) in regulating alternative splicing to produce isoforms in sheep are largely unknown. The current project studied six sheep from two crossbred populations (Dorper × Hu × Hu, DHH and Dorper × Dorper × Hu, DDH) with divergent meat tenderness. Pooled Iso-seq data were used to annotate the sheep genomes. Then, the updated genome annotation and six RNA-seq data were combined to identify differentially expressed isoforms (DEIs) in muscles between DHH and DDH. These data were also combined with peaks detected from CTCF ChIP-seq data to investigate the regulatory role of CTCF for the alternative splicing. As a result, a total of 624 DEIs were identified between DDH and DHH. For example, isoform 7.524.18 transcribed from CAPN3 may be associated with meat tenderness. In addition, a total of 86 genes were overlapped between genes with transcribed DEIs and genes in differential peaks identified by CTCF ChIP-seq. Among these overlapped genes, ANKRD23 produces different isoforms which may be regulated by CTCF via methylation. As preliminary research, our results identified novel isoforms associated with meat tenderness and revealed the possible regulating mechanisms of alternative splicing to produce isoforms.

1. Introduction

Meat tenderness is determined by several intrinsic factors (e.g., proteolytic activity, amount of glycogen, fiber type and connective tissue) and extrinsic factors (e.g., animal breed and nutritional conditions) [1]. Among these factors, the calpain system plays an essential role in meat tenderness [2]. In beef, the most important factor that determines tenderness is the proteolytic activity of the calpain system. Calpain3 (CAPN3), a member of the calpain system, is mainly expressed in skeletal muscle, especially in type II fibers [3,4,5], and it may play an important role in the calpain system to regulate meat tenderness [6]. A previous study suggested that the expression level of CAPN3 protein and mRNA was significantly associated with tenderness [7,8]. Another member of the calpain system is calpastatin (CAST). Accumulated evidence suggests that different isoforms transcribed from a gene due to alternative splicing may play different biological roles. For example, distinct isoforms due to alternative splicing [9] in CAST have different biological functions in beef meat tenderness.
PacBio long-reads isoform sequencing (Iso-seq), a third-generation sequencing technology, can directly obtain full-length isoforms [10]. By combing with traditional RNA-sequencing (RNA-seq), Iso-seq can contribute to identify novel isoform associated with complex traits [11,12,13]. Although differentially expressed gene expression studies have been implemented for meat quality traits in sheep [14,15], studies focusing on the investigation of novel isoforms in sheep muscle are still scarce. Alternative splicing is the main mechanism to diversify isoforms and is regulated by numerous interacting components [16]. For example, the CCCTC-binding factor (CTCF) has been identified as a modulator of the alternative splicing process [17]. However, the role of CTCF in regulating alternative splicing to produce isoforms in sheep is largely unknown. In this context, using Iso-seq, RNA-seq and CTCF chromatin immunoprecipitation sequencing (ChIP-seq) data in sheep muscle can help to identify key isoforms associated with meat tenderness and to reveal the biological mechanisms of alternative splicing regulation to produce isoforms.
Tenderness is an important indicator of meat quality with the greatest consumer appreciation [18]. Crossbreeding is a common method for increasing sheep production. When using crossbreeding methods to increase sheep meat production, some meat quality traits, e.g., tenderness, show a significant difference across different cross-breed populations. Unveiling the genetic factors associated with tenderness is essential for guiding cross-breeding efficiently and profitably. Therefore, revealing the genetic processes that regulate meat quality between different cross-breeds is of great interest to the sheep industry.
The goals of the current study were to identify key isoforms associated with meat tenderness and to reveal the biological mechanisms of alternative splicing regulations to produce isoforms. To achieve these goals, differentially expressed isoforms (DEIs) were identified between two cross-bred populations with divergent meat tenderness by integrating Iso-seq and RNA-seq data. Then, we combined peaks detected from CTCF chromatin immunoprecipitation sequencing (ChIP-seq) data with DEIs to investigate the possible mechanism with which CTCF regulates the alternative splicing process. Our results have the potential to reveal novel isoforms associated with meat tenderness and uncover the possible regulating mechanism of alternative splicing to produce isoforms.

2. Materials and Methods

The animal experiment of this study was approved by the Experimental Animal Ethical Committee of Yangzhou University (File NO. 202103294). The overall experimental design of the current study is shown in Figure 1.

2.1. Animal Samples and Meat Traits

In the current study, a total of six unrelated 6-month-old male sheep, including three Dorper × Hu × Hu (DHH) and three Dorper × Dorper × Hu (DDH), were selected from a large population for sampling and measuring meat tenderness. Six male sheep were fed in the same environment from 2 months old to 6 months old. Six sheep were selected with similar live weights (DHH = 40.80 ± 4.57 kg, DDH = 41.50 ± 0.95 kg, Table 1) to rule out the effect of live weight on meat tenderness. Live weight was measured before slaughter (24 h in advance). After slaughter, longissimus dorsi from between the 12th and 13th ribs from six sheep were sampled within 30 min. Each muscle sample was packed into a 1.5 mL cryotube with triplicate. All the muscle samples were frozen in liquid and stored in a −80 °C refrigerator. After slaughter, carcass weight was measured by an electronic scale. Then, 45 min after slaughter, longissimus dorsi were cut after the 12th rib and the samples were cooked for 35 min at 71 °C [19]. Then, six 1 cm2 cylinders were cut and all samples were tested using a meat tenderness analyzer (CL3, NanJing, CN). The difference of meat tenderness between DDH and DHH was tested by the t-test function in the basic environment of R 4.1.1 [20]. A p value smaller than 0.05 denotes significance.

2.2. RNA-Extraction, RNA-seq and Iso-seq

A Trizol reagent Kit (TaKaRa, Kusatsu, Shiga, Japan) was used to extract the total RNA from six muscle tissues. Then, Nanodrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA) and a 2100 Bioanalyzer (Agilent Technologies, Waldronn, Germany) were used to evaluate the RNA Integrity Number (RIN >7 passed quality control) and 28S/18S ratio (ratio > 1.0 passed quality control). After quality control of all six RNA samples, they were used for library construction. Following library quality control, six libraries were sequenced on the MGISEQ-2000 (BIG) platform [13].
Pooled RNA samples of six muscle tissues were used for Iso-seq. Briefly, cDNA was transcribed from a pooled RNA sample using the SMARTerTM PCR cDNA synthesis kit (Takara Biotechnology, Dalian, China). Following PCR amplification, PCR product purification, size selection (>1 kb), SMRTbell library construction, and library quality control, the library was sequenced on the PacBio sequencing platform. The raw Iso-seq data were processed using the SMRT Link v8.0 pipeline. Briefly, the BAM file was processed to obtain the circular consensus sequence (CCS). Full-length CCS reads with 5′ and 3′ cDNA primers and polyA were defined as full-length non-concatemer (FLNC) reads. FLNC reads were corrected by six clean RNA-seq data using LoRDEC v0.9 software [21]. The quality control of six RNA-seq data was implemented by SOAPnuke v2.1.0 [22].

2.3. Update Sheep Reference Annotation File

The pooled Iso-seq data were used to update the sheep reference annotation file. After correction, FLNC reads were mapped to the reference genome (http://ftp.ensembl.org/pub/release-102/fasta/ovis_aries_rambouillet/dna/, accessed on 29 March 2022) using GMAP [23]. Isoforms that met at least one of the following criteria were kept: first, an isoform was identified by at least two FLNC reads; second, an isoform was identified by only one FLNC read with a percentage identity (PID) of greater than 99%; third, all splicing sites in an isoform were identified by RNA-seq data; fourth, all alternative splicing events in a identified isoform were also annotated by the reference genome annotation file. Isoforms that overlapped less than 20% of their length on the same strand were identified as distinct isoforms. A novel gene was defined as a gene locus overlapping less than 20% of their length with known genes. An isoform with a new emerged intron (exon) or a final splice site of 3′ changed ends was defined as a novel isoform. In addition, the alternative splicing events were detected by the AStalavista v3.2 software [24].

2.4. Differentially Expressed Isoform Detection and Functional Analyses

Bowtie2 was used to align short reads to novel annotation files [25]. The DESeq2 R Bioconductor package was used to identify DEIs between two groups [26]. An isoform with a fold change (FC) greater than two and a false discovery rate (FDR) smaller than 0.05 was defined as a significant DEI. All identified DEIs were blasted to GO database using Diamond [27] and blasted to the KEGG database using KOBAS [28] to extract their potential biological functions. The GOplot 1.0.2 R package was used to visualize GO enrichment results [29]. GO terms and KEGG pathways with an FDR smaller than 0.05 were regarded as having significant enrichment.

2.5. Validation of Differentially Expressed Isoforms

The isoform abundance of five DEIs were quantified for validating the result of identified DEIs by quantitative real-time PCR (qRT-PCR). Five isoform-specific paired premiers were designed by SnapGene® v2.3.2 software (from Insightful Science; available at snapgene.com) and were synthesized by Tsingke Biotechnology Co., Ltd. (Nanjing, China, Table S1). cDNA was synthesized using a FastKing gDNA Dispelling RT Super Mix (TIANGEN, Beijing, China). The qRT-PCR was performed on a CFX96 Connect™ Real-Time System (BIO-RAD, Hercules, CA, USA) using a 20 μL reaction volume, including 1 μL of cDNA in 10 μL of 2× TSINGKE Master qPCR Mix (SYBR Green I) (Tsingke, Nanjing, China), 0.8 μL (10 μm/μL) each of the forward and reverse primers, and 7.4 μL of distilled water. The abundance of the GAPDH was used as the control. Each biological sample was implemented in triplicate, and the 2−ΔΔCt method was used for calculating the relative expression level of isoforms.

2.6. CTCF ChIP-seq

We aimed to explore the possible regulating mechanism of CTCF in alternative splicing in sheep muscle. Two muscle samples with extreme meat tenderness were selected to implement CTCF ChIP-seq. CTCF ChIP-seq was implemented by a commercial sequencing provider (igenebook Technology Co., Ltd., Wuhan, China). Briefly, chromatin was crosslinked by formaldehyde, following nuclear processing, chromatin digestion, DNA-protein compound capture, decrosslinking of DNA-protein compounds, and the purification of DNA. Finally, the input and ChIP DNA samples were sequenced on the Illumina Hiseq X ten platform.

2.7. Bioinformatics Analysis of CTCF ChIP-seq

The quality of CTCF ChIP-seq data was evaluated by fastqc v0.11.5 [30]. The quality control of the raw CTCF ChIP-seq data was implemented by Trimmomatic v0.36 [31]. After quality control, clean reads were aligned to the sheep reference genome by BWA v0.7.1 [32]. The read distribution in different genomic regions was investigated using the ChIPseeker R Bioconductor package [33]. An Upset plot was plotted by the UpSetR R package [34]. Potential peaks were called by MACS v2.1.1 [35]. Differential peaks were detected by DiffBind v1.16.3 (https://bioconductor.org/packages/release/bioc/html/DiffBind.html, accessed on 29 March 2022). Peaks with an FDR < 0.05 and a Fold value >0 were defined as significant differential peaks. Motifs in significant differential peaks were predicted by HOMER v3 [36].

2.8. Overlapping between DEIs and Differential Peaks Called from ChIP-seq Data

The overlapping between genes with transcribed DEIs and genes located in differential peaks were investigated to explore the regulating role of CTCF in the alternative splicing process. The significance of overlapping was tested by the GeneOverlap R Bioconductor package (https://bioconductor.org/packages/release/bioc/html/GeneOverlap.html, accessed on 29 March 2022).

3. Results

3.1. Meat Traits

In the current study, a total of three meat traits were measured after slaughter. The mean value and standard deviation of three meat traits are documented in Table 1. Among these traits, shear force in DDH was significantly (p = 0.0406) higher than that in DHH. Other meat traits did not show a significant difference between DDH and DHH.

3.2. Update of Reference Genome Annotation File

To update the sheep genome annotation (Generic Feature Format, GFF), which is essential for accurately quantifying the abundance of isoforms, six RNA samples were pooled to implement Iso-seq. In total, 442,966 polymerase reads were produced by pooled Iso-seq (Table S2). After pre-processing, a total of 247,201 FLNC reads were kept for further analysis (Table S2). Novel genes and isoforms were identified according to the genomic position of each FLNC. In total, 18,959 gene loci (9184 known and 9775 novel) were identified (Figure S1a). Among these detected genes, 5104 (26.92%) loci had two or more transcripts (Figure S1b). In sheep muscle, 20,205 novel isoforms (57.51%) were transcribed from annotated genes (Figure S1c), followed by those transcribed from novel genes. Identified potential novel gene loci and novel isoforms as well as sheep reference annotation files (GFF file) were merged to obtain an updated genome annotation file (Table S3). Alternative splicing is the main mechanism to diversify isoforms. In this study, a total of 38,070 alternative splicing events were detected by pooled Iso-seq data (Figure S2). For example, seven novel isoforms in six pooled samples were identified in CAST due to alternative splicing (Figure 2).

3.3. Differentially Expressed Isoforms and Functional Analysis

Six RNA-seq data were aligned to the updated genome annotation file (Table S3) to identify the DEIs between DDH and DHH. As a result, a total of 624 DEIs were identified (Figure 3, Table S4). These 624 DEIs were transcribed from 492 genes, suggesting that some genes could produce more than one DEI. The most significant DEIs were X.351.3 (FHL1, FDR = 8.21 × 10−15) and 11.673.91 (MYH2, FDR = 8.21 × 10−15).
All detected DEIs were significantly (FDR < 0.05) enriched in 280 GO terms (Figure 4a, Table S5). The most significant GO term was contractile fiber (GO:0043292, FDR = 5.38 × 10−9, Figure 4a, Table S5). All detected DEIs were significantly (FDR < 0.05) enriched in 17 KEGG pathways (Figure 4b, Table S6). The most significant pathway was the glucagon signaling pathway (ko049222).

3.4. Validation for Target Isoforms

To validate the result of the identified DEIs, the qRT-PCR of five DEIs was implemented. Among these five isoforms, three DEIs identified from the transcriptome analysis were significantly expressed (Figure 5). The expression trends of these five isoforms, determined by qRT-PCR, were consistent with the transcriptome data.

3.5. CTCF ChIP-seq

Two muscle samples with divergent meat tenderness were selected for CTCF ChIP-seq. After quality control, more than 40 M clean reads in each sample were obtained (Table S7). All clean reads were mapped to the sheep reference genome with a mapping rate greater than 96.46% in each sample (Table S7). Above half of the mapped reads were located in the intergenic region (Figure 6), followed by an intron, promoter, exon, 3′ untranslated regions (UTR) and 5′ UTR.
A total of 4388 differential peaks were detected between two samples with divergent meat tenderness (Table S8). The greatest number of differential peaks were distributed in the intergenic genome region, followed by the intron region, promoter region, exon region, 5′ UTR and 3′ UTR (Figure 7). The motifs in the differential peaks were predicted. A total of 66 motifs were predicted in differential peaks (Table S9).

3.6. Overlap between DEIs and Differential Peaks

To investigate the potential role of CTCF in regulating alternative splicing in sheep muscle, the overlapping analysis between genes with transcribed DEIs and genes located in differential peaks was implemented. A total of 86 overlapped genes were found (Figure 8). The overlapping p-value was smaller than 0.05, indicating that the number of overlapped genes was too great to have been by chance.

4. Discussion

4.1. Update Sheep Reference Genome Annotation File

It has been reported that one possible method to quantify transcripts is mapping the reads to the annotated transcriptome [37]. In this study, pooled Iso-seq data were used to update the reference genome annotation file. As a result, 9775 possible novel genes and 30,513 novel isoforms were identified (Figure S1) in sheep muscle and were added to the sheep reference genome annotation file (Table S3). In a pig study, a total of 10,465 novel genes were identified by integrating Iso-seq and RNA-seq data [38]. Similarly, in sheep tail fat, a total of 9001 novel genes and 36,667 novel isoforms were detected using pooled Iso-seq data [13]. Our results are in line with these published works, which suggests that Iso-seq is a useful method to improve genome annotation in sheep.
Alternative splicing is the main mechanism to diversify isoforms. Here, we highlight an example of an alternative splicing event that may regulate meat tenderness. CAST is linked with meat tenderness across many farm animals [39]. In beef, it has been reported that alternative splicing [9] in CAST could regulate meat tenderness. In sheep, a previous study suggested that SNPs in CAST may regulate exon excision events in CAST (Oar.3.1, Chr5:93439378-93444596), which was identified by RNA-seq [40]. In the current study, this exon excision event was further validated by Figure 2. Taken together, all of these pieces of evidence indicate that alternative splicing in CAST regulated by SNPs may play an important role in the meat tenderness of sheep.

4.2. Novel Isoforms Linked with Meat Tenderness

The calpain system plays an essential role in meat tenderness [2]. Calpain3 (CAPN3), a member of the calpain system, is mainly expressed in skeletal muscle and may play an important role in the calpain system to regulate meat tenderness [6]. In the current study, a novel transcript, 7.524.18, transcribed from CAPN3, was identified as a DEI (Table S4), suggesting its potential role in meat tenderness.
Muscle fiber type also plays an important role in meat tenderness [41,42]. Muscle fiber type is related to the abundance of myosin type [43]. The ratio of slow and fast heavy chain myosin is linked to the ratio of type I and type II fibers [44]. In the current study, 12 DEIs transcribed from four myosin heavy chain genes were identified (11.673.36, 11.673.38, 11.673.91 and 11.673.53 transcribed from MYH2; ENSOART00020028698 transcribed from MYH4; 7.469.116, 7.469.121 and 7.469.129 transcribed from MYH7; 14.723.9, 14.723.27, 14.723.58 and 14.723.65 transcribed from MYBPC2; Table S4). Previous studies have reported that MYH2, MYH4, MYH7 and MYBPC2 play an important role in degerming muscle fiber types [45,46]. Overall, the myosin heavy chain-related isoforms were differentially expressed between DDH and DHH, suggesting that these isoforms may regulate meat tenderness directly or indirectly.
In addition, the amount of glycogen could affect meat tenderness [47]. In the current study, seven DEIs transcribed by two glycogen-related genes were detected (14.271.7 transcribed from GYS1; 21.138.57, 21.138.27, 21.138.13, 21.138.58, 21.138.9 and ENSOART00020017812 transcribed from PYGM). A previous study suggested that PYGM is related to shear force in cattle [48]. These DEIs might regulate tenderness by controlling glycogen content in muscles.
GO and KEGG enrichment analyses were implemented to extract the potential function of DEIs. In pigs, differentially expressed genes between fast and slow muscle were relevant to myofibril and contractile fiber GO terms. In the current study, DEIs were most significantly enriched in contractile fiber (GO:0043292) and myofibril (GO:0030016) GO terms. Our results provide further evidence for a relationship between muscle fiber type transformation, meat tenderness and DEIs. For the KEGG pathway enrichment analysis, a previous study suggested that a gene in the calcium signaling pathway (ko04020) is related to stiffness and affects the speed of fiber degradation during the meat aging process [49]. In the current study, 23 DEIs were significantly enriched in calcium signaling pathways, indicating their potential role in the meat aging process.

4.3. CTCF Might Regulate Alternative Splicing in Sheep Muscle

In the current study, CTCF ChIP-seq was implanted to investigate its potential role in regulating alternative splicing in sheep muscle with divergent meat tenderness. The main function of alternative splicing processes are to diversify isoforms [50]. In recent years, CTCF has been identified as a modulator of the alternative splicing process [17]. In this study, 86 overlapped genes were found between genes transcribing DEIs and genes in CTCF peaks, and the number of overlapped genes was too great to be accidental (Figure 7). These results suggest that CTCF may regulate alternative splicing in sheep muscle.
The CTCF-regulated alternative splicing mechanism can be divided into co-transcriptional, genomic and epigenetic mechanisms [17]. Here, our study example suggests that CTCF may regulate alternative splicing through the CTCF-mediated DNA methylation process in sheep muscle. In the current study, two DEIs, 3.1299.2 and 26.152.10, were transcribed from ANKRD23 (Table S4). In a cattle study, the ANKRD23-202 mRNA isoform, which is a splice form of the ANKRD, was defined as a DEI between the tough group and the tender group [49]. The result of a rabbit study suggested that ANKRD23 methylated in promoter and gene body regions is associated with exon skipping alternative splicing in the skeletal muscle [51]. Previous results suggest that CTCF could regulate alternative splicing by controlling DNA methylation [52,53,54]. In the current study, a CTCF peak was observed in ANKRD23 (Table S8). Taken together, CTCF may regulate ANKRD23 to produce different isoforms mediated by methylation.
As the first analysis of lamb meat tenderness based on the study of novel isoforms and alternative splicing regulation pathways, our study has its limitations. In our study, the sample size was relatively small, which may reduce the power to detect DEIs. In addition, tenderness is a parameter linked to DNA, exposition, individual characteristics, and a multitude of other factors that genetics cannot fully explain. Thus, in the future, we should increase sample size and consider more factors to investigate lamb meat tenderness. Nevertheless, our results provide the first insights about lamb meat tenderness based on the study of novel isoforms and alternative splicing regulation.

5. Conclusions

In this study, a total of 624 DEIs were identified between DDH and DHH. For example, isoform 7.524.18 transcribed from CAPN3 may be associated with meat tenderness. In addition, a total of 86 overlapped genes were found between genes with transcribed DEIs and genes in differential peaks identified by CTCF ChIP-seq. Among overlapped genes, ANKRD23 produces different isoforms that may be regulated by CTCF mediated by methylation. As preliminary research, our results identified novel isoforms associated with meat tenderness and revealed the possible regulating mechanism of alternative splicing to produce isoforms using bioinformatic analysis. In the future, more samples should be collected and more molecular experiments, e.g., TA clone, ChIP-qPCR, etc., should be implemented to improve these preliminary results.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods11081068/s1: Figure S1: Identified novel genes and isoforms: (a) Venn plot of identified genes and annotated genes; (b) frequency of the number of transcripts transcribed by a gene; (c) classification of identified isoforms; Figure S2: Classification of alternative splicing events; Table S1: Designed isoform-specific primers; Table S2: Data information of pooled isoform sequencing (Iso-seq) in sheep muscle; Table S3: Updated sheep genome annotation file; Table S4: Differentially expressed isoforms; Table S5: Gene ontology (GO) enrichment analysis of differentially expressed isoforms; Table S6: Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed isoforms; Table S7: Summarized information of chromatin immunoprecipitation sequencing (ChIP-seq); Table S8: Differential peaks between two muscle samples with divergent tenderness; Table S9: Predicated motifs in differential peaks.

Author Contributions

Conceptualization, W.S. and Z.Y.; methodology, Z.Y.; software, Z.Y.; validation, Z.Y., L.G., W.Z. and S.W.; formal analysis, Z.Y.; data curation, L.G., W.Z., S.W., X.C., X.L. and Z.Y.; writing—original draft preparation, Z.Y.; writing—review and editing, L.G., W.Z., S.W., X.L., X.C. and W.S.; visualization, Z.Y.; supervision, W.S.; project administration, W.S.; funding acquisition, W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Major New Varieties of Agricultural Projects in Jiangsu Province (PZCZ201739), Natural Science Foundation of Jiangsu Province (BK20210811), National Natural Science Foundation of China (31872333, 32172689), National Natural Science Foundation of China-CGIAR (32061143036), The Projects of Domesticated Animals Platform of the Ministry of Science, Key Research and Development Plan (modern agriculture) in Jiangsu Province (BE2018354) and Jiangsu Agricultural Science and Technology Innovation Fund (CX(18)2003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All of the raw sequencing data will be available through NCBI SRA: PRJNA745517.

Acknowledgments

We thank Mohammed Gagaoua (Teagasc Food Research Centre) for his valuable comments on the manuscript and Ruidong Xiang (The University of Melbourne) for his work in revising English phrasing.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Raza, S.H.A.; Kaster, N.; Khan, R.; Abdelnour, S.A.; El-Hack, M.E.A.; Khafaga, A.F.; Taha, A.; Ohran, H.; Swelum, A.A.; Schreurs, N.M.; et al. The role of MicroRNAs in muscle tissue development in beef cattle. Genes 2020, 11, 295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bhat, Z.F.; Morton, J.D.; Mason, S.L.; Bekhit, A.E.-D.A. Role of calpain system in meat tenderness: A review. Food Sci. Hum. Wellness 2018, 7, 196–204. [Google Scholar] [CrossRef]
  3. Jones, S.; Parr, T.; Sensky, P.; Scothern, G.; Bardsley, R.; Buttery, P. Fibre type-specific expression of p94, a skeletal muscle-specific calpain. J. Muscle Res. Cell Motil. 1999, 20, 417–424. [Google Scholar] [CrossRef]
  4. Ojima, K.; Kawabata, Y.; Nakao, H.; Nakao, K.; Doi, N.; Kitamura, F.; Ono, Y.; Hata, S.; Suzuki, H.; Kawahara, H. Dynamic distribution of muscle-specific calpain in mice has a key role in physical-stress adaptation and is impaired in muscular dystrophy. J. Clin. Investig. 2010, 120, 2672–2683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sorimachi, H.; Imajoh-Ohmi, S.; Emori, Y.; Kawasaki, H.; Ohno, S.; Minami, Y.; Suzuki, K. Molecular cloning of a novel mammalian calcium-dependent protease distinct from both m-and μ-types: Specific expression of the mRNA in skeletal muscle. J. Biol. Chem. 1989, 264, 20106–20111. [Google Scholar] [CrossRef]
  6. Basson, A.; Strydom, P.E.; van Marle-Köster, E.; Webb, E.C.; Frylinck, L. Sustained effects of muscle calpain system genotypes on tenderness pPhenotypes of South African beef Bulls during ageing up to 20 days. Animals 2022, 12, 686. [Google Scholar] [CrossRef] [PubMed]
  7. Ilian, M.A.; Morton, J.D.; Bekhit, A.E.-D.; Roberts, N.; Palmer, B.; Sorimachi, H.; Bickerstaffe, R. Effect of preslaughter feed withdrawal period on longissimus tenderness and the expression of calpains in the ovine. J. Agric. Food Chem. 2001, 49, 1990–1998. [Google Scholar] [CrossRef] [PubMed]
  8. Ilian, M.A.; Bekhit, A.E.-D.A.; Stevenson, B.; Morton, J.D.; Isherwood, P.; Bickerstaffe, R. Up-and down-regulation of longissimus tenderness parallels changes in the myofibril-bound calpain 3 protein. Meat Sci. 2004, 67, 433–445. [Google Scholar] [CrossRef]
  9. Motter, M.M.; Corva, P.M.; Soria, L.A. Expression of calpastatin isoforms in three skeletal muscles of Angus steers and their association with fiber type composition and proteolytic potential. Meat Sci. 2021, 171, 108267. [Google Scholar] [CrossRef]
  10. Rhoads, A.; Au, K.F. PacBio Sequencing and Its Applications. Genom. Proteom. Bioinform. 2015, 13, 278–289. [Google Scholar] [CrossRef] [Green Version]
  11. Chao, Q.; Gao, Z.F.; Zhang, D.; Zhao, B.G.; Dong, F.Q.; Fu, C.X.; Liu, L.J.; Wang, B.C. The developmental dynamics of the Populus stem transcriptome. Plant Biotechnol. J. 2019, 17, 206–219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Li, Y.; Dai, C.; Hu, C.; Liu, Z.; Kang, C. Global identification of alternative splicing via comparative analysis of SMRT- and Illumina-based RNA-seq in strawberry. Plant J. 2017, 90, 164–176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Yuan, Z.; Ge, L.; Sun, J.; Zhang, W.; Wang, S.; Cao, X.; Sun, W. Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat. PeerJ 2021, 9, e12454. [Google Scholar] [CrossRef]
  14. Cheng, S.; Wang, X.; Wang, Q.; Yang, L.; Shi, J.; Zhang, Q. Comparative analysis of Longissimus dorsi tissue from two sheep groups identifies differentially expressed genes related to growth, development and meat quality. Genomics 2020, 112, 3322–3330. [Google Scholar] [CrossRef]
  15. Cheng, S.; Wang, X.; Zhang, Q.; He, Y.; Zhang, X.; Yang, L.; Shi, J. Comparative transcriptome analysis identifying the different molecular genetic markers related to production performance and meat quality in longissimus dorsi tissues of MG x STH and STH Sheep. Genes 2020, 11, 183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wang, Y.; Liu, J.; Huang, B.O.; Xu, Y.M.; Li, J.; Huang, L.F.; Lin, J.; Zhang, J.; Min, Q.H.; Yang, W.M.; et al. Mechanism of alternative splicing and its regulation. Biomed. Rep. 2015, 3, 152–158. [Google Scholar] [CrossRef] [Green Version]
  17. Alharbi, A.B.; Schmitz, U.; Bailey, C.G.; Rasko, J.E.J. CTCF as a regulator of alternative splicing: New tricks for an old player. Nucleic Acids Res. 2021, 49, 7825–7838. [Google Scholar] [CrossRef]
  18. Hwang, I.H.; Devine, C.E.; Hopkins, D.L. The biochemical and physical effects of electrical stimulation on beef and sheep meat tenderness. Meat Sci. 2003, 65, 677–691. [Google Scholar] [CrossRef]
  19. Hopkins, D.L.; Toohey, E.S.; Warner, R.D.; Kerr, M.J.; van de Ven, R. Measuring the shear force of lamb meat cooked from frozen samples: Comparison of two laboratories. Anim. Prod. Sci. 2010, 50, 382–385. [Google Scholar] [CrossRef]
  20. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; Available online: https://www.R-project.org (accessed on 29 March 2022).
  21. Salmela, L.; Rivals, E. LoRDEC: Accurate and efficient long read error correction. Bioinformatics 2014, 30, 3506–3514. [Google Scholar] [CrossRef]
  22. Chen, Y.; Chen, Y.; Shi, C.; Huang, Z.; Zhang, Y.; Li, S.; Li, Y.; Ye, J.; Yu, C.; Li, Z.; et al. SOAPnuke: A MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience 2018, 7, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wu, T.D.; Watanabe, C.K. GMAP: A genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 2005, 21, 1859–1875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Foissac, S.; Sammeth, M. ASTALAVISTA: Dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 2007, 35, W297–W299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, e550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef]
  28. Xie, C.; Mao, X.; Huang, J.; Ding, Y.; Wu, J.; Dong, S.; Kong, L.; Gao, G.; Li, C.Y.; Wei, L. KOBAS 2.0: A web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011, 39, W316–W322. [Google Scholar] [CrossRef] [Green Version]
  29. Walter, W.; Sanchez-Cabo, F.; Ricote, M. GOplot: An R package for visually combining expression data with functional analysis. Bioinformatics 2015, 31, 2912–2914. [Google Scholar] [CrossRef]
  30. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: https://www.bioinformat-ics.babraham.ac.uk/projects/fastqc (accessed on 7 January 2022).
  31. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  32. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  33. Yu, G.; Wang, L.G.; He, Q.Y. ChIPseeker: An R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics 2015, 31, 2382–2383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Conway, J.R.; Lex, A.; Gehlenborg, N. UpSetR: An R package for the visualization of intersecting sets and their properties. Bioinformatics 2017, 33, 2938–2940. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhang, Y.; Liu, T.; Meyer, C.A.; Eeckhoute, J.; Johnson, D.S.; Bernstein, B.E.; Nusbaum, C.; Myers, R.M.; Brown, M.; Li, W.; et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9, R137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Heinz, S.; Benner, C.; Spann, N.; Bertolino, E.; Lin, Y.C.; Laslo, P.; Cheng, J.X.; Murre, C.; Singh, H.; Glass, C.K. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 2010, 38, 576–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Conesa, A.; Madrigal, P.; Tarazona, S.; Gomez-Cabrero, D.; Cervera, A.; McPherson, A.; Szczesniak, M.W.; Gaffney, D.J.; Elo, L.L.; Zhang, X.; et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016, 17, e13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Beiki, H.; Liu, H.; Huang, J.; Manchanda, N.; Nonneman, D.; Smith, T.P.L.; Reecy, J.M.; Tuggle, C.K. Improved annotation of the domestic pig genome through integration of Iso-Seq and RNA-seq data. BMC Genom. 2019, 20, e344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Palma, G.A.; Carranza, P.G.; Coria, M.S. El sistema proteolítico calpaina en la tenderización de la carne: Un enfoque molecular. Rev. MVZ Córdoba 2018, 23, 6523–6536. [Google Scholar] [CrossRef] [Green Version]
  40. Yuan, Z.; Sunduimijid, B.; Xiang, R.; Behrendt, R.; Knight, M.I.; Mason, B.A.; Reich, C.M.; Prowse-Wilkins, C.; Vander Jagt, C.J.; Chamberlain, A.J.; et al. Expression quantitative trait loci in sheep liver and muscle contribute to variations in meat traits. Genet. Sel. Evol. 2021, 53, e8. [Google Scholar] [CrossRef]
  41. Şirin, E.; Aksoy, Y.; Uğurlu, M.; Çiçek, Ü.; Önenç, A.; Ulutaş, Z.; Şen, U.; Kuran, M. The relationship between muscle fiber characteristics and some meat quality parameters in Turkish native sheep breeds. Small Rumin. Res. 2017, 150, 46–51. [Google Scholar] [CrossRef]
  42. Picard, B.; Gagaoua, M. Muscle Fiber Properties in Cattle and Their Relationships with Meat Qualities: An Overview. J. Agric. Food Chem. 2020, 68, 6021–6039. [Google Scholar] [CrossRef]
  43. Duris, M.-P.; Picard, B.; Geay, Y. Specificity of different anti-myosin heavy chain antibodies in bovine muscle. Meat Sci. 2000, 55, 67–78. [Google Scholar] [CrossRef]
  44. Sazili, A.Q.; Parr, T.; Sensky, P.L.; Jones, S.W.; Bardsley, R.G.; Buttery, P.J. The relationship between slow and fast myosin heavy chain content, calpastatin and meat tenderness in different ovine skeletal muscles. Meat Sci. 2005, 69, 17–25. [Google Scholar] [CrossRef] [PubMed]
  45. Yu, H.; Waddell, J.N.; Kuang, S.; Tellam, R.L.; Cockett, N.E.; Bidwell, C.A. Identification of genes directly responding to DLK1 signaling in Callipyge sheep. BMC Genom. 2018, 19, 283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Vuocolo, T.; Byrne, K.; White, J.; McWilliam, S.; Reverter, A.; Cockett, N.E.; Tellam, R.L. Identification of a gene network contributing to hypertrophy in callipyge skeletal muscle. Physiol. Genom. 2007, 28, 253–272. [Google Scholar] [CrossRef] [Green Version]
  47. Onopiuk, A.; Półtorak, A.; Wierzbicka, A. Influence of post-mortem muscle glycogen content on the quality of beef during aging. J. Vet. Res. 2016, 60, 301–307. [Google Scholar] [CrossRef] [Green Version]
  48. Gagaoua, M.; Terlouw, E.M.C.; Mullen, A.M.; Franco, D.; Warner, R.D.; Lorenzo, J.M.; Purslow, P.P.; Gerrard, D.; Hopkins, D.L.; Troy, D.; et al. Molecular signatures of beef tenderness: Underlying mechanisms based on integromics of protein biomarkers from multi-platform proteomics studies. Meat Sci. 2021, 172, 108311. [Google Scholar] [CrossRef]
  49. Muniz, M.M.M.; Fonseca, L.F.S.; Dos Santos Silva, D.B.; de Oliveira, H.R.; Baldi, F.; Chardulo, A.L.; Ferro, J.A.; Canovas, A.; de Albuquerque, L.G. Identification of novel mRNA isoforms associated with meat tenderness using RNA sequencing data in beef cattle. Meat Sci. 2021, 173, 108378. [Google Scholar] [CrossRef] [PubMed]
  50. Keren, H.; Lev-Maor, G.; Ast, G. Alternative splicing and evolution: Diversification, exon definition and function. Nat. Rev. Genet. 2010, 11, 345–355. [Google Scholar] [CrossRef] [PubMed]
  51. Li, Y.; Wang, J.; Elzo, M.A.; Fan, H.; Du, K.; Xia, S.; Shao, J.; Lai, T.; Hu, S.; Jia, X.; et al. Molecular profiling of DNA methylation and alternative splicing of genes in skeletal muscle of obese rabbits. Curr. Issues Mol. Biol. 2021, 43, 1558–1575. [Google Scholar] [CrossRef] [PubMed]
  52. Guastafierro, T.; Cecchinelli, B.; Zampieri, M.; Reale, A.; Riggio, G.; Sthandier, O.; Zupi, G.; Calabrese, L.; Caiafa, P. CCCTC-binding factor activates PARP-1 affecting DNA methylation machinery. J. Biol. Chem. 2008, 283, 21873–21880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zampieri, M.; Guastafierro, T.; Calabrese, R.; Ciccarone, F.; Bacalini, M.G.; Reale, A.; Perilli, M.; Passananti, C.; Caiafa, P. ADP-ribose polymers localized on Ctcf-Parp1-Dnmt1 complex prevent methylation of Ctcf target sites. Biochem. J. 2012, 441, 645–652. [Google Scholar] [CrossRef] [PubMed]
  54. Dubois-Chevalier, J.; Oger, F.; Dehondt, H.; Firmin, F.F.; Gheeraert, C.; Staels, B.; Lefebvre, P.; Eeckhoute, J. A dynamic CTCF chromatin binding landscape promotes DNA hydroxymethylation and transcriptional induction of adipocyte differentiation. Nucleic Acids Res. 2014, 42, 10943–10959. [Google Scholar] [CrossRef] [PubMed]
Foods 11 01068 g001 550
Figure 1. The overall experimental design. DEIs, differentially expressed isoforms.
Figure 1. The overall experimental design. DEIs, differentially expressed isoforms.
Foods 11 01068 g001
Foods 11 01068 g002 550
Figure 2. Identified isoforms in CAST. The numbers at the top of the figure represent the genome coordinates. Each yellow horizontal line represents an isoform. Each yellow vertical dash represents an exon. Stars on the right denote the novel isoforms identified in the current study. Blue arrows denote the direction of transcription. Stars in right denote novel isoforms.
Figure 2. Identified isoforms in CAST. The numbers at the top of the figure represent the genome coordinates. Each yellow horizontal line represents an isoform. Each yellow vertical dash represents an exon. Stars on the right denote the novel isoforms identified in the current study. Blue arrows denote the direction of transcription. Stars in right denote novel isoforms.
Foods 11 01068 g002
Foods 11 01068 g003 550
Figure 3. Volcano plot of differentially expressed isoforms (DEIs). X-axis denotes log2 (fold change). Y-axis denotes −log10 (adjusted p value). Red dots denote the up-regulated isoforms. Blue dots denote down-regulated isoforms.
Figure 3. Volcano plot of differentially expressed isoforms (DEIs). X-axis denotes log2 (fold change). Y-axis denotes −log10 (adjusted p value). Red dots denote the up-regulated isoforms. Blue dots denote down-regulated isoforms.
Foods 11 01068 g003
Foods 11 01068 g004 550
Figure 4. Functional annotation of differentially expressed isoforms (DEIs); (a) gene ontology (GO) enrichment analysis; X-axis denotes z-score; Y-axis denotes -log10(adjusted p value); bubble size denotes the enriched gene numbers; (b) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis; X-axis denotes −log10 (adjusted p value); Y-axis denotes enriched pathway; bubble size denotes the enriched gene numbers; bubble color denotes the adjusted p value.
Figure 4. Functional annotation of differentially expressed isoforms (DEIs); (a) gene ontology (GO) enrichment analysis; X-axis denotes z-score; Y-axis denotes -log10(adjusted p value); bubble size denotes the enriched gene numbers; (b) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis; X-axis denotes −log10 (adjusted p value); Y-axis denotes enriched pathway; bubble size denotes the enriched gene numbers; bubble color denotes the adjusted p value.
Foods 11 01068 g004
Foods 11 01068 g005 550
Figure 5. Validation of the differentially expressed isoforms by qRT-PCR; Y-axis denotes isoforms; X-axis denotes isoform relative expression level; DDH denotes Dorper × Dorper × Hu; DHH denotes Dorper × Hu × Hu. Results were presented as the mean ± SEM, * p < 0.05, significant difference.
Figure 5. Validation of the differentially expressed isoforms by qRT-PCR; Y-axis denotes isoforms; X-axis denotes isoform relative expression level; DDH denotes Dorper × Dorper × Hu; DHH denotes Dorper × Hu × Hu. Results were presented as the mean ± SEM, * p < 0.05, significant difference.
Foods 11 01068 g005
Foods 11 01068 g006 550
Figure 6. Read distribution of CTCF chromatin immunoprecipitation sequencing (ChIP-seq) among sheep genomes. According to the genome coordinates of reads, reads can be divided into six categories, including in intergenic region, intron region, promoter region, exon region, 3′ untranslated region (UTR) and 5′ untranslated region.
Figure 6. Read distribution of CTCF chromatin immunoprecipitation sequencing (ChIP-seq) among sheep genomes. According to the genome coordinates of reads, reads can be divided into six categories, including in intergenic region, intron region, promoter region, exon region, 3′ untranslated region (UTR) and 5′ untranslated region.
Foods 11 01068 g006
Foods 11 01068 g007 550
Figure 7. Distribution of differential peaks among the sheep genomes. (a) Up-regulated peaks; (b) down-regulated peaks. According to the genome coordinates of detected peaks, peaks can be divided into six categories, including in intergenic region, intron region, promoter region, exon region, 3′ untranslated region (UTR) and 5′ untranslated region. Dots in the bottom denote the peak categories. Black lines connecting the dots show peak numbers shared among several categories. The number above each bar shows the number of peaks for each category or shared category.
Figure 7. Distribution of differential peaks among the sheep genomes. (a) Up-regulated peaks; (b) down-regulated peaks. According to the genome coordinates of detected peaks, peaks can be divided into six categories, including in intergenic region, intron region, promoter region, exon region, 3′ untranslated region (UTR) and 5′ untranslated region. Dots in the bottom denote the peak categories. Black lines connecting the dots show peak numbers shared among several categories. The number above each bar shows the number of peaks for each category or shared category.
Foods 11 01068 g007
Foods 11 01068 g008 550
Figure 8. Veen plot of genes in differential peaks (left) and genes with transcribed differentially expressed isoforms (DEIs, right).
Figure 8. Veen plot of genes in differential peaks (left) and genes with transcribed differentially expressed isoforms (DEIs, right).
Foods 11 01068 g008
Table
Table 1. Meat traits of Dorper × Hu × Hu (DHH) and Dorper × Dorper × Hu (DDH) sheep.
Table 1. Meat traits of Dorper × Hu × Hu (DHH) and Dorper × Dorper × Hu (DDH) sheep.
Meat Trait (Unit)Crossing TypeStatistic (t)p Value
DDH 1 (n = 3)DHH 2 (n = 3)
Live weight (kg)40.80 ± 4.5741.50 ± 0.95−0.27200.8094
Carcass weight (kg) 322.00 ± 2.2723.70 ± 0.64−1.27170.3161
Shear force (N) 383.20 ± 13.6053.30 ± 8.393.25580.0406
1 DDH denotes Dorper × Dorper × Hu. 2 DHH denotes Dorper × Hu × Hu. 3 These traits were measured on hot carcass.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yuan, Z.; Ge, L.; Zhang, W.; Lv, X.; Wang, S.; Cao, X.; Sun, W. Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data. Foods 2022, 11, 1068. https://doi.org/10.3390/foods11081068

AMA Style

Yuan Z, Ge L, Zhang W, Lv X, Wang S, Cao X, Sun W. Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data. Foods. 2022; 11(8):1068. https://doi.org/10.3390/foods11081068

Chicago/Turabian Style

Yuan, Zehu, Ling Ge, Weibo Zhang, Xiaoyang Lv, Shanhe Wang, Xiukai Cao, and Wei Sun. 2022. "Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data" Foods 11, no. 8: 1068. https://doi.org/10.3390/foods11081068

APA Style

Yuan, Z., Ge, L., Zhang, W., Lv, X., Wang, S., Cao, X., & Sun, W. (2022). Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data. Foods, 11(8), 1068. https://doi.org/10.3390/foods11081068

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop