Novel phenotypes of feeding and social behaviour and their relationship with individual rabbit growth and feed efficiency

Due to the lack of a recording system for individual consumption of group-housed rabbits, published studies about feeding behaviour are based on information recorded at the group- and not at the individual level and periods covering only a few days or, in some cases, only part of a day. Such information could be used to inform rabbit management systems but cannot be used for genetic selection. We aimed to generate and use information from a novel automated feeder for group-housed rabbits to identify new phenotypes for individual animals that could be incorporated into breeding programs to improve feed efficiency and social behaviour under different feeding regimens. At 39 d of age, rabbits from 15 batches were placed in cages and fed ad libitum to become used to the electronic feeder. From 42 to 58-59 d, one group of 1 086 rabbits was fed ad libitum (AL), while another group of 1 134 rabbits was fed on a restricted feeding schedule (R) by limiting the feeding time to the period between 1800 and 0600 h of the following day. We implemented a reliable multivariate method to remove anomalous feeding behaviour records. We then defined novel traits for feeding behaviour that apply to both types of feeding regimes, and for social behaviour that indicates an animal's rank within the cage hierarchy. We based these traits on feeder records and a biologically sound definition of a meal. Finally, we estimated the phenotypic correlations of those traits with growth and feed efficiency traits. Our findings demonstrate that variables about resource distribution among cage mates and an animal's priority for feed access were found to be good indicators of an animal's dominant or subordinate status within the cage. Based on results obtained in R animals (results were similar in AL animals), the most efficient animals were those that ate less frequently (phenotypic correlation with feed conversion ratio, rho = 0.6), and consumed smaller amounts per meal (rho = 0.7), spent less time at the feeder (rho = 0.4), and appeared to be subordinate, as they did not have priority access to the feeder (rho = -0.3), and had the smallest share of resources (range of rho = 0.2-0.6). We conclude that quantifying feeding and social behaviour traits can enhance the understanding of the mechanisms through which individuals exert their effects on the performance of their cage mates.

Survival analysis of productive life in Florida dairy goats using a Cox proportional hazards model

Longevity is an economically important trait, since extending the functional life of a doe would allow us to keep the most productive females in the herd as long as possible, and this could result in the increased profitability of dairy farms. Thus, the objectives of this study were to determine the most important factors that influence the length of productive life (LPL) of female Florida goats and to estimate its genetic additive variance using a Cox proportional hazards model. The data consisted of 70,695 productive life records from 25,722 Florida females kidding between 2006 and 2020. A total of 19,495 does had completed their productive life while 6227 (24.2%) does had censored information. The pedigree contained information on 56,901 animals. The average censoring age and average failure age after first kidding for LPL were 36 and 47 months respectively. The model included, as time-independent effects, the age at first kidding and the interaction between herd, year and season of birth of the doe, and as time-dependent effects, the age at kidding, the interaction between herd, year and season of kidding, the within-herd class of milk production deviation, and the interaction between the lactation number and the stage of lactation. All fixed effects had a significant effect on LPL (p < 0.05). Does with older ages at the first kidding and an earlier age at kidding were at higher risk of being culled. A large difference among herds was observed in terms of culling risk, which highlighted the importance of adequate management practices. Also, high-producing does were less likely to be culled. The estimate of the additive genetic variance was 1.844 (in genetic standard deviation), with a heritability estimate of 0.58 ± 0.012. The results of this study are expected to contribute to the development of a genetic model for genetic evaluation of the length of the productive life of Spanish dairy goat breeds.

Molecular detection of Rickettsia in ectoparasites (Siphonaptera and Phthiraptera) of domestic and feral pigs from Argentina

Rickettsioses are distributed among a variety of hematophagous arthropods, and represent an emergent threat. The presence of rickettsial bacteria in ectoparasites collected from pigs from Argentina is still unknown. This study investigated the presence and identity of Rickettsia spp. in fleas, Pulex irritans, and sucking lice, Haematopinus suis, of domestic and feral pigs, Sus scrofa, from Central-Northern Argentina, through the genes gltA and ompB. Rickettsial bacteria were detected in 50% of fleas and 24% of lice. The BLASTn analysis of the ompB gene fragments in P. irritans samples showed identities 99% and 100% with R. felis. Positive samples of H. suis were 99% similar with species from the spotted fever group, future amplifications of a more polymorphic fragment of the ompB gene will allow to corroborate the identity of the Rickettsia species present in these lice samples. The Rickettsia spp. reported in the present study are having eventually been associated with cases of human diseases, and the circulation of these agents in arthropods has already been reported in several countries. Therefore, the identification of circulating pathogenic agents, such as reported in this study, is crucial for development of preventive measures for the control of ectoparasite-borne rickettsiosis diseases. Further studies, using serology techniques, will be allow to explore the ability of pigs as a possible Rickettsia reservoir and its role as part of transmission cycle of Rickettsia spp. in the studied scenarios.

Mutations in HID1 Cause Syndromic Infantile Encephalopathy and Hypopituitarism

