Rachel Rupp

Detection of genes influencing economic traits in three French dairy cattle breeds

A project of QTL detection was carried out in the French Holstein, Normande, and Montbéliarde dairy cattle breeds. This granddaughter design included 1 548 artificial insemination bulls distributed in 14 sire families and evaluated after a progeny-test for 24 traits (production, milk composition, persistency, type, fertility, mastitis resistance, and milking ease). These bulls were also genotyped for 169 genetic markers, mostly microsatellites. The QTL were analysed by within-sire linear regression of daughter yield deviations or deregressed proofs on the probability that the son receives one or the other paternal QTL allele, given the marker information. QTL were detected for all traits, including those with a low heritability. One hundred and twenty QTL with a chromosome-wise significance lower than 3% were tabulated. This threshold corresponded to a 15% false discovery rate. Amongst them, 32 were genome-wise significant. Estimates of their contribution to genetic variance ranged from 6 to 40%. Most substitution effects ranged from 0.6 to 1.0 genetic standard deviation. For a given QTL, only 1 to 5 families out of 14 were informative. The confidence intervals of the QTL locations were large and always greater than 20 cM. This experiment confirmed several already published QTL but most of them were original, particularly for non-production traits.

Differential response of bovine mammary epithelial cells to Staphylococcus aureus or Escherichia coli agonists of the innate immune system

Mastitis caused by Escherichia coli and Staphylococcus aureus is a major pathology of dairy cows. To better understand the differential response of the mammary gland to these two pathogens, we stimulated bovine mammary epithelial cells (bMEC) with either E. coli crude lipopolysaccharide (LPS) or with S. aureus culture supernatant (SaS) to compare the transcriptomic profiles of the initial bMEC response. By using HEK 293 reporter cells for pattern recognition receptors, the LPS preparation was found to stimulate TLR2 and TLR4 but not TLR5, Nod1 or Nod2, whereas SaS stimulated TLR2. Biochemical analysis revealed that lipoteichoic acid, protein A and α-hemolysin were all present in SaS, and bMEC were found to be responsive to each of these molecules. Transcriptome profiling revealed a core innate immune response partly shared by LPS and SaS. However, LPS induced expression of a significant higher number of genes and the fold changes were of greater magnitude than those induced by SaS. Microarray data analysis suggests that the activation pathways and the early chemokine and cytokine production preceded the defense and stress responses. A major differential response was the activation of the type I IFN pathway by LPS but not by SaS. The higher upregulation of chemokines (Cxcl10, Ccl2, Ccl5 and Ccl20) that target mononuclear leucocytes by LPS than by SaS is likely to be related to the differential activation of the type I IFN pathway, and could induce a different profile of the initial recruitment of leucocytes. The MEC responses to the two stimuli were different, as LPS was associated with NF-κB and Fas signaling pathways, whereas SaS was associated with AP-1 and IL-17A signaling pathways. It is noteworthy that at the protein level secretion of TNF-α and IL-1β was not induced by either stimulus. These results suggest that the response of MEC to diffusible stimuli from E. coli and S. aureus contributes to the onset of the response with differential leucocyte recruitment and distinct inflammatory and innate immune reactions of the mammary gland to infection.

Design and Characterization of a 52K SNP Chip for Goats

The success of Genome Wide Association Studies in the discovery of sequence variation linked to complex traits in humans has increased interest in high throughput SNP genotyping assays in livestock species. Primary goals are QTL detection and genomic selection. The purpose here was design of a 50–60,000 SNP chip for goats. The success of a moderate density SNP assay depends on reliable bioinformatic SNP detection procedures, the technological success rate of the SNP design, even spacing of SNPs on the genome and selection of Minor Allele Frequencies (MAF) suitable to use in diverse breeds. Through the federation of three SNP discovery projects consolidated as the International Goat Genome Consortium, we have identified approximately twelve million high quality SNP variants in the goat genome stored in a database together with their biological and technical characteristics. These SNPs were identified within and between six breeds (meat, milk and mixed): Alpine, Boer, Creole, Katjang, Saanen and Savanna, comprising a total of 97 animals. Whole genome and Reduced Representation Library sequences were aligned on >10 kb scaffolds of the de novo goat genome assembly. The 60,000 selected SNPs, evenly spaced on the goat genome, were submitted for oligo manufacturing (Illumina, Inc) and published in dbSNP along with flanking sequences and map position on goat assemblies (i.e. scaffolds and pseudo-chromosomes), sheep genome V2 and cattle UMD3.1 assembly. Ten breeds were then used to validate the SNP content and 52,295 loci could be successfully genotyped and used to generate a final cluster file. The combined strategy of using mainly whole genome Next Generation Sequencing and mapping on a contig genome assembly, complemented with Illumina design tools proved to be efficient in producing this GoatSNP50 chip. Advances in use of molecular markers are expected to accelerate goat genomic studies in coming years.

