Denis Laloë

A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis

During the last 3 years, a number of approaches for the normalization of RNA sequencing data have emerged in the literature, differing both in the type of bias adjustment and in the statistical strategy adopted. However, as data continue to accumulate, there has been no clear consensus on the appropriate normalization method to be used or the impact of a chosen method on the downstream analysis. In this work, we focus on a comprehensive comparison of seven recently proposed normalization methods for the differential analysis of RNA-seq data, with an emphasis on the use of varied real and simulated datasets involving different species and experimental designs to represent data characteristics commonly observed in practice. Based on this comparison study, we propose practical recommendations on the appropriate normalization method to be used and its impact on the differential analysis of RNA-seq data.

Genetic diversity measures of local European beef cattle breeds for conservation purposes

This study was undertaken to determine the genetic structure, evolutionary relationships, and the genetic diversity among 18 local cattle breeds from Spain, Portugal, and France using 16 microsatellites. Heterozygosities, estimates of Fst, genetic distances, multivariate and diversity analyses, and assignment tests were performed. Heterozygosities ranged from 0.54 in the Pirenaica breed to 0.72 in the Barrosã breed. Seven percent of the total genetic variability can be attributed to differences among breeds (mean Fst = 0.07; P < 0.01). Five different genetic distances were computed and compared with no correlation found to be significantly different from 0 between distances based on the effective size of the population and those which use the size of the alleles. The Weitzman recursive approach and a multivariate analysis were used to measure the contribution of the breeds diversity. The Weitzman approach suggests that the most important breeds to be preserved are those grouped into two clusters: the cluster formed by the Mirandesa and Alistana breeds and that of the Sayaguesa and Tudanca breeds. The hypothetical extinction of one of those clusters represents a 17% loss of diversity. A correspondence analysis not only distinguished four breed groups but also confirmed results of previous studies classifying the important breeds contributing to diversity. In addition, the variation between breeds was sufficiently high so as to allow individuals to be assigned to their breed of origin with a probability of 99% for simulated samples.

Maximizing the Reliability of Genomic Selection by Optimizing the Calibration Set of Reference Individuals: Comparison of Methods in Two Diverse Groups of Maize Inbreds ( Zea mays L.)

Genomic selection refers to the use of genotypic information for predicting breeding values of selection candidates. A prediction formula is calibrated with the genotypes and phenotypes of reference individuals constituting the calibration set. The size and the composition of this set are essential parameters affecting the prediction reliabilities. The objective of this study was to maximize reliabilities by optimizing the calibration set. Different criteria based on the diversity or on the prediction error variance (PEV) derived from the realized additive relationship matrix–best linear unbiased predictions model (RA–BLUP) were used to select the reference individuals. For the latter, we considered the mean of the PEV of the contrasts between each selection candidate and the mean of the population (PEVmean) and the mean of the expected reliabilities of the same contrasts (CDmean). These criteria were tested with phenotypic data collected on two diversity panels of maize (Zea mays L.) genotyped with a 50k SNPs array. In the two panels, samples chosen based on CDmean gave higher reliabilities than random samples for various calibration set sizes. CDmean also appeared superior to PEVmean, which can be explained by the fact that it takes into account the reduction of variance due to the relatedness between individuals. Selected samples were close to optimality for a wide range of trait heritabilities, which suggests that the strategy presented here can efficiently sample subsets in panels of inbred lines. A script to optimize reference samples based on CDmean is available on request.

A whole genome Bayesian scan for adaptive genetic divergence in West African cattle

