The taxonomic status of the endangered thin-spined porcupine, Chaetomys subspinosus (Olfers, 1818), based on molecular and karyologic data

Background: The thin-spined porcupine, also known as the bristle-spined rat, Chaetomys subspinosus (Olfers, 1818), the only member of its genus, figures among Brazilian endangered species. In addition to being threatened, it is poorly known, and even its taxonomic status at the family level has long been controversial. The genus Chaetomys was originally regarded as a porcupine in the family Erethizontidae, but some authors classified it as a spiny-rat in the family Echimyidae. Although the dispute seems to be settled in favor of the erethizontid advocates, further discussion of its affinities should be based on a phylogenetic framework. In the present study, we used nucleotide-sequence data from the complete mitochondrial cytochrome b gene and karyotypic information to address this issue. Our molecular analyses included one individual of Chaetomys subspinosus from the state of Bahia in northeastern Brazil, and other hystricognaths. Results: All topologies recovered in our molecular phylogenetic analyses strongly supported Chaetomys subspinosus as a sister clade of the erethizontids. Cytogenetically, Chaetomys subspinosus showed 2n = 52 and FN = 76. Although the sexual pair could not be identified, we assumed that the X chromosome is biarmed. The karyotype included 13 large to medium metacentric and submetacentric chromosome pairs, one small subtelocentric pair, and 12 small acrocentric pairs. The subtelocentric pair 14 had a terminal secondary constriction in the short arm, corresponding to the nucleolar organizer region (Ag-NOR), similar to the erethizontid Sphiggurus villosus, 2n = 42 and FN = 76, and different from the echimyids, in which the secondary constriction is interstitial. Conclusion: Both molecular phylogenies and karyotypical evidence indicated that Chaetomys is closely related to the Erethizontidae rather than to the Echimyidae, although in a basal position relative to the rest of the Erethizontidae. The high levels of molecular and morphological divergence suggest that Chaetomys belongs to an early radiation of the Erethizontidae that may have occurred in the Early Miocene, and should be assigned to its own subfamily, the Chaetomyinae.

Deep sequencing of New World screw-worm transcripts to discover genes involved in insecticide resistance

Background: The New World screw-worm (NWS), Cochliomyia hominivorax, is one of the most important myiasis-causing flies, causing severe losses to the livestock industry. In its current geographical distribution, this species has been controlled by the application of insecticides, mainly organophosphate (OP) compounds, but a number of lineages have been identified that are resistant to such chemicals. Despite its economic importance, only limited genetic information is available for the NWS. Here, as a part of an effort to characterize the C. hominivorax genome and identify putative genes involved in insecticide resistance, we sampled its transcriptome by deep sequencing of polyadenylated transcripts using the 454 sequencing technology. Results: Deep sequencing on the 454 platform of three normalized libraries (larval, adult male and adult female) generated a total of 548,940 reads. Eighteen candidate genes coding for three metabolic detoxification enzyme families, cytochrome P450 monooxygenases, glutathione S transferases and carboxyl/cholinesterases were selected and gene expression levels were measured using quantitative real-time polymerase chain reaction (qRT-PCR). Of the investigated candidates, only one gene was expressed differently between control and resistant larvae with, at least, a 10-fold down-regulation in the resistant larvae. The presence of mutations in the acetylcholinesterase (target site) and carboxylesterase E3 genes was investigated and all of the resistant flies presented E3 mutations previously associated with insecticide resistance. Conclusions: Here, we provided the largest database of NWS expressed sequence tags that is an important resource, not only for further studies on the molecular basis of the OP resistance in NWS fly, but also for functional and comparative studies among Calliphoridae flies. Among our candidates, only one gene was found differentially expressed in resistant individuals, and its role on insecticide resistance should be further investigated. Furthermore, the absence of mutations in the OP target site and the high frequency of mutant carboxylesterase E3 indicate that metabolic resistance mechanisms have evolved predominantly in this species.

Detection of polymorphisms of the mtDNA control region of Caretta caretta (Testudines: Cheloniidae) by PCR-SSCP