Genetic analysis for mastitis resistance and milk somatic cell score in French Lacaune dairy sheep

Genetic analysis for mastitis resistance was studied from two data sets. Firstly, risk factors for different mastitis traits, i.e. culling due to clinical or chronic mastitis and subclinical mastitis predicted from somatic cell count (SCC), were explored using data from 957 first lactation Lacaune ewes of an experimental INRA flock composed of two divergent lines for milk yield. Secondly, genetic parameters for SCC were estimated from 5 272 first lactation Lacaune ewes recorded among 38 flocks, using an animal model. In the experimental flock, the frequency of culling due to clinical mastitis (5%) was lower than that of subclinical mastitis (10%) predicted from SCC. Predicted subclinical mastitis was unfavourably associated with the milk yield level. Such an antagonism was not detected for clinical mastitis, which could result, to some extent, from its low frequency or from the limited amount of data. In practice, however, selection for mastitis resistance could be limited in a first approach to selection against subclinical mastitis using SCC. The heritability estimate of SCC was 0.15 for the lactation mean trait and varied from 0.04 to 0.12 from the first to the fifth test-day. The genetic correlation between lactation SCC and milk yield was slightly positive (0.15) but showed a strong evolution during lactation, i.e. from favourable (-0.48) to antagonistic (0.27). On a lactation basis, our results suggest that selection for mastitis resistance based on SCC is feasible. Patterns for genetic parameters within first lactation, however, require further confirmation and investigation.

Identification of seven haplotypes of the caprine PrP gene at codons 127, 142, 154, 211, 222 and 240 in French Alpine and Saanen breeds and their association with classical scrapie

In sheep, susceptibility to scrapie is mainly influenced by polymorphisms of the PrP gene. In goats, there are to date few data related to scrapie susceptibility association with PrP gene polymorphisms. In this study, we first investigated PrP gene polymorphisms of the French Alpine and Saanen breeds. Based on PrP gene open reading frame sequencing of artificial insemination bucks (n=404), six encoding mutations were identified at codons 127, 142, 154, 211, 222 and 240. However, only seven haplotypes could be detected: four (GIH(154)RQS, GIRQ(211)QS, GIRRK(222)S and GIRRQP(240)) derived from the wild-type allele (G(127)I(142)R(154)R(211)Q(222)S(240)) by a single-codon mutation, and two (S(127)IRRQP(240) and GM(142)RRQP(240)) by a double-codon mutation. A case-control study was then implemented in a highly affected Alpine and Saanen breed herd (90 cases/164 controls). Mutations at codon 142 (I/M), 154 (R/H), 211 (R/Q) and 222 (Q/K) were found to induce a significant degree of protection towards natural scrapie infection. Compared with the baseline homozygote wild-type genotype I(142)R(154)R(211)Q(222)/IRRQ goats, the odds of scrapie cases in IRQ(211)Q/IRRQ and IRRK(222)/IRRQ heterozygous animals were significantly lower [odds ratio (OR)=0.133, P<0.0001; and OR=0.048, P<0.0001, respectively]. The heterozygote M(142)RRQ/IRRQ genotype was only protective (OR=0.243, P=0.0186) in goats also PP(240) homozygous at codon 240. However, mutated allele frequencies in French Alpine and Saanen breeds were low (0.5-18.5 %), which prevent us from assessing the influence of all the possible genotypes in natural exposure conditions.

Response to somatic cell count-based selection for mastitis resistance in a divergent selection experiment in sheep

