Multivariate Analysis for Animal Selection in Experimental Research

Background: Several researchers seek methods for the selection of homogeneous groups of animals in experimental studies, a fact justified because homogeneity is an indispensable prerequisite for casualization of treatments. The lack of robust methods that comply with statistical and biological principles is the reason why researchers use empirical or subjective methods, influencing their results. Objective: To develop a multivariate statistical model for the selection of a homogeneous group of animals for experimental research and to elaborate a computational package to use it. Methods: The set of echocardiographic data of 115 male Wistar rats with supravalvular aortic stenosis (AoS) was used as an example of model development. Initially, the data were standardized, and became dimensionless. Then, the variance matrix of the set was submitted to principal components analysis (PCA), aiming at reducing the parametric space and at retaining the relevant variability. That technique established a new Cartesian system into which the animals were allocated, and finally the confidence region (ellipsoid) was built for the profile of the animals’ homogeneous responses. The animals located inside the ellipsoid were considered as belonging to the homogeneous batch; those outside the ellipsoid were considered spurious. Results: The PCA established eight descriptive axes that represented the accumulated variance of the data set in 88.71%. The allocation of the animals in the new system and the construction of the confidence region revealed six spurious animals as compared to the homogeneous batch of 109 animals. Conclusion: The biometric criterion presented proved to be effective, because it considers the animal as a whole, analyzing jointly all parameters measured, in addition to having a small discard rate.

Multivariate analysis for animal selection in experimental research

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Characterization of mitochondrial genotypes in the foundation herd of the Canchim beef cattle breed

The Canchim (5/8 Charolais + 3/8 Zebu) beef cattle breed was developed at Southeast-Embrapa Cattle to take advantage of hybrid vigor and to combine the higher growth rate and beef quality of Charolais with tropical adaptations of Zebu. The development of three lineages (old, new, and crossbred) has increased its genetic basis. The genotypic origin (Bos taurus or Bos indicus) of the mitochondrial DNA (mtDNA) of the Canchim breed was unknown. We characterized the mtDNA genotype of this founder herd by allele-specific polymerase chain reaction. The 173 founder Zebu females (62 Indubrasil, 3 Guzerat, and 108 Nellore) and their 6749 offspring were identified. The frequency of B. indicus mtDNA ranged from 1.15 to 2.05% among the descendants (N = 6404) of each maternal line with available DNA, and among animals that were alive (N = 689) in December 2007 among the three lineages. Though mtDNA characterization can be used to direct animal selection, the low frequency of B. indicus mtDNA impairs the evaluation of its effects on production traits in these animals. The high prevalence of B. taurus mtDNA in Canchim proves that the founder Zebu females from the Indubrasil, Guzerat and Nellore breeds were obtained from crosses of Zebu sires with local B. taurus dams.

Rats with high left ventricular end-diastolic pressure can be identified by Doppler echocardiography one week after myocardial infarction

The severity of left ventricular (LV) dysfunction in rats with myocardial infarction (MI) varies widely. Because homogeneity in baseline parameters is essential for experimental investigations, a study was conducted to establish whether Doppler echocardiography (DE) could accurately identify animals with high LV end-diastolic pressure as a marker of LV dysfunction soon after MI. Direct measurements of LV end-diastolic pressure were made and DE was performed simultaneously 1 week after surgically induced MI (N = 16) or sham-operation (N = 17) in female Wistar rats (200 to 250 g). The ratio of peak early (E) to late (A) diastolic LV filling velocities and the ratio of E velocity to peak early (Em) diastolic myocardial velocity were the best predictors of high LV end-diastolic pressure (>12 mmHg) soon after MI. Cut-off values of 1.77 for the E/A ratio (P = 0.001) identified rats with elevated LV end-diastolic pressure with 90% sensitivity and 80% specificity. Cut-off values of 20.4 for the E/Em ratio (P = 0.0001) identified rats with elevated LV end-diastolic pressure with 81.8% sensitivity and 80% specificity. Moreover, E/A and E/Em ratios were the only echocardiographic parameters independently associated with LV end-diastolic pressure in multiple linear regression analysis. Therefore, DE identifies rats with high LV end-diastolic pressure soon after MI. These findings have implications for using serial DE in animal selection and in the assessment of their response to experimental therapies.

Efeito de polimorfismos no gene LEP na expressão da leptina em adipócitos de bovinos de corte