OBJECTIVE

Precursors of peptide hormones undergo posttranslational modifications within the trans-Golgi network (TGN). Dysfunction of proteins involved at different steps of this process cause several complex syndromes affecting the central nervous system (CNS). We aimed to clarify the genetic cause in a group of patients characterized by hypopituitarism in combination with brain atrophy, thin corpus callosum, severe developmental delay, visual impairment, and epilepsy.

METHODS

Whole exome sequencing was performed in seven individuals of six unrelated families with these features. Postmortem histopathological and HID1 expression analysis of brain tissue and pituitary gland were conducted in one patient. Functional consequences of the homozygous HID1 variant p.R433W were investigated by Seahorse XF Assay in fibroblasts of two patients.

RESULTS

Bi-allelic variants in the gene HID1 domain-containing protein 1 (HID1) were identified in all patients. Postmortem examination confirmed cerebral atrophy with enlarged lateral ventricles. Markedly reduced expression of pituitary hormones was found in pituitary gland tissue. Colocalization of HID1 protein with the TGN was not altered in fibroblasts of patients compared to controls, while the extracellular acidification rate upon stimulation with potassium chloride was significantly reduced in patient fibroblasts compared to controls.

INTERPRETATION

Our findings indicate that mutations in HID1 cause an early infantile encephalopathy with hypopituitarism as the leading presentation, and expand the list of syndromic CNS diseases caused by interference of TGN function. ANN NEUROL 2021;90:149-164.

Anticipating the Impact of the USMLE Step 1 Pass/Fail Scoring Decision on Underrepresented-in-Medicine Students

Three-digit United States Medical Licensing Examination (USMLE) Step 1 scores have assumed an outsized role in residency selection decisions, creating intense pressure for medical students to obtain a high score on this exam. In February 2020, the Federation of State Medical Boards and the National Board of Medical Examiners announced that Step 1 would transition to pass/fail scoring beginning in 2022.The authors discuss the potential advantages and disadvantages of the pass/fail scoring change for underrepresented-in-medicine (UiM) trainees. UiM students may benefit from this change because it reduces the effect of an inequitable exam; helps correct for students who attend medical schools with a curriculum heavier on nontested formative elements; and decreases stress, improves quality of life, and undermines imposter syndrome. However, this change may also precipitate unforeseen challenges, such as increased discrimination toward UiM trainees, an increase in high-stakes test failures due to a reduced focus on preparing for standardized exams, or the development of new (e.g., subject exams) or overreliance on existing (e.g., school ranking) metrics that would make UiM residency candidates less competitive.To enhance UiM representation in the future health care workforce, it is imperative that national organizations (e.g., accrediting, licensing, regulatory, professional, honor, student, and faculty), hospitals, residency programs, and patient advocacy groups undertake a shared, rigorous approach in assessing the impact of the pass/fail scoring change on UiM applicants' selection to specialty and subspecialty residencies.

Genome-wide association study for feed efficiency in collective cage-raised rabbits under full and restricted feeding

Feed efficiency (FE) is one of the most economically and environmentally relevant traits in the animal production sector. The objective of this study was to gain knowledge about the genetic control of FE in rabbits. To this end, GWASs were conducted for individual growth under two feeding regimes (full feeding and restricted) and FE traits collected from cage groups, using 114 604 autosome SNPs segregating in 438 rabbits. Two different models were implemented: (1) an animal model with a linear regression on each SNP allele for growth trait; and (2) a two‐trait animal model, jointly fitting the performance trait and each SNP allele content, for FE traits. This last modeling strategy is a new tool applied to GWAS and allows information to be considered from non‐genotyped individuals whose contribution is relevant in the group average traits. A total of 189 SNPs in 17 chromosomal regions were declared to be significantly associated with any of the five analyzed traits at a chromosome‐wide level. In 12 of these regions, 20 candidate genes were proposed to explain the variation of the analyzed traits, including genes such as FTO, NDUFAF6 and CEBPA previously associated with growth and FE traits in monogastric species. Candidate genes associated with behavioral patterns were also identified. Overall, our results can be considered as the foundation for future functional research to unravel the actual causal mutations regulating growth and FE in rabbits.

Adherence and long-term outcomes of growth hormone therapy with easypod™ in pediatric subjects: Spanish ECOS study

BACKGROUND

Non-adherence to r-hGH treatments occurs in a variable percentage of subjects. One problem found when evaluating adherence is the great variability in methods of detection and definitions utilized in studies. This study assessed the level of adherence in subjects receiving r-hGH with the easypod™ electronic device.

METHODS

National, multicenter, prospective and observational study involving 238 subjects (144 with GH deficiency (GHD), and 86 with small for gestational age (SGA), 8 with Turner Syndrome), who received r-hGH with easypod™ for at least 3 months before inclusion. The follow-up period was 4 years.

RESULTS