A divergent selection experiment in sheep was implemented to study the consequences of log-transformed somatic cell score (SCS)-based selection on resistance to natural intramammary infections. Using dams and progeny-tested rams selected for extreme breeding values for SCS, we created 2 groups of ewes with a strong divergence in SCS of approximately 3 genetic standard deviations. A survey of 84 first-lactation ewes of both the High and Low SCS lines indicated favorable responses to SCS-based selection on resistance to both clinical and subclinical mastitis. All clinical cases (n = 5) occurred in the High SCS line. Additionally, the frequency of chronic clinical mastitis, as detected by the presence of parenchymal abscesses, was much greater in the High SCS line (n = 21) than in the Low SCS line (n = 1). According to monthly milk bacteriological examinations of udder halves, the prevalence of infection was significantly greater (odds ratio = 3.1) in the High SCS line than in the Low SCS line, with predicted probabilities of 37 and 16%, respectively. The most frequently isolated bacteria responsible for mastitis were staphylococci: Staphylococcus auricularis (42.6% of positive samples), Staphylococcus simulans, Staphylococcus haemoliticus, Staphylococcus xylosus, Staphylococcus chromogenes, Staphylococcus lentus, Staphylococcus warneri, and Staphylococcus aureus. The incidence of positive bacteriology was greater in the High SCS line (39%) than in the Low SCS line (12%) at lambing, indicating that High SCS line ewes were especially susceptible to postpartum subclinical mastitis. Negativation of bacteriological results from one sampling time point to the next was markedly different between lines after weaning (e.g., 41 and 84% in the High and Low SCS lines, respectively). This result was consistent with differences in the duration of infection, which was much greater in the High SCS line compared with the Low SCS line. Finally, ewes from the High SCS line consistently had greater SCS in positive milk samples than did ewes from the Low SCS line (+2.04 SCS, on average), with an especially large difference between lines during the suckling period (+3.42 SCS). Altogether, the preliminary results suggest that the better resistance of Low SCS line ewes, compared with High SCS line ewes, was principally characterized by a better ability to limit infections during the peripartum period, to eliminate infections during lactation, and quantitatively to limit the inflammation process and its clinical consequences.

A genome scan for QTL affecting resistance to Haemonchus contortus in sheep.

Gastrointestinal nematodes are one of the main health issues in sheep breeding. To identify loci affecting the resistance to Haemonchus contortus, a genome scan was carried out using 1,275 Romane × Martinik Black Belly backcross lambs. The entire population was challenged with Haemonchus contortus in 2 consecutive experimental infections, and fecal egg counts (FEC) and packed cell volumes were measured. A subgroup of 332 lambs with extreme FEC was necropsied to determine the total worm burden, length of female worms, sex ratio in the worm population, abomasal pH, and serum and mucosal G immunoglobulins (IgG) responses. Pepsinogen concentration was measured in another subset of 229 lambs. For QTL detection, 160 microsatellite markers were used as well as the Illumina OvineSNP50 BeadChip that provided 42,469 SNP markers after quality control. Linkage, association, and joint linkage and association analyses were performed with the QTLMAP software. Linkage disequilibrium (LD) was estimated within each pure breed, and association analyses were carried out either considering or not the breed origin of the haplotypes. Four QTL regions on sheep chromosomes (OAR)5, 12, 13, and 21 were identified as key players among many other QTL with small to moderate effects. A QTL on OAR21 affecting pepsinogen concentration exactly matched the pepsinogen (PGA5) locus. A 10-Mbp region affecting FEC after the 1st and 2nd infections was found on OAR12. The SNP markers outperformed microsatellites in the linkage analysis. Taking advantage of the LD helped to refine the locations of the QTL mapped on OAR5 and 13.

Transcriptomic analysis of milk somatic cells in mastitis resistant and susceptible sheep upon challenge with Staphylococcus epidermidis and Staphylococcus aureus.