Marine turtles are increasingly being threatened worldwide by anthropogenic activities. Better understanding of their life cycle, behavior and population structure is imperative for the design of adequate conservation strategies. The mtDNA control region is a fast-evolving matrilineal marker that has been employed in the study of marine turtle populations. We developed and tested a simple molecular tracing system for Caretta caretta mtDNA haplotypes by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Using this technique, we were able to distinguish the SSCP patterns of 18 individuals of the haplotypes CC-A4, CC-A24 and CCxLO, which are commonly found in turtles sampled on the Brazilian coast. When we analyzed 15 turtles with previously unknown sequences, we detected two other haplotypes, in addition to the other four. Based on DNA sequencing, they were identified as the CC-A17 and CC-A1 haplotypes. Further analyses were made with the sea turtles, Chelonia mydas (N = 8), Lepidochelys olivacea (N = 3) and Eretmochelys imbricata (N = 1), demonstrating that the PCR-SSCP technique is able to distinguish intra-and interspecific variation in the family Cheloniidae. We found that this technique can be useful for identifying sea turtle mtDNA haplotypes, reducing the need for sequencing.

Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

Control of the Intracellular Redox State by Glucose Participates in the Insulin Secretion Mechanism

Background: Production of reactive oxygen species (ROS) due to chronic exposure to glucose has been associated with impaired beta cell function and diabetes. However, physiologically, beta cells are well equipped to deal with episodic glucose loads, to which they respond with a fine tuned glucose-stimulated insulin secretion (GSIS). In the present study, a systematic investigation in rat pancreatic islets about the changes in the redox environment induced by acute exposure to glucose was carried out. Methodology/Principal Findings: Short term incubations were performed in isolated rat pancreatic islets. Glucose dose- and time-dependently reduced the intracellular ROS content in pancreatic islets as assayed by fluorescence in a confocal microscope. This decrease was due to activation of pentose-phosphate pathway (PPP). Inhibition of PPP blunted the redox control as well as GSIS in a dose-dependent manner. The addition of low doses of ROS scavengers at high glucose concentration acutely improved beta cell function. The ROS scavenger N-acetyl-L-cysteine increased the intracellular calcium response to glucose that was associated with a small decrease in ROS content. Additionally, the presence of the hydrogen peroxide-specific scavenger catalase, in its membrane-permeable form, nearly doubled glucose metabolism. Interestingly, though an increase in GSIS was also observed, this did not match the effect on glucose metabolism. Conclusions: The control of ROS content via PPP activation by glucose importantly contributes to the mechanisms that couple the glucose stimulus to insulin secretion. Moreover, we identified intracellular hydrogen peroxide as an inhibitor of glucose metabolism intrinsic to rat pancreatic islets. These findings suggest that the intracellular adjustment of the redox environment by glucose plays an important role in the mechanism of GSIS.

SCFAs Induce Mouse Neutrophil Chemotaxis through the GPR43 Receptor

Short chain fatty acids (SCFAs) have recently attracted attention as potential mediators of the effects of gut microbiota on intestinal inflammation. Some of these effects have been suggested to occur through the direct actions of SCFAs on the GPR43 receptor in neutrophils, though the precise role of this receptor in neutrophil activation is still unclear. We show that mouse bone marrow derived neutrophils (BMNs) can chemotax effectively through polycarbonate filters towards a source of acetate, propionate or butyrate. Moreover, we show that BMNs move with good speed and directionality towards a source of propionate in an EZ-Taxiscan chamber coated with fibrinogen. These effects of SCFAs were mimicked by low concentrations of the synthetic GPR43 agonist phenylacetamide-1 and were abolished in GPR43(-/-) BMNs. SCFAs and phenylacetamide-1 also elicited GPR43-dependent activation of PKB, p38 and ERK and these responses were sensitive to pertussis toxin, indicating a role for Gi proteins. Phenylacetamide-1 also elicited rapid and transient activation of Rac1/2 GTPases and phosphorylation of ribosomal protein S6. Genetic and pharmacological intervention identified important roles for PI3K gamma, Rac2, p38 and ERK, but not mTOR, in GPR43-dependent chemotaxis. These results identify GPR43 as a bona fide chemotactic receptor for neutrophils in vitro and start to define important elements in its signal transduction pathways.

Normal and defective mariner-like elements in Rhynchosciara species (Sciaridae, Diptera)

Mariner-like elements are widely present in diverse organisms. These elements constitute a large fraction of the eukaryotic genome; they transpose by a ""cut-and-paste"" mechanism with their own transposase protein. We found two groups of mobile elements in the genus Rhynchosciara. PCR using primers designed from R. americana transposons (Ramar1 and Ramar2) were the starting point for this comparative study. Genomic DNA templates of four species: R. hollaenderi, R. millerii, R. baschanti, and Rhynchosciara sp were used and genomic sequences were amplified, sequenced and the molecular structures of the elements characterized as being putative mariner-like elements. The first group included the putative full-length elements. The second group was composed of defective mariner elements that contain a deletion overlapping most of the internal region of the transposase open reading frame. They were named Rmar1 (type 1) and Rmar2 (type 2), respectively. Many conserved amino acid blocks were identified, as well as a specific D,D(34) D signature motif that was defective in some elements. Based on predicted transposase sequences, these elements encode truncated proteins and are phylogenetically very close to mariner-like elements of the mauritiana subfamily. The inverted terminal repeat sequences that flanked the mariner-like elements are responsible for their mobility. These inverted terminal repeat sequences were identified by inverse PCR.