Os caracteres produtivos são normalmente influenciados por muitos fatores, sendo difícil determinar todos os locos envolvidos em um fenótipo específico. Por isso, a seleção animal tem se baseado principalmente em uma estimativa direta ou indireta do fenótipo. A leptina é um importante regulador do metabolismo energético, da adiposidade e da reprodução. E por desempenhar diferentes funções, pode ser considerado um bom gene candidato para o estudo de associações entre marcadores moleculares e a eficiência reprodutiva ou ganho de peso. Em várias espécies, têm sido descritos diversos polimorfismos no gene da leptina, influenciando o ganho de peso, a reprodução, e outras características produtivas. Em bovinos, o STR IDVGA51 e o SNP LEPSau3A1, foram descritos por afetarem a performance reprodutiva, os alelos IDVGA51*181 e LEPSau3A1*2 estando associados a um aumento no intervalo entre partos de 79 e 81 dias, respectivamente, e os STRs BMS1074 e BM1500 afetam o ganho de peso, em vacas, no pós-parto: os alelos BMS1074*151 e BM1500*135 reduzindo e aumentando, respectivamente, o ganho de peso diário. Para confirmar o efeito ou não destes alelos na expressão do gene da leptina, este trabalho comparou os níveis de mRNA de leptina em animais portadores e não portadores dos alelos IDVGA51*181, LEPSau3A1*2, BMS1074*151 e BM1500*135, com os objetivos de: 1. Verificar as distribuições genotípicas e alélicas dos marcadores IDVGA51, BMS1074, BM1500 e LEPSau3A1, em uma amostra de 137 bovinos da raça Brangus-Ibagé, provenientes da Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA-Pecuária Sul, Bagé-RS). 2. Identificar, no rebanho, animais portadores e não portadores dos alelos IDVGA51*181, LEPSau3A1*2, BMS1074*151 ou BM1500*135, previamente descritos como associados à eficiência reprodutiva e ganho de peso. 3. Avaliar os animais, criados em campo nativo, com pastagem natural, quanto à sua condição corporal. 4. Realizar procedimento cirúrgico, para obtenção de amostras de tecido adiposo subcutâneo e omental. 5. Dosar os níveis de mRNA de leptina nos adipócitos destes dois tecidos, através do método quantitativo em tempo real (Real-Time RT-PCR), comparando a expressão do gene LEP entre indivíduos portadores e não portadores dos alelos de predisposição a ganho de peso e desempenho reprodutivo.

Correlações genéticas entre escores visuais e características reprodutivas em bovinos Nelore usando inferência bayesiana

O objetivo deste trabalho foi estimar a herdabilidade e as correlações genéticas entre escores visuais e características reprodutivas de animais da raça Nelore. As características avaliadas foram: precocidade, musculatura, e escores de conformação à desmama (PD, MD e CD, respectivamente) e ao sobreano (PS, MS e CS, respectivamente); idade ao primeiro parto (IPP); e perímetro escrotal (PE). Foram utilizadas informações de 66.244 animais, nascidos entre 1990 e 2006. Os parâmetros genéticos foram estimados em análises bicaracterísticas, com inferência bayesiana. Foi utilizado um modelo linear para IPP e PE, e um modelo não linear (threshold) para os escores visuais. As herdabilidades estimadas foram: CD, 0,19±0,02; PD, 0,23±0,02; MD, 0,20±0,02; CS, 0,26±0,01; PS, 0,33±0,02; MS, 0,32±0,02; IPP, 0,16±0,03; e PE, 0,36±0,02. As correlações genéticas estimadas entre os escores visuais e IPP foram negativas, de -0,18±0,03 a -0,29±0,02. Correlações genéticas positivas foram obtidas entre os escores visuais e o PE, de 0,19±0,01 a 0,31±0,01. A seleção de animais com os maiores escores visuais, principalmente ao sobreano, permite melhorar o desempenho reprodutivo dos rebanhos

Mecanismos celulares e moleculares que controlam o desenvolvimento e o crescimento muscular