BackgroundThe existence of a genetic basis for host responses to bacterial intramammary infections has been widely documented, but the underlying mechanisms and the genes are still largely unknown. Previously, two divergent lines of sheep selected for high/low milk somatic cell scores have been shown to be respectively susceptible and resistant to intramammary infections by Staphylococcus spp. Transcriptional profiling with an 15K ovine-specific microarray of the milk somatic cells of susceptible and resistant sheep infected successively by S. epidermidis and S. aureus was performed in order to enhance our understanding of the molecular and cellular events associated with mastitis resistance.ResultsThe bacteriological titre was lower in the resistant than in the susceptible animals in the 48 hours following inoculation, although milk somatic cell concentration was similar. Gene expression was analysed in milk somatic cells, mainly represented by neutrophils, collected 12 hours post-challenge. A high number of differentially expressed genes between the two challenges indicated that more T cells are recruited upon inoculation by S. aureus than S. epidermidis. A total of 52 genes were significantly differentially expressed between the resistant and susceptible animals. Further Gene Ontology analysis indicated that differentially expressed genes were associated with immune and inflammatory responses, leukocyte adhesion, cell migration, and signal transduction. Close biological relationships could be established between most genes using gene network analysis. Furthermore, gene expression suggests that the cell turn-over, as a consequence of apoptosis/granulopoiesis, may be enhanced in the resistant line when compared to the susceptible line.ConclusionsGene profiling in resistant and susceptible lines has provided good candidates for mapping the biological pathways and genes underlying genetically determined resistance and susceptibility towards Staphylococcus infections, and opens new fields for further investigation.

Strengthening insights into host responses to mastitis infection in ruminants by combining heterogeneous microarray data sources

BackgroundGene expression profiling studies of mastitis in ruminants have provided key but fragmented knowledge for the understanding of the disease. A systematic combination of different expression profiling studies via meta-analysis techniques has the potential to test the extensibility of conclusions based on single studies. Using the program Pointillist, we performed meta-analysis of transcription-profiling data from six independent studies of infections with mammary gland pathogens, including samples from cattle challenged in vivo with S. aureus, E. coli, and S. uberis, samples from goats challenged in vivo with S. aureus, as well as cattle macrophages and ovine dendritic cells infected in vitro with S. aureus. We combined different time points from those studies, testing different responses to mastitis infection: overall (common signature), early stage, late stage, and cattle-specific.ResultsIngenuity Pathway Analysis of affected genes showed that the four meta-analysis combinations share biological functions and pathways (e.g. protein ubiquitination and polyamine regulation) which are intrinsic to the general disease response. In the overall response, pathways related to immune response and inflammation, as well as biological functions related to lipid metabolism were altered. This latter observation is consistent with the milk fat content depression commonly observed during mastitis infection. Complementarities between early and late stage responses were found, with a prominence of metabolic and stress signals in the early stage and of the immune response related to the lipid metabolism in the late stage; both mechanisms apparently modulated by few genes, including XBP1 and SREBF1.The cattle-specific response was characterized by alteration of the immune response and by modification of lipid metabolism. Comparison of E. coli and S. aureus infections in cattle in vivo revealed that affected genes showing opposite regulation had the same altered biological functions and provided evidence that E. coli caused a stronger host response.ConclusionsThis meta-analysis approach reinforces previous findings but also reveals several novel themes, including the involvement of genes, biological functions, and pathways that were not identified in individual studies. As such, it provides an interesting proof of principle for future studies combining information from diverse heterogeneous sources.

A first step toward genomic selection in the multi-breed French dairy goat population.

The objectives of this study were to describe, using the goat SNP50 BeadChip (Illumina Inc., San Diego, CA), molecular data for the French dairy goat population and compare the effect of using genomic information on breeding value accuracy in different reference populations. Several multi-breed (Alpine and Saanen) reference population sizes, including or excluding female genotypes (from 67 males to 677 males, and 1,985 females), were used. Genomic evaluations were performed using genomic best linear unbiased predictor for milk production traits, somatic cell score, and some udder type traits. At a marker distance of 50kb, the average r(2) (squared correlation coefficient) value of linkage disequilibrium was 0.14, and persistence of linkage disequilibrium as correlation of r-values among Saanen and Alpine breeds was 0.56. Genomic evaluation accuracies obtained from cross validation ranged from 36 to 53%. Biases of these estimations assessed by regression coefficients (from 0.73 to 0.98) of phenotypes on genomic breeding values were higher for traits such as protein yield than for udder type traits. Using the reference population that included all males and females, accuracies of genomic breeding values derived from prediction error variances (model accuracy) obtained for young buck candidates without phenotypes ranged from 52 to 56%. This was lower than the average pedigree-derived breeding value accuracies obtained at birth for these males from the official genetic evaluation (62%). Adding females to the reference population of 677 males improved accuracy by 5 to 9% depending on the trait considered. Gains in model accuracies of genomic breeding values ranged from 1 to 7%, lower than reported in other studies. The gains in breeding value accuracy obtained using genomic information were not as good as expected because of the limited size (at most 677 males and 1,985 females) and the structure of the reference population.