Abstract
More than 1000 humans have acquired the febrile disease tularemia in Spain since the first notification of human cases in 1997. We here aimed to study the recent molecular evolution of the causative bacterium Francisella tularensis during disease establishment in Spain. Single-nucleotide polymorphisms (SNPs) and variable-number tandem repeats (VNTRs) were analyzed in whole-genome sequences (WGS) of F. tularensis. Short-read WGS data for 20 F. tularensis strains from humans infected in the periods 2014–2015 and 2018–2020 in Spain were generated. These data were combined with WGS data of 25 Spanish strains from 1998 to 2008 and two reference strains. Capillary electrophoresis data of VNTR genetic regions were generated and compared with the WGS data for the 11 strains from 2014 to 2015. Evolutionary relationships among strains were analyzed by phylogenetic methods. We identified 117 informative SNPs in a 1,577,289-nucleotide WGS alignment of 47 F. tularensis genomes. Forty-five strains from Spain formed a star-like SNP phylogeny with six branches emerging from a basal common node. The most recently evolved genomes formed four additional star-like structures that were derived from four branches of the basal common node. VNTR copy number variation was detected in two out of 10 VNTR regions examined. Genetic clustering of strains by VNTRs agreed with the clustering by SNPs. The SNP data provided higher resolution among strains than the VNTRs data in all but one cases. There was an excellent correlation between VNTR marker sizing by capillary electrophoresis and prediction from WGS data. The genetic data strongly support that tularemia, indeed, emerged recently in Spain. Distinct genetic patterns of local F. tularensis population expansions imply that the pathogen has colonized a previously disease-free geographical area. We also found that genome-wide SNPs provide higher genetic resolution among F. tularensis genomes than the use of VNTRs, and that VNTR copy numbers can be accurately predicted using short-read WGS data.
Highlights
Transfer-RNA-derived fragments are small non-coding RNAs with emerging functions in regulating gene expression [1,2,3,4,5]
In the MINTbase v2.0 database, expression of 10,158 Transfer-RNA-derived fragments (tRFs) were recorded for this cohort, representing 43% of all tRFs reported from all
Levels of tRFs ranged from 1.002 RPM to 61,097.254 RPM
Summary
Transfer-RNA-derived fragments (tRFs) are small non-coding RNAs with emerging functions in regulating gene expression [1,2,3,4,5]. TRFs are produced constitutively or in response to a variety of stress stimuli [6,7,8]. Based on the cleavage position of tRNAs, tRFs are generally classified into six subtypes, with lengths of 14–50 nucleotides [7,9,10]. The endoribonuclease RNase Z processes a precursor tRNA to generate tRF-1 ( known as tsRNA), and angiogenin cleaves mature tRNAs within the anticodon loop to generate 50 -tRNA halves and 30 -tRNA halves. TRNA halves are typically 33 nucleotides in length, whereas other tRFs are of variable length [11,12]. Genes 2020, 11, 1344; doi:10.3390/genes11111344 www.mdpi.com/journal/genes tRFs were first observed in 1977 in human tumor tissue [13] and suggested as markers of cancer [14] Other tRFs generated from the 50 or 30 ends or the internal region of mature tRNA, probably by the action of angiogenin or Dicer, are called 50 -tRFs, 30 -tRFs or i-tRFs, respectively [7,8,10]. tRNA halves are typically 33 nucleotides in length, whereas other tRFs are of variable length [11,12].
Sign-in/Register to view highlights & summary Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
