- Research article
- Open Access
Influence of endogamy and mitochondrial DNA on immunological parameters in cattle
BMC Veterinary Researchvolume 10, Article number: 79 (2014)
Endogamy increases the risk of manifestation of deleterious recessive genes. Mitochondrial DNA allows the separation of American Zebu (Bos indicus and Bos taurus) and evaluate the effect of mitochondrial DNA on productive traits of cattle. However, the effect of endogamy and mitochondrial DNA (mtDNA) on the immune system remains unclear. The aim of this study was to evaluate the association between endogamy, mtDNA and immune parameters.
A total of 86 cattle (43 cows and 43 calves) were used in this study. Age, endogamy, milk yield, and origin of mtDNA were measured and their influence on immunological parameters was evaluated. Older cows had increased CD4+ T cells, decreased CD21+ and γδhigh T cells as well as increased CD4+/CD8+ and T/B ratio. Multiple regression analysis indicated that endogamy in calves was associated with increased CD8+ T and CD21+ B lymphocytes, and decreased γδhigh T cells in peripheral blood. Cows with medium and lower endogamy had a lower percentage of B lymphocytes and γδlow T cells and cows with lower endogamy had higher levels of γδ T cells and γδhigh T cells, as well as the CD4+/CD48+ cell ratio. Calves with higher endogamy had higher levels of CD8+ T lymphocytes, whereas calves with lower endogamy had lower levels of γδlow T cells.
These results demonstrated for the first time that endogamy influences the immune system of cattle.
The degree of endogamy or consanguinity in the offspring increases proportionally to the genetic similarity of its ancestors, in contrast to individuals that are randomly selected from a given population. The endogamy coefficient of an individual indicates the probability that a pair of randomly selected genes of the same locus is completely identical. Therefore, endogamic individuals can only result from mating between closely related ancestors [1, 2].
The study of endogamy in cattle is extremely important, particularly due to the widespread use of artificial insemination, which allows a high selection pressure and rapid genetic gain, resulting in thousands of descendants from one single sire. Despite the unquestionable genetic gain and consequent improvement in productive traits, extensive use of a few bulls may lead to reduced genetic variability and potential deleterious effects caused by endogamy [1, 3].
Endogamy decreases the frequency of heterozygotes in the population, increasing the risk of deleterious effects of recessive genes, which may negatively influence the average productive performance in the population, a phenomenon known as endogamy depression [4, 5]. Therefore, the primary effect of endogamy is an increase the frequency of homozygous genes, favoring manifestation of several recessive genes, which usually cause impairment on the average individual phenotypic merit. Possible causes for declining of phenotypic value as a consequence of endogamy includes the fact that favorable genes tend to be dominant or partially dominant [6–8]. Deleterious effects of endogamy on the reproductive ability of cows and heifers have been reported [1, 2, 9, 10]. However, whether endogamy influences the immune system remains to be investigated. Importantly, some studies report an association between low genetic diversity and increased susceptibility to diseases [8, 11].
Evaluation of mitocondrial DNA allows separation of American Zebu cattle, according to its maternal lineage ancestry, into two groups: one with only Bos taurus indicus mtDNA (Bi-mtDNA) and other with Bos taurus taurus mtDNA (Bt-mtDNA) . Studies to evaluate the effect of mitochondrial DNA on productive and reproductive traits in Gir and Guzerat have been conducted and show that the mitochondrial origin (taurus or indicus) significantly affects the age at first calving  and it does not seem to have significant effects on lactation milk yield, days in milk, or calving interval .
Considering the scarcity of studies regarding the effect of endogamy and mtDNA on immunological parameters, the goal of this study was to evaluate the relationship association between endogamy and the immune system. Our focus was to understand variations of immune parameters, including CD4+ T cells, CD8+ T cells, γδ T cells, and B cells (CD21+) in peripheral blood, that may be associated to endogamy in Guzerat (Zebu) cattle.
Immunological parameters are influenced by endogamy
Multiple regression models including milk yield, age, mtDNA, and endogamy coefficient for cows, and mtDNA and endogamy coefficient for calves in relation to levels of lymphocyte subsets in the peripheral blood were performed and the results are summarized in Tables 1 and 2. The endogamy coefficient was significantly associated with increased CD21+ levels (p < 0.05; Table 1). In calves, multiple regression analysis indicated that endogamy was significantly associated with increased CD8+ T and CD21+ cells, and with significantly decreased γδhigh T cells and T/B cell ratio (p < 0.05; Table 2). Importantly, the age of cows was significantly associated with increased CD4+ T cells, CD4+/CD8+ and T/B cell ratio, and with significantly decreased CD21+ and γδhigh T cells (p < 0.05; Table 1).
Regression analysis was also performed to verify possible associations between the endogamy coefficient and immunological parameters. In cows, there were no statistically significant correlations. However, in calves there was a positive correlation between endogamy and CD8+ T lymphocytes (p < 0.05, r = 0.3080) and CD21+ B cells (p < 0.05, r = 0.3530), and a negative correlation with γδhigh T cells (p < 0.05, r = –0.3113) and T/B cell ratio (p < 0.05, r = –0.3334).
Endogamy alters the immune profile of Guzerat cattle
Based on endogamy coefficients, cows and calves were grouped into high, medium, and low endogamy as described in Methods. The immune profiles of these cattle were then analyzed by flow cytometry (Figures 1 and 2). These results indicated that cows with medium and low endogamy had a decreased percentage of B lymphocytes when compared to cows with high endogamy (p < 0.05; Figure 1C). Furthermore, cows with lower endogamy had a higher percentage of γδlow T cells when compared with cows with medium endogamy (p < 0.05; Figure 1F). There was also a clear influence of endogamy on the immunological profile of calves. Calves with high endogamy had an increased percentage of CD8+ T lymphocytes when compared to calves with medium and low endogamy (p < 0.05, Figure 2B). In contrast, calves with low endogamy had an increased percentage of γδT cells (p < 0.05; Figure 2D) and γδhigh T cells (p < 0.05; Figure 2E) as well as the CD4+/CD8+ T cell ratio (p < 0.05; Figure 2G) and T/B ratio (p < 0.05; Figure 2H) when compared to calves with high endogamy (p < 0.05). Calves with low endogamy had a decreased percentage of γδlow T cells only in comparison to those with medium endogamy (p < 0.05; Figure 2F).
European mtDNA in Guzerat is associated with low endogamy and differences in immunological parameters
All cattle used in this study were mtDNA genotyped. This analysis allowed us to assess whether individual Zebu cattle (Bos taurus indicus) remained with no introduction of maternal DNA of B. taurus taurus origin, thus tending to have higher levels of endogamy . Our results demonstrated that 42% (18/43) of the cows had mtDNA from B. taurus taurus, and 58% (25/43) had only B. taurus indicus (Zebu) mtDNA. The same result was observed in calves, which was expected since these calves were the offspring of those same cows mentioned above. After grouping cows and calves according to the presence or absence of B. taurus taurus (European) mtDNA, which were classified into “Bt-mtDNA” or “Bi-mtDNA”, respectively, the profile of leukocytes in peripheral blood was analyzed by flow cytometry. In cows, the percentages of B lymphocytes, γδT cells, and γδhigh T cells were decreased in individuals with Bi-mtDNA when compared to individuals with Bt-mtDNA (p < 0.05; Figure 3C, D, and E, respectively). In calves, individuals with Bi-mtDNA had a decreased percentage of γδT cells and γδhigh T cells when compared to individuals with Bt-mtDNA (p < 0.05; Figure 4D and E, respectively). As expected, the endogamy coefficient was lower in calves with Bt-mtDNAwhen compared to calves with Bi-mtDNA (p < 0.05; Figure 5).
Here we demonstrated a clear effect of endogamy on immune parameters of Guzerat cattle. This is the first evidence that endogamy may affect the immunological parameters in cattle. The main consequences of endogamy increase are the reduction of genetic variability within inbred lines and increased frequency of homozygosity, which may be desirable for alleles with favorable effects. However, a small fraction of undesirable deleterious recessive alleles may have an enhancement of their manifestation under these conditions. Therefore, endogamy tends to increase expression of undesirable recessive alleles .
Deleterious effects of endogamy have been demonstrated, especially on reproductive traits [1, 2, 9, 10]. A study of inbreed and outbreed populations of several different species of captive ungulates indicates that the rate of neonatal mortality is higher in endogamic when compared to non endogamic animals . In humans, consanguinity is an important risk factor for susceptibility to infectious diseases . However, direct associations of endogamy with the immune system have not been previously investigated in animals.
Higher levels of endogamy clearly affected some of the immunological parameters, particularly γδ T cells. These cells are the major circulating lymphocyte subset population in ruminants, especially in calves. These cells play a role in the transition between innate and adaptive immune . γδ T lymphocytes can recognize and be activated by pathogen-associated molecular patterns (PAMPs) or danger-associated molecular patterns (DAMPs) in absence of other stimuli [18, 19], supporting the notion of an innate immune function of these cells. However, other studies demonstrated an important role of γδ T lymphocytes in adaptive response since they recognize antigens presented on the surface of antigen-presenting cells (APCs), and they have cytotoxic effect mediated by expression of granulysin  and perforin [19, 21]. These cells play an important role during infection with Mycobacterium bovis, the causative agent of bovine tuberculosis [21–23]. Recent studies have also demonstrated the role of γδ T lymphocytes as APCs [23–25] so they may play a role in developing protective immune responses following vaccination , and are important for host defense during Brucella abortus infection, when these cells enhance macrophage function via IFN-γ .
Although resistance or susceptibility to disease involves multifactorial mechanisms, it has been reported that Zebu cattle tend to have a higher degree of natural resistance to diseases when compared to European cattle [27–31]. Analyzing the profile of peripheral blood leukocytes of different bovine breeds, we recently demonstrated that Guzerat cattle comparatively has increased levels of circulating leukocytes, especially monocytes, eosinophils, T lymphocytes (particularly CD4+ T cells) and non-B/non-T (NTNB) lymphocytes, with a decreased percentage of B lymphocytes, and consequent increase in the T/B and CD4+/CD8+ ratio . In another study, we compared the macrophage function of Zebu cattle with that of macrophages from European cattle. Macrophages from Zebu are more efficient in controlling intracellular replication of Brucella abortus. Furthermore, Zebu macrophages generate higher levels of nitric oxide and cytokines such as IL-12 and TNF-α . Therefore, increasing lines of evidences support the notion that resistance of Zebu cattle may be due to quantitative and qualitative differences in immune cells. In this study, we addressed for the first time whether endogamy have any influence on the profile of immune cells in the peripheral blood of Zebu cattle.
Here we demonstrated that age affected most of the immunological parameters of cows. The influence of age on the immune system has already been demonstrated in both animals and humans, affecting both innate and adaptive immune systems [32–35]. Age-related decrease in B cell activity, potentially related to the decline in helper T cell activity, and decreased cytolytic and delayed-type hypersensitivity responses have been described [32, 36]. The age-related alteration in helper T cell activity appears to be due to intrinsic or functional changes within the cells, as well as to a shift in subset proportions of CD4+ and CD8+ T cell populations, with increases in the memory phenotype and corresponding decreases in the naive phenotype .
In conclusion, our data indicate that endogamy influences the immune system of cattle since it is associated with significant changes in the profile of immune cells in the peripheral blood, particularly decreased levels of γδ T lymphocytes, an important subpopulation of lymphocytes in cattle, especially in calves.
Cattle and samples
A total of 43 Guzerat (Bos taurus indicus) cows with an average of 8 (3–17) years of age and their progeny (a total of 43 male or female calves) with an average of 12 (2.4–25) month-old were included in this study. The average daily milk yield of cows was 18.4 (4.4–31) liters per day. These cattle were maintained under regular management conditions in their farm. They were free of ticks and hemoparasites. Blood samples were collected within 5 mL tubes containing ethylene diamine tetra acetic acid (EDTA) for ex vivo cellular evaluation, and stored at room temperature until further processing. This study was approved by the Ethical Committee for the use of Experimental Animals of the Universidade Federal de Minas Gerais, Brazil (CETEA), under protocol 139/2010.
Analysis of endogamy coefficient
All cattle included in this study were subjected to artificial insemination with a highly stringent breeding control, and have registered Bos indicus breed Guzerat pedigrees since 1895 (over 100 years) that guarantee a very low risk of pedigree error. The individual inbreeding coefficients were calculated by inverting the diagonal of the inverse matrix of kinship using the set of multiple-trait derivative-free - MTDFREML software package . The software calculates the inverse of the relationship matrix directly from a list of animals and their parents, providing individual identification for matching phenotypic records to individuals and calculating the inbreeding coefficients as well as the logarithm of the determinant of the relationship matrix needed to calculate the logarithm of the likelihood function. The program also prepares the coefficients for the mixed model equations based on the statistical model for single and multiple trait analyses. Then, the program solves the mixed model equations and finds variance component estimates that maximize the restricted likelihood given the phenotypic data. Based on the inbreeding coefficients, the animals were grouped at low, medium and high endogamy. Endogamy coefficients were first ordered from the highest to the lowest and the upper third was considered high endogamy (coefficients ranging from 26.4 to 16.3, and 26.1 to 18.1, for cows and calves, respectively), the mid third was considered medium endogamy (coefficients ranging from 16.2 to 14.2, and 18.0 to 15.4, for cows and calves, respectively), and the lower third was considered low endogamy (coefficients ranging from 14.0 to 9.8, and 15.3 to 0.0, for cows and calves, respectively).
Mitochondrial DNA analysis
Analysis of mitochondrial DNA (mtDNA) was performed as previously described . Polymorphisms for B. indicus and B. taurus mtDNA were identified using two other RFLP markers. Briefly, total DNA samples from several animals were re-amplified by PCR using primers 5’-CCCAACGAGGAAAATATACC-3’ (BosmtF1) and 5’-AACCGCAAACAACCTCTTCC-3’ (BosmtR1) targeting a region of the ND5 gene of the mitochondrial genome. Amplified mtDNA was then digested HindIII for 1 h at 37°C, electrophoresed on 1.5% agarose gel, stained with ethidium bromide and evaluated with a Fuji Fla 3000G Laser Scanner (Fuji Film Co., Tokyo, Japan), using Image Gauge version 3.12 software. The presence of the HindIII restriction site within the amplified region is indicative of B. taurus taurus mtDNA. Efficiency of HindIII digestion was assured by the complete absence of partially digested fragments.
Immunophenotyping of peripheral lymphocytes
Immunophenotyping was performed using specific monoclonal anti-bovine cell receptor antibodies as previously described . Briefly, these antibodies were tittered, and 1:15, 1:30, 1:60, and 1:120 dilutions were used for anti-CD4, anti-CD8, γδ T and anti-CD21 antibodies, respectively. A single cell type was analyzed in each tube. After incubation with primary antibodies, erythrocytes were lysed by incubating with 2 mL lysis solution (FACS™ Lysing Solution – BD Biosciences, San Jose, CA, USA) for 8 min at room temperature. After lysis, leukocytes were washed in PBS and fixed with 200 μL Max FaxFix. Acquisition, storage, and analysis of the data were performed using a FACScalibur flow cytometer (BD Biosciences, San Jose, CA, USA) with CellQuest software (BD Biosciences, San Jose, CA, USA). Leukocyte subpopulations were quantified based on size (laser forward scatter; FSC) and granularity. Lymphocyte subsets were analyzed considering cells within the selected lymphocyte region based on the relative fluorescence intensity observed in the FL1 versus FL2 dot plot distributions. Quadrant statistical analysis was applied to quantify the fluorescent positive lymphocyte subset within the lower right quadrant. Calculations of the CD4+/CD8+ T cell and T/B cell ratio and percentage of γδhigh T andγδlow T cells (based on the expression profile among the total population of γδ T cells) were also performed.
Several parameters including age, milk yield, mtDNA, and the endogamy coefficient were considered in these analyses. Multiple regression analyzes and linear correlation were performed using the Graphpad Instat3 software (GraphPad Software, San Diego, CA, USA).
Flow cytometric data were submitted to analysis of variance and means were compared by the Student-Newman-Keuls test (SNK) using the Prisma Graphpad 5.0 software (GraphPad Software, San Diego, CA, USA). Data were considered statistically different when p < 0.05.
Young CW: Inbreeding and the gene pool. J Dairy Sci. 1984, 67: 472-477. 10.3168/jds.S0022-0302(84)81326-0.
Hermas SA, Young CW, Rust JW: Effects of mild inbreeding on productive and reproductive performance of Guernsey cattle. J Dairy Sci. 1987, 70: 712-715. 10.3168/jds.S0022-0302(87)80062-0.
Hudson GFS, Van Vleck LD: Inbreeding of artificially bred dairy cattle in the northeastern United States. J Dairy Sci. 1984, 67: 161-170. 10.3168/jds.S0022-0302(84)81280-1.
Falconer DS, Mackay TFC: Introduction to Quantitative Genetics. Harlow: Longman Group Limited; 1996.
Curik I, Sölkner J, Stipic N: The influence of selection and epistasis on inbreeding depression estimates. J Anim Breed Genet. 2001, 118: 247-262. 10.1046/j.1439-0388.2001.00284.x.
Hansson B, Westerberg L: On the correlation between heterozygosity and fitness in natural populations. Mol Ecol. 2002, 11: 2467-2474. 10.1046/j.1365-294X.2002.01644.x.
Kaeuffer R, Réale D, Pontier D, Chapuis JL, Coltman DW: Local effects of inbreeding on embryo number and consequences for genetic diversity in Kerguelen mouflon. Biol Lett. 2008, 4: 504-507. 10.1098/rsbl.2008.0222.
Acevedo-Whitehouse K, Spraker TR, Lyons E, Melin SR, Gulland F, DeLong RL, Amos W: Contrasting effects of heterozygosity on survival and hookworm resistance in California sealion pups. Mol Ecol. 2006, 15: 1973-1982. 10.1111/j.1365-294X.2006.02903.x.
MacNeil MD, Urick JJ, Newman S, Knapp BW: Selection for postweaning growth in inbred Hereford cattle: The Fort Keogh, Montana line I example. J Anim Sci. 1992, 70: 723-733.
Weigel KA: Controlling inbreeding in modern breeding programs. J Dairy Sci. 2001, 84: 177-184. 10.3168/jds.S0022-0302(01)74467-0.
Coltman DW, Pilkington JG, Smith JA, Pemberton JM: Parasite-mediated selection against inbred Soay sheep in a free-living island population. Evolution. 1999, 53: 1259-1267. 10.2307/2640828.
Paneto JCC, Ferraz JBS, Balieiro JCC, Bittar JFF, Ferreira MBD, Leite MB, Merighe GKF, Meirelles FV: Bos indicus or Bos taurus mitochondrial DNA-comparison of productive and reproductive breeding values in a Guzerat dairy herd. Genet Mol Res. 2008, 7: 592-602. 10.4238/vol7-3gmr449.
Ribeiro SHA, Pereira JCC, Verneque RS, Silva MA, Bergmann JAG, Ledic IL, Morais OR: Efeitos da origem e da linhagem do DNA mitocondrial sobre características produtivas e reprodutivas de bovinos leiteiros da raça Gir. Arq Bras Med Vet Zootec. 2009, 61: 232-242. 10.1590/S0102-09352009000100033.
Smith LA, Cassell BG, Pearson RE: The effects of inbreeding on the lifetime performance of dairy cattle. J Dairy Sci. 1998, 81: 2729-2737. 10.3168/jds.S0022-0302(98)75830-8.
Ralls K, Brugger K, Ballou J: Inbreeding and juvenile mortality in small populations of ungulates. Science. 1979, 206: 1101-1103. 10.1126/science.493997.
Lyons EJ, Frodsham AJ, Zhang L, Hill AVS, Amos W: Consanguinity and susceptibility to infectious diseases in humans. Biol Lett. 2009, 5: 5745-5776.
Guzman E, Price S, Poulsom H, Hope J: Bovine γδ T cells: cells with multiple functions and important roles in immunity. Vet Immunol Immunopathol. 2012, 148: 161-167. 10.1016/j.vetimm.2011.03.013.
Jutila MA, Holderness J, Graff JC, Hedges JF: Antigen independent priming: a transitional response of bovine gamma delta T-cells to infection. Anim Health Res Rev. 2008, 9: 47-57. 10.1017/S1466252307001363.
Guzman E, Birch JR, Ellis SA: Cattle MIC is a ligand for the activating NK cell receptor NKG2D. Vet Immunol Immunopathol. 2010, 136: 227-234. 10.1016/j.vetimm.2010.03.012.
Endsley JJ, Furrer JL, Endsley MA, McIntosh MA, Maue AC, Waters WR, Lee DR, Estes DM: Characterization of bovine homologues of granulysin and NK-lysin. J Immunol. 2004, 173: 2607-2614.
Alvarez AJ, Endsley JJ, Werling D, Mark Estes D: WC1(+) gammadelta T cells indirectly regulate chemokine production during Mycobacterium bovis infection in SCID-bo mice. Transbound Emerg Dis. 2009, 56: 275-284. 10.1111/j.1865-1682.2009.01081.x.
Welsh MD, Kennedy HE, Smyth AJ, Girvin RM, Andersen P, Pollock JM: Responses of bovine WC1(+) gamma delta T cells to protein and non protein antigens of Mycobacterium bovis. Infect Immun. 2002, 70: 6114-6120. 10.1128/IAI.70.11.6114-6120.2002.
Takamatsu HH, Denyer MS, Wileman TE: A subpopulation of circulating porcine gamma delta T cells can act as professional antigen presenting cells. Vet Immunol Immunopathol. 2002, 87: 223-224. 10.1016/S0165-2427(02)00083-1.
Cheng L, Cui Y, Shao H, Han G, Zhu L, Huang Y, O’Brien RL, Born WK, Kaplan HJ, Sun D: Mouse gamma delta T cells are capable of expressing MHC class II molecules, and of functioning as antigen-presenting cells. J Neuroimmunol. 2008, 203: 3-11. 10.1016/j.jneuroim.2008.06.007.
Brandes M, Willimann K, Bioley G, Levy N, Eberl M, Luo M, Tampe R, Levy F, Romero P, Moser B: Cross-presenting human gammadelta T cells induce robust CD8+alphabeta T cell responses. Proc Natl Acad Sci U S A. 2009, 106: 2307-2312. 10.1073/pnas.0810059106.
Skyberg JA, Thornburg T, Rollins M, Huarte E, Jutila MA, Pascual DW: Murine and bovine γδ T cells enhance innate immunity against Brucella abortus infections. PLoS One. 2011, 6: e21978-10.1371/journal.pone.0021978.
Bock RE, de Vos AJ, Kingston TG, McLellan DJ: Effect of breed of cattle on innate resistance to infection with Babesia bovis, B. bigemina and Anaplasma marginale. Aust Vet J. 1997, 75: 337-340. 10.1111/j.1751-0813.1997.tb15706.x.
Naessens J, Teale AJ, Sileghem M: Identification of mechanisms of natural resistance to African trypanosomiasis in cattle. Vet Immunol Immunopathol. 2002, 87: 187-194. 10.1016/S0165-2427(02)00070-3.
Macedo AA, Marciano APV, Rocha LM, Alves-Júnior JRF, Faria AMC, Bittar JFF, Araújo MSS, Santos RL, Martins-Filho OA: Comparative phenotypic profile of subpopulations of peripheral blood leukocytes in European (Bos taurus taurus) and Zebu cattle (Bos taurus indicus). Genet Mol Res. 2013, 12: 6838-6849. 10.4238/2013.December.19.2.
Macedo AA, Costa EA, Silva APC, Paixão TA, Santos RL: Monocyte-derived macrophages from Zebu (Bos taurus indicus) are more efficient to control Brucella abortus intracellular survival than macrophages from European cattle (Bos taurus taurus). Vet Immunol Immunopathol. 2013, 151: 294-302. 10.1016/j.vetimm.2012.12.001.
Martínez R, Dunner S, Toro R, Tobón J, Gallego J, Cañón J: Effect of polymorphisms in the Slc11a1 coding region on resistance to brucellosis by macrophages in vitro and after challenge in two Bos breeds (Blanco Orejinegro and Zebu). Genet Mol Biol. 2010, 33: 463-470. 10.1590/S1415-47572010000300014.
Greeley EH, Kealy RD, Ballam JM, Lawler DF, Segre M: The influence of age on the canine immune system. Vet Immunol Immunopathol. 1996, 5: 1-10.
Greeley EH, Ballam JM, Harrison JM, Kealy RD, Lawler DF, Segre M: The influence of age and gender on the immune system: a longitudinal study in Labrador Retriever dogs. Vet Immunol Immunopathol. 2001, 82: 57-71. 10.1016/S0165-2427(01)00336-1.
Naylor K, Li G, Vallejo AN, Lee WW, Koetz K, Bryl E, Witkowski J, Fulbright J, Weyand CM, Goronzy JJ: The influence of age on T cell generation. J Immunol. 2005, 174: 7446-7452.
May L, van den Biggelaar AH, van Bodegom D, Meij HJ, de Craen AJM, Amankwa J, Frölich M, Kuningas M, Westendorp RGJ: Adverse environmental conditions influence age-related innate immune responsiveness. Immun Ageing. 2009, 30: 6-7.
Paganelli R, Quinti I, Fagialo U, Cassarizza A, Ortolani C, Guerra E, Sansoni P, Pucillo LP, Scala E, Cozzi E, Bertollo L, Monti D, Franceschi C: Changes in circulating B cells and Ig classes and subclasses in a healthy aged population. Clin Exp Immunol. 1992, 90: 351-354.
Witkowski JM, Miller RA: Increased function of P-glycoprotein in T lymphocyte subsets of aging mice. J Immunol. 1993, 150: 1296-1306.
Boldman KG, Kriese LA, VanVleck LD, Van Tassell CP, Kachman SD: A Manual for Use of MTDFREML. A Set of Programs To Obtain Estimates of Variances and Covariances [DRAFT]. Lincoln: U.S. Department of Agriculture, Agricultural Research Service; 1995.
Work in OAMF and RLS labs is supported by CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazil), FAPEMIG (Fundação de Amparo a Pesquisa do Estado de Minas Gerais, Brazil) and PAPE-UNIUBE (Programa de Apoio a Pesquisa da Universidade de Uberaba, Minas Gerais, Brazil). OAMF and RLS are thankful to the CNPq for the PQ fellowship program. The authors thank the Program for Technological Development in Tools for Health - PDTIS - FIOCRUZ for use of its facilities.
The authors confirm that they have no conflicts of interest in this work.
AAM and PBB designed and prepared the manuscript. JFFB, JBR, ERB and JCCP performed the experiments. MSSA, RLS and OAMF supervised the study. All authors have read and approved the final manuscript.
Auricélio A Macedo, Joely F F Bittar contributed equally to this work.