Overall adherence was 94.5%; 97.5% after 6 months, 95.3% after 1 year, 93.7% after 2, 94.4% after 3 and 95.5% after 4 years of treatment. No differences in adherence were observed between prepubertal and pubertal groups and GHD and SGA groups. Change in height after 1 and 2 years, change in height SDS after 1 and 2 years, HV after 1 year, HV SDS after at 1 and 4 years, change in BMI after 1 year and change in BMI SDS at 1 and 2 years showed significant correlation with adherence. No significant differences in adherence according to IGF-I levels were found in follow-up visits or between groups.

CONCLUSIONS

The easypod™ electronic device, apart from being a precise and objective measure of adherence to r-hGH treatment, allows high compliance rates to be achieved over long periods of time. Adherence significantly impacts growth outcomes associated with r-hGH treatment.

Safety Outcomes During Pediatric GH Therapy: Final Results From the Prospective GeNeSIS Observational Program

Context

Safety concerns have been raised regarding premature mortality, diabetes, neoplasia, and cerebrovascular disease in association with GH therapy.

Objective

To assess incidence of key safety outcomes.

Design

Prospective, multinational, observational study (1999 to 2015).

Setting

A total of 22,311 GH-treated children from 827 investigative sites in 30 countries.

Patients

Children with growth disorders.

Interventions

GH treatment.

Main outcome measures

Standardized mortality ratio (SMR) and standardized incidence ratio (SIR) with 95% CIs for mortality, diabetes, and primary cancer using general population registries.

Results

Predominant short stature diagnoses were GH deficiency (63%), idiopathic short stature (13%), and Turner syndrome (8%), with mean ± SD follow-up of 4.2 ± 3.2 years (∼92,000 person-years [PY]). Forty-two deaths occurred in patients with follow-up, with an SMR (95% CI) of 0.61 (0.44, 0.82); the SMR was elevated for patients with cancer-related organic GH deficiency [5.87 (3.21, 9.85)]. Based on 18 cases, type 2 diabetes mellitus (T2DM) risk was elevated [SIR: 3.77 (2.24, 5.96)], but 72% had risk factors. In patients without cancer history, 14 primary cancers were observed [SIR: 0.71 (0.39, 1.20)]. Second neoplasms occurred in 31 of 622 cancer survivors [5.0%; 10.7 (7.5, 15.2) cases/1000 PY] and intracranial tumor recurrences in 67 of 823 tumor survivors [8.1%; 16.9 (13.3, 21.5) cases/1000 PY]. All three hemorrhagic stroke cases had risk factors.

Conclusions

GeNeSIS (Genetics and Neuroendocrinology of Short Stature International Study) data support the favorable safety profile of pediatric GH treatment. Overall risk of death or primary cancer was not elevated in GH-treated children, and no hemorrhagic strokes occurred in patients without risk factors. T2DM incidence was elevated compared with the general population, but most cases had diabetes risk factors.

Pedigree-based estimation of covariance between dominance deviations and additive genetic effects in closed rabbit lines considering inbreeding and using a computationally simpler equivalent model

Inbreeding generates covariances between additive and dominance effects (breeding values and dominance deviations). In this work, we developed and applied models for estimation of dominance and additive genetic variances and their covariance, a model that we call "full dominance," from pedigree and phenotypic data. Estimates with this model such as presented here are very scarce both in livestock and in wild genetics. First, we estimated pedigree-based condensed probabilities of identity using recursion. Second, we developed an equivalent linear model in which variance components can be estimated using closed-form algorithms such as REML or Gibbs sampling and existing software. Third, we present a new method to refer the estimated variance components to meaningful parameters in a particular population, i.e., final partially inbred generations as opposed to outbred base populations. We applied these developments to three closed rabbit lines (A, V and H) selected for number of weaned at the Polytechnic University of Valencia. Pedigree and phenotypes are complete and span 43, 39 and 14 generations, respectively. Estimates of broad-sense heritability are 0.07, 0.07 and 0.05 at the base versus 0.07, 0.07 and 0.09 in the final generations. Narrow-sense heritability estimates are 0.06, 0.06 and 0.02 at the base versus 0.04, 0.04 and 0.01 at the final generations. There is also a reduction in the genotypic variance due to the negative additive-dominance correlation. Thus, the contribution of dominance variation is fairly large and increases with inbreeding and (over)compensates for the loss in additive variation. In addition, estimates of the additive-dominance correlation are -0.37, -0.31 and 0.00, in agreement with the few published estimates and theoretical considerations.

Role of inbreeding depression, non-inbred dominance deviations and random year-season effect in genetic trends for prolificacy in closed rabbit lines

