<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="brief-report">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Vet. Sci.</journal-id>
<journal-title>Frontiers in Veterinary Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Vet. Sci.</abbrev-journal-title>
<issn pub-type="epub">2297-1769</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fvets.2022.1050561</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Veterinary Science</subject>
<subj-group>
<subject>Brief Research Report</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Molecular epidemiological study of <italic>Trichomonas gallinae</italic> focusing on central and southeastern Europe</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Tuska-Szalay</surname> <given-names>Barbara</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2017753/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Sipos</surname> <given-names>G&#x000E1;bor</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Tak&#x000E1;cs</surname> <given-names>N&#x000F3;ra</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Kontsch&#x000E1;n</surname> <given-names>Jen&#x00151;</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>S&#x000E1;ndor</surname> <given-names>Attila D.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/91784/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>P&#x000E9;ter</surname> <given-names>&#x000C1;ron</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Berta</surname> <given-names>Kriszti&#x000E1;n</given-names></name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Kerek</surname> <given-names>&#x000C1;d&#x000E1;m</given-names></name>
<xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1761070/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Jerzsele</surname> <given-names>&#x000C1;kos</given-names></name>
<xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1947546/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Vot&#x000FD;pka</surname> <given-names>Jan</given-names></name>
<xref ref-type="aff" rid="aff7"><sup>7</sup></xref>
<xref ref-type="aff" rid="aff8"><sup>8</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/664702/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Hornok</surname> <given-names>S&#x000E1;ndor</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1451868/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Parasitology and Zoology, University of Veterinary Medicine</institution>, <addr-line>Budapest</addr-line>, <country>Hungary</country></aff>
<aff id="aff2"><sup>2</sup><institution>ELKH-&#x000C1;TE Climate Change: New Blood-Sucking Parasites and Vector-Borne Pathogens Research Group</institution>, <addr-line>Budapest</addr-line>, <country>Hungary</country></aff>
<aff id="aff3"><sup>3</sup><institution>Centre for Agricultural Research, Plant Protection Institute, ELKH</institution>, <addr-line>Budapest</addr-line>, <country>Hungary</country></aff>
<aff id="aff4"><sup>4</sup><institution>Department of Parasitology and Parasitic Diseases, University of Agricultural Sciences and Veterinary Medicine</institution>, <addr-line>Cluj-Napoca</addr-line>, <country>Romania</country></aff>
<aff id="aff5"><sup>5</sup><institution>Szent Bern&#x000E1;t Small Animal Ambulance</institution>, <addr-line>Budapest</addr-line>, <country>Hungary</country></aff>
<aff id="aff6"><sup>6</sup><institution>Department of Pharmacology and Toxicology, University of Veterinary Medicine</institution>, <addr-line>Budapest</addr-line>, <country>Hungary</country></aff>
<aff id="aff7"><sup>7</sup><institution>Biology Centre, Institute of Parasitology, Czech Academy of Sciences</institution>, <addr-line>&#x0010C;esk&#x000E9; Bud&#x0011B;jovice</addr-line>, <country>Czechia</country></aff>
<aff id="aff8"><sup>8</sup><institution>Department of Parasitology, Faculty of Science, Charles University</institution>, <addr-line>Prague</addr-line>, <country>Czechia</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Mar&#x000ED;a Teresa G&#x000F3;mez-Mu&#x000F1;oz, Complutense University of Madrid, Spain</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Francisco Ponce-Gordo, Universidad Complutense, Spain; Ivana Bilic, University of Veterinary Medicine Vienna, Austria</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Barbara Tuska-Szalay <email>tuska-szalay.barbara&#x00040;univet.hu</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Parasitology, a section of the journal Frontiers in Veterinary Science</p></fn></author-notes>
<pub-date pub-type="epub">
<day>15</day>
<month>12</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>9</volume>
<elocation-id>1050561</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>09</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>18</day>
<month>11</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2022 Tuska-Szalay, Sipos, Tak&#x000E1;cs, Kontsch&#x000E1;n, S&#x000E1;ndor, P&#x000E9;ter, Berta, Kerek, Jerzsele, Vot&#x000FD;pka and Hornok.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Tuska-Szalay, Sipos, Tak&#x000E1;cs, Kontsch&#x000E1;n, S&#x000E1;ndor, P&#x000E9;ter, Berta, Kerek, Jerzsele, Vot&#x000FD;pka and Hornok</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p><italic>Trichomonas gallinae</italic> is a geographically widespread protozoan parasite of birds. In this study, oropharyngeal swab samples were collected in Hungary and Romania from 99 columbiform birds, including 76 feral pigeons (<italic>Columba livia domestica</italic>: 42 kept for racing, 32 with urban and two with rural habitat), four common wood pigeons (<italic>C. palumbus</italic>), 16 ring doves (<italic>Streptopelia risoria</italic>) and three Eurasian collared doves (<italic>S. decaocto</italic>). These samples were analyzed for the presence of <italic>T. gallinae</italic> using molecular methods. Racing feral pigeons had significantly higher prevalence of <italic>T. gallinae</italic> infection than urban feral pigeons. The rate of PCR-positivity was the highest among wood pigeons and ring doves. Based on 18S rRNA gene, <italic>T. gallinae</italic> was the most heterogenous among racing feral pigeons sampled in a trading-breeding place. Clinical signs were associated with only one 18S rRNA gene subtype. The most divergent 18S rRNA gene subtype, <italic>Trichomonas</italic> sp. Hu-TG37 clustered with <italic>T. canistomae</italic> and <italic>T. tenax</italic> and represents probably a new species. To our knowledge, this is the first report on the genetic diversity of <italic>T. gallinae</italic> in the southeastern European region. The results suggest that most detected <italic>T. gallinae</italic> 18S rRNA gene subtypes are not host-specific and do not cause clinical signs. The highest number of 18S rRNA gene subtypes was demonstrated among racing feral pigeons. Significantly more captive than free-living columbiform birds had <italic>T. gallinae</italic> infection. These data highlight the importance of epizootic monitoring of the genetic diversity and presence of <italic>T. gallinae</italic> in trading-breeding places of pigeons and doves.</p>
</abstract>
<kwd-group>
<kwd>Trichomonadea</kwd>
<kwd>18S rRNA gene</kwd>
<kwd>alpha-tubulin gene</kwd>
<kwd>Columbiformes</kwd>
<kwd>pigeon</kwd>
<kwd>dove</kwd>
</kwd-group>
<contract-num rid="cn001">TKP2020-NKA-01</contract-num>
<contract-sponsor id="cn001">Nemzeti Kutat&#x000E1;si, Fejleszt&#x000E9;si &#x000E9;s Innov&#x000E1;ci&#x000F3;s Alap<named-content content-type="fundref-id">10.13039/501100012550</named-content></contract-sponsor>
<counts>
<fig-count count="2"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="25"/>
<page-count count="8"/>
<word-count count="5253"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p><italic>Trichomonas gallinae</italic> (Parabasalia: Trichomonadida) is a widespread flagellated protozoan parasite of birds from various orders, including Columbiformes, Accipitriformes, Strigiformes, Psittaciformes, Falconiformes and Passeriformes (<xref ref-type="bibr" rid="B1">1</xref>&#x02013;<xref ref-type="bibr" rid="B3">3</xref>). Among columbiform birds, the rock pigeon (<italic>Columba livia</italic>) is the main reservoir of this parasite. The most important route of trophozoite transmission between birds is oral by saliva, through shared water and food sources (<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B5">5</xref>). In columbids the predominant way of spreading to nestlings is <italic>via</italic> crop milk (<xref ref-type="bibr" rid="B1">1</xref>). Additionally, predatory birds can become infected by consuming a carrier prey item, since trichomonas trophozoites may survive in carcasses for at least 48 h (<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B5">5</xref>). <italic>Trichomonas gallinae</italic> can persist for up to 1 h in various water sources, e.g., in gutters and drinkers (<xref ref-type="bibr" rid="B6">6</xref>), but higher temperatures (30&#x02013;35&#x000B0;C) can further prolong its survival (<xref ref-type="bibr" rid="B7">7</xref>). Although it is able to form pseudocyst in unfavorable conditions, the moist environment is essential to maintain its infectivity (<xref ref-type="bibr" rid="B2">2</xref>).</p>
<p><italic>Trichomonas gallinae</italic> is often considered a normal inhabitant (commensal) of the mucosal surface in the upper gastrointestinal tract (<xref ref-type="bibr" rid="B2">2</xref>). However, by eliciting inflammation in the underlying tissues or when entering more distally the digestive tract of birds, this protozoan parasite might cause mild to severe lesions depending on strain virulence and host susceptibility. Infection with highly pathogenic strains may lead to death. However, columbiform birds may also be asymptomatic carriers of <italic>T. gallinae</italic>, ensuring the carefree spread of this protozoon. In addition, if protective immunity develops, affected birds become resistant to a new infection (<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B9">9</xref>).</p>
<p>Regarding pathogenesis, <italic>T. gallinae</italic> trophozoites establish preferentially in the upper gastrointestinal tract (the oropharynx, esophagus and crop), where they can cause yellowish necrotic lesions. In severe cases trichomonosis can lead to starvation and suffocation. Furthermore, <italic>T. gallinae</italic> can spread to tissues of the cranium, thorax, and abdomen, as well as of the liver and air sacs causing similar lesions, with deepening tissue involvement, referred to as canker (<xref ref-type="bibr" rid="B1">1</xref>). Recently, massive death of passeriform birds (the greenfinch, <italic>Chloris chloris</italic>, and the goldfinch, <italic>Carduelis carduelis</italic>) due to infection with a specific lineage of <italic>T. gallinae</italic> has been reported throughout Europe (<xref ref-type="bibr" rid="B3">3</xref>).</p>
<p>Avian trichomonosis has been reported to have a worldwide occurrence (<xref ref-type="bibr" rid="B2">2</xref>). In Europe, as on other continents, columbiform birds play the most significant role in the maintenance of <italic>T. gallinae</italic> (<xref ref-type="bibr" rid="B2">2</xref>). The prevalence in western and southern Europe is high among wild columbids (74%) (<xref ref-type="bibr" rid="B9">9</xref>), unlike in northern central Europe where only a little more than one third of racing pigeons proved to be PCR positive (<xref ref-type="bibr" rid="B10">10</xref>).</p>
<p>The main objective of this study was to investigate the prevalence, genetic diversity and phylogenetic relationships of avian <italic>Trichomonas</italic> species in Hungary and Romania, where no similar data are available. Sampled hosts included racing and urban feral pigeons as well as other highly urbanized or pet columbiform bird species (Eurasian collared doves: <italic>Streptopelia decaocto</italic>, common wood pigeons: <italic>C. palumbus</italic> and ring doves: <italic>S. risoria</italic>, respectively). The nomenclature used in the context of feral pigeons complies with genetic and ornithological studies on this species (<xref ref-type="bibr" rid="B11">11</xref>, <xref ref-type="bibr" rid="B12">12</xref>).</p>
</sec>
<sec sec-type="methods" id="s2">
<title>Methods</title>
<sec>
<title>Sample collection</title>
<p>Oropharyngeal mucosal samples were collected with sterile cotton swab applicators randomly from 99 columbiform birds that underwent routine veterinary examination in Hungary (<italic>n</italic> = 77) and Romania (<italic>n</italic> = 22) between May and August, 2021. Four avian host species were sampled, each bird on one occasion, including feral pigeons (<italic>Columba livia domestica</italic>: 42 racing, 32 urban and 2 rural), wood pigeons (<italic>n</italic> = 4), ring doves (<italic>n</italic> = 16) and Eurasian collared doves (<italic>n</italic> = 3) (<xref ref-type="supplementary-material" rid="SM3">Supplementary Table 1</xref>). Racing pigeons were sampled at a trading-breeding place in Csepel (Budapest, Hungary). To assess the necessity of culturing <italic>T. gallinae</italic> prior to DNA extraction, the swab sampling was performed in duplicates from 20 racing pigeons, and one swab sample was placed into 8 ml CPLM culture medium with <italic>Trichomonas</italic> selective supplement (Biolab Diagnostics Laboratory Inc., Budapest Hungary), containing streptomycin, penicillin and sterile inactivated horse serum (pH adjusted to 6). These cultures were kept at 37&#x000B0;C for 2 days. All other swab samples were placed in 2 ml sterile Sarstedt tubes and frozen at &#x02212;20&#x000B0;C.</p>
</sec>
<sec>
<title>DNA extraction and PCR methods</title>
<p>DNA was extracted with the QIAamp DNA Mini Kit (Qiagen, Hilden Germany) according to the manufacturer&#x00027;s blood or tissue protocol, with slight modification. In particular, DNA extraction was performed from 200 &#x003BC;l of culture medium in duplicates, after adding 200 &#x003BC;l AL buffer and continuing with the blood DNA extraction protocol. On the other hand, thawed swabs were overlaid with 200 &#x003BC;l AL buffer and 200 &#x003BC;l sterile PBS, incubated for 10 min at 56&#x000B0;C prior to removal of cotton swab from the fluid, followed by adding proteinase-K and continuing the procedure according to the tissue protocol. In each group of 23 samples an extraction control (180 &#x003BC;l tissue lysis buffer) was included to monitor cross-contamination.</p>
<p>All DNA extracts and extraction controls were analyzed with three conventional PCRs: first with a screening assay amplifying a short, approx. 500-bp-long fragment of the 18S rRNA gene to detect the presence of Trichomonadea, followed by a primary and a secondary assay for sequencing approx. 1,550&#x02013;1,600 and 1,200-bp-long parts of two genetic markers (18S rRNA gene and alpha-tubulin genes, respectively). The reasons for selecting these two genetic markers were to include a conserved gene (18S rRNA) that is widely used in molecular characterization of <italic>T. gallinae</italic>, and a protein encoding gene (alpha-tubulin) for which sequences corresponding to 18S rRNA gene subtypes are available from North America (<xref ref-type="bibr" rid="B8">8</xref>) but not from Europe. The suitability of the screening assay was checked by sequencing PCR products of 13 samples which verified the presence of <italic>T. gallinae</italic> in all cases. Primers and cycling conditions of PCRs are summarized in the <xref ref-type="supplementary-material" rid="SM2">Technical Appendix</xref>. In these PCRs 5 &#x003BC;l of extracted DNA was added to 20 &#x003BC;l of reaction mixture containing 1.0 U HotStar Taq Plus DNA Polymerase (5 U/&#x003BC;l) (Qiagen, Hilden, Germany), 0.5 &#x003BC;l dNTP Mix (10 mM), 0.5 &#x003BC;l of each primer (50 &#x003BC;M), 2.5 &#x003BC;l of 10&#x000D7; Coral Load PCR buffer (15 mM MgCl<sub>2</sub> included), 1 &#x003BC;l extra MgCl<sub>2</sub> (25 mM) and 14.8 &#x003BC;l distilled water. Except for alpha-tubulin PCR where 15.8 &#x003BC;l distilled water was added without extra MgCl<sub>2.</sub> Sequence-verified <italic>T. vaginalis</italic> served as a positive control.</p>
</sec>
<sec>
<title>Sequencing and phylogenetic analyses</title>
<p>The purification and sequencing of the PCR products were performed at Biomi Ltd. (G&#x000F6;d&#x000F6;ll&#x00151;, Hungary). The newly generated sequences were submitted to GenBank under accession numbers <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631556">ON631556</ext-link>-<ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631566">ON631566</ext-link> (18S rRNA gene, long fragment) and <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON808545">ON808545</ext-link>-<ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON808550">ON808550</ext-link> (alpha-tubulin gene) (<xref ref-type="supplementary-material" rid="SM3">Supplementary Table 1</xref>). The 18S rRNA gene subtypes are designated as A to E here independently from other studies, for instance where ITS was used for this purpose (<xref ref-type="bibr" rid="B13">13</xref>).</p>
<p>Obtained sequences were compared to GenBank data using the nucleotide BLASTn program (<ext-link ext-link-type="uri" xlink:href="https://blast.ncbi.nlm.nih.gov">https://blast.ncbi.nlm.nih.gov</ext-link>). All sequences retrieved from GenBank and included in the phylogenetic analysis had nearly or exactly 100% coverage with sequences from this study. The dataset was resampled 1,000 times to generate bootstrap values. Phylogenetic analysis was conducted by using the maximum-likelihood method and the Jukes-Cantor model according to the best-fit selection with the program MEGA 7.0 (<xref ref-type="bibr" rid="B14">14</xref>). Prevalence data were compared with Fisher exact test. Differences were regarded significant when <italic>P</italic> &#x0003C; 0.05.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<p>According to the preliminary comparison on the efficacy of molecular detection of <italic>T. gallinae</italic>, 19 out of 20 swab samples but only 18 out of 20 culture medium samples were PCR positive (i.e., one pigeon was positive only by its culture medium sample vs. two pigeons diagnosed as infected only from their swab samples) (<xref ref-type="supplementary-material" rid="SM3">Supplementary Table 1</xref>: 7th column, samples TG1-50). Therefore, swab samples were used in the remaining part of the study.</p>
<p><italic>Trichomonas gallinae</italic> was detected in all four studied bird species, with an overall prevalence of 73% (72 out of 99) (<xref ref-type="table" rid="T1">Table 1</xref>). Based on the screening assay, racing feral pigeons had a significantly (<italic>P</italic> &#x0003C; 0.0001) higher prevalence of <italic>T. gallinae</italic> infection (95%: 40 of 42) than urban feral pigeons (34%: 11 of 32). Among other columbiform bird species, the rate of PCR-positivity was the highest, 100% among wood pigeons (4 out of 4), followed in decreasing order by ring doves (94%: 15 out of 16) and collared doves (33%: 1 out of 3) (<xref ref-type="table" rid="T1">Table 1</xref>; <xref ref-type="supplementary-material" rid="SM3">Supplementary Table 1</xref>). These results also imply that the association of <italic>T. gallinae</italic> infection with artificially bred, captive columbiform birds (racing feral pigeons and ring doves) was highly significant (<italic>P</italic> &#x0003C; 0.0001) in comparison with free-living columbiform birds (urban and rural feral pigeons, wood pigeons and collared doves) (55/58 vs. 17/41, respectively).</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Results of molecular analyses for the 18S rRNA gene of Trichomonadea according to sample types of columbiform birds.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Species or type of sample source</bold></th>
<th valign="top" align="left"><bold>Country of origin</bold></th>
<th valign="top" align="center"><bold>Prevalence<xref ref-type="table-fn" rid="TN2"><sup>&#x0002A;&#x0002A;</sup></xref> (positive/all)</bold></th>
<th valign="top" align="left"><bold>Long 18S rRNA gene subtype (<italic>n</italic>)</bold></th>
<th valign="top" align="left"><bold>GenBank accession number</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Racing feral pigeon (<italic>Columba livia</italic>)</td>
<td valign="top" align="left">Hungary</td>
<td valign="top" align="center">97% (29/30)</td>
<td valign="top" align="left"><bold>B</bold> (16), <bold>C</bold> (3),<break/> <italic>Trichomonas</italic> sp.<break/> Hu-TG37 (1)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631556">ON631556</ext-link>, <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631557">ON631557</ext-link>,<break/> <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631566">ON631566</ext-link></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Germany<xref ref-type="table-fn" rid="TN1"><sup>&#x0002A;</sup></xref></td>
<td valign="top" align="center">89% (8/9)</td>
<td valign="top" align="left"><bold>D</bold> (3)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631558">ON631558</ext-link></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Denmark<xref ref-type="table-fn" rid="TN1"><sup>&#x0002A;</sup></xref></td>
<td valign="top" align="center">100% (2/2)</td>
<td valign="top" align="left"><bold>B</bold> (1)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631556">ON631556</ext-link></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Belgium<xref ref-type="table-fn" rid="TN1"><sup>&#x0002A;</sup></xref></td>
<td valign="top" align="center">&#x02013;(1/1)</td>
<td valign="top" align="left"><bold>B</bold> (1)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631556">ON631556</ext-link></td>
</tr>
<tr>
<td valign="top" align="left">Urban feral pigeon (<italic>Columba livia</italic>)</td>
<td valign="top" align="left">Hungary</td>
<td valign="top" align="center">40% (4/10)</td>
<td valign="top" align="left"><bold>B</bold> (1), <bold>D</bold> (1)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631559">ON631559</ext-link>, <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631560">ON631560</ext-link></td>
</tr>
<tr>
<td/>
<td valign="top" align="left">Romania</td>
<td valign="top" align="center">32% (7/22)</td>
<td valign="top" align="left"><bold>A</bold> (1), <bold>D</bold> (2)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631561">ON631561</ext-link>, <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631562">ON631562</ext-link></td>
</tr>
<tr>
<td valign="top" align="left">Rural feral pigeon (<italic>Columba livia</italic>)</td>
<td valign="top" align="left">Hungary</td>
<td valign="top" align="center">50% (1/2)</td>
<td valign="top" align="left">&#x02013;</td>
<td valign="top" align="left">&#x02013;</td>
</tr>
<tr>
<td valign="top" align="left">Common wood pigeon (<italic>Columba palumbus</italic>)</td>
<td valign="top" align="left">Hungary</td>
<td valign="top" align="center">100% (4/4)</td>
<td valign="top" align="left"><bold>A</bold> (1)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631563">ON631563</ext-link></td>
</tr>
<tr>
<td valign="top" align="left">Ring dove (<italic>Streptopelia risoria</italic>)</td>
<td valign="top" align="left">Hungary</td>
<td valign="top" align="center">94% (15/16)</td>
<td valign="top" align="left"><bold>B</bold> (4), <bold>E</bold> (2)</td>
<td valign="top" align="left"><ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631564">ON631564</ext-link>, <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="ON631565">ON631565</ext-link></td>
</tr>
<tr>
<td valign="top" align="left">Eurasian collared dove (<italic>Streptopelia decaocto</italic>)</td>
<td valign="top" align="left">Hungary</td>
<td valign="top" align="center">33% (1/3)</td>
<td valign="top" align="left">&#x02013;</td>
<td valign="top" align="left">&#x02013;</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="TN1"><label>&#x0002A;</label><p>Kept separately after arrival, but contact with Hungarian birds cannot be excluded;</p></fn>
<fn id="TN2"><label>&#x0002A;&#x0002A;</label><p>based on the short 18S rRNA gene screening PCR.</p></fn>
</table-wrap-foot>
</table-wrap>
<p>Based on the long fragment of the 18S rRNA gene, six <italic>Trichomonas</italic> subtypes were detected in 37 columbiform birds (<xref ref-type="table" rid="T1">Table 1</xref>). Compared to a reference sequence from North America (GenBank: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="EU215373">EU215373</ext-link>), these had four to eight nucleotide differences (<xref ref-type="table" rid="T2">Table 2</xref>), i.e., 99.5&#x02013;99.7% (1,471&#x02013;1,475/1,479 bp) identity. In a smaller geographical context, the 18S rRNA gene sequences obtained in this study from central and southeastern Europe had higher, up to 11 nucleotide differences with conspecific sequences from western and southwestern Europe: only 99.3&#x02013;99.7% (1,468&#x02013;1,474/1,479 bp) sequence identity with <italic>T. gallinae</italic> from passeriform birds sampled in 2019 in France (e.g., MK172846), whereas 99.5&#x02013;100% (1,443/1,450&#x02013;1,449/1,449 bp) sequence identity with samples collected from columbiform birds in 2019 in Portugal (e.g., MK932772).</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Site-specific variations among 18S rRNA gene subtypes of <italic>Trichomonas gallinae</italic> and <italic>Trichomonas</italic> sp. Hu-TG37 compared to positions in EU215373 used as a reference sequence (identical with genotype A).</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>18S rRNA</bold></th>
<th valign="top" align="center" colspan="15" style="border-bottom: thin solid #000000;"><bold>Position (in which the nucleotide is indicated below)</bold></th>
</tr>
<tr>
<th valign="top" align="left"><bold>gene subtype</bold></th>
<th valign="top" align="center"><bold>33</bold></th>
<th valign="top" align="center"><bold>215</bold></th>
<th valign="top" align="center"><bold>384</bold></th>
<th valign="top" align="center"><bold>386</bold></th>
<th valign="top" align="center"><bold>408</bold></th>
<th valign="top" align="center"><bold>418</bold></th>
<th valign="top" align="center"><bold>590</bold></th>
<th valign="top" align="center"><bold>592</bold></th>
<th valign="top" align="center"><bold>648</bold></th>
<th valign="top" align="center"><bold>825</bold></th>
<th valign="top" align="center"><bold>837</bold></th>
<th valign="top" align="center"><bold>850</bold></th>
<th valign="top" align="center"><bold>1120</bold></th>
<th valign="top" align="center"><bold>1143</bold></th>
<th valign="top" align="center"><bold>1281</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><bold>A</bold></td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">C</td>
</tr>
<tr>
<td valign="top" align="left"><bold>B</bold></td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">&#x02013;</td>
</tr>
<tr>
<td valign="top" align="left"><bold>C</bold></td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">&#x02013;</td>
</tr>
<tr>
<td valign="top" align="left"><bold>D</bold></td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">&#x02013;</td>
</tr>
<tr>
<td valign="top" align="left"><bold>E</bold></td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">C</td>
<td valign="top" align="center">G</td>
<td valign="top" align="center">&#x02013;</td>
</tr>
<tr>
<td valign="top" align="left"><bold>Hu-TG37</bold></td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">A</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">&#x02013;</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">T</td>
<td valign="top" align="center">T</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Nucleotides that are identical to the reference are indicated with &#x0201C;&#x02013;&#x0201D;. Note that only the last three positions shown in this table are covered by the screening (short 18S rRNA gene) PCR, i.e., genotyping was possible from 37 of 72 <italic>Trichomonas</italic>-positive specimens yielding sequenceable product in the long 18S rRNA PCR.</p>
</table-wrap-foot>
</table-wrap>
<p>Only two genetic variants were detected among sympatric urban feral pigeons (subtypes B, D in Hungary; A, D in Romania), while four subtypes (B, C, D and the most divergent Hu-TG37) occurred in racing feral pigeons kept in one trading-breeding place (<xref ref-type="table" rid="T1">Table 1</xref>). Among racing feral pigeons, 18S rRNA gene subtype D was only found in birds that originated from Germany, while in Hungary subtypes A and E were exclusively found in wood pigeons and ring doves, respectively (<xref ref-type="table" rid="T1">Table 1</xref>). In the latter species, the occurrence of 18S rRNA gene subtypes was related to the origin (breeding place) of captive birds: subtype E was detected in two ring doves of one sampling locality, whereas subtype B in four birds at the other sampling locality (<xref ref-type="supplementary-material" rid="SM3">Supplementary Table 1</xref>). Clinical signs (lesions in the oropharyngeal cavity or the eyes, conjunctiva) were only associated with subtype D (<xref ref-type="supplementary-material" rid="SM3">Supplementary Table 1</xref>; <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1</xref>).</p>
<p>Amplification and sequencing of part of the alpha-tubulin gene were successful from at least one sample representing each 18S subtype. Alpha-tubulin gene sequences obtained in this study had a lower (98.5&#x02013;98.6%: 1,008&#x02013;1,009/1,023 bp) or a higher level (99.6&#x02013;99.7%: 1,019&#x02013;1,020/1,023 bp) of sequence identity to that of an isolate (EU215382) used as a reference. Based on the corresponding amino acid sequences, most mutations in this protein encoding gene were synonymous, but subtypes E and D had a single amino acid difference compared to the reference isolate (valine instead of alanine at position 149, and isoleucine instead of valine at position 34, respectively).</p>
<p>Considering the results of phylogenetic analyses, the topology of the 18S rRNA phylogenetic tree (<xref ref-type="fig" rid="F1">Figure 1</xref>) did not show clear clustering of <italic>T. gallinae</italic> subtypes obtained in this study according to host species, living place (or keeping modes) of columbiform birds. However, subtype A (identified both in Hungary and Romania) formed a sister group to all other 18S rRNA gene subtypes (B, C, D, E): although with only moderate (65%) support (<xref ref-type="fig" rid="F1">Figure 1</xref>) but confirmed by the alpha-tubulin phylogenetic tree (<xref ref-type="fig" rid="F2">Figure 2</xref>). More importantly, <italic>Trichomonas</italic> sp. Hu-TG37 (detected in a racing feral pigeon from south Hungary) belonged to the phylogenetic group of <italic>T. canistomae</italic> and <italic>T. tenax</italic> with moderate (61%) support, implying, that this is a separate species (in other words: if this isolate would belong to <italic>T. gallinae</italic>, this species would not be monophyletic). The separation of <italic>Trichomonas</italic> sp. Hu-TG37 from <italic>T. gallinae</italic> was also confirmed by the alpha-tubulin phylogenetic tree (<xref ref-type="fig" rid="F2">Figure 2</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Phylogenetic tree of Trichomonadea based on the 18S rRNA gene, made with the maximum likelihood method and the Jukes-Cantor model. In each row, after the species or genus name, the isolation source, for <italic>Trichomonas gallinae</italic> the country of origin and the GenBank accession number are shown. Sequences obtained in this study and representing each subtype are in red and bold accession numbers. The analysis involved 44 nucleotide sequences and 1,000 bootstrap replications. There were a total of 1,121 positions in the final dataset. The scale-bar indicates the number of substitutions per site.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fvets-09-1050561-g0001.tif"/>
</fig>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Phylogenetic tree of Trichomonadea based on the alpha-tubulin gene, made with the maximum likelihood method and the Jukes-Cantor model. In each row, after the species or genus name, the isolation source, for species closely related to <italic>Trichomonas gallinae</italic> the country of origin and the GenBank accession number are shown. Sequences obtained in this study and representing each subtype are in red and bold accession numbers. The analysis involved 20 nucleotide sequences and 1,000 bootstrap replications. There were a total of 977 positions in the final dataset. The scale-bar indicates the number of substitutions per site.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fvets-09-1050561-g0002.tif"/>
</fig>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>To our knowledge, this is the first study on the genetic diversity of <italic>T. gallinae</italic> in Hungary, Romania and the whole southeastern European region, complementing previous reports from western, central and southern Europe (see below). Most infected birds in this study did not show clinical signs of trichomonosis except five pigeons (6.9%). A lower prevalence of clinical trichomonosis (0.37%) has been reported in a study involving 612 wild and domestic pigeons (<xref ref-type="bibr" rid="B13">13</xref>). Since <italic>T. gallinae</italic> can cause the death of infected hosts (<xref ref-type="bibr" rid="B15">15</xref>), the rarity of symptoms can in part be explained by the death of severely affected birds, which can die before they are examined (<xref ref-type="bibr" rid="B13">13</xref>). Furthermore, from an epidemiological point of view, subclinical cases might ensure easier spread of these protozoa. Therefore, screening pigeons for <italic>Trichomonas</italic> spp. should become an integral part of veterinary practice.</p>
<p>In this study, <italic>T. gallinae</italic> was detected in all four columbiform bird species examined. The infection rate was 73% which is similar to the 74% overall prevalence reported in wild columbids from western and southern Europe (<xref ref-type="bibr" rid="B9">9</xref>). Within the Mediterranean region, in the Iberian Peninsula 44.8% of wild and domestic pigeons were shown to carry <italic>T. gallinae</italic> (<xref ref-type="bibr" rid="B13">13</xref>), but in another study from Spain the prevalence was much higher (79.4%) (<xref ref-type="bibr" rid="B16">16</xref>). Furthermore, in the UK other columbiform species including <italic>C. palumbus</italic> and <italic>S. decaocto</italic> were also examined and a 60% incidence was reported (<xref ref-type="bibr" rid="B17">17</xref>). In Germany, four species (<italic>C. livia, C. oenas, C. palumbus</italic>, and <italic>S. decaocto</italic>) were screened and 50% of the birds were infected with <italic>T. gallinae</italic> (<xref ref-type="bibr" rid="B18">18</xref>). The prevalence was 37% among racing pigeons in Poland (<xref ref-type="bibr" rid="B10">10</xref>).</p>
<p>In Hungary, based on our screening assay, racing feral pigeons had a significantly higher prevalence of <italic>T. gallinae</italic> infection than urban feral pigeons (95 vs. 33%). The possible reason for this difference might be that at trading-breeding places birds of different origin are housed close to each other, and the chances for infection are higher if naive racing pigeons can get into contact with carrier birds not only outside, but inside such enclosures, e.g., by sharing food and drinking water or by kissing. Therefore, we propose that similar places play a crucial role as hotspots in the transmission of <italic>T. gallinae</italic>. To our knowledge, there is no similar study in Europe that compares <italic>T. gallinae</italic> according to the place where pigeons are kept.</p>
<p>All four examined wood pigeons were infected. Although the sample size was limited, this apparently high (100%) infection rate is similar to what was reported among wood pigeons in Germany (70%) (<xref ref-type="bibr" rid="B9">9</xref>) and the Iberian Peninsula (83.3%) (<xref ref-type="bibr" rid="B19">19</xref>). This high prevalence is likely associated with urbanization of wood pigeons in Hungary, resulting in close contact of birds (e.g., <italic>via</italic> sharing common drinking sources) in green areas where the number of wood pigeons has recently increased significantly (<xref ref-type="bibr" rid="B20">20</xref>).</p>
<p>Concerning other studied bird species, the prevalence was also high (94%) among ring doves in Hungary. There is a lack of literature data on trichomonosis of both collared dove and ring dove in Europe, despite the fact that ring doves were found to be susceptible to <italic>T. gallinae</italic> during experimental infection (<xref ref-type="bibr" rid="B21">21</xref>). In the Caribbean an outbreak was reported (<xref ref-type="bibr" rid="B22">22</xref>), drawing the attention of veterinarians to the necessity to monitor <italic>T. gallinae</italic> in this bird species.</p>
<p>Based on the long fragment of the 18S rRNA gene, six <italic>Trichomonas</italic> subtypes were detected in columbiform birds in Hungary and Romania. Compared to a reference sequence, they had up to eight nucleotide differences, meaning that the maximum genetic difference was low (0.5%) compared to what was reported from North America (3.4%) (<xref ref-type="bibr" rid="B8">8</xref>) and even from Austria, a country neighboring Hungary (2.9%) (<xref ref-type="bibr" rid="B23">23</xref>). In Hungary, two genetic variants were detected among urban feral pigeons, vs. four 18S rRNA gene subtypes occurred in racing feral pigeons kept in the same trading-breeding place, highlighting the epidemiological importance of similar facilities in general.</p>
<p>It is noteworthy that each bird from which <italic>T. gallinae</italic> was sequenced, carried a single 18S rRNA gene subtype, as also demonstrated repeatedly with different modes of detection (<xref ref-type="supplementary-material" rid="SM3">Supplementary Table 1</xref>, 8th column: samples TG1-50). Since it was reported that <italic>T. gallinae</italic> triggers premunition (<xref ref-type="bibr" rid="B24">24</xref>) and immunity lasts until the loss of infection (<xref ref-type="bibr" rid="B1">1</xref>), the <italic>a priori</italic> presence of any variant probably protected the relevant birds during a heterologous challenge, which is likely to occur in the environment of a pigeon trading-breeding place. It is also highly relevant to note that among racing pigeons, 18S rRNA gene subtype D was only found in pigeons that originated from Germany. We suspect that subtype D was already present in these birds when imported to Hungary and probably premunition protected them from becoming infected in the pigeon trading place by other variants.</p>
<p>The clinical signs relevant to trichomonosis were only associated with 18S rRNA gene subtype D in Hungary. There are several reports indicating that certain haplotypes are highly correlated with more severe lesions in various bird species (<xref ref-type="bibr" rid="B13">13</xref>, <xref ref-type="bibr" rid="B16">16</xref>, <xref ref-type="bibr" rid="B25">25</xref>). However, no close correlation has been found in terms of pathogenicity and geographical distribution when columbids were examined across Europe (<xref ref-type="bibr" rid="B9">9</xref>).</p>
<p>Considering the results in a phylogenetic context, the topology of the 18S rRNA phylogenetic tree did not show clear clustering of <italic>T. gallinae</italic> subtypes from this study according to host species, living place or keeping mode of columbiform birds. Clustering of <italic>Trichomonas</italic> sp. Hu-TG37 to the phylogenetic group of <italic>T. canistomae</italic> and <italic>T. tenax</italic> suggests that it might represent a new species. This separate position of <italic>Trichomonas</italic> sp. Hu-TG37 was confirmed also by the alpha-tubulin phylogenetic tree. It is known that some <italic>Trichomonas</italic> spp. show high genetic diversity depending on bird species, and few of these infections are caused by variants/species closely related to <italic>T. vaginalis, T. tenax</italic> (<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B23">23</xref>) or <italic>T. canistomae</italic> (<xref ref-type="bibr" rid="B16">16</xref>). <italic>Trichomonas</italic> sp. Hu-TG37 identified for the first time in the present study belongs to the phylogenetic group of <italic>Trichomonas</italic> spp. infecting (among the others) domestic carnivores. This calls for further epidemiological studies on the possible contact between dogs and pigeons (e.g., <italic>via</italic> water deposited in gardens in drinking bowls) and its role in the transmission of these protozoan parasites.</p>
<p>In conclusion, this is the first report in Hungary and Romania on the prevalence and 18S rRNA gene subtypes of <italic>T. gallinae</italic> in various columbiform birds using molecular methods. The results suggest that most of these variants are not host-specific and do not cause clinical signs. The highest degree of genetic diversity and high prevalence of infection was observed among racing pigeons and captive ring doves, thereby highlighting the epidemiological importance of pigeon/dove trading-breeding places.</p>
</sec>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s10">Supplementary material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec sec-type="ethics-statement" id="s6">
<title>Ethics statement</title>
<p>Ethical review and approval was not required for the animal study because all columbiform birds in this study were handled and sampled during regular veterinary care, therefore no ethical permission was needed. Written informed consent was obtained from the owners for the participation of their animals in this study.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>BT-S: study design, DNA extraction, data analysis, and manuscript writing. GS, AS, &#x000C1;P, and KB: sample collection. NT: PCR tests and sequencing. JK: phylogenetic analyses. &#x000C1;K and &#x000C1;J: contribution to sample collection. JV: supervision and contribution to molecular analyses. SH: conceptualization, study design, primer design, GenBank processing, and manuscript writing. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>This study was funded by project no. TKP2020-NKA-01 implemented with the support provided from the National Research, Development and Innovation Fund of Hungary, financed under the T&#x000E9;mater&#x000FC;leti Kiv&#x000E1;l&#x000F3;s&#x000E1;gi Program 2020 (2020-4.1.1-TKP2020) funding scheme.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="s10">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fvets.2022.1050561/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fvets.2022.1050561/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Image_1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Figure 1</label>
<caption><p>Lesions associated with trichomonosis in two racing pigeons: <bold>(A)</bold> yellowish debris in the corner of the beaks and close to the eye; <bold>(B)</bold> small necrotic-inflammatory foci on the palate.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Data_Sheet_1.pdf" id="SM2" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Data_Sheet_2.pdf" id="SM3" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stabler</surname> <given-names>RM</given-names></name></person-group>. <article-title><italic>Trichomonas gallinae</italic>: A review</article-title>. <source>Exp Parasitol.</source> (<year>1954</year>) <volume>3</volume>:<fpage>368</fpage>&#x02013;<lpage>402</lpage>. <pub-id pub-id-type="doi">10.1016/0014-4894(54)90035-1</pub-id><pub-id pub-id-type="pmid">13183096</pub-id></citation></ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Amin</surname> <given-names>A</given-names></name> <name><surname>Bilic</surname> <given-names>I</given-names></name> <name><surname>Liebhart</surname> <given-names>D</given-names></name> <name><surname>Hess</surname> <given-names>M</given-names></name></person-group>. <article-title>Trichomonads in birds&#x02013;a review</article-title>. <source>Parasitology.</source> (<year>2014</year>) <volume>141</volume>:<fpage>733</fpage>&#x02013;<lpage>47</lpage>. <pub-id pub-id-type="doi">10.1017/S0031182013002096</pub-id><pub-id pub-id-type="pmid">24476968</pub-id></citation></ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chavatte</surname> <given-names>J-M</given-names></name> <name><surname>Giraud</surname> <given-names>P</given-names></name> <name><surname>Esperet</surname> <given-names>D</given-names></name> <name><surname>Place</surname> <given-names>G</given-names></name> <name><surname>Cavalier</surname> <given-names>F</given-names></name> <name><surname>Landau</surname> <given-names>I</given-names></name></person-group>. <article-title>An outbreak of trichomonosis in European greenfinches chloris chloris and European goldfinches carduelis carduelis wintering in Northern France</article-title>. <source>Parasite.</source> (<year>2019</year>) <volume>26</volume>:<fpage>21</fpage>. <pub-id pub-id-type="doi">10.1051/parasite/2019022</pub-id><pub-id pub-id-type="pmid">30957740</pub-id></citation></ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gerhold</surname> <given-names>RW</given-names></name> <name><surname>Tate</surname> <given-names>CM</given-names></name> <name><surname>Gibbs</surname> <given-names>SE</given-names></name> <name><surname>Mead</surname> <given-names>DG</given-names></name> <name><surname>Allison</surname> <given-names>AB</given-names></name> <name><surname>Fischer</surname> <given-names>JR</given-names></name></person-group>. <article-title>Necropsy findings and arbovirus surveillance in mourning doves from the southeastern United States</article-title>. <source>J Wildl Dis.</source> (<year>2007</year>) <volume>43</volume>:<fpage>129</fpage>&#x02013;<lpage>35</lpage>. <pub-id pub-id-type="doi">10.7589/0090-3558-43.1.129</pub-id><pub-id pub-id-type="pmid">17347403</pub-id></citation></ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Villan&#x000FA;a</surname> <given-names>D</given-names></name> <name><surname>H&#x000F6;fle</surname> <given-names>U</given-names></name> <name><surname>P&#x000E9;rez-Rodr&#x000ED;guez</surname> <given-names>L</given-names></name> <name><surname>Gort&#x000E1;zar</surname> <given-names>C</given-names></name></person-group>. <article-title><italic>Trichomonas gallinae</italic> in wintering common wood pigeons <italic>Columba palumbus</italic> in Spain</article-title>. <source>Ibis.</source> (<year>2006</year>) <volume>148</volume>:<fpage>641</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1111/j.1474-919X.2006.00561.x</pub-id></citation>
</ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Purple</surname> <given-names>KE</given-names></name> <name><surname>Humm</surname> <given-names>JM</given-names></name> <name><surname>Kirby</surname> <given-names>RB</given-names></name> <name><surname>Saidak</surname> <given-names>CG</given-names></name> <name><surname>Gerhold</surname> <given-names>R</given-names></name></person-group>. <article-title><italic>Trichomonas gallinae</italic> persistence in four water treatments</article-title>. <source>J Wildl Dis.</source> (<year>2015</year>) <volume>51</volume>:<fpage>739</fpage>&#x02013;<lpage>42</lpage>. <pub-id pub-id-type="doi">10.7589/2014-05-137</pub-id><pub-id pub-id-type="pmid">25919469</pub-id></citation></ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kocan</surname> <given-names>RM</given-names></name></person-group>. <article-title>Various grains and liquid as potential vehicles of transmission for <italic>Trichomonas gallinae</italic></article-title>. <source>Wildl Dis</source>. (<year>1969</year>) <volume>5</volume>:<fpage>148</fpage>&#x02013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.7589/0090-3558-5.3.148</pub-id><pub-id pub-id-type="pmid">5817771</pub-id></citation></ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gerhold</surname> <given-names>RW</given-names></name> <name><surname>Yabsley</surname> <given-names>MJ</given-names></name> <name><surname>Smith</surname> <given-names>AJ</given-names></name> <name><surname>Ostergaard</surname> <given-names>E</given-names></name> <name><surname>Mannan</surname> <given-names>W</given-names></name> <name><surname>Cann</surname> <given-names>JD</given-names></name> <etal/></person-group>. <article-title>Molecular characterization of the <italic>Trichomonas gallinae</italic> morphologic complex in the United States</article-title>. <source>J Parasitol.</source> (<year>2008</year>) <volume>94</volume>:<fpage>1335</fpage>&#x02013;<lpage>41</lpage>. <pub-id pub-id-type="doi">10.1645/GE-1585.1</pub-id><pub-id pub-id-type="pmid">18576862</pub-id></citation></ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marx</surname> <given-names>M</given-names></name> <name><surname>Reiner</surname> <given-names>G</given-names></name> <name><surname>Willems</surname> <given-names>H</given-names></name> <name><surname>Rocha</surname> <given-names>G</given-names></name> <name><surname>Hillerich</surname> <given-names>K</given-names></name> <name><surname>Masello</surname> <given-names>JF</given-names></name> <etal/></person-group>. <article-title>High prevalence of <italic>Trichomonas gallinae</italic> in wild columbids across western and southern Europe</article-title>. <source>Parasit Vect.</source> (<year>2017</year>) <volume>10</volume>:<fpage>242</fpage>. <pub-id pub-id-type="doi">10.1186/s13071-017-2170-0</pub-id><pub-id pub-id-type="pmid">28521843</pub-id></citation></ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bobrek</surname> <given-names>K</given-names></name> <name><surname>Urbanowicz</surname> <given-names>J</given-names></name> <name><surname>Chorbi&#x00144;ski</surname> <given-names>P</given-names></name> <name><surname>Gawe&#x00142;</surname> <given-names>A</given-names></name></person-group>. <article-title>Molecular analysis of <italic>Trichomonas gallinae</italic> in racing pigeons from upper Silesia, Poland</article-title>. <source>Pol J Vet Sci.</source> (<year>2017</year>) <volume>20</volume>:<fpage>185</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1515/pjvs-2017-0023</pub-id><pub-id pub-id-type="pmid">28525339</pub-id></citation></ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Haag-Wackernagel</surname> <given-names>D</given-names></name> <name><surname>Heeb</surname> <given-names>P</given-names></name> <name><surname>Leiss</surname> <given-names>A</given-names></name></person-group>. <article-title>Phenotype-dependent selection of juvenile urban feral pigeons <italic>Columba livia</italic></article-title>. <source>Bird Study</source>. (<year>2006</year>) <volume>53</volume>:<fpage>163</fpage>&#x02013;<lpage>70</lpage>. <pub-id pub-id-type="doi">10.1080/00063650609461429</pub-id></citation>
</ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>WJ</given-names></name> <name><surname>Sendell-Price</surname> <given-names>AT</given-names></name> <name><surname>Fayet</surname> <given-names>AL</given-names></name> <name><surname>Schweizer</surname> <given-names>TM</given-names></name> <name><surname>Jezierski</surname> <given-names>MT</given-names></name> <name><surname>van de Kerkhof</surname> <given-names>C</given-names></name> <etal/></person-group>. <article-title>Limited domestic introgression in a final refuge of the wild pigeoni</article-title>. <source>Science.</source> (<year>2022</year>) <volume>25</volume>:<fpage>104620</fpage>. <pub-id pub-id-type="doi">10.1016/j.isci.2022.104620</pub-id><pub-id pub-id-type="pmid">35880028</pub-id></citation></ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sansano-Maestre</surname> <given-names>J</given-names></name> <name><surname>Garijo-Toledo</surname> <given-names>MM</given-names></name> <name><surname>G&#x000F3;mez-Mu&#x000F1;oz</surname> <given-names>MT</given-names></name></person-group>. <article-title>Prevalence and genotyping of <italic>Trichomonas gallinae</italic> in pigeons and birds of prey</article-title>. <source>Avian Pathol J WVPA.</source> (<year>2009</year>) <volume>38</volume>:<fpage>201</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1080/03079450902912135</pub-id><pub-id pub-id-type="pmid">19468936</pub-id></citation></ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname> <given-names>S</given-names></name> <name><surname>Stecher</surname> <given-names>G</given-names></name></person-group>. <article-title>Tamura K. MEGA7: molecular evolutionary genetics analysis version 70 for bigger datasets</article-title>. <source>Mol Biol Evol.</source> (<year>2016</year>) <volume>33</volume>:<fpage>1870</fpage>&#x02013;<lpage>4</lpage>. <pub-id pub-id-type="doi">10.1093/molbev/msw054</pub-id><pub-id pub-id-type="pmid">27004904</pub-id></citation></ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stockdale</surname> <given-names>JE</given-names></name> <name><surname>Dunn</surname> <given-names>JC</given-names></name> <name><surname>Goodman</surname> <given-names>SJ</given-names></name> <name><surname>Morris</surname> <given-names>AJ</given-names></name> <name><surname>Sheehan</surname> <given-names>DK</given-names></name> <name><surname>Grice</surname> <given-names>PV</given-names></name> <etal/></person-group>. <article-title>The protozoan parasite <italic>Trichomonas gallinae</italic> causes adult and nestling mortality in a declining population of European turtle doves, <italic>Streptopelia turtur</italic></article-title>. <source>Parasitology.</source> (<year>2015</year>) <volume>142</volume>:<fpage>490</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1017/S0031182014001474</pub-id><pub-id pub-id-type="pmid">25212735</pub-id></citation></ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mart&#x000ED;nez-Herrero</surname> <given-names>MC</given-names></name> <name><surname>Sansano-Maestre</surname> <given-names>J</given-names></name> <name><surname>L&#x000F3;pez M&#x000E1;rquez</surname> <given-names>I</given-names></name> <name><surname>Ob&#x000F3;n</surname> <given-names>E</given-names></name> <name><surname>Ponce</surname> <given-names>C</given-names></name> <name><surname>Gonz&#x000E1;lez</surname> <given-names>J</given-names></name> <etal/></person-group>. <article-title>Genetic characterization of oropharyngeal trichomonad isolates from wild birds indicates that genotype is associated with host species, diet and presence of pathognomonic lesions</article-title>. <source>Avian Pathol J WVPA.</source> (<year>2014</year>) <volume>43</volume>:<fpage>535</fpage>&#x02013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.1080/03079457.2014.967660</pub-id><pub-id pub-id-type="pmid">25262786</pub-id></citation></ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lennon</surname> <given-names>RJ</given-names></name> <name><surname>Dunn</surname> <given-names>JC</given-names></name> <name><surname>Stockdale</surname> <given-names>JE</given-names></name> <name><surname>Goodman</surname> <given-names>SJ</given-names></name> <name><surname>Morris</surname> <given-names>AJ</given-names></name> <name><surname>Hamer</surname> <given-names>KC</given-names></name></person-group>. <article-title>Trichomonad parasite infection in four species of Columbidae in the UK</article-title>. <source>Parasitology.</source> (<year>2013</year>) <volume>140</volume>:<fpage>1368</fpage>&#x02013;<lpage>76</lpage>. <pub-id pub-id-type="doi">10.1017/S0031182013000887</pub-id><pub-id pub-id-type="pmid">23866933</pub-id></citation></ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quillfeldt</surname> <given-names>P</given-names></name> <name><surname>Schumm</surname> <given-names>YR</given-names></name> <name><surname>Marek</surname> <given-names>C</given-names></name> <name><surname>Mader</surname> <given-names>V</given-names></name> <name><surname>Fischer</surname> <given-names>D</given-names></name> <name><surname>Marx</surname> <given-names>M</given-names></name></person-group>. <article-title>Prevalence and genotyping of <italic>Trichomonas</italic> infections in wild birds in central Germany</article-title>. <source>PLoS ONE.</source> (<year>2018</year>) <volume>13</volume>:<fpage>e0200798</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0200798</pub-id><pub-id pub-id-type="pmid">30092001</pub-id></citation></ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Santos</surname> <given-names>N</given-names></name> <name><surname>Jambas</surname> <given-names>J</given-names></name> <name><surname>Monteiro</surname> <given-names>A</given-names></name> <name><surname>Amaral</surname> <given-names>J</given-names></name> <name><surname>Martins</surname> <given-names>N</given-names></name> <name><surname>Garcia</surname> <given-names>J</given-names></name> <etal/></person-group>. <article-title><italic>Trichomonas</italic> infection in a community of free-ranging domestic and wild Columbiformes and bonelli&#x00027;s eagle (<italic>Aquila fasciata</italic>)</article-title>. <source>Front Vet Sci.</source> (<year>2019</year>) <volume>6</volume>:<fpage>148</fpage>. <pub-id pub-id-type="doi">10.3389/fvets.2019.00148</pub-id><pub-id pub-id-type="pmid">31192234</pub-id></citation></ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Bankovics</surname> <given-names>A</given-names></name></person-group>. <article-title>Status of wild pigeons and doves in Hungary</article-title>. <source>Naturzale.</source> (<year>2001</year>) <volume>16</volume>:<fpage>61</fpage>&#x02013;<lpage>70</lpage>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://core.ac.uk/download/pdf/11502607.pdf">https://core.ac.uk/download/pdf/11502607.pdf</ext-link></citation>
</ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Powell</surname> <given-names>EC</given-names></name> <name><surname>Hollander</surname> <given-names>WF</given-names></name></person-group>. <article-title><italic>Trichomonas gallinae</italic> infections in the ringdove (<italic>Streptopelia risoria</italic>)</article-title>. <source>J Wildl Dis.</source> (<year>1982</year>) <volume>18</volume>:<fpage>89</fpage>&#x02013;<lpage>90</lpage>. <pub-id pub-id-type="doi">10.7589/0090-3558-18.1.89</pub-id><pub-id pub-id-type="pmid">7097875</pub-id></citation></ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stimmelmayr</surname> <given-names>R</given-names></name> <name><surname>Stefani</surname> <given-names>LM</given-names></name> <name><surname>Thrall</surname> <given-names>MA</given-names></name> <name><surname>Landers</surname> <given-names>K</given-names></name> <name><surname>Revan</surname> <given-names>F</given-names></name> <name><surname>Miller</surname> <given-names>A</given-names></name> <etal/></person-group>. <article-title>Trichomonosis in free-ranging Eurasian collared doves (<italic>Streptopelia decaocto</italic>) and African collared dove hybrids (<italic>Streptopelia risoria</italic>) in the Caribbean and description of ITS-1 region genotypes</article-title>. <source>Avian Dis.</source> (<year>2012</year>) <volume>56</volume>:<fpage>441</fpage>&#x02013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.1637/9905-082311-Case.1</pub-id><pub-id pub-id-type="pmid">22856210</pub-id></citation></ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grabensteiner</surname> <given-names>E</given-names></name> <name><surname>Bilic</surname> <given-names>I</given-names></name> <name><surname>Kolbe</surname> <given-names>T</given-names></name> <name><surname>Hess</surname> <given-names>M</given-names></name></person-group>. <article-title>Molecular analysis of clonal trichomonad isolates indicate the existence of heterogenic species present in different birds and within the same host</article-title>. <source>Vet Parasitol.</source> (<year>2010</year>) <volume>172</volume>:<fpage>53</fpage>&#x02013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1016/j.vetpar.2010.04.015</pub-id><pub-id pub-id-type="pmid">20471174</pub-id></citation></ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kocan</surname> <given-names>RM</given-names></name> <name><surname>Knisley</surname> <given-names>JO</given-names></name></person-group>. <article-title>Challenge infection as a means of determining the rate of disease resistant <italic>Trichomonas gallinae</italic>-free birds in a population</article-title>. <source>J Wildl Dis.</source> (<year>1970</year>) <volume>6</volume>:<fpage>13</fpage>&#x02013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.7589/0090-3558-6.1.13</pub-id><pub-id pub-id-type="pmid">16509124</pub-id></citation></ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mart&#x000ED;nez-Herrero</surname> <given-names>MC</given-names></name> <name><surname>Sansano-Maestre</surname> <given-names>J</given-names></name> <name><surname>Azami-Conesa</surname> <given-names>I</given-names></name> <name><surname>Gonz&#x000E1;lez-Gonz&#x000E1;lez</surname> <given-names>F</given-names></name> <name><surname>Su&#x000E1;rez Regalado</surname> <given-names>L</given-names></name> <name><surname>Garijo-Toledo</surname> <given-names>MM</given-names></name> <etal/></person-group>. <article-title>Sequence subtyping of <italic>Trichomonas gallinae</italic> from Bonelli&#x00027;s eagle (<italic>Aquila fasciata</italic>) during four years (2014&#x02013;2017) reveals that MLS type is associated with lesions</article-title>. <source>Avian Pathol.</source> (<year>2021</year>) <volume>50</volume>:<fpage>339</fpage>&#x02013;<lpage>49</lpage>. <pub-id pub-id-type="doi">10.1080/03079457.2021.1940099</pub-id><pub-id pub-id-type="pmid">34126819</pub-id></citation></ref>
</ref-list> 
</back>
</article>