Hormone-regulated expression and distribution of versican in mouse uterine tissues

Background: Remodeling of the extracellular matrix is one of the most striking features observed in the uterus during the estrous cycle and after hormone replacement. Versican (VER) is a hyaluronan-binding proteoglycan that undergoes RNA alternative splicing, generating four distinct isoforms. This study analyzed the synthesis and distribution of VER in mouse uterine tissues during the estrous cycle, in ovariectomized (OVX) animals and after 17beta-estradiol (E2) and medroxyprogesterone (MPA) treatments, either alone or in combination. Methods: Uteri from mice in all phases of the estrous cycle, and animals subjected to ovariectomy and hormone replacement were collected for immunoperoxidase staining for versican, as well as PCR and quantitative Real Time PCR. Results: In diestrus and proestrus, VER was exclusively expressed in the endometrial stroma. In estrus and metaestrus, VER was present in both endometrial stroma and myometrium. In OVX mice, VER immunoreaction was abolished in all uterine tissues. VER expression was restored by E2, MPA and E2+MPA treatments. Real Time PCR analysis showed that VER expression increases considerably in the MPA-treated group. Analysis of mRNA identified isoforms V0, V1 and V3 in the mouse uterus. Conclusion: These results show that the expression of versican in uterine tissues is modulated by ovarian steroid hormones, in a tissue-specific manner. VER is induced in the myometrium exclusively by E2, whereas MPA induces VER deposition only in the endometrial stroma.

A systematic review with procedural assessments and meta-analysis of Low Level Laser Therapy in lateral elbow tendinopathy (tennis elbow)

Background: Recent reviews have indicated that low level level laser therapy (LLLT) is ineffective in lateral elbow tendinopathy (LET) without assessing validity of treatment procedures and doses or the influence of prior steroid injections. Methods: Systematic review with meta-analysis, with primary outcome measures of pain relief and/or global improvement and subgroup analyses of methodological quality, wavelengths and treatment procedures. Results: 18 randomised placebo-controlled trials (RCTs) were identified with 13 RCTs (730 patients) meeting the criteria for meta-analysis. 12 RCTs satisfied half or more of the methodological criteria. Publication bias was detected by Egger's graphical test, which showed a negative direction of bias. Ten of the trials included patients with poor prognosis caused by failed steroid injections or other treatment failures, or long symptom duration or severe baseline pain. The weighted mean difference (WMD) for pain relief was 10.2 mm [95% CI: 3.0 to 17.5] and the RR for global improvement was 1.36 [1.16 to 1.60]. Trials which targeted acupuncture points reported negative results, as did trials with wavelengths 820, 830 and 1064 nm. In a subgroup of five trials with 904 nm lasers and one trial with 632 nm wavelength where the lateral elbow tendon insertions were directly irradiated, WMD for pain relief was 17.2 mm [95% CI: 8.5 to 25.9] and 14.0 mm [95% CI: 7.4 to 20.6] respectively, while RR for global pain improvement was only reported for 904 nm at 1.53 [95% CI: 1.28 to 1.83]. LLLT doses in this subgroup ranged between 0.5 and 7.2 Joules. Secondary outcome measures of painfree grip strength, pain pressure threshold, sick leave and follow-up data from 3 to 8 weeks after the end of treatment, showed consistently significant results in favour of the same LLLT subgroup (p < 0.02). No serious side-effects were reported. Conclusion: LLLT administered with optimal doses of 904 nm and possibly 632 nm wavelengths directly to the lateral elbow tendon insertions, seem to offer short-term pain relief and less disability in LET, both alone and in conjunction with an exercise regimen. This finding contradicts the conclusions of previous reviews which failed to assess treatment procedures, wavelengths and optimal doses.

Genes up- and down-regulated by dermcidin in breast cancer: a microarray analysis