In closed rabbit lines selected for prolificacy at the Polytechnic University of Valencia, genetic responses are predicted using BLUP. With a standard additive BLUP model and year-season (YS) effects fitted as fixed, genetic trends were overestimated compared to responses estimated using control populations obtained from frozen embryos. In these lines, there is a confounding between genetic trend, YS effects and inbreeding, and the role of dominance is uncertain. This is a common situation in data from reproductively closed selection lines. This paper fits different genetic evaluation models to data of these lines, aiming to identify the source of these biases: dominance, inbreeding depression and/or an ill-conditioned model due to the strong collinearity between YS, inbreeding and genetic trend. The study involved three maternal lines (A, V and H) and analysed two traits, total born (TB) and the number of kits at weaning (NW). Models fitting YS effect as fixed or random were implemented, in addition to additive genetic, permanent environment effects and non-inbred dominance deviations effects. When YS was fitted as a fixed effect, the genetic trends were overestimated compared to control populations, inbreeding had an apparent positive effect on litter size and the environmental trends were negative. When YS was fitted as random, the genetic trends were compatible with control populations results, inbreeding had a negative effect (lower prolificacy) and environmental trends were flat. The model fitting random YS, inbreeding and non-inbred dominance deviations yielded the following ratios of additive and dominance variances to total variance for NW: 0.06 and 0.01 for line A, 0.06 and 0.00 for line V and 0.01 and 0.08 for line H. Except for line H, dominance deviations seem to be of low relevance. When it is confounded with inbreeding as in these lines, fitting YS effect as random allows correct estimation of genetic trends.

Genetic analysis of meat quality traits in maternal lines of rabbit and their diallel cross

Young rabbits, the dams of which came from a full diallel cross among four maternal lines (A, V, H and LP) and the sires from a single paternal line (R), that produce sixteen genetic groups, was carried out to evaluate the genetic groups and to estimate the crossbreeding genetic parameters of meat quality. The meat quality traits were recorded by NIRS from a sample of 285 longissimus lumborum muscles. Crossbreeding parameters were estimated according to Dickerson model. No differences in protein were found. The line A had significant differences with V line for intramuscular fat, and fatty acids groups. Significant differences for these traits appeared between the crossbred AH and VV (in favor of AH). As conclusion, the significant contrasts between genetic types for chemical composition of the meat are mainly consequence of direct-maternal genetic effects, having grandmaternal and maternal heterosis effects a less relevant role.

Expected responses to different strategies of selection to increase heat tolerance assessed by changes in litter size in rabbit

Response on litter size (LS) to selection based on EBV of the own trait and several indexes including EBV of 1 or more parameters of a reaction norm model (RNM) was evaluated by simulation. The RNM described animal's performance under different environmental temperatures assuming that this has an animal-specific effect on LS beyond an animal-specific temperature threshold (i.e., it had 3 parameters: intercept [Int], slope [Slp], and threshold [Thr]). Litter size data were generated assuming this model under 2 environmental conditions (comfort [CC] and comfort-to-heat [CtHC]). Variance components for Int, Slp, Thr, and LS were taken from literature. The initial base population consisted of 125 females and 25 males. Ten generations of selection were conducted keeping constant the population size. Eight different selection criteria were considered, depending on both the evaluation model and the index combining EBV for 1 or more parameters of the RNM. In 1 case selection was based on EBV of the own trait predicted by using a repeatability animal model. In the other 7 cases the genetic evaluation was conducted using the same model as that used in the simulation. For each scenario 25 replicates were conducted and response to selection was assessed within replicate as the difference between generations in the average of trait. Results indicate that, under the studied conditions and for the used genetic parameters, selection based on the observed trait seems to have the same effect as selection based on some index, including EBV of the RNM parameters. In addition, response to selection could be greater under CtHC rather than CC. Animals selected exclusively for EBV of the Slp and Thr are not expected to have good performance under CC. Under CtHC, selection for LS has a major response on those parameters, whereas response on Int was very small. Under CC, response to selection on LS is mainly determined by a change in Int, whereas Slp does not change and Thr slightly increases. Selection based on EBV of Slp, Thr, or on an index including both seems to have the same effect on the trait and it was doubly effective in modifying the shape of the RNM under CtHC than under CC. Selection based on EBV of Thr and Slp does not seem to lead to any response in LS. Selection based on EBV of Int seems to have no effect on the trait or on animal's tolerance to heat under CtHC, but it would lead to a positive response in LS under CC.

Effective population size and inbreeding depression on litter size in rabbits. A case study

The purpose of this study is to use demographic and litter size data on four Spanish maternal lines of rabbits (A, V, H and LP), as a case study, in order to: (i) estimate the effective population size of the lines, as a measure of the rate of increase of inbreeding, and (ii) study whether the inbreeding effect on litter size traits depends on the pattern of its accumulation over time. The lines are being selected for litter size at weaning and are kept closed at the same selection nucleus under the same selection and management programme. The study considered 47,794 l and a pedigree of 14,622 animals. Some practices in mating and selection management allow an increase of the inbreeding coefficient lower than 0.01 per generation in these lines of around 25 males and 125 females. Their effective population size (Ne) was around 57.3, showing that the effect of selection, increasing the inbreeding, was counterbalanced by the management practices, intended to reduce the rate of inbreeding increase. The inbreeding of each individual was broken down into three components: old, intermediate and new inbreeding. The coefficients of regression of the old, intermediate and new inbreeding on total born (TB), number born alive (NBA) and number weaned (NW) per litter showed a decreasing trend from positive to negative values. Regression coefficients significantly different from zero were those for the old inbreeding on TB (6.79 ± 2.37) and NBA (5.92 ± 2.37). The contrast between the coefficients of regression between the old and new inbreeding were significant for the three litter size traits: 7.57 ± 1.72 for TB; 6.66 ± 1.73 for NBA and 5.13 ± 1.67 for NW. These results have been interpreted as the combined action of purging unfavourable genes and artificial selection favoured by the inbreeding throughout the generations of selection.

