Abstract: Severe equine asthma is a chronic inflammatory disease of the lung in horses similar to low-Th2 late-onset asthma in humans. This study aimed to determine the utility of RNA-Seq to call gene sequence variants, and to identify sequence variants of potential relevance to the pathogenesis of asthma. Methods: RNA-Seq data were generated from endobronchial biopsies collected from six asthmatic and seven non-asthmatic horses before and after challenge (26 samples total). Sequences were aligned to the equine genome with Spliced Transcripts Alignment to Reference software. Read preparation for sequence variant calling was performed with Picard tools and Genome Analysis Toolkit (GATK). Sequence variants were called and filtered using GATK and Ensembl Variant Effect Predictor (VEP) tools, and two RNA-Seq predicted sequence variants were investigated with both PCR and Sanger sequencing. Supplementary analysis of novel sequence variant selection with VEP was based on a score of <0.01 predicted with Sorting Intolerant from Tolerant software, missense nature, location within the protein coding sequence and presence in all asthmatic individuals. For select variants, effect on protein function was assessed with Polymorphism Phenotyping 2 and screening for non-acceptable polymorphism 2 software. Sequences were aligned and 3D protein structures predicted with Geneious software. Difference in allele frequency between the groups was assessed using a Pearson's Chi-squared test with Yates' continuity correction, and difference in genotype frequency was calculated using the Fisher's exact test for count data. Results: RNA-Seq variant calling and filtering correctly identified substitution variants in and . Sanger sequencing confirmed that the substitution was appropriately identified in all 26 samples while the substitution was identified correctly in 24 of 26 samples. These variants of uncertain significance had substitutions that were predicted to result in loss of function and to be non-neutral. Amino acid substitutions projected no change of hydrophobicity and isoelectric point in PACRG, and a change in both for RTTN. For , no difference in allele frequency between the two groups was detected but a higher proportion of asthmatic horses had the altered allele compared to non-asthmatic animals. Conclusions: RNA-Seq was sensitive and specific for calling gene sequence variants in this disease model. Even moderate coverage (<10-20 counts per million) yielded correct identification in 92% of samples, suggesting RNA-Seq may be suitable to detect sequence variants in low coverage samples. The impact of amino acid alterations in PACRG and RTTN proteins, and possible association of the sequence variants with asthma, is of uncertain significance, but their role in ciliary function may be of future interest.
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.
This research aimed to identify genetic variants potentially associated with severe equine asthma using RNA-Seq, a technology that involves analyzing sequences of RNA. This type of asthma in horses is a chronic lung disease similar to a type of late-onset asthma in humans. The study targeted variant sequence identification to shed light on the origins of asthma.
Methodology
The study processed RNA-Seq data from lung biopsies of both asthmatic and non-asthmatic horses, both before and after exposure to asthma triggers.
The sequences were then matched to the equine genome using a Spliced Transcripts Alignment to Reference software.
Read preparation for sequence variant calling was then performed using Picard tools and Genome Analysis Toolkit (GATK).
The sequence variants were identified and filtered using GATK and Ensembl Variant Effect Predictor (VEP) tools.
Two RNA-Seq predicted sequence variants were further investigated utilizing both PCR and Sanger sequencing.
Additional analysis of novel sequence variant selection was conducted with VEP, relying on a score of under 0.01 predicted with Sorting Intolerant from Tolerant software, and considering other factors like missense nature, location within the protein coding sequence, and presence in all asthmatic individuals.
The consequences of select variants in protein functionality were assessed with Polymorphism Phenotyping 2 and screening for non-acceptable polymorphism 2 software.
Sequences were aligned and 3D protein structures predicted with Geneious software.
Difference of frequency in alleles between the groups was calculated using a Pearson’s Chi-squared test with Yates’ continuity correction. The difference of genotype frequency was calculated using the Fisher’s exact test for count data.
Results
The RNA-Seq variant calling and filtering method successfully identified two substitution variants.
These sequences were then confirmed useful by Sanger sequencing. Specifically, one variant was established in all samples, and the other was set in 24 out of 26 samples.
The identified variants are believed to cause loss of function and not be neutral, expected to impact the development of severe equine asthma.
Amino acid substitutions predicted no change for PACRG protein, but a change for RTTN.
One of the variants did not display any difference between the two groups, but a higher proportion of asthmatic horses carried the altered allele.
Conclusions
RNA-Seq proved both sensitive and specific for calling gene sequence variants in this disease model.
The methodology had a 92% correct identification rate even with moderate coverage. This suggests the method’s possible suitability to detect sequence variants in samples with low coverage.
Although the effects of the amino acid changes in the PACRG and RTTN proteins are still uncertain, their connection with ciliary function may be important for future research.
Cite This Article
APA
Tessier L, Côté O, Bienzle D.
(2018).
Sequence variant analysis of RNA sequences in severe equine asthma.
PeerJ, 6, e5759.
https://doi.org/10.7717/peerj.5759
Department of Pathobiology, University of Guelph, Guelph, ON, Canada.
BenchSci, Toronto, ON, Canada.
Côté, Olivier
Department of Pathobiology, University of Guelph, Guelph, ON, Canada.
BioAssay Works, Ijamsville, MD, USA.
Bienzle, Dorothee
Department of Pathobiology, University of Guelph, Guelph, ON, Canada.
Conflict of Interest Statement
Olivier Côté is employed by BioAssay Works LLC, Ijamsville, MD, USA. Laurence Tessier became an employee of BenchSci, Toronto, ON after completion of this study. Dorothee Bienzle is an Academic Editor for PeerJ.
References
This article includes 94 references
Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2.. Curr Protoc Hum Genet 2013 Jan;Chapter 7:Unit7.20.
Betts MJ, Russell RB. Amino acid properties and consequences of substitutions. In: Barnes MR, Gray IC, editors. Bioinformatics for Geneticists. Hoboken: John Wiley & Sons; 2003. pp. 289–316.
Bisgaard H, Bønnelykke K, Sleiman PM, Brasholt M, Chawes B, Kreiner-Møller E, Stage M, Kim C, Tavendale R, Baty F, Pipper CB, Palmer CN, Hakonarsson H. Chromosome 17q21 gene variants are associated with asthma and exacerbations but not atopy in early childhood.. Am J Respir Crit Care Med 2009 Feb 1;179(3):179-85.
Chen HY, Wu CT, Tang CC, Lin YN, Wang WJ, Tang TK. Human microcephaly protein RTTN interacts with STIL and is required to build full-length centrioles.. Nat Commun 2017 Aug 15;8(1):247.
Choudhry S, Taub M, Mei R, Rodriguez-Santana J, Rodriguez-Cintron W, Shriver MD, Ziv E, Risch NJ, Burchard EG. Genome-wide screen for asthma in Puerto Ricans: evidence for association with 5q23 region.. Hum Genet 2008 Jun;123(5):455-68.
Dawe HR, Farr H, Portman N, Shaw MK, Gull K. The Parkin co-regulated gene product, PACRG, is an evolutionarily conserved axonemal protein that functions in outer-doublet microtubule morphogenesis.. J Cell Sci 2005 Dec 1;118(Pt 23):5421-30.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner.. Bioinformatics 2013 Jan 1;29(1):15-21.
Dong C, Wei P, Jian X, Gibbs R, Boerwinkle E, Wang K, Liu X. Comparison and integration of deleteriousness prediction methods for nonsynonymous SNVs in whole exome sequencing studies.. Hum Mol Genet 2015 Apr 15;24(8):2125-37.
Faisst AM, Alvarez-Bolado G, Treichel D, Gruss P. Rotatin is a novel gene required for axial rotation and left-right specification in mouse embryos.. Mech Dev 2002 Apr;113(1):15-28.
Flanagan SE, Patch AM, Ellard S. Using SIFT and PolyPhen to predict loss-of-function and gain-of-function mutations.. Genet Test Mol Biomarkers 2010 Aug;14(4):533-7.
Galanter J, Choudhry S, Eng C, Nazario S, Rodríguez-Santana JR, Casal J, Torres-Palacios A, Salas J, Chapela R, Watson HG, Meade K, LeNoir M, Rodríguez-Cintrón W, Avila PC, Burchard EG. ORMDL3 gene is associated with asthma in three ethnically diverse populations.. Am J Respir Crit Care Med 2008 Jun 1;177(11):1194-200.
Gavino C, Richard S. Patched1 haploinsufficiency impairs ependymal cilia function of the quaking viable mice, leading to fatal hydrocephalus.. Mol Cell Neurosci 2011 Jun;47(2):100-7.
Gheber L, Priel Z. On metachronism in ciliary systems: a model describing the dependence of the metachronal wave properties on the intrinsic ciliary parameters.. Cell Motil Cytoskeleton 1990;16(3):167-81.
Guo Y, Zhao S, Sheng Q, Samuels DC, Shyr Y. The discrepancy among single nucleotide variants detected by DNA and RNA high throughput sequencing data.. BMC Genomics 2017 Oct 3;18(Suppl 6):690.
Hancock DB, Romieu I, Shi M, Sienra-Monge JJ, Wu H, Chiu GY, Li H, del Rio-Navarro BE, Willis-Owen SA, Weiss ST, Raby BA, Gao H, Eng C, Chapela R, Burchard EG, Tang H, Sullivan PF, London SJ. Genome-wide association study implicates chromosome 9q21.31 as a susceptibility locus for asthma in mexican children.. PLoS Genet 2009 Aug;5(8):e1000623.
Huangfu D, Liu A, Rakeman AS, Murcia NS, Niswander L, Anderson KV. Hedgehog signalling in the mouse requires intraflagellar transport proteins.. Nature 2003 Nov 6;426(6962):83-7.
Ikeda K, Ikeda T, Morikawa K, Kamiya R. Axonemal localization of Chlamydomonas PACRG, a homologue of the human Parkin-coregulated gene product.. Cell Motil Cytoskeleton 2007 Nov;64(11):814-21.
Imai Y, Soda M, Murakami T, Shoji M, Abe K, Takahashi R. A product of the human gene adjacent to parkin is a component of Lewy bodies and suppresses Pael receptor-induced cell death.. J Biol Chem 2003 Dec 19;278(51):51901-10.
Kazemi-Esfarjani P, Beitel LK, Trifiro M, Kaufman M, Rennie P, Sheppard P, Matusik R, Pinsky L. Substitution of valine-865 by methionine or leucine in the human androgen receptor causes complete or partial androgen insensitivity, respectively with distinct androgen receptor phenotypes.. Mol Endocrinol 1993 Jan;7(1):37-46.
Kheradmand Kia S, Verbeek E, Engelen E, Schot R, Poot RA, de Coo IF, Lequin MH, Poulton CJ, Pourfarzad F, Grosveld FG, Brehm A, de Wit MC, Oegema R, Dobyns WB, Verheijen FW, Mancini GM. RTTN mutations link primary cilia function to organization of the human cerebral cortex.. Am J Hum Genet 2012 Sep 7;91(3):533-40.
Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm.. Nat Protoc 2009;4(7):1073-81.
Lechtreck KF, Luro S, Awata J, Witman GB. HA-tagging of putative flagellar proteins in Chlamydomonas reinhardtii identifies a novel protein of intraflagellar transport complex B.. Cell Motil Cytoskeleton 2009 Aug;66(8):469-82.
Li X, Howard TD, Moore WC, Ampleford EJ, Li H, Busse WW, Calhoun WJ, Castro M, Chung KF, Erzurum SC, Fitzpatrick AM, Gaston B, Israel E, Jarjour NN, Teague WG, Wenzel SE, Peters SP, Hawkins GA, Bleecker ER, Meyers DA. Importance of hedgehog interacting protein and other lung function genes in asthma.. J Allergy Clin Immunol 2011 Jun;127(6):1457-65.
Li X, Howard TD, Zheng SL, Haselkorn T, Peters SP, Meyers DA, Bleecker ER. Genome-wide association study of asthma identifies RAD50-IL13 and HLA-DR/DQ regions.. J Allergy Clin Immunol 2010 Feb;125(2):328-335.e11.
Li Y, Xiao X, Ji X, Liu B, Amos CI. RNA-seq analysis of lung adenocarcinomas reveals different gene expression profiles between smoking and nonsmoking patients.. Tumour Biol 2015 Nov;36(11):8993-9003.
Madore AM, Tremblay K, Hudson TJ, Laprise C. Replication of an association between 17q21 SNPs and asthma in a French-Canadian familial collection.. Hum Genet 2008 Feb;123(1):93-5.
Majewski J, Ott J. Amino acid substitutions in the human genome: evolutionary implications of single nucleotide polymorphisms.. Gene 2003 Feb 27;305(2):167-73.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data.. Genome Res 2010 Sep;20(9):1297-303.
McLaren W, Pritchard B, Rios D, Chen Y, Flicek P, Cunningham F. Deriving the consequences of genomic variants with the Ensembl API and SNP Effect Predictor.. Bioinformatics 2010 Aug 15;26(16):2069-70.
Mizuno K, Dymek EE, Smith EF. Microtubule binding protein PACRG plays a role in regulating specific ciliary dyneins during microtubule sliding.. Cytoskeleton (Hoboken) 2016 Dec;73(12):703-711.
Moffatt MF, Gut IG, Demenais F, Strachan DP, Bouzigon E, Heath S, von Mutius E, Farrall M, Lathrop M, Cookson WOCM. A large-scale, consortium-based genomewide association study of asthma.. N Engl J Med 2010 Sep 23;363(13):1211-1221.
Moffatt MF, Kabesch M, Liang L, Dixon AL, Strachan D, Heath S, Depner M, von Berg A, Bufe A, Rietschel E, Heinzmann A, Simma B, Frischer T, Willis-Owen SA, Wong KC, Illig T, Vogelberg C, Weiland SK, von Mutius E, Abecasis GR, Farrall M, Gut IG, Lathrop GM, Cookson WO. Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthma.. Nature 2007 Jul 26;448(7152):470-3.
Quinn EM, Cormican P, Kenny EM, Hill M, Anney R, Gill M, Corvin AP, Morris DW. Development of strategies for SNP detection in RNA-seq data: application to lymphoblastoid cell lines and evaluation using 1000 Genomes data.. PLoS One 2013;8(3):e58815.
Scholtens S, Postma DS, Moffatt MF, Panasevich S, Granell R, Henderson AJ, Melén E, Nyberg F, Pershagen G, Jarvis D, Ramasamy A, Wjst M, Svanes C, Bouzigon E, Demenais F, Kauffmann F, Siroux V, von Mutius E, Ege MJ, Braun-Fahrländer C, Genuneit J, Brunekreef B, Smit HA, Wijga AH, Kerkhof M, Curjuric I, Imboden M, Thun GA, Probst-Hensch N, Freidin MB, Bragina EIu, Deev IA, Puzyrev VP, Daley D, Park J, Becker A, Chan-Yeung M, Kozyrskyj AL, Pare P, Marenholz I, Lau S, Keil T, Lee YA, Kabesch M, Wijmenga C, Franke L, Nolte IM, Vonk J, Kumar A, Farrall M, Cookson WO, Strachan DP, Koppelman GH, Boezen HM. Novel childhood asthma genes interact with in utero and early-life tobacco smoke exposure.. J Allergy Clin Immunol 2014 Mar;133(3):885-8.
Sheng Q, Zhao S, Li CI, Shyr Y, Guo Y. Practicability of detecting somatic point mutation from RNA high throughput sequencing data.. Genomics 2016 May;107(5):163-9.
Sleiman PM, Flory J, Imielinski M, Bradfield JP, Annaiah K, Willis-Owen SA, Wang K, Rafaels NM, Michel S, Bonnelykke K, Zhang H, Kim CE, Frackelton EC, Glessner JT, Hou C, Otieno FG, Santa E, Thomas K, Smith RM, Glaberson WR, Garris M, Chiavacci RM, Beaty TH, Ruczinski I, Orange JS, Allen J, Spergel JM, Grundmeier R, Mathias RA, Christie JD, von Mutius E, Cookson WO, Kabesch M, Moffatt MF, Grunstein MM, Barnes KC, Devoto M, Magnusson M, Li H, Grant SF, Bisgaard H, Hakonarson H. Variants of DENND1B associated with asthma in children.. N Engl J Med 2010 Jan 7;362(1):36-44.
Tavendale R, Macgregor DF, Mukhopadhyay S, Palmer CN. A polymorphism controlling ORMDL3 expression is associated with asthma that is poorly controlled by current medications.. J Allergy Clin Immunol 2008 Apr;121(4):860-3.
Thumberger T, Hagenlocher C, Tisler M, Beyer T, Tietze N, Schweickert A, Feistel K, Blum M. Ciliary and non-ciliary expression and function of PACRG during vertebrate development.. Cilia 2012 Aug 1;1(1):13.
Thusberg J, Vihinen M. Pathogenic or not? And if so, then how? Studying the effects of missense mutations using bioinformatics methods.. Hum Mutat 2009 May;30(5):703-14.
Torgerson DG, Ampleford EJ, Chiu GY, Gauderman WJ, Gignoux CR, Graves PE, Himes BE, Levin AM, Mathias RA, Hancock DB, Baurley JW, Eng C, Stern DA, Celedón JC, Rafaels N, Capurso D, Conti DV, Roth LA, Soto-Quiros M, Togias A, Li X, Myers RA, Romieu I, Van Den Berg DJ, Hu D, Hansel NN, Hernandez RD, Israel E, Salam MT, Galanter J, Avila PC, Avila L, Rodriquez-Santana JR, Chapela R, Rodriguez-Cintron W, Diette GB, Adkinson NF, Abel RA, Ross KD, Shi M, Faruque MU, Dunston GM, Watson HR, Mantese VJ, Ezurum SC, Liang L, Ruczinski I, Ford JG, Huntsman S, Chung KF, Vora H, Li X, Calhoun WJ, Castro M, Sienra-Monge JJ, del Rio-Navarro B, Deichmann KA, Heinzmann A, Wenzel SE, Busse WW, Gern JE, Lemanske RF Jr, Beaty TH, Bleecker ER, Raby BA, Meyers DA, London SJ, Gilliland FD, Burchard EG, Martinez FD, Weiss ST, Williams LK, Barnes KC, Ober C, Nicolae DL. Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations.. Nat Genet 2011 Jul 31;43(9):887-92.
Udina IG, Kordicheva SIu, Pospelov LE, Malenko AF, Gergert VIa, Pospelov AL, Matrashkin AG, Kyzyl-Ool MM, Nachin AA, Zhivotovskiĭ LA. [Study of the polymorphic markers--the PARK2 and PACRG genes due to the incidence of pulmonary tuberculosis in two districts of the Republic of Tyva].. Probl Tuberk Bolezn Legk 2007;(7):27-9.
Valley CC, Cembran A, Perlmutter JD, Lewis AK, Labello NP, Gao J, Sachs JN. The methionine-aromatic motif plays a unique role in stabilizing protein structure.. J Biol Chem 2012 Oct 12;287(42):34979-34991.
Vargas A, Roux-Dalvai F, Droit A, Lavoie JP. Neutrophil-Derived Exosomes: A New Mechanism Contributing to Airway Smooth Muscle Remodeling.. Am J Respir Cell Mol Biol 2016 Sep;55(3):450-61.
West AB, Lockhart PJ, O'Farell C, Farrer MJ. Identification of a novel gene linked to parkin via a bi-directional promoter.. J Mol Biol 2003 Feb 7;326(1):11-9.
Wilson GR, Wang HX, Egan GF, Robinson PJ, Delatycki MB, O'Bryan MK, Lockhart PJ. Deletion of the Parkin co-regulated gene causes defects in ependymal ciliary motility and hydrocephalus in the quakingviable mutant mouse.. Hum Mol Genet 2010 Apr 15;19(8):1593-602.
Wolfenden R. Experimental measures of amino acid hydrophobicity and the topology of transmembrane and globular proteins.. J Gen Physiol 2007 May;129(5):357-62.
Wu H, Romieu I, Sienra-Monge JJ, Li H, del Rio-Navarro BE, London SJ. Genetic variation in ORM1-like 3 (ORMDL3) and gasdermin-like (GSDML) and childhood asthma.. Allergy 2009 Apr;64(4):629-35.
Ribeiro G, Baldi F, Cesar ASM, Alexandre PA, Peripolli E, Ferraz JBS, Fukumasu H. Detection of potential functional variants based on systems-biology: the case of feed efficiency in beef cattle. BMC Genomics 2022 Nov 25;23(1):774.
Xiang F, Sheng J, Li G, Ma J, Wang X, Jiang C, Zhang Z. Black soldier fly larvae vermicompost alters soil biochemistry and bacterial community composition. Appl Microbiol Biotechnol 2022 Jun;106(11):4315-4328.
Jiao W, Sun H, Zhang Z, Xiao Z, Song H, Liu J, Xu X, Wang J, Wang G, Zhang J, Wang C, Li L, Chen L. Construction of a Heterotrophic Nitrification-Aerobic Denitrification Composite Microbial Consortium and Its Bioaugmentation Role in Wastewater Treatment. Biology (Basel) 2025 Dec 4;14(12).
Li W, Chen W, Wang Y, Wang Q, Yang H, Wang Q, Wang B. APOBEC1-Dependent RNA Eiting of TNF Signaling Orchestrates Ileal Villus Morphogenesis in Pigs: Integrative Transcriptomic and Editomic Insights. Animals (Basel) 2025 Aug 18;15(16).
Chen L, Wang G, Song H, Yang Q, Fu J, Liu J, Sun H, Wang Y, Tian Q, Sun Y, Sun L, Xin H, Xiao Z, Wang G, Zhang Z, Zhao Y, Yang H, Li L. Adding Fruit Fermentation Liquid Improves the Efficiency of the Black Soldier Fly in Converting Chicken Manure and Reshapes the Structure of Its Intestinal Microbial Community. Insects 2025 Apr 29;16(5).
Chen L, Li L, Wang G, Xu M, Xin Y, Song H, Liu J, Fu J, Yang Q, Tian Q, Wang Y, Sun H, Lin J, Chen L, Zhang J, Lin J. Insights into a Novel and Efficient Microbial Nest System for Treating Pig Farm Wastewater. Microorganisms 2025 Mar 19;13(3).