Dermcidin (DCD) is a human gene mapped to chromosome 12q13 region, which is co-amplified with multiple oncogenes with a well-established role in the growth, survival and progression of breast cancers. Here, we present a summary of a DNA microarray-based study that identified the genes that are up- and down-regulated in a human MDA-361 pLKO control clone and three clones expressing short hairpin RNA against three different regions of DCD mRNA. A list of 235 genes was differentially expressed among independent clones (> 3-fold change and P < 0.005). The gene expression of 208 was reduced and of 27 was increased in the three DCD-RNAi clones compared to pLKO control clone. The expression of 77 genes (37%) encoding for enzymes involved in amino acid metabolism, glucose metabolism and oxidoreductase activity and several genes required for cell survival and DNA repair were decreased. The expression of EGFR/ErbB-1 gene, an important predictor of outcome in breast cancer, was reduced together with the genes for betacellulin and amphiregulin, two known ligands of EGFR/ErbB receptors. Many of the 27 genes up-regulated by DCD-RNAi expression have not yet been fully characterized; among those with known function, we identified the calcium-calmodulin-dependent protein kinase-II delta and calcineurin A alpha. We compared 132 up-regulated and 12 down-regulated genes in our dataset with those genes up- and down-regulated by inhibitors targeting various signaling pathway components. The analysis showed that the genes in the DCD pathway are aligned with those functionally influenced by the drugs sirolimus, LY-294002 and wortmannin. Therefore, DCD may exert its function by activating the PI3K/AKT/mTOR signaling pathway. Together, these bioinformatic approaches suggest the involvement of DCD in the regulation of genes for breast cancer cell metabolism, proliferation and survival.

The case of thyroid hormones: how to learn physiology by solving a detective case

Lellis-Santos C, Giannocco G, Nunes MT. The case of thyroid hormones: how to learn physiology by solving a detective case. Adv Physiol Educ 35: 219-226, 2011; doi:10.1152/advan.00135.2010.Thyroid diseases are prevalent among endocrine disorders, and careful evaluation of patients' symptoms is a very important part in their diagnosis. Developing new pedagogical strategies, such as problem-based learning (PBL), is extremely important to stimulate and encourage medical and biomedical students to learn thyroid physiology and identify the signs and symptoms of thyroid dysfunction. The present study aimed to create a new pedagogical approach to build deep knowledge about hypo-/hyperthyroidism by proposing a hands-on activity based on a detective case, using alternative materials in place of laboratory animals. After receiving a description of a criminal story involving changes in thyroid hormone economy, students collected data from clues, such as body weight, mesenteric vascularization, visceral fat, heart and thyroid size, heart rate, and thyroid-stimulating hormone serum concentration to solve the case. Nevertheless, there was one missing clue for each panel of data. Four different materials were proposed to perform the same practical lesson. Animals, pictures, small stuffed toy rats, and illustrations were all effective to promote learning, and the detective case context was considered by students as inviting and stimulating. The activity can be easily performed independently of the institution's purchasing power. The practical lesson stimulated the scientific method of data collection and organization, discussion, and review of thyroid hormone actions to solve the case. Hence, this activity provides a new strategy and alternative materials to teach without animal euthanization.

Complete Genome Sequence of Mycoplasma suis and Insights into Its Biology and Adaption to an Erythrocyte Niche

Mycoplasma suis, the causative agent of porcine infectious anemia, has never been cultured in vitro and mechanisms by which it causes disease are poorly understood. Thus, the objective herein was to use whole genome sequencing and analysis of M. suis to define pathogenicity mechanisms and biochemical pathways. M. suis was harvested from the blood of an experimentally infected pig. Following DNA extraction and construction of a paired end library, whole-genome sequencing was performed using GS-FLX (454) and Titanium chemistry. Reads on paired-end constructs were assembled using GS De Novo Assembler and gaps closed by primer walking; assembly was validated by PFGE. Glimmer and Manatee Annotation Engine were used to predict and annotate protein-coding sequences (CDS). The M. suis genome consists of a single, 742,431 bp chromosome with low G+C content of 31.1%. A total of 844 CDS, 3 single copies, unlinked rRNA genes and 32 tRNAs were identified. Gene homologies and GC skew graph show that M. suis has a typical Mollicutes oriC. The predicted metabolic pathway is concise, showing evidence of adaptation to blood environment. M. suis is a glycolytic species, obtaining energy through sugars fermentation and ATP-synthase. The pentose-phosphate pathway, metabolism of cofactors and vitamins, pyruvate dehydrogenase and NAD(+) kinase are missing. Thus, ribose, NADH, NADPH and coenzyme A are possibly essential for its growth. M. suis can generate purines from hypoxanthine, which is secreted by RBCs, and cytidine nucleotides from uracil. Toxins orthologs were not identified. We suggest that M. suis may cause disease by scavenging and competing for host nutrients, leading to decreased life-span of RBCs. In summary, genome analysis shows that M. suis is dependent on host cell metabolism and this characteristic is likely to be linked to its pathogenicity. The prediction of essential nutrients will aid the development of in vitro cultivation systems.

