2.1 Tick specimens
R. microplus tick specimens of colonies established in the National Laboratory of Parasitology (LNP) from isolations of a small key in the north Cuba named Cayo Coco (22°30'32.45"N 78°24'25.1"W), and the Manga Larga locality (22°04′30″N 78°21′1″W) in the Ciego de Ávila province in Cuba were used to perform the present studies. These tick strains were named Cayo Coco (CC) and Manga Larga (ML), respectively. Tick specimens of a colony kept at the LNP from the Media Joya (MJ) strain of Mexican R. microplus was used as reference for phylogenetic analyses. This Mexican reference tick strain was kindly provided by CENAPA, Mexico and it was originally established at the National Institute of Forestry, Agriculture and Livestock (INIFAP) in Jiutepec, Morelos in 2001 from cattle infested with R. microplus ticks in the municipality of Tapalpa in the state of Jalisco, Mexico (19° 57' 0" N, 103° 46' 0"W) [14].
2.2 Life cycle of R. microplus strains under laboratory conditions
All procedures involving animals and samplings were carried out in accordance with the Guide for the Care and Use of Laboratory Animals [15] and were approved by the Ethic Committee of the LNP.
Tick-free cattle of the Cuban Siboney breed (5/8 Holstein and 3/8 Cebu) from uninfested pastures neither exposed to chemical treatments, aged between 1 and 2 years were individually housed in stalls at the LNP to study the parasitic phase of the CC and ML strains. During the experiment, animals were fed with forage and water ad libitum and supplemented with a pellet diet (produced by CENPALAB, Havana, Cuba). Forty thousand larvae of each tick strain were delivered on the back of bovines to allow free infestations. Three bovines were used for each tick strain. Random tick samples were removed at daily intervals from each bovine. The number of unfed larvae, fed larvae, unfed nymphs, fed nymphs, unfed adults, partially fed adults and fully engorged females were recorded in the daily batches.
In order to study non parasitic stage, twenty fully engorged females collected from each tick strain were placed into individual glass vials and maintained during oviposition period at 28°C and relative humidity of 80%. Egg masses were individually weighted and placed in the same incubation conditions. After hatching, unhatched eggs and eggshells from each egg mass were counted under stereoscope.
2.3 Morphological analysis of CC and ML tick strains
Morphological characterization was carried out by using 40 unfed adult specimens (20 female and 20 male ticks) from CC and ML strains. Unfed adults were collected from bovines around day 15 after larvae infestation, cleaned and fixated with glutaraldehyde 3% during 20 minutes at 4°C. After, they were washed three times with phosphate-buffered saline (PBS -135mM NaCl, 8mM Na2HPO4, 3mM KCL, 1.5mM KH2PO4, and pH 7.2) and finally submerged in 1% of osmium tetroxide (OsO4) during 20 minutes at 4°C. Posteriorly, ticks were again washed three times with PBS and dehydrated in a bath series of increasing concentrations of ethanol from 50% to 99% during 10 min each, at 4°C. At the end, ticks were dried and lyophilized (Freezer Dryer Model FD- 10V) to -62°C and 1.2Pa during 18 hours and coated with a 10 nm gold layer with a Desk Sputter Coater DSR1. Main structures with taxonomic value were characterized by using a TESCAN MIRA-3 FE- Scanning Electron Microscope (SEM) operated at 5kV at the Center for Advanced Studies in Cuba. These structures were analysed and measured by using the Digital Micrograph ™ program (Version 2.32.888.0). Measure averages of these structures from female and male ticks were compared by using a Student's t test performed on Prism (version 6.0 for Windows; GraphPad Software).
2.4 DNA extraction and Polymerase Chain Reaction (PCR)
Egg masses from the CC and ML strains were crushed in a mortar with liquid Nitrogen. Genomic DNA was extracted by using the DNeasy Blood & Tissue Kit (QIAGEN, Germany) according to the manufacturer’s recommendations. The DNA integrity was verified by agarose gel electrophoresis. Primers were designed for specific amplification of 12S rRNA, 16S rRNA, COXI and ITS2 sequences. PCR reactions were performed in a thermal cycler (MinicyclerTM, MJ Research, Inc., USA) by using the System GoTaq® Green Master Mix (Promega, USA). Primers and conditions used in the PCR reactions were described in the Table 1 [7, 9]. The amplified DNA fragments were extracted from agarose gels using the QIAquick Gel Extraction kit (Qiagen, Germany) and sequenced through the services of the Microsynth SeqLab GmbH (Germany).
2.5 Phylogenetic analyses
Sequence alignments were performed by using Clustalw Omega [16]. Phylogenetic trees were constructed with the “Maximum Likelihood (ML)” method based on the Tamura 3-parameter model [17] with 1000 replicates by using the MEGA 7 software [18]. The GenBank accession numbers of sequences included in the analyses were: R. microplus clade A (12S- EU921764, EU921761; 16S- EU918184, EU918183, EU918181, EU918180, EU918187, L34310, EU918178, EU918177, EU918176.1, EU918182; COXI- KP143546, KF200106, KP226177, KP226172, KP226169, KP226160, KP226174; ITS2- KC503273, U97715, EU520392), R. microplus clade B (12S- JF906025, JF906024, JF906019, JF906028, JF906029, JF906021; 16S- JX051068, JF979381, JX051072; COXI- JF758636, JQ625683, JF758630; ITS2- JF758640, JF758642, JQ625705, JQ625709, KC503274, KC203364), R. microplus clade C (12S- MG459958, MG459960, MG459959; 16S- HM536972, HM536977, EU918188, HM536976; COXI- MG459964, KP698516, KP698515, MG459961, MG459962, KP792580, KP318133, KP792579, KP792578, KP792586, MG459963, KP792587, KP792588, KP792589, KM246869, KM246868; ITS2- KC503264, KC503265, MG459967, MG459966, MG459965, KC503276, KC503272, JX974346), R. annulatus (12S- EU921773, KF219716, AF133058, KT335263; 16S- Z97877, L34311; COXI- KM235716, KM235715, KM494917, KF219739, AF132825; ITS2- KJ410770, KC503267) and R. australis (12S- EU921769, EU921767; 16S- EU918185, EU918186, EU918192, EU918190, EU918191; COXI- KC503255, AF132827; ITS2- U97712, KC503268).These sequences were selected only if a species’ validation was available according to the recent definition of the R. microplus complex [7-9]. Sequences from the Amblyomma mixtum tick species were used as outgroups (12S- KF527329; 16S- KT820359; COXI- KT820364 and ITS2- JN866886). The estimates of genetic distances among sequences of each gene were conducted by using the Tamura 3-parameter model [17] performed on Mega7 [18].