BackgroundThe recent settlement of cattle in West Africa after several waves of migration from remote centres of domestication has imposed dramatic changes in their environmental conditions, in particular through exposure to new pathogens. West African cattle populations thus represent an appealing model to unravel the genome response to adaptation to tropical conditions. The purpose of this study was to identify footprints of adaptive selection at the whole genome level in a newly collected data set comprising 36,320 SNPs genotyped in 9 West African cattle populations.ResultsAfter a detailed analysis of population structure, we performed a scan for SNP differentiation via a previously proposed Bayesian procedure including extensions to improve the detection of loci under selection. Based on these results we identified 53 genomic regions and 42 strong candidate genes. Their physiological functions were mainly related to immune response (MHC region which was found under strong balancing selection, CD79A, CXCR4, DLK1, RFX3, SEMA4A, TICAM1 and TRIM21), nervous system (NEUROD6, OLFM2, MAGI1, SEMA4A and HTR4) and skin and hair properties (EDNRB, TRSP1 and KRTAP8-1).ConclusionThe main possible underlying selective pressures may be related to climatic conditions but also to the host response to pathogens such as Trypanosoma(sp). Overall, these results might open the way towards the identification of important variants involved in adaptation to tropical conditions and in particular to resistance to tropical infectious diseases.

Insights into the Genetic History of French Cattle from Dense SNP Data on 47 Worldwide Breeds

Background Modern cattle originate from populations of the wild extinct aurochs through a few domestication events which occurred about 8,000 years ago. Newly domesticated populations subsequently spread worldwide following breeder migration routes. The resulting complex historical origins associated with both natural and artificial selection have led to the differentiation of numerous different cattle breeds displaying a broad phenotypic variety over a short period of time. Methodology/Principal Findings This study gives a detailed assessment of cattle genetic diversity based on 1,121 individuals sampled in 47 populations from different parts of the world (with a special focus on French cattle) genotyped for 44,706 autosomal SNPs. The analyzed data set consisted of new genotypes for 296 individuals representing 14 French cattle breeds which were combined to those available from three previously published studies. After characterizing SNP polymorphism in the different populations, we performed a detailed analysis of genetic structure at both the individual and population levels. We further searched for spatial patterns of genetic diversity among 23 European populations, most of them being of French origin, under the recently developed spatial Principal Component analysis framework. Conclusions/Significance Overall, such high throughput genotyping data confirmed a clear partitioning of the cattle genetic diversity into distinct breeds. In addition, patterns of differentiation among the three main groups of populations—the African taurine, the European taurine and zebus—may provide some additional support for three distinct domestication centres. Finally, among the European cattle breeds investigated, spatial patterns of genetic diversity were found in good agreement with the two main migration routes towards France, initially postulated based on archeological evidence.

MICROSATELLITE-BASED PHYLOGENY OF INDIAN DOMESTIC GOATS

BackgroundThe domestic goat is one of the important livestock species of India. In the present study we assess genetic diversity of Indian goats using 17 microsatellite markers. Breeds were sampled from their natural habitat, covering different agroclimatic zones.ResultsThe mean number of alleles per locus (NA) ranged from 8.1 in Barbari to 9.7 in Jakhrana goats. The mean expected heterozygosity (He) ranged from 0.739 in Barbari to 0.783 in Jakhrana goats. Deviations from Hardy-Weinberg Equilibrium (HWE) were statistically significant (P < 0.05) for 5 loci breed combinations. The DA measure of genetic distance between pairs of breeds indicated that the lowest distance was between Marwari and Sirohi (0.135). The highest distance was between Pashmina and Black Bengal. An analysis of molecular variance indicated that 6.59% of variance exists among the Indian goat breeds. Both a phylogenetic tree and Principal Component Analysis showed the distribution of breeds in two major clusters with respect to their geographic distribution.ConclusionOur study concludes that Indian goat populations can be classified into distinct genetic groups or breeds based on the microsatellites as well as mtDNA information.

Using molecular markers and multivariate methods to study the genetic diversity of local European and Asian chicken breeds