Genome Size, Karyotype Polymorphism and Chromosomal Evolution in Trypanosoma cruzi

Background: The Trypanosoma cruzi genome was sequenced from a hybrid strain (CL Brener). However, high allelic variation and the repetitive nature of the genome have prevented the complete linear sequence of chromosomes being determined. Determining the full complement of chromosomes and establishing syntenic groups will be important in defining the structure of T. cruzi chromosomes. A large amount of information is now available for T. cruzi and Trypanosoma brucei, providing the opportunity to compare and describe the overall patterns of chromosomal evolution in these parasites. Methodology/Principal Findings: The genome sizes, repetitive DNA contents, and the numbers and sizes of chromosomes of nine strains of T. cruzi from four lineages (TcI, TcII, TcV and TcVI) were determined. The genome of the TcI group was statistically smaller than other lineages, with the exception of the TcI isolate Tc1161 (Jose-IMT). Satellite DNA content was correlated with genome size for all isolates, but this was not accompanied by simultaneous amplification of retrotransposons. Regardless of chromosomal polymorphism, large syntenic groups are conserved among T. cruzi lineages. Duplicated chromosome-sized regions were identified and could be retained as paralogous loci, increasing the dosage of several genes. By comparing T. cruzi and T. brucei chromosomes, homologous chromosomal regions in T. brucei were identified. Chromosomes Tb9 and Tb11 of T. brucei share regions of syntenic homology with three and six T. cruzi chromosomal bands, respectively. Conclusions: Despite genome size variation and karyotype polymorphism, T. cruzi lineages exhibit conservation of chromosome structure. Several syntenic groups are conserved among all isolates analyzed in this study. The syntenic regions are larger than expected if rearrangements occur randomly, suggesting that they are conserved owing to positive selection. Mapping of the syntenic regions on T. cruzi chromosomal bands provides evidence for the occurrence of fusion and split events involving T. brucei and T. cruzi chromosomes.

Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus)microplus in response to infection with Anaplasma marginale

Background: Bovine anaplasmosis, caused by the rickettsial tick-borne pathogen Anaplasma marginale (Rickettsiales: Anaplasmataceae), is vectored by Rhipicephalus (Boophilus) microplus in many tropical and subtropical regions of the world. A. marginale undergoes a complex developmental cycle in ticks which results in infection of salivary glands from where the pathogen is transmitted to cattle. In previous studies, we reported modification of gene expression in Dermacentor variabilis and cultured Ixodes scapularis tick cells in response to infection with A. marginale. In these studies, we extended these findings by use of a functional genomics approach to identify genes differentially expressed in R. microplus male salivary glands in response to A. marginale infection. Additionally, a R. microplus-derived cell line, BME26, was used for the first time to also study tick cell gene expression in response to A. marginale infection. Results: Suppression subtractive hybridization libraries were constructed from infected and uninfected ticks and used to identify genes differentially expressed in male R. microplus salivary glands infected with A. marginale. A total of 279 ESTs were identified as candidate differentially expressed genes. Of these, five genes encoding for putative histamine-binding protein (22Hbp), von Willebrand factor (94Will), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz) and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were confirmed by real-time RT-PCR to be down-regulated in tick salivary glands infected with A. marginale. The impact of selected tick genes on A. marginale infections in tick salivary glands and BME26 cells was characterized by RNA interference. Silencing of the gene encoding for putative flagelliform silk protein (100Silk) resulted in reduced A. marginale infection in both tick salivary glands and cultured BME26 cells, while silencing of the gene encoding for subolesin (4D8) significantly reduced infection only in cultured BME26 cells. The knockdown of the gene encoding for putative metallothionein (93 Meth), significantly up-regulated in infected cultured BME26 cells, resulted in higher A. marginale infection levels in tick cells. Conclusions: Characterization of differential gene expression in salivary glands of R. microplus in response to A. marginale infection expands our understanding of the molecular mechanisms at the tick-pathogen interface. Functional studies suggested that differentially expressed genes encoding for subolesin, putative von Willebrand factor and flagelliform silk protein could play a role in A. marginale infection and multiplication in ticks. These tick genes found to be functionally relevant for tick-pathogen interactions will likely be candidates for development of vaccines designed for control of both ticks and tick-borne pathogens.