O músculo estriado esquelético é formado pela associação de fibras musculares com a matriz extracelular. Esse tecido possui alta plasticidade e o conhecimento das características morfológicas, da miogênese, e da dinâmica do crescimento é importante para o entendimento da morfofisiologia bem como para a seleção de animais visando a melhoria na produção de carne. A maioria dos músculos estriados originam-se de células precursoras do mesoderma a partir dos somitos do embrião e o controle da diferenciação ocorre pela ação de fatores indutores ou inibidores. Um grupo de fatores transcricionais, pertencentes à família MyoD tem um papel central na diferenciação muscular. Coletivamente chamados de Fatores de Regulação Miogênica (MRFs), são conhecidos quatro tipos: MyoD, myf-5, miogenina e MRF4. Esses fatores ligam-se à seqüências de DNA conhecidas como Ebox (CANNTG) na região promotora de vários genes músculo-específicos, levando à expressão dos mesmos. As células embrionárias com potencial para diferenciação em células musculares (células precursoras miogênicas) expressam MyoD e Myf-5 e são denominadas de mioblastos. Essas células proliferam, saem do ciclo celular, expressam miogenina e MRF4, que regulam a fusão e a diferenciação da fibra muscular. Uma população de mioblastos que se diferencia mais tardiamente, as células miossatélites, são responsáveis pelo crescimento muscular no período pós natal, que pode ocorrer por hiperplasia e hipertrofia das fibras. As células satélites quiescentes não expressam os MRFs, porém, sob a ação de estímulos como fatores de crescimento ou citocinas, ocorre a ativação desse tipo celular que prolifera e expressa os MRFs de maneira similar ao que ocorre com as células precursoras miogênicas durante a miogênese. Os mecanismos de crescimento muscular são regulados pela expressão temporal dos (MRFs), que controlam a expressão dos genes relacionados com o crescimento muscular.

Degradabilidade in situ do Capim Panicum maximum cv. Tanzânia Incubado Cortado ou na Forma de Extrusa

The aim of this study was to compare the chemical composition, in vitro dry matter disappearance (IVDMD) and ruminal degradation of Panicum maximum, J. cv. Tanzania samples obtained by clipping (square method) or extrusa collection (animal selection). In the in situ trial, three ruminal fistulated dry crossbred cows, with 499 kg LW, were used in a completely randomized block design with split-plot arrangement design. Five grams of clipped (+/- 2 cm) grass or extrusa samples were placed in nylon bags (7 x 14 cm) and rumen incubated during 3, 6, 12, 24, 48, 96 and 120 hours. The IVDMD and the CP, NDF and ADF content were, respectively, 55.8, 7.6, 81.9 e 43.6%, for the clipped grass and .66.5, 12.1, 78.8 e 39.5%, for the extrusa samples. The potential degradability of DM, C P, NDF and ADF were 62.59, 80.88, 50.73 and 46.65%, for clipped grass; and 79.53, 90.97, 71.21 and 65.68%, for extrusa samples. The quality of the selected animal diet (extrusa) was better than the available forage in terms of IVDMD and chemical composition (high protein and low fiber content). In situ degradability trials carried out with clipped samples, and non selected by animal, could not supply reliable results closed to the animal diet.

Association of the Leptin Gene with Carcass Characteristics in Nellore Cattle

Advances in DNA technology have created biotechnological tools that can be used in animal selection and new strategies for increasing herd productivity and quality. The objective of the present work was to associate the genotypes of leptin gene exon 2 polymorphisms with productive traits in Nellore cattle. Blood was collected from Nellore males and PCR-RFLP reactions were performed with the restriction enzymes ClaI and Kpn2I. The gene frequencies resulting from digestion by ClaI were 0.60 and 0.40 for allele A and T, respectively; the genotypic frequencies were AA = 0.20 and AT = 0.80. The gene frequencies from digestion by Kpn2I were 0.81 for allele C and 0.194 for allele T; the genotypic frequencies were CC = 0.62 and CT = 0.38. The populations in both cases were not in Hardy-Weinberg equilibrium (p > 0.05), and the TT genotype was not found. Significant associations were noted between leptin gene exon 2 polymorphisms and five productive traits in Nellore cattle: carcass fat distribution, the intensity of red muscle coloration, pH, marbling, and post-slaughter fat thickness. © 2013 Copyright Taylor and Francis Group, LLC.

Estimativa de componentes de (co)variância de características de crescimento na raça Brahman utilizando inferência bayesiana

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Componentes principais na elaboração de índices de provas de desempenho em bovinos

