{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,20]],"date-time":"2026-05-20T17:17:11Z","timestamp":1779297431067,"version":"3.51.4"},"reference-count":55,"publisher":"MDPI AG","issue":"3","license":[{"start":{"date-parts":[[2022,3,11]],"date-time":"2022-03-11T00:00:00Z","timestamp":1646956800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100001809","name":"National Natural Science Foundation of China","doi-asserted-by":"publisher","award":["31472065"],"award-info":[{"award-number":["31472065"]}],"id":[{"id":"10.13039\/501100001809","id-type":"DOI","asserted-by":"publisher"}]},{"name":"the National key Research and Develop-ment Project of China","award":["2021YFF1000700"],"award-info":[{"award-number":["2021YFF1000700"]}]},{"name":"the Program for Changjiang Scholar and Innovation Research Team in University","award":["IRT_15R62"],"award-info":[{"award-number":["IRT_15R62"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Genes"],"abstract":"<jats:p>Paratuberculosis is a major endemic disease caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection and leads to huge economic loss in the dairy sector worldwide. Alternative splicing (AS) events, playing indispensable regulatory roles in many protein functions and biological pathways, are shown to be associated with complex traits and diseases. In this study, by integrating the RNA sequencing (RNA-seq) data of 24 samples from three tissues (peripheral blood, jejunum and salivary gland) of Holstein cows, we obtained 2,706,541,696 uniquely mapped reads in total that represented 12,870 expressed genes, and detected 4285 differentially expressed genes (DEGs) between MAP-infected and healthy cows (p &lt; 0.05). Of them, 92 differentially expressed splicing factors (DESFs) were included. Further, 119, 150 and 68 differential alternative splicing (DAS) events between MAP-infected and healthy cows were identified in peripheral blood, jejunum and salivary glands, respectively. Of note, six DAS events were highly and significantly correlated with the DESFs (R2 &gt; 0.9; p &lt; 0.01), and their corresponding genes (COPI coat complex subunit gamma 2gene (COPG2), kinesin family member 2C gene (KIF2C), exocyst complex component 7 (EXOC7), Rab9 effector protein with kelch motifs gene (RABEPK), deoxyribonuclease 1 gene (DNASE1) and early endosome antigen 1gene (EEA1)) were significantly enriched in immune response such as vesicle-mediated transport, regulation of acute inflammatory response and tuberculosis through gene ontology (GO) and KEGG analysis. KS test showed that the DAS events in the EXOC7 and KIF2C genes indeed displayed differences between MAP-infected cows and healthy cows. The DAS in EXOC7 might produce a new protein sequence with lack of 23 amino acids, and the DAS in KIF2C induced a stop codon of premature occurrence and resulted in a lack of functional domain. In summary, this study identified the DAS events and corresponding genes related to MAP-infection base on the RNA-seq data from multiple tissues of Holstein cows, providing novel insights into the regulatory mechanisms underpinning paratuberculosis in dairy cattle.<\/jats:p>","DOI":"10.3390\/genes13030497","type":"journal-article","created":{"date-parts":[[2022,3,11]],"date-time":"2022-03-11T12:58:36Z","timestamp":1647003516000},"page":"497","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":11,"title":["Identification of Alternative Splicing Events Associated with Paratuberculosis in Dairy Cattle Using Multi-Tissue RNA Sequencing Data"],"prefix":"10.3390","volume":"13","author":[{"given":"Houcheng","family":"Li","sequence":"first","affiliation":[{"name":"National Engineering Laboratory of Animal Breeding, Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jinfeng","family":"Huang","sequence":"additional","affiliation":[{"name":"National Engineering Laboratory of Animal Breeding, Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Junnan","family":"Zhang","sequence":"additional","affiliation":[{"name":"National Engineering Laboratory of Animal Breeding, Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yahui","family":"Gao","sequence":"additional","affiliation":[{"name":"National Engineering Laboratory of Animal Breeding, Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Bo","family":"Han","sequence":"additional","affiliation":[{"name":"National Engineering Laboratory of Animal Breeding, Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7512-5946","authenticated-orcid":false,"given":"Dongxiao","family":"Sun","sequence":"additional","affiliation":[{"name":"National Engineering Laboratory of Animal Breeding, Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2022,3,11]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"217","DOI":"10.1016\/S0021-9975(97)80001-1","article-title":"The pathology and pathogenesis of paratuberculosis in ruminants and other species","volume":"116","author":"Clarke","year":"1997","journal-title":"J. Comp. Pathol."},{"key":"ref_2","doi-asserted-by":"crossref","unstructured":"Stevenson, K., Alvarez, J., Bakker, D., Biet, F., de Juan, L., Denham, S., Dimareli, Z., Dohmann, K., Gerlach, G.F., and Heron, I. (2009). Occurrence of Mycobacterium avium subspecies paratuberculosis across host species and European countries with evidence for transmission between wildlife and domestic ruminants. BMC Microbiol., 9.","DOI":"10.1186\/1471-2180-9-212"},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"5019","DOI":"10.3168\/jds.2014-9241","article-title":"Invited review: The economic impact and control of paratuberculosis in cattle","volume":"98","author":"Garcia","year":"2015","journal-title":"J. Dairy Sci."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"3513","DOI":"10.3168\/jds.2009-2742","article-title":"Effect of Johne\u2019s disease status on reproduction and culling in dairy cattle","volume":"93","author":"Smith","year":"2010","journal-title":"J. Dairy Sci."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"389","DOI":"10.1111\/j.1365-2052.2007.01627.x","article-title":"Identification of a QTL on BTA20 affecting susceptibility to Mycobacterium avium ssp. paratuberculosis infection in US Holsteins","volume":"38","author":"Gonda","year":"2007","journal-title":"Anim. Genet."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"176","DOI":"10.1016\/j.ygeno.2010.01.001","article-title":"A principal component regression based genome wide analysis approach reveals the presence of a novel QTL on BTA7 for MAP resistance in holstein cattle","volume":"95","author":"Pant","year":"2010","journal-title":"Genomics"},{"key":"ref_7","first-page":"1080","article-title":"Identification of loci associated with susceptibility to subspecies tissue infection in cattle","volume":"95","author":"Kiser","year":"2017","journal-title":"J. Anim. Sci."},{"key":"ref_8","doi-asserted-by":"crossref","unstructured":"Gao, Y., Jiang, J., Yang, S., Cao, J., Han, B., Wang, Y., Zhang, Y., Yu, Y., Zhang, S., and Zhang, Q. (2018). Genome-wide association study of Mycobacterium avium subspecies Paratuberculosis infection in Chinese Holstein. BMC Genom., 19.","DOI":"10.1186\/s12864-018-5385-3"},{"key":"ref_9","unstructured":"Gao, Y. (2017). Identification of Susceptibility\/Resistance and Milk Component Trait Functional Genes of Para-Tuberculosis in Dairy Cows Based on Omics Technique, China Agricultural University."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"359","DOI":"10.1146\/annurev-biochem-091719-064225","article-title":"RNA Splicing by the Spliceosome","volume":"89","author":"Wilkinson","year":"2020","journal-title":"Annu. Rev. Biochem."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"741","DOI":"10.1038\/nrm2777","article-title":"Mechanisms of alternative splicing regulation: Insights from molecular and genomics approaches","volume":"10","author":"Chen","year":"2009","journal-title":"Nat. Rev. Mol. Cell Biol."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"437","DOI":"10.1038\/nrm.2017.27","article-title":"Alternative splicing as a regulator of development and tissue identity","volume":"18","author":"Baralle","year":"2017","journal-title":"Nat. Rev. Mol. Cell Biol."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"D1055","DOI":"10.1093\/nar\/gkr1171","article-title":"SpliceDisease database: Linking RNA splicing and disease","volume":"40","author":"Wang","year":"2012","journal-title":"Nucleic Acids Res."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"68","DOI":"10.1016\/j.tig.2018.10.002","article-title":"RNA Splicing and Disease: Animal Models to Therapies","volume":"35","author":"Montes","year":"2019","journal-title":"Trends Genet. TIG"},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"43198","DOI":"10.1038\/srep43198","article-title":"Increased 4R tau expression and behavioural changes in a novel MAPT-N296H genomic mouse model of tauopathy","volume":"7","author":"Wobst","year":"2017","journal-title":"Sci. Rep."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"981","DOI":"10.1172\/JCI72723","article-title":"ALS-associated mutation FUS-R521C causes DNA damage and RNA splicing defects","volume":"124","author":"Qiu","year":"2014","journal-title":"J. Clin. Investig."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"155","DOI":"10.4049\/jimmunol.1102042","article-title":"Overexpression of the Ctla-4 isoform lacking exons 2 and 3 causes autoimmunity","volume":"188","author":"Liu","year":"2012","journal-title":"J. Immunol."},{"key":"ref_18","doi-asserted-by":"crossref","unstructured":"Wang, X.G., Ju, Z.H., Hou, M.H., Jiang, Q., Yang, C.H., Zhang, Y., Sun, Y., Li, R.L., Wang, C.F., and Zhong, J.F. (2016). Deciphering Transcriptome and Complex Alternative Splicing Transcripts in Mammary Gland Tissues from Cows Naturally Infected with Staphylococcus aureus Mastitis. PLoS ONE, 11.","DOI":"10.1371\/journal.pone.0159719"},{"key":"ref_19","first-page":"4159","article-title":"Production of an aberrant splice variant of CCL5 is not caused by genetic mutation in the mammary glands of mastitis-infected Holstein cows","volume":"19","author":"Yang","year":"2019","journal-title":"Mol. Med. Rep."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"24964","DOI":"10.1038\/srep24964","article-title":"Altered microRNA expression and pre-mRNA splicing events reveal new mechanisms associated with early stage Mycobacterium avium subspecies paratuberculosis infection","volume":"6","author":"Liang","year":"2016","journal-title":"Sci. Rep."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"3445","DOI":"10.1093\/nar\/22.17.3445","article-title":"The European Bioinformatics Institute (EBI) databases","volume":"22","author":"Emmert","year":"1994","journal-title":"Nucleic Acids Res."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"7274","DOI":"10.3168\/jds.2017-13264","article-title":"Short communication: Heritability estimates for susceptibility to Mycobacterium avium ssp. paratuberculosis infection in Chinese Holstein cattle","volume":"101","author":"Gao","year":"2018","journal-title":"J. Dairy Sci."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"14845","DOI":"10.1038\/s41598-019-51328-0","article-title":"RNA-Seq analysis of ileocecal valve and peripheral blood from Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis revealed dysregulation of the CXCL8\/IL8 signaling pathway","volume":"9","author":"Canive","year":"2019","journal-title":"Sci. Rep."},{"key":"ref_24","doi-asserted-by":"crossref","unstructured":"Mallikarjunappa, S., Adnane, M., Cormican, P., Karrow, N.A., and Meade, K.G. (2019). Characterization of the bovine salivary gland transcriptome associated with Mycobacterium avium subsp. paratuberculosis experimental challenge. BMC Genom., 20.","DOI":"10.1186\/s12864-019-5845-4"},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"295","DOI":"10.1007\/978-1-4939-7286-9_23","article-title":"RNA-Seq Data Analysis from Raw Data Quality Control to Differential Expression Analysis","volume":"1669","author":"Qi","year":"2017","journal-title":"Methods Mol. Biol."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"3","DOI":"10.1093\/gigascience\/giaa021","article-title":"De novo assembly of the cattle reference genome with single-molecule sequencing","volume":"9","author":"Rosen","year":"2020","journal-title":"GigaScience"},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"D884","DOI":"10.1093\/nar\/gkaa942","article-title":"Ensembl 2021","volume":"49","author":"Howe","year":"2021","journal-title":"Nucleic Acids Res."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1093\/bioinformatics\/bts635","article-title":"STAR: Ultrafast universal RNA-seq aligner","volume":"29","author":"Dobin","year":"2013","journal-title":"Bioinformatics"},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"923","DOI":"10.1093\/bioinformatics\/btt656","article-title":"featureCounts: An efficient general purpose program for assigning sequence reads to genomic features","volume":"30","author":"Liao","year":"2014","journal-title":"Bioinformatics"},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol."},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"1523","DOI":"10.1038\/s41467-019-09234-6","article-title":"Metascape provides a biologist-oriented resource for the analysis of systems-level datasets","volume":"10","author":"Zhou","year":"2019","journal-title":"Nat. Commun."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1089\/omi.2011.0118","article-title":"clusterProfiler: An R package for comparing biological themes among gene clusters","volume":"16","author":"Yu","year":"2012","journal-title":"Omics"},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"E5593","DOI":"10.1073\/pnas.1419161111","article-title":"rMATS: Robust and flexible detection of differential alternative splicing from replicate RNA-Seq data","volume":"111","author":"Shen","year":"2014","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"151","DOI":"10.1038\/s41588-017-0004-9","article-title":"Annotation-free quantification of RNA splicing using LeafCutter","volume":"50","author":"Li","year":"2018","journal-title":"Nat. Genet."},{"key":"ref_35","doi-asserted-by":"crossref","unstructured":"Jin, Y.J., Byun, S., Han, S., Chamberlin, J., Kim, D., Kim, M.J., and Lee, Y. (2019). Differential alternative splicing regulation among hepatocellular carcinoma with different risk factors. BMC Med. Genom., 12.","DOI":"10.1186\/s12920-019-0635-z"},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"83","DOI":"10.1016\/j.canlet.2020.11.043","article-title":"Splicing factors: Insights into their regulatory network in alternative splicing in cancer","volume":"501","author":"Du","year":"2021","journal-title":"Cancer Lett."},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"D125","DOI":"10.1093\/nar\/gks997","article-title":"SpliceAid-F: A database of human splicing factors and their RNA-binding sites","volume":"41","author":"Giulietti","year":"2013","journal-title":"Nucleic Acids Res."},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"2498","DOI":"10.1101\/gr.1239303","article-title":"Cytoscape: A software environment for integrated models of biomolecular interaction networks","volume":"13","author":"Shannon","year":"2003","journal-title":"Genome Res."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"431","DOI":"10.1093\/bioinformatics\/btq675","article-title":"Cytoscape 2.8: New features for data integration and network visualization","volume":"27","author":"Smoot","year":"2011","journal-title":"Bioinformatics"},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"68","DOI":"10.1080\/01621459.1951.10500769","article-title":"The Kolmogorov-Smirnov test for goodness of fit","volume":"46","author":"Massey","year":"1951","journal-title":"J. Am. Stat. Assoc."},{"key":"ref_41","unstructured":"(2018). NCBI Resource Coordinators: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res., 46(D1), D8\u2013D13."},{"key":"ref_42","doi-asserted-by":"crossref","first-page":"231","DOI":"10.1093\/nar\/28.1.231","article-title":"SMART: A web-based tool for the study of genetically mobile domains","volume":"28","author":"Schultz","year":"2000","journal-title":"Nucleic Acids Res."},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"3960","DOI":"10.1128\/IAI.00575-13","article-title":"No holes barred: Invasion of the intestinal mucosa by Mycobacterium avium subsp. paratuberculosis","volume":"81","author":"Bannantine","year":"2013","journal-title":"Infect. Immun."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"54","DOI":"10.1186\/1297-9716-45-54","article-title":"From mouth to macrophage: Mechanisms of innate immune subversion by Mycobacterium avium subsp. paratuberculosis","volume":"45","author":"Arsenault","year":"2014","journal-title":"Vet. Res."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"2861","DOI":"10.1128\/IAI.00339-13","article-title":"Divergent immune responses to Mycobacterium avium subsp. paratuberculosis infection correlate with kinome responses at the site of intestinal infection","volume":"81","author":"Trost","year":"2013","journal-title":"Infect. Immun."},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"167","DOI":"10.1016\/S0034-5288(18)34970-1","article-title":"The Pathogenesis of Experimental Johne\u2019s Disease in Calves","volume":"2","author":"Payne","year":"1961","journal-title":"Res. Vet. Sci."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"253","DOI":"10.1016\/j.rvsc.2005.07.007","article-title":"Tissue predilection sites and effect of dose on Mycobacterium avium subs. paratuberculosis organism recovery in a short-term bovine experimental oral infection model","volume":"80","author":"Sweeney","year":"2006","journal-title":"Res. Vet. Sci."},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"4295","DOI":"10.3390\/ijms13044295","article-title":"Salivary defense proteins: Their network and role in innate and acquired oral immunity","volume":"13","author":"Hermann","year":"2012","journal-title":"Int. J. Mol. Sci."},{"key":"ref_49","doi-asserted-by":"crossref","first-page":"2740","DOI":"10.1128\/IAI.67.6.2740-2745.1999","article-title":"Production of \u03b2-defensin antimicrobial peptides by the oral mucosa and salivary glands","volume":"67","author":"Mathews","year":"1999","journal-title":"Infect. Immun."},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"R42","DOI":"10.1186\/gb-2009-10-4-r42","article-title":"A whole-genome assembly of the domestic cow, Bos taurus","volume":"10","author":"Zimin","year":"2009","journal-title":"Genome Biol."},{"key":"ref_51","doi-asserted-by":"crossref","first-page":"19988","DOI":"10.1074\/jbc.RA119.010821","article-title":"Inducible Exoc7\/Exo70 knockout reveals a critical role of the exocyst in insulin-regulated GLUT4 exocytosis","volume":"294","author":"Wang","year":"2019","journal-title":"J. Biol. Chem."},{"key":"ref_52","doi-asserted-by":"crossref","first-page":"85","DOI":"10.1016\/j.ccell.2018.06.007","article-title":"PTBP1-Mediated Alternative Splicing Regulates the Inflammatory Secretome and the Pro-tumorigenic Effects of Senescent Cells","volume":"34","author":"Georgilis","year":"2018","journal-title":"Cancer Cell"},{"key":"ref_53","doi-asserted-by":"crossref","first-page":"342","DOI":"10.3389\/fgene.2020.00342","article-title":"Identification of Crucial Genes Associated With Immune Cell Infiltration in Hepatocellular Carcinoma by Weighted Gene Co-expression Network Analysis","volume":"11","author":"Wang","year":"2020","journal-title":"Front. Genet."},{"key":"ref_54","doi-asserted-by":"crossref","first-page":"1434856","DOI":"10.1155\/2021\/1434856","article-title":"KIF2C Is a Novel Prognostic Biomarker and Correlated with Immune Infiltration in Endometrial Cancer","volume":"2021","author":"An","year":"2021","journal-title":"Stem Cells Int."},{"key":"ref_55","doi-asserted-by":"crossref","first-page":"2052","DOI":"10.1093\/bib\/bbz126","article-title":"Systematic evaluation of differential splicing tools for RNA-seq studies","volume":"21","author":"Mehmood","year":"2020","journal-title":"Brief. Bioinform."}],"container-title":["Genes"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2073-4425\/13\/3\/497\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,10]],"date-time":"2025-10-10T22:34:54Z","timestamp":1760135694000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2073-4425\/13\/3\/497"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2022,3,11]]},"references-count":55,"journal-issue":{"issue":"3","published-online":{"date-parts":[[2022,3]]}},"alternative-id":["genes13030497"],"URL":"https:\/\/doi.org\/10.3390\/genes13030497","relation":{},"ISSN":["2073-4425"],"issn-type":[{"value":"2073-4425","type":"electronic"}],"subject":[],"published":{"date-parts":[[2022,3,11]]}}}