Racial and ethnic minority medical students' perceptions of and interest in careers in academic medicine

PURPOSE

To describe diverse medical students' perceptions of and interest in careers in academic medicine.

METHOD

In 2010, the authors invited students attending three national medical student conferences to respond to a survey and participate in six focus groups. The authors identified trends in data through bivariate analyses of the quantitative dataset and using a grounded theory approach in their analysis of focus group transcripts.

RESULTS

The 601 survey respondents represented 103 U.S. medical schools. The majority (72%) were in their first or second year; 34% were black and 17% were Hispanic. Many respondents (64%) expressed interest in careers in academic medicine; teaching and research were viewed as positive influences on that interest. However, black and Hispanic respondents felt they would have a harder time succeeding in academia. The 73 focus group participants (25% black, 29% Hispanic) described individual- and institutional-level challenges to academic medicine careers and offered recommendations. They desired deliberate and coordinated exposure to academic career paths, research training, clarification of the promotion process, mentorship, protected time for faculty to provide teaching and research training, and an enhanced infrastructure to support diversity and inclusion.

CONCLUSIONS

Medical students expressed an early interest in academic medicine but lacked clarity about the career path. Black and Hispanic students' perceptions of having greater difficulty succeeding in academia may be an obstacle to engaging them in the prospective pool of academicians. Strategic and dedicated institutional resources are needed to encourage racial and ethnic minority medical students to explore careers in academic medicine.

Comparison of immune response to lipopolysaccharide of rabbit does selected for litter size at weaning or founded for reproductive longevity

To evaluate differences in maternal lines to the immune response of reproductive rabbit does, a total of 64 animals of two different lines: (1) founded for hyper-longevity and litter size criteria (LP) and (2) selected for litter size at weaning (V) were used. Females were subjected to three different reproductive efforts: post-partum (PP) mating at first lactation and 9 kits during the second; post-weaning (PW) mating at first lactation and 9 kits during the second; and PW mating at first lactation and 5 kits during the second. At second weaning (30 days PP), an acute response was induced by intravenous infusion of lipopolysaccharide (LPS). LP females seemed to be lower affected during the hyper-acute phase than V females, showing lower plasma glucose content at 1.5 h post infusion (pi) and rectal temperature at 6 h pi; and showed higher ulterior immune response, with higher levels of C-reactive protein at 48 h pi and haptoglobin in plasma from 24 h pi. Survival test conferred a higher risk of culling for V than for LP females during the first hours after challenge. These results may suggest that, regarding immune response to LPS challenge, foundation by hyper-longevity productive criteria lead to obtain a more robust population of rabbit does, characterized by improved response ability.

Sources of individual variation to heat tolerance in a rabbit line

The objective of the present study was to quantify the extent of the interaction between temperature and animal-specific effects on prolificacy and average weaning weight in rabbit does and to assess to what extent this interaction has a genetic determination. For this purpose, 18,491 total born (TB) records from 5,929 Caldes does, as well as 16,868 weaned kits (NW) and average weaning weight (AvgWW) records from 5,353 does of the same line were used. Four models were considered, ranging from the simplest, which assumed no effect of temperature on the traits studied, to the most complex, which allowed for different slopes of response and thresholds for heat tolerance for each animal. Deviance information criteria (DIC) were adopted as model choice criteria. The model postulating that each animal is subjected to a linear change in its performance when the temperature reaches a certain animal-specific threshold was preferred for all the considered traits. This interaction has the greatest effect for AvgWW; the ratio between the variance associated to the intercept and the total variance declined from 29.1% at 19°C to 19.4% at 25°C. This decline results from the increase in the variation associated with the interaction of the individual effects with temperature and also from the covariance between this interaction term and the intercepts. An important part of this interaction has a genetic origin, particularly for AvgWW and NW, for which the estimated heritability increased from 0.11 to 0.24 and 0.06 to 0.10, respectively, from 19 to 25°C. Similarly, it was found that the proportion of the total variance for the individual effects involved in the interaction with temperature which have a genetic origin is high for all the traits, around one-half for the variation on the linear changes in the performance and more than one-third for the variation on the onsets for heat stress. However, the error associated with these estimates is large and results must be considered with caution. Thus, in conclusion the extent of the interaction between individual variation and in farm temperature seems to be important and has a genetic origin. Nevertheless, the accuracy of the estimated genetic parameters is low.