The principal component analysis assists the producers in making decision of which evaluated features must be maintained in performance tests indexes, according to the variation present in these animals evaluated. The objective in this study was to evaluate a set of characteristics measured in a performance test in semifeedlot cattle of the Simmental and Angus breeds, by means principal component analysis (PC), aim to identify the features that represent most of the phenotypic variation for preparation of indexes. It was used data from 39 Angus and 38 Simmental bulls from the Santa Éster farm, located in Silvianópolis - MG. The performance test period was from october 2014 to february 2015. The features evaluated in the test were: final weight (FW), average daily gain weight (GW), respiratory rate (RR), haircoat temperature (HT) and rectal (RT), hair number (HN), hair length (HL), hair thickness (HT), muscularity (MUSC), racial characteristics, angulation, reproductive and balance (BAL), height of the front and back, width and length of croup, body length, depth and heart girth, subcutaneous fat thickness and rump (FTR), loin eye area and marbling (MAR). It was used PRINCOMP from SAS program for procedure the PC analysis. It was found that of the 27 features evaluated, the first four PC for Simmental breed explained 74% total variation data. The four PC selected with the corresponding weighting coefficients formed the following index: (0.27 * FW) + (0.47 * MUSC) + (0.50 * HL) + (0.39 * HT). Since the characteristics related to the adaptability of great importance for the studied breed, it was decided to keep the index of evidence for the Angus breed, the feature hair number, because there is a feature that presented a great variability and occupied one of the first principal component. Thus, the Angus index was composed by five features, with 79% total variation data, resulting in the following formula: (0.26 * FW) + (0.33 * BAL) + (0.58 * MAR) - (0.43 * FTR) – (0.38 * HN). By the principal component analysis it was possible to minimize the features number to be evaluated on performance tests from that farm, making the animal selection rapidly and accurate.

Estudo genômico da produção de leite e seus constituintes em bovinos da raça Guzerá

Pós-graduação em Genética e Melhoramento Animal - FCAV

Navigating uncertain waters: a critical review of inferring foraging behaviour from location and dive data in pinnipeds

In the last thirty years, the emergence and progression of biologging technology has led to great advances in marine predator ecology. Large databases of location and dive observations from biologging devices have been compiled for an increasing number of diving predator species (such as pinnipeds, sea turtles, seabirds and cetaceans), enabling complex questions about animal activity budgets and habitat use to be addressed. Central to answering these questions is our ability to correctly identify and quantify the frequency of essential behaviours, such as foraging. Despite technological advances that have increased the quality and resolution of location and dive data, accurately interpreting behaviour from such data remains a challenge, and analytical methods are only beginning to unlock the full potential of existing datasets. This review evaluates both traditional and emerging methods and presents a starting platform of options for future studies of marine predator foraging ecology, particularly from location and two-dimensional (time-depth) dive data. We outline the different devices and data types available, discuss the limitations and advantages of commonly-used analytical techniques, and highlight key areas for future research. We focus our review on pinnipeds - one of the most studied taxa of marine predators - but offer insights that will be applicable to other air-breathing marine predator tracking studies. We highlight that traditionally-used methods for inferring foraging from location and dive data, such as first-passage time and dive shape analysis, have important caveats and limitations depending on the nature of the data and the research question. We suggest that more holistic statistical techniques, such as state-space models, which can synthesise multiple track, dive and environmental metrics whilst simultaneously accounting for measurement error, offer more robust alternatives. Finally, we identify a need for more research to elucidate the role of physical oceanography, device effects, study animal selection, and developmental stages in predator behaviour and data interpretation.

Bioinformatic study of selection in animal genomes

Les gènes orthologues divergent sur plusieurs aspects durant l'évolution. Après une revue de la littérature cherchant à montrer de la divergence entre les orthologues de l'humain et de la souris, j'ai souligné les différentes causes de cette divergence. En comparant les gènes qui divergent en fonction, je n'ai pas trouvé de lien avec la divergence des séquences, pour cette raison je me suis penché sur l'étude de l'expression. Notamment, j'ai étudié le niveau, la spécificité ainsi que la présence/absence d'expression des orthologues humain-souris liés aux maladies Mendéliennes. Malgré les similarités trouvées entre l'humain et la souris, j'ai détecté une différence d'expression spécifique à une des deux espèces liée a un phénotype précis (gène essentiel/non-essentiel). Cela m'a permis de conclure que la différence sur le plan phénotypique entre l'humain et la souris est mieux expliquée par les patrons d'expression plutôt que le niveau d'expression ou la sélection. J'ai été également intéressé par l'évolution des séquences d'ADN codantes pour des protéines, en particulier sur le rôle de la sélection. J'ai commencé par une étude sur la fiabilité de détection de la sélection positive en comparant des séquences divergentes. J'ai trouvé, en utilisant le model de branche-site que la sélection peut être détectée sur des séquences qui ont divergé il y a plus de 500 millions d'années. J'ai analysé le biais de GC entres les séquences sans trouver une influence sur l'estimation de la sélection positive. Finalement, Je crois que ce travail est une première étape dans l'établissement d'un lien entre la sélection et les patrons d'expression des gènes chez les vertébrés.

Strategies to improve the efficiency of genomic selection in animal breeding programs

Pós-graduação em Zootecnia - FCAV