French and Asian subsets of chicken breeds were first analysed using 22 microsatellites and then compared to the AVIANDIV European set using 14 loci. Positive correlations were observed between F(IT) or F(ST) and typological values or variance of markers using the multivariate analysis mcoa. The first axis of the multivariate representation separated Asian from European breeds, revealing breeds with Asian ancestor. Using all or 14 loci, correct assignation rate was always higher than 93%. The Weitzman index and the aggregate diversity D were calculated using 22 loci within French and Asian breeds. The French breed Coucou de Rennes and the Hua-Tung breed seemed to contribute the most to the global diversity of each subset. This approach on French-only breeds and then on French with AVIANDIV domestic breeds (14 loci) showed that the Marans breed contributed the most. The AVIANDIV framework could be useful to evaluate the genetic diversity of local breeds and to help in connecting national and regional conservation policies.

Immunity Traits in Pigs: Substantial Genetic Variation and Limited Covariation

Background Increasing robustness via improvement of resistance to pathogens is a major selection objective in livestock breeding. As resistance traits are difficult or impossible to measure directly, potential indirect criteria are measures of immune traits (ITs). Our underlying hypothesis is that levels of ITs with no focus on specific pathogens define an individuals immunocompetence and thus predict response to pathogens in general. Since variation in ITs depends on genetic, environmental and probably epigenetic factors, our aim was to estimate the relative importance of genetics. In this report, we present a large genetic survey of innate and adaptive ITs in pig families bred in the same environment. Methodology/Principal Findings Fifty four ITs were studied on 443 Large White pigs vaccinated against Mycoplasma hyopneumoniae and analyzed by combining a principal component analysis (PCA) and genetic parameter estimation. ITs include specific and non specific antibodies, seric inflammatory proteins, cell subsets by hemogram and flow cytometry, ex vivo production of cytokines (IFNα, TNFα, IL6, IL8, IL12, IFNγ, IL2, IL4, IL10), phagocytosis and lymphocyte proliferation. While six ITs had heritabilities that were weak or not significantly different from zero, 18 and 30 ITs had moderate (0.10.4) heritability values, respectively. Phenotypic and genetic correlations between ITs were weak except for a few traits that mostly include cell subsets. PCA revealed no cluster of innate or adaptive ITs. Conclusions/Significance Our results demonstrate that variation in many innate and adaptive ITs is genetically controlled in swine, as already reported for a smaller number of traits by other laboratories. A limited redundancy of the traits was also observed confirming the high degree of complementarity between innate and adaptive ITs. Our data provide a genetic framework for choosing ITs to be included as selection criteria in multitrait selection programmes that aim to improve both production and health traits.

Precision and information in linear models of genetic evaluation

Some criteria for measuring the overall precision of a genetic evaluation using linear mixed-model methodology are presented. They are derived via an extension of the coefficient of determination to linear combinations of estimates and via the use of the Kullback information. A parallel is drawn between inestimability of fixed-effects contrasts and the zero coefficient of determination for contrasts of random effects. The procedure is illustrated with 2 minor hypothetical examples of genetic evaluation based on an animal model and on a sire model.

Overexpression of miR-30b in the Developing Mouse Mammary Gland Causes a Lactation Defect and Delays Involution

Background MicroRNA (miRNA) are negative regulators of gene expression, capable of exerting pronounced influences upon the translation and stability of mRNA. They are potential regulators of normal mammary gland development and of the maintenance of mammary epithelial progenitor cells. This study was undertaken to determine the role of miR-30b on the establishment of a functional mouse mammary gland. miR-30b is a member of the miR-30 family, composed of 6 miRNA that are highly conserved in vertebrates. It has been suggested to play a role in the differentiation of several cell types. Methodology/Principal Findings The expression of miR-30b was found to be regulated during mammary gland development. Transgenic mice overexpressing miR-30b in mammary epithelial cells were used to investigate its role. During lactation, mammary histological analysis of the transgenic mice showed a reduction in the size of alveolar lumen, a defect of the lipid droplets and a growth defect of pups fed by transgenic females. Moreover some mammary epithelial differentiated structures persisted during involution, suggesting a delay in the process. The genes whose expression was affected by the overexpression of miR-30b were characterized by microarray analysis. Conclusion/Significance Our data suggests that miR-30b is important for the biology of the mammary gland and demonstrates that the deregulation of only one miRNA could affect lactation and involution.