Health and body condition of lactating females on rabbit farms

The aim was to study relationships between morbidity, defined through prevalence of diseases of rabbit females, and BCS, and to assess the effect of several risk factors on both morbidity and BCS. Our study was based on individual examinations of 18,510 does in lactation on 103 farms. We evaluated BCS on a linear scale from 1 to 9, with 5 being the optimum. Prevalence of major diseases were 22.7%, 4.0%, and 6.4% for coryza, mastitis, and ulcerative pododermatitis, respectively. In addition, prevalence was 3.0% for diseases of minor presentation, including mange, which had a prevalence of 1.9%. The BCS of the R line (selected for growth) was 5.55 ± 0.14, whereas for the A line (selected for litter size) it was 4.40 ± 0.11. Females with more than 20 kindlings had on average a BCS 0.25 ± 0.07 units less than those in the 12th lactation (P = 0.0002). Optimal BCS 4.60 ± 0.11 was reached during the third lactation week. Sick females had a BCS of 0.6 ± 0.01 units less than healthy females. Females with a footrest had on average a BCS 0.19 ± 0.05 units greater than those without. The absence of footrests was an enabling risk factor for ulcerative pododermatitis, the prevalence of which increased by 53%. Ulcerative pododermatitis was associated (P = 0.045) with diet; females consuming a rich energy diet were prone to having this disorder; 1 SD increase in DE (0.32 MJ) determined an increase in ulcerative pododermatitis prevalence of 0.8 percentage points. Diet was not an enabling risk factor for the other diseases. The genetic type to which a female belongs is a predisposing risk factor of disease; P, V and H were also maternal lines, while S group was exclusively formed by maternal lines. With regard to coryza, the S group had the greatest prevalence (44.0%), followed by A, P, R (19.0 to 21.0%); the V line, selected for prolificacy, showed the least prevalence (12.0%). For the case of mastitis although significant (P < 0.05), the magnitude of the differences between disease prevalence was less; R line had a mastitis prevalence of 11.0% while the least prevalence was observed for V does (4.0%). Simultaneous evaluation of both BCS and morbidity on the rabbit farm is recommended for the right assessment of welfare conditions. In this study, the relationships between both variables have been shown, as well as how other intrinsic and extrinsic risk factors modulate these variables; and thus, these factors should be considered during a welfare assessment.

Replication and refinement of a quantitative trait locus influencing milk protein percentage on ovine chromosome 3

A previous genome scan that was conducted in Spanish Churra sheep identified a significant quantitative trait locus (QTL) for milk protein percentage (PP) on chromosome 3 (OAR3), between markers KD103 and OARVH34. The aim of this study was to replicate these results and to refine the mapped position of this QTL. To accomplish this goal, we analysed 14 new half-sib families of Spanish Churra sheep including 1661 ewes from 29 different flocks. These animals were genotyped for 21 microsatellite markers mapping to OAR3. In addition to a classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) was performed with the aim of enhancing the resolution of the QTL mapping. The LA that was performed in this sheep population identified the presence of a highly significant QTL for PP near marker KD103 (P(c) < 0.001; P(exp) < 0.001). The phenotypic variance that was owing to the QTL was 2.74%. Two segregating families for the target QTL were identified in this population with QTL effect estimates of 0.47 and 0.95 SD. The LDLA identified the same QTL as the previous analyses with a high level of statistical significance (P = 9.184 E-11) and narrowed the confidence interval (CI) to a 13 cM region. These results confirm the segregation of the previously identified OAR3 QTL that influences PP in Spanish Churra sheep. Future research will aim to increase the marker density across the refined CI and to analyse the corresponding candidate genes to identify the allelic variant or variants that underlie this genetic effect.

The role of bovine causal genes underlying dairy traits in Spanish Churra sheep

In dairy cattle, quantitative trait nucleotides (QTNs) underlying quantitative trait loci (QTL) for milk production traits have been identified in bovine DGAT1, GHR and ABCG2 genes. The SPP1 gene has also been proposed to be a regulator of lactation. In sheep, QTL underlying milk production traits have been reported only recently, and no proven QTN has been identified. Taking into account the close phylogenetic relationship between sheep and cattle, this study examined the possible effects of the aforementioned genes on sheep milk production traits. We first studied the genetic variability of the DGAT1, GHR, ABCG2 and SPP1 genes in 15 rams of the Spanish Churra dairy sheep breed. Second, we performed an association analysis between SNPs identified in these genes and three milk production traits recorded in a commercial population of Churra sheep. This analysis revealed only three significant associations at the nominal level (P-value <0.05) involving allelic variants of the ABCG2 gene, whereas no significant association was found for the DGAT1, GHR and SPP1 genes. When the Bonferroni correction was applied to take into account the multiple tests performed, none of the associations identified at the nominal level remained significant. Nevertheless, taking into account the high level of false-negative findings that can arise when applying the stringent Bonferroni correction, we think that our results provide a valuable primary assessment of strong candidate genes for milk traits in sheep.

Detection of quantitative trait loci affecting the milk fatty acid profile on sheep chromosome 22: role of the stearoyl-CoA desaturase gene in Spanish Churra sheep

Sheep milk fat contains several components that may provide human health benefits, such as monounsaturated fatty acids and conjugated linoleic acid (CLA). Most of the CLA in ruminant milk is synthesized in the mammary gland by the action of the enzyme stearoyl-CoA desaturase (SCD) on circulating vaccenic acid (trans-11 C18:2; VA). Previous studies have found significant associations between polymorphisms in the SCD gene and the fatty acid composition of ruminant products, including sheep milk. Based on this, we performed a quantitative trait loci (QTL) analysis of an ovine chromosome (22) that harbors the SCD gene for effects on milk fatty acid composition traits and classical milk production traits. We identified a suggestive QTL influencing the CLA/VA ratio with the maximum statistic at position 26 cM of the studied chromosome, whereas the SCD gene has been mapped to position 41.6 cM. The individual introduction of 4 SCD single nucleotide polymorphisms in the QTL model did not cause a reduction of the variance explained by the QTL, which suggests that the SCD gene is not directly responsible for the detected effect in the Churra population studied herein. This conclusion was supported by the lack of any significant association identified between the 4 SCD single nucleotide polymorphisms and the CLA/VA ratio. This association analysis suggested a possible effect of the SCD gene on milk fat percentage in Churra sheep. An independent confirmation of these primary results will be required before attempting its practical implementation in selection programs.

The identification of QTL that affect the fatty acid composition of milk on sheep chromosome 11

In this work, we analysed 11 genetic markers localized on OAR11 in a commercial population of Spanish Churra sheep to detect QTL that underlie milk fatty acid (FA) composition traits. Following a daughter design, we analysed 799 ewes distributed in 15 half-sib families. Eight microsatellite markers and three novel SNPs identified in two genes related to fatty acid metabolism, acetyl-CoA carboxylase alpha (ACACA) and fatty acid synthase (FASN), were genotyped in the whole population under study. The phenotypic traits considered in the study included 22 measurements related to the FA composition of the milk and three other milk production traits (milk protein percentage, milk fat percentage and milk yield). Across-family regression analysis revealed four significant QTL at the 5% chromosome-wise level influencing contents of capric acid (C10:0), lauric acid (C12:0), linoleic conjugated acid (CLA) and polyunsaturated fatty acids (PUFA) respectively. The peaks of the QTL affecting C10:0 and PUFA contents in milk map close to the FASN gene, which has been evaluated as a putative positional candidate for these QTL. The QTL influencing C12:0 content reaches its maximum significance at 58 cM, close to the gene coding for the glucose-dependent insulinotropic polypeptide. We were not able to find any candidate genes related to fat metabolism at the QTL influencing CLA content, which is located at the proximal end of the chromosome. Further research efforts will be needed to confirm and refine the QTL locations reported here.

Genetic determination of the onset of heat stress on daily milk production in the US Holstein cattle

Existence of individual variation in the onset of heat stress for daily milk yield of dairy cows was assessed. Data included 353,376 test-day records of 38,383 first-parity Holsteins from a random sample of US herds. Three hierarchical models were investigated. Model 1 inferred the value of a temperature-humidity index (THI) at which mean yield began to decline as well as the extent of that decline. Model 2 assumed individual variation in yield decline beyond a common THI threshold. Model 3 additionally assumed individual variation for the onset of heat stress. Deviance information criteria indicated the superiority of model 3 over model 2. For model 2, genetic correlation between milk yield in the absence of heat stress and the THI threshold for heat stress was -0.4 (0.11) [marginal posterior mean (marginal posterior standard deviation)]. For model 3, genetic correlations were -0.53 (0.05) between milk yield and THI threshold and -0.62 (0.08) between milk yield and yield decay beyond the THI threshold. Total standard deviation (sum of additive genetic and permanent environmental standard deviations) for the THI threshold was 3.95 (0.06), and more than half of that variation had an additive genetic origin [56% (5%)]. Because of the high genetic correlation [0.95 (0.03)] between yield decay and THI threshold with model 3, using only one of them as a selection criterion for heat tolerance would modify the other in the desired direction.

Multi-breed genetic evaluation in a Gelbvieh population

A multi-breed model was presented for the genetic evaluation of growth traits in beef cattle. In addition to the fixed effects, random direct and maternal genetic effects, and random maternal permanent environmental effects are considered; the model also fits direct and maternal heterosis and direct and maternal breed-of-founder (BOF) x generation group effects using a Bayesian approach that weights prior literature estimates relative to information supplied by the dataset to which the model will be applied. The multi-breed evaluation procedures also allow the inclusion of external evaluations for animals of other breeds. The multi-breed model was applied to a dataset provided by the American Gelbvieh Association. Different analyses were conducted by varying the weights given to the prior literature relative to the information provided by the dataset. Large differences were observed for the heterosis estimates, the BOF x generation group effect estimates, and the predicted breeding values across breeds due to the weights posed on prior literature estimates versus estimates derived directly from data. However, the rankings within breed were observed to be relatively robust to the different weights on prior information.

Genetic and environmental correlations between longevity and litter size in rabbits

An estimation of the genetic and environmental correlations between prolificacy (kits born alive and number of kits at weaning) and functional longevity in a population selected for litter size at weaning since 1992 was carried out. The method is only an approximation, but suggests that genetically longevity and litter size are not antagonistic objectives in breeding programmes, because the estimated genetic correlations are not significantly different from zero.

Crossbreeding parameter estimation for functional longevity in rabbits using survival analysis methodology

A complete diallel cross involving 3 maternal lines of rabbit was performed to estimate cross-breeding parameters for functional longevity. This trait was defined as the ability to delay involuntary culling. The lines considered, A, V, and Prat, had all been selected by litter size at weaning for a long period. Data were related to a total of 653 does belonging to the 9 genetic types from the diallel cross; does were reared and bred on the same commercial farm. Survival analysis was performed using a Cox proportional hazard model. The model incorporated time-dependent factors, such as year-season, litter size, and the interaction between cycle and physiological status of the female; time-independent factors, such as the genetic type of the doe; and sire and dam random factors. Crossbreeding parameters were estimated from the solutions obtained for the type of doe and its estimated variance-covariance matrix, using a generalized least squares procedure. The estimated parameters were the differences between lines in direct genetic effects and maternal genetic effects and individual heterosis. Relevant differences were observed in direct genetic effects between lines A and Prat but not in any maternal effects. Heterosis was found to be significant and favorable between lines A and Prat, and between the lines V and Prat. The magnitude of this effect was variable but important, especially in the first cross. Results stress the importance of using crosses between specialized lines to produce does for intensive meat rabbit production.

Effect of changing female stimulus on intensive semen collection in young Murciano-Granadina male goats1

The aim of this research was to study the effect of changing female stimulus on libido and semen characteristics from young Murciano-granadina male goats submitted to intensive semen collection using females not in estrus as teasers. Males were submitted to two different sexual stimulation procedures. In the first procedure, the same doe was used as the female stimulus for three consecutive presentations. In the second, the doe was replaced after the second presentation by a new female. Semen volume, concentration, forward progressive motility, and live spermatozoa were scored. To analyze reaction time (RT), three types of analysis were performed. In the first one, RT was analyzed by multifactor ANOVA, taking as a missing value 300 s when a buck did not ejaculate. In the second, RT also was analyzed by multifactor ANOVA, but data from males that did not ejaculate were removed. In the third, a Cox Survival analysis was carried out by censoring data when a buck did not ejaculate within 5 min of entering the test arena. A decrease in semen volume and sperm concentration in the successive ejaculations was observed, being highly marked in the third ejaculation independent of the stimulation procedure (0.62 vs. 0.38 and 0.43 mL, and 2,828 vs. 2,183 and 2,223 million spermatozoa/mL to the first and third ejaculation respectively; P < 0.05). No significant differences were observed either in forward progressive motility or live sperm rate. Changing the female stimulus in the third presentation had no significant effect on any seminal characteristic. Regarding libido and mounting behavior variables, there was a substantial decrease in RT in the third service when the female was changed (with both types of ANOVA). When censored data were taken into account, the relative risk showed that the probability of a male ejaculating in the third presentation increased almost fourfold when the female was replaced than when the female was the same in all services (P < 0.05). In conclusion, young Murciano-granadina bucks can be used as semen donors because none of the most important semen variables used to reject or accept an ejaculate before freezing process decreased after intensive semen collection. We also recommend changing the female stimulus to make the semen collection procedure more efficient and using survival analysis methodology to analyze time data, mainly when a high rate of censored data are scored.

[Biotechnological aspects in "loco" larvae]

The biology of planktotrophic larvae of Concholepas concholepas is the main bottleneck towards developing biotechnologies to rear this muricid. Data concerning planktonic larvae development, diets and environmental signals triggering larval settlement and recruitment is scarce. We have begun the study of the molecular and cell biology of embryos, larvae and recruits having as a final goal, the development of appropriate biotechnologies to rear this gastropod. First, an inverse ratio between BuChE and AChE enzyme activities was established. This ratio may be a precise developmental marker for this species. Second, for the first time a phosphoinositide related regulatory pathway is reported in a muricid, opening a new approach to the biotechnological management of larvae. Third, the relation between sulfate in sea water and larval motility was studied. Concentrations below 125 microM sulfate decreases larval motility. The sulfate is incorporated in proteoglycans which participate in different developmental phenomena. Lastly, a genomic Concholepas concholepas DNA sequence, similar to that of a human growth hormone probe was detected. This is very interesting since growth factors are key molecules during development, growth and are involved in food conversion rates in fish and also, in a variety of marine invertebrates.