Cell specification and functional interactions in the pig blastocyst inferred from single-cell transcriptomics and uterine fluids proteomics

The embryonic development of the pig comprises a long in utero pre- and peri-implantation development, which dramatically differs from mice and humans. During this peri-implantation period, a complex series of paracrine signals establishes an intimate dialogue between the embryo and the uterus. To better understand the biology of the pig blastocyst during this period, we generated a large dataset of single-cell RNAseq from early and hatched blastocysts, spheroid and ovoid conceptus and proteomic datasets from corresponding uterine fluids. Our results confirm the molecular specificity and functionality of the three main cell populations. We also discovered two previously unknown subpopulations of the trophectoderm, one characterised by the expression of LRP2, which could represent progenitor cells, and the other, expressing pro-apoptotic markers, which could correspond to the Rauber's layer. Our work provides new insights into the biology of these populations, their reciprocal functional interactions, and the molecular dialogue with the maternal uterine environment.

MiRNome variations in milk fractions during feed restrictions of different intensities in dairy cows

BACKGROUND

In dairy cows, diet is one factor that can affect their milk production and composition. However, the effect of feed restriction on milk miRNome has not yet been described. Indeed, milk is the body fluid with the highest RNA concentration, which includes numerous microRNA. Its presence in the four different milk fractions, whole milk, fat globules, mammary epithelial cells and extracellular vesicles, is still poorly documented. This study aimed to describe the effects of different feed restrictions on the miRNome composition of different milk fractions.

RESULTS

Two feed restrictions were applied to lactating dairy cows, one of high intensity and one of moderate intensity. 2,896 mature microRNA were identified in the different milk fractions studied, including 1,493 that were already known in the bovine species. Among the 1,096 microRNA that were sufficiently abundant to be informative, the abundance of 1,027 of them varied between fractions: 36 of those were exclusive to one milk fraction. Feed restriction affected the abundance of 155 microRNA, with whole milk and milk extracellular vesicles being the most affected, whereas milk fat globules and exfoliated mammary epithelial cells were little or not affected at all. The high intensity feed restriction led to more microRNA variations in milk than moderate restriction. The target prediction of known microRNA that varied under feed restriction suggested the modification of some key pathways for lactation related to milk fat and protein metabolisms, cell cycle, and stress responses.

CONCLUSIONS

This study highlighted that the miRNome of each milk fraction is specific, with mostly the same microRNA composition but with variations in abundance between fractions. These specific miRNomes were affected differently by feed restrictions, the intensity of which appeared to be a major factor modulating milk miRNomes. These findings offer opportunities for future research on the use of milk miRNA as biomarkers of energy status in dairy cows, which is affected by feed restrictions.

Specific Milk Composition of miR-30b Transgenic Mice Associated with Early Duodenum Maturation in Offspring with Lasting Consequences for Growth

BACKGROUND

Milk composition is complex and includes numerous components essential for offspring growth and development. In addition to the high abundance of miR-30b microRNA, milk produced by the transgenic mouse model of miR-30b-mammary deregulation displays a significantly altered fatty acid profile. Moreover, wild-type adopted pups fed miR-30b milk present an early growth defect.

OBJECTIVE

This study aimed to investigate the consequences of miR-30b milk feeding on the duodenal development of wild-type neonates, a prime target of suckled milk, along with comprehensive milk phenotyping.

METHODS

The duodenums of wild-type pups fed miR-30b milk were extensively characterized at postnatal day (PND)-5, PND-6, and PND-15 using histological, transcriptomic, proteomic, and duodenal permeability analyses and compared with those of pups fed wild-type milk. Milk of miR-30b foster dams collected at mid-lactation was extensively analyzed using proteomic, metabolomic, and lipidomic approaches and hormonal immunoassays.

RESULTS

At PND-5, wild-type pups fed miR-30b milk showed maturation of their duodenum with 1.5-fold (P < 0.05) and 1.3-fold (P < 0.10) increased expression of Claudin-3 and Claudin-4, respectively, and changes in 8 duodenal proteins (P < 0.10), with an earlier reduction in paracellular and transcellular permeability (183 ng/mL fluorescein sulfonic acid [FSA] and 12 ng/mL horseradish peroxidase [HRP], respectively, compared with 5700 ng/mL FSA and 90 ng/mL HRP in wild-type; P < 0.001). Compared with wild-type milk, miR-30b milk displayed an increase in total lipid (219 g/L compared with 151 g/L; P < 0.05), ceramide (17.6 μM compared with 6.9 μM; P < 0.05), and sphingomyelin concentrations (163.7 μM compared with 76.3 μM; P < 0.05); overexpression of 9 proteins involved in the gut barrier (P < 0.1); and higher insulin and leptin concentrations (1.88 ng/mL and 2.04 ng/mL, respectively, compared with 0.79 ng/mL and 1.06 ng/mL; P < 0.01).

CONCLUSIONS

miR-30b milk displays significant changes in bioactive components associated with neonatal duodenal integrity and maturation, which could be involved in the earlier intestinal closure phenotype of the wild-type pups associated with a lower growth rate.

Potential genetic robustness of Prnp and Sprn double knockout mouse embryos towards ShRNA-lentiviral inoculation

The Shadoo and PrP prion protein family members are thought to be functionally related, but previous knockdown/knockout experiments in early mouse embryogenesis have provided seemingly contradictory results. In particular, Shadoo was found to be indispensable in the absence of PrP in knockdown analyses, but a double-knockout of the two had little phenotypic impact. We investigated this apparent discrepancy by comparing transcriptomes of WT, Prnp^0/0 and Prnp^0/0Sprn^0/0 E6.5 mouse embryos following inoculation by Sprn- or Prnp-ShRNA lentiviral vectors. Our results suggest the possibility of genetic adaptation in Prnp^0/0Sprn^0/0 mice, thus providing a potential explanation for their previously observed resilience.

A unified framework for the integration of multiple hierarchical clusterings or networks from multi-source data

Background

Integrating data from different sources is a recurring question in computational biology. Much effort has been devoted to the integration of data sets of the same type, typically multiple numerical data tables. However, data types are generally heterogeneous: it is a common place to gather data in the form of trees, networks or factorial maps, as these representations all have an appealing visual interpretation that helps to study grouping patterns and interactions between entities. The question we aim to answer in this paper is that of the integration of such representations.

Results

To this end, we provide a simple procedure to compare data with various types, in particular trees or networks, that relies essentially on two steps: the first step projects the representations into a common coordinate system; the second step then uses a multi-table integration approach to compare the projected data. We rely on efficient and well-known methodologies for each step: the projection step is achieved by retrieving a distance matrix for each representation form and then applying multidimensional scaling to provide a new set of coordinates from all the pairwise distances. The integration step is then achieved by applying a multiple factor analysis to the multiple tables of the new coordinates. This procedure provides tools to integrate and compare data available, for instance, as tree or network structures. Our approach is complementary to kernel methods, traditionally used to answer the same question.

Conclusion

Our approach is evaluated on simulation and used to analyze two real-world data sets: first, we compare several clusterings for different cell-types obtained from a transcriptomics single-cell data set in mouse embryos; second, we use our procedure to aggregate a multi-table data set from the TCGA breast cancer database, in order to compare several protein networks inferred for different breast cancer subtypes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04303-4.

Assessment of trade-offs between feed efficiency, growth-related traits, and immune activity in experimental lines of layer chickens

Background

In all organisms, life-history traits are constrained by trade-offs, which may represent physiological limitations or be related to energy resource management. To detect trade-offs within a population, one promising approach is the use of artificial selection, because intensive selection on one trait can induce unplanned changes in others. In chickens, the breeding industry has achieved remarkable genetic progress in production and feed efficiency over the last 60 years. However, this may have been accomplished at the expense of other important biological functions, such as immunity. In the present study, we used three experimental lines of layer chicken—two that have been divergently selected for feed efficiency and one that has been selected for increased antibody response to inactivated Newcastle disease virus (ND3)—to explore the impact of improved feed efficiency on animals’ immunocompetence and, vice versa, the impact of improved antibody response on animals’ growth and feed efficiency.

Results

There were detectable differences between the low (R+) and high (R−) feed-efficiency lines with respect to vaccine-specific antibody responses and counts of monocytes, heterophils, and/or T cell population. The ND3 line presented reduced body weight and feed intake compared to the control line. ND3 chickens also demonstrated an improved antibody response against a set of commercial viral vaccines, but lower blood leucocyte counts.

Conclusions

This study demonstrates the value of using experimental chicken lines that are divergently selected for RFI or for a high antibody production, to investigate the modulation of immune parameters in relation to growth and feed efficiency. Our results provide further evidence that long-term selection for the improvement of one trait may have consequences on other important biological functions. Hence, strategies to ensure optimal trade-offs among competing functions will ultimately be required in multi-trait selection programs in livestock.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00636-z.

PiQSARS: A pipeline for quantitative and statistical analyses of ratiometric fluorescent biosensors

Genetically encoded ratiometric fluorescent probes are cutting-edge tools in biology. They allow precise and dynamic measurement of various physiological parameters within cell compartments. Because data extraction and analysis are time consuming and may lead to inconsistencies between results, we describe here a standardized pipeline for•Semi-automated treatment of time-lapse fluorescence microscopy images.•Quantification of individual cell signal.•Statistical analysis of the data.First, a dedicated macro was developed using the FIJI software to reproducibly quantify the fluorescence ratio as a function of time. Raw data are then exported and analyzed using R and MATLAB softwares. Calculation and statistical analysis of selected graphic parameters are performed. In addition, a functional principal component analysis allows summarizing the dataset. Finally, a principal component analysis is performed to check consistency and final analysis is presented as a visual diagram. The method is adapted to any ratiometric fluorescent probe. As an example, the analysis of the cytoplasmic HyPer probe in response to an acute cell treatment with increasing amounts of hydrogen peroxide is shown. In conclusion, the pipeline allows to save time and analyze a larger amount of samples while reducing manual interventions and consequently increasing the robustness of the analysis.

Individualized multi-omic pathway deviation scores using multiple factor analysis

Malignant progression of normal tissue is typically driven by complex networks of somatic changes, including genetic mutations, copy number aberrations, epigenetic changes, and transcriptional reprogramming. To delineate aberrant multi-omic tumor features that correlate with clinical outcomes, we present a novel pathway-centric tool based on the multiple factor analysis framework called padma. Using a multi-omic consensus representation, padma quantifies and characterizes individualized pathway-specific multi-omic deviations and their underlying drivers, with respect to the sampled population. We demonstrate the utility of padma to correlate patient outcomes with complex genetic, epigenetic, and transcriptomic perturbations in clinically actionable pathways in breast and lung cancer.

Characterization of Holstein and Normande whole milk miRNomes highlights breed specificities

The concept of milk as a healthy food has opened the way for studies on milk components, from nutrients to microRNAs, molecules with broad regulatory properties present in large quantities in milk. Characterization of these components has been performed in several species, such as humans and bovine, depending on the stages of lactation. Here, we have studied the variation in milk microRNA composition according to genetic background. Using high throughput sequencing, we have characterized and compared the milk miRNomes of Holstein and Normande cattle, dairy breeds with distinct milk production features, in order to highlight microRNAs that are essential for regulation of the lactation process. In Holstein and Normande milk, 2,038 and 2,030 microRNAs were identified, respectively, with 1,771 common microRNAs, of which 1,049 were annotated and 722 were predicted. The comparison of the milk miRNomes of two breeds allowed to highlight 182 microRNAs displaying significant differences in the abundance. They are involved in the regulation of lipid metabolism and mammary morphogenesis and development, which affects lactation. Our results provide new insights into the regulation of molecular mechanisms involved in milk production.

A genomic map of climate adaptation in Mediterranean cattle breeds

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short-term evolutionary response to climate pressure induced by their post-domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome-wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.

The impact of genomic selection on genetic diversity and genetic gain in three French dairy cattle breeds

Background

In France, implementation of genomic evaluations in dairy cattle breeds started in 2009 and this has modified the breeding schemes drastically. In this context, the goal of our study was to understand the impact of genomic selection on the genetic diversity of bulls from three French dairy cattle breeds born between 2005 and 2015 (Montbéliarde, Normande and Holstein) and the factors that are involved.

Methods

We compared annual genetic gains, inbreeding rates based on runs of homozygosity (ROH) and pedigree data, and mean ROH length within breeds, before and after the implementation of genomic selection.

Results

Genomic selection induced an increase in mean annual genetic gains of 50, 71 and 33% for Montbéliarde, Normande and Holstein bulls, respectively, and in parallel, the generation intervals were reduced by a factor of 1.7, 1.9 and 2, respectively. We found no significant change in inbreeding rate for the two national breeds, Montbéliarde and Normande, and a significant increase in inbreeding rate for the Holstein international breed, which is now as high as 0.55% per year based on ROH and 0.49% per year based on pedigree data (equivalent to a rate of 1.36 and 1.39% per generation, respectively). The mean ROH length was longer for bulls from the Holstein breed than for those from the other two breeds.

Conclusions

With the implementation of genomic selection, the annual genetic gain increased for bulls from the three major French dairy cattle breeds. At the same time, the annual loss of genetic diversity increased for Holstein bulls, possibly because of the massive use of a few elite bulls in this breed, but not for Montbéliarde and Normande bulls. The increase in mean ROH length in Holstein may reflect the occurrence of recent inbreeding. New strategies in breeding schemes, such as female donor stations and embryo transfer, and recent implementation of genomic evaluations in small regional breeds should be studied carefully in order to ensure the sustainability of breeding schemes in the future.

Genetic improvement of canine hip dysplasia through sire selection across countries

Breeding against canine hip dysplasia (HD) may benefit from the importation of foreign sires. When foreign sires are evaluated on a different HD scale, this may diminish the efficacy. Using stochastic simulations, we evaluated genetic change and inbreeding levels for different scenarios of importing sires with high genetic merit for HD. Population size and genetic parameters (e.g. heritability, accuracy of selection, genetic correlation) were based on actual data for HD in Golden retrievers and Labrador retrievers in the UK and Sweden. For countries with different HD scales and an estimated breeding value (EBV) evaluation in place, the importation was useful if imported sires had EBV rankings in the top 50% and if genetic correlations between EBV systems were above 0.85. When importing sires with EBV rankings in the top 10%, moderate accuracies of EBVs (>0.40) and moderately strong genetic correlations (>0.70) were needed. Selection against HD without the importation of sires may increase inbreeding levels, while the importation of sires can decrease inbreeding levels. For national genetic evaluation and selection programmes, importing sires with high genetic merit can be an effective breeding strategy, but care is needed to estimate reliable EBVs.

Domestic Fowl Breed Variation in Egg White Protein Expression: Application of Proteomics and Transcriptomics

Avian egg white is essential for protecting and nourishing bird embryos during their development. Being produced in the female magnum, variability in hen oviduct gene expression may affect egg white composition in domestic chickens. Since traditional poultry breeds may represent a source of variation, in the present study we describe the egg white proteome (mass spectrometry) and corresponding magnum transcriptome (high-throughput sequencing) for 20 hens from five domestic fowl breeds (large breeds: Araucana, Czech golden pencilled, Minorca; and small breeds: Booted bantam, Rosecomb bantam). In total, we identified 189 egg white proteins and 16391 magnum-expressed genes. The majority of egg white protein content comprised proteins with an antimicrobial function. Despite general similarity, Between-class Principal Component Analysis revealed significant breed-specific variability in protein abundances, differentiating especially small and large breeds. Though we found strong association between magnum mRNA expression and egg white protein abundance across genes, coinertia analysis revealed no transcriptome/proteome costructure at the individual level. Our study is the first to show variation in protein abundances in egg white across chicken breeds with potential effects on egg quality, biosafety, and chick development. The observed interindividual variation probably results from post-transcriptional regulation creating a discrepancy between proteomic and transcriptomic data.

Detection of new pathways involved in the acceptance and the utilisation of a plant-based diet in isogenic lines of rainbow trout fry

To meet the growing demand of fish feed for aquaculture, an increasing proportion of marine ingredients are being replaced by blends of plant products. However, the total replacement of marine ingredients in salmonid diets impairs fish performance. This is particularly true during the early fry stage and this stage is therefore considered of particular importance. In rainbow trout (RBT), the existence of a genetic variability to survive and grow with plant-based diets devoid of marine ingredients has now been proved, but the mechanisms behind are little studied especially at early stage. To investigate these, we analysed the whole transcriptome of three isogenic lines of RBT fry, which have similar growth when fed a marine resources-based diet (M diet) but which highly differ in their responses to a plant-based diet (V diet). Analysis of transcriptomes profiles revealed 1740, 1834 and 246 probes differentially expressed among the three genotypes when fed the V diet. The use of these lines led to the discovery of potential molecular markers linked to plant-based diet utilisation, some of them belonging to new pathways, never described before. An important number of genes was related to immunity, but further investigations are needed to better understand the difference between the genotypes in their immune status response to V diet exposure. Finally, differences in expression of genes related to feed intake and sensory perception among genotypes suggested that the mechanisms underlying the differences in growth on plant-based diet are closely linked to diet acceptance. Research on plants components affecting feed intake should be thus further explored.

Genetic and transcriptomic analyses provide new insights on the early antiviral response to VHSV in resistant and susceptible rainbow trout

BACKGROUND

The viral hemorrhagic septicemia virus (VHSV) is a major threat for salmonid farming and for wild fish populations worldwide. Previous studies have highlighted the importance of innate factors regulated by a major quantitative trait locus (QTL) for the natural resistance to waterborne VHSV infection in rainbow trout. The aim of this study was to analyze the early transcriptomic response to VHSV inoculation in cell lines derived from previously described resistant and susceptible homozygous isogenic lines of rainbow trout to obtain insights into the molecular mechanisms responsible for the resistance to the viral infection.

RESULTS

We first confirmed the presence of the major QTL in a backcross involving a highly resistant fish isogenic line (B57) and a highly susceptible one (A22), and were able to define the confidence interval of the QTL and to identify its precise position. We extended the definition of the QTL since it controls not only resistance to waterborne infection but also the kinetics of mortality after intra-peritoneal injection. Deep sequencing of the transcriptome of B57 and A22 derived cell lines exposed to inactivated VHSV showed a stronger response to virus inoculation in the resistant background. In line with our previous observations, an early and strong induction of interferon and interferon-stimulated genes was correlated with the resistance to VHSV, highlighting the major role of innate immune factors in natural trout resistance to the virus. Interestingly, major factors of the antiviral innate immunity were much more expressed in naive B57 cells compared to naive A22 cells, which likely contributes to the ability of B57 to mount a fast antiviral response after viral infection. These observations were further extended by the identification of several innate immune-related genes localized close to the QTL area on the rainbow trout genome.

CONCLUSIONS

Taken together, our results improve our knowledge in virus-host interactions in vertebrates and provide novel insights in the molecular mechanisms explaining the resistance to VHSV in rainbow trout. Our data also provide a collection of potential markers for resistance and susceptibility of rainbow trout to VHSV infection.

Sunflower oil supplementation affects the expression of miR-20a-5p and miR-142-5p in the lactating bovine mammary gland

Oil supplementation in dairy cattle diets is used to modulate milk fat composition, as well as the expression of mammary lipogenic genes, whose regulation remains unclear. MiRNAs are small non-coding RNA considered as crucial regulators of gene expression, offering clues to explain the mechanism underlying gene nutriregulation. The present study was designed to identify miRNAs whose expression in the cow mammary gland is modulated by sunflower oil supplementation. MiRNomes were obtained using RNAseq technology from the mammary gland of lactating cows receiving a low forage diet, supplemented or not with 4% sunflower oil. Among the 272 miRNAs characterized, eight were selected for RT-qPCR validations, showing the significant down-regulation of miR-142-5p and miR-20a-5p by sunflower supplementation. These two miRNAs are predicted to target genes whose expression was reported as differentially expressed by sunflower supplementation. Among their putative targets, ELOVL6 gene involved in lipid metabolism has been studied. However, a first analysis did not show its significant down-regulation, in response to the over-expression of miR-142-5p, of miR-20a-5p, or both, in a bovine mammary epithelial cell line. However, a clearer understanding of the miRNA expression by lipid supplementation would help to decipher the regulation of lactating cow mammary gland in response to nutrition.

Epithelial response to a high-protein diet in rat colon

Background

High-protein diets (HPD) alter the large intestine microbiota composition in association with a metabolic shift towards protein degradation. Some amino acid-derived metabolites produced by the colon bacteria are beneficial for the mucosa while others are deleterious at high concentrations. The aim of the present work was to define the colonic epithelial response to an HPD. Transcriptome profiling was performed on colonocytes of rats fed an HPD or an isocaloric normal-protein diet (NPD) for 2 weeks.

Results

The HPD downregulated the expression of genes notably implicated in pathways related to cellular metabolism, NF-κB signaling, DNA repair, glutathione metabolism and cellular adhesion in colonocytes. In contrast, the HPD upregulated the expression of genes related to cell proliferation and chemical barrier function. These changes at the mRNA level in colonocytes were not associated with detrimental effects of the HPD on DNA integrity (comet assay), epithelium renewal (quantification of proliferation and apoptosis markers by immunohistochemistry and western blot) and colonic barrier integrity (Ussing chamber experiments).

Conclusion

The modifications of the luminal environment after an HPD were associated with maintenance of the colonic homeostasis that might be the result of adaptive processes in the epithelium related to the observed transcriptional regulations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3514-z) contains supplementary material, which is available to authorized users.

Differences during the first lactation between cows cloned by somatic cell nuclear transfer and noncloned cows

Lactation performance is dependent on both the genetic characteristics and the environmental conditions surrounding lactating cows. However, individual variations can still be observed within a given breed under similar environmental conditions. The role of the environment between birth and lactation could be better appreciated in cloned cows, which are presumed to be genetically identical, but differences in lactation performance between cloned and noncloned cows first need to be clearly evaluated. Conflicting results have been described in the literature, so our aim was to clarify this situation. Nine cloned Prim' Holstein cows were produced by the transfer of nuclei from a single fibroblast cell line after cell fusion with enucleated oocytes. The cloned cows and 9 noncloned counterparts were raised under similar conditions. Milk production and composition were recorded monthly from calving until 200d in milk. At 67d in milk, biopsies were sampled from the rear quarter of the udder, their mammary epithelial cell content was evaluated, and mammary cell renewal, RNA, and DNA were then analyzed in relevant samples. The results showed that milk production did not differ significantly between cloned and noncloned cows, but milk protein and fat contents were less variable in cloned cows. Furthermore, milk fat yield and contents were lower in cloned cows during early lactation. At around 67 DIM, milk fat and protein yields, as well as milk fat, protein, and lactose contents, were also lower in cloned cows. These lower yields could be linked to the higher apoptotic rate observed in cloned cows. Apoptosis is triggered by insulin-like factor growth binding protein 5 (IGFBP5) and plasminogen activator inhibitor (PAI), which both interact with CSN1S2. During our experiments, CSN1S2 transcript levels were lower in the mammary gland of cloned cows. The mammary cell apoptotic rate observed in cloned cows may have been related to the higher levels of DNA (cytosine-5-)-methyltransferase 1 (DNMT1) transcripts, coding for products that maintain the epigenetic status of cells. We conclude, therefore, that milk production in cloned cows differs slightly from that of noncloned cows. These differences may be due, in part, to a higher incidence of subclinical mastitis. They were associated with differences in cell apoptosis and linked to variations in DNMT1 mRNA. However, milk protein and fat contents were more similar among cloned cows than among noncloned cows.

Milk from dams fed an obesogenic diet combined with a high-fat/high-sugar diet induces long-term abnormal mammary gland development in the rabbit

Alterations to the metabolic endocrine environment during early life are crucial to mammary gland development. Among these environmental parameters, the initial nutritional event after birth is the consumption of milk, which represents the first maternal support provided to mammalian newborns. Milk is a complex fluid that exerts effects far beyond its immediate nutritional value. The present study, therefore, aimed to determine the effect of the nutritional changes during the neonatal and prepubertal periods on the adult mammary phenotype. Newborn rabbits were suckled by dams fed a high-fat/high-sugar obesogenic (OD) or a control (CON) diet and then subsequently fed either the OD or CON diets from the onset of puberty and throughout early pregnancy. Mammary glands were collected during early pregnancy (Day 8 of pregnancy). Rabbits fed with OD milk and then subjected to an OD diet displayed an abnormal development of the mammary gland: the mammary ducts were markedly enlarged (P < 0.05) and filled with abundant secretory products. Moreover, the alveolar secretory structures were disorganized, with an abnormal aspect characterized by large lumina. Mammary epithelial cells contained numerous large lipid droplets and exhibited fingering of the apical membrane and abnormally enlarged intercellular spaces filled with casein micelles. Leptin has been shown to be involved in modulating several developmental processes. We therefore analyzed its expression in the mammary gland. Mammary leptin mRNA was strongly expressed in rabbits fed with OD milk and subjected to an OD diet by comparison with the CON rabbits. Leptin transcripts and protein were localized in the epithelial cells, indicating that the increase in leptin synthesis occurs in this compartment. Taken together, these findings suggest that early-life nutritional history, in particular through the milking period, can determine subsequent mammary gland development. Moreover, they highlight the potentially important regulatory role that leptin may play during critical early-life nutritional windows with respect to long-term growth and mammary function.

WIDDE: a Web-Interfaced next generation database for genetic diversity exploration, with a first application in cattle

BACKGROUND

The advent and democratization of next generation sequencing and genotyping technologies lead to a huge amount of data for the characterization of population genetic diversity in model and non model-species. However, efficient storage, management, cross-analyzing and exploration of such dense genotyping datasets remain challenging. This is particularly true for the bovine species where many SNP datasets have been generated in various cattle populations with different genotyping tools.

DESCRIPTION

We developed WIDDE, a Web-Interfaced Next Generation Database that stands as a generic tool applicable to a wide range of species and marker types ( http://widde.toulouse.inra.fr). As a first illustration, we hereby describe its first version dedicated to cattle biodiversity, which includes a large and evolving cattle genotyping dataset for over 750,000 SNPs available on 129 (89 public) different cattle populations representative of the world-wide bovine genetic diversity and on 7 outgroup bovid species. This version proposes an optional marker and individual filtering step, an export of genotyping data in different popular formats, and an exploration of genetic diversity through a principal component analysis. Users can also explore their own genotyping data together with data from WIDDE, assign their samples to WIDDE populations based on distance assignment method and supervised clustering, and estimate their ancestry composition relative to the populations represented in the database.

CONCLUSION

The cattle version of WIDDE represents to our knowledge the first database dedicated to cattle biodiversity and SNP genotyping data that will be very useful for researchers interested in this field. As a generic tool applicable to a wide range of marker types, WIDDE is overall intended to the genetic diversity exploration of any species and will be extended to other species shortly. The structure makes it easy to include additional output formats and new tools dedicated to genetic diversity exploration.

DNA methylation and transcription in a distal region upstream from the bovine AlphaS1 casein gene after once or twice daily milking

Once daily milking (ODM) induces a reduction in milk production when compared to twice daily milking (TDM). Unilateral ODM of one udder half and TDM of the other half, enables the study of underlying mechanisms independently of inter-individual variability (same genetic background) and of environmental factors. Our results show that in first-calf heifers three CpG, located 10 kb upstream from the CSN1S1 gene were methylated to 33, 34 and 28%, respectively, after TDM but these levels were higher after ODM, 38, 38 and 33%, respectively. These methylation levels were much lower than those observed in the mammary gland during pregnancy (57, 59 and 50%, respectively) or in the liver (74, 78 and 61%, respectively). The methylation level of a fourth CpG (CpG4), located close by (29% during TDM) was not altered after ODM. CpG4 methylation reached 39.7% and 59.5%, during pregnancy or in the liver, respectively. CpG4 is located within a weak STAT5 binding element, arranged in tandem with a second high affinity STAT5 element. STAT5 binding is only marginally modulated by CpG4 methylation, but it may be altered by the methylation levels of the three other CpG nearby. Our results therefore shed light on mechanisms that help to explain how milk production is almost, but not fully, restored when TDM is resumed (15.1±0.2 kg/day instead of 16.2±0.2 kg/day, p<0.01). The STAT5 elements are 100 bp away from a region transcribed in the antisense orientation, in the mammary gland during lactation, but not during pregnancy or in other reproductive organs (ovary or testes). We now need to clarify whether the transcription of this novel RNA is a consequence of STAT5 interacting with the CSN1S1 distal region, or whether it plays a role in the chromatin structure of this region.

Stable methylation at promoters distinguishes epiblast stem cells from embryonic stem cells and the in vivo epiblasts

Embryonic Stem Cells (ESCs) and Epiblast Stem Cells (EpiSCs) are the in vitro representatives of naïve and primed pluripotency, respectively. It is currently unclear how their epigenomes underpin the phenotypic and molecular characteristics of these distinct pluripotent states. Here, we performed a genome-wide comparison of DNA methylation between ESCs and EpiSCs by MethylCap-Seq. We observe that promoters are preferential targets for methylation in EpiSC compared to ESCs, in particular high CpG island promoters. This is in line with upregulation of the de novo methyltransferases Dnmt3a1 and Dnmt3b in EpiSC, and downregulation of the demethylases Tet1 and Tet2. Remarkably, the observed DNA methylation signature is specific to EpiSCs and differs from that of their in vivo counterpart, the postimplantation epiblast. Using a subset of promoters that are differentially methylated, we show that DNA methylation is established within a few days during in vitro outgrowth of the epiblast, and also occurs when ESCs are converted to EpiSCs in vitro. Once established, this methylation is stable, as ES-like cells obtained by in vitro reversion of EpiSCs display an epigenetic memory that only extensive passaging and sub-cloning are able to almost completely erase.

Characterisation and comparison of lactating mouse and bovine mammary gland miRNomes

BACKGROUND

The mammary gland is a dynamic organ that undergoes important physiological changes during reproductive cycles. Until now, data regarding the characterisation of miRNA in the mammary gland have been scarce and mainly focused on their abnormal expression in breast cancer. Our goal was to characterise the microRNA (miRNA) involved in mechanisms regulating the mammary function, with particular focus on the lactation stage.

METHODOLOGY/PRINCIPAL FINDINGS

Using high-throughput sequencing technology, the exhaustive repertoires of miRNA expressed (miRNome) in mouse and bovine mammary glands during established lactation were identified, characterized and compared. Furthermore, in order to obtain more information on miRNA loading in the RNA-induced silencing complex (RISC), the miRNome was compared with that obtained from RNA associated with the AGO2 protein (AGO2-miRNome) in mouse lactating mammary gland. This study enabled the identification of 164 and 167 miRNA in mouse and bovine, respectively. Among the 30 miRNA most highly expressed in each species, 24 were common to both species and six of them were preferentially highly expressed in lactating than non-lactating mammary gland. The potential functional roles of these 24 miRNA were deduced using DIANA-miRPath software, based on miRNA/mRNA interactions. Moreover, seven putative novel miRNA were identified. Using DAVID analysis, it was concluded that the predicted targets of two of these putative novel miRNA are involved in mammary gland morphogenesis.

CONCLUSION/SIGNIFICANCE

Our study provides an overview of the characteristics of lactating mouse and bovine mammary gland miRNA expression profiles. Moreover, species-conserved miRNA involved in this fundamental biological function were identified. These miRNomes will now be used as references for further studies during which the impact of animal breeding on the miRNA expression will be analysed.

Out of America: tracing the genetic footprints of the global diffusion of maize

Maize was first domesticated in a restricted valley in south-central Mexico. It was diffused throughout the Americas over thousands of years, and following the discovery of the New World by Columbus, was introduced into Europe. Trade and colonization introduced it further into all parts of the world to which it could adapt. Repeated introductions, local selection and adaptation, a highly diverse gene pool and outcrossing nature, and global trade in maize led to difficulty understanding exactly where the diversity of many of the local maize landraces originated. This is particularly true in Africa and Asia, where historical accounts are scarce or contradictory. Knowledge of post-domestication movements of maize around the world would assist in germplasm conservation and plant breeding efforts. To this end, we used SSR markers to genotype multiple individuals from hundreds of representative landraces from around the world. Applying a multidisciplinary approach combining genetic, linguistic, and historical data, we reconstructed possible patterns of maize diffusion throughout the world from American "contribution" centers, which we propose reflect the origins of maize worldwide. These results shed new light on introductions of maize into Africa and Asia. By providing a first globally comprehensive genetic characterization of landraces using markers appropriate to this evolutionary time frame, we explore the post-domestication evolutionary history of maize and highlight original diversity sources that may be tapped for plant improvement in different regions of the world.

Molecular fingerprinting of the podocyte reveals novel gene and protein regulatory networks

A thorough characterization of the transcriptome and proteome of endogenous podocytes has been hampered by low cell yields during isolation. Here we describe a double fluorescent reporter mouse model combined with an optimized bead perfusion protocol and efficient single cell dissociation to yield more than 500,000 podocytes per mouse allowing for global, unbiased downstream applications. Combining mRNA and miRNA transcriptional profiling with quantitative proteomic analyses revealed programs of highly specific gene regulation tightly controlling cytoskeleton, cell differentiation, endosomal transport, and peroxisome function in podocytes. Strikingly, the analyses further predict that these podocyte-specific gene regulatory networks are accompanied by alternative splicing of respective genes. Thus, our 'omics' approach will facilitate the discovery and integration of novel gene, protein, and organelle regulatory networks that deepen our systematic understanding of podocyte biology.

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.

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.

Trends of the genetic connectedness measures among Nelore beef cattle herds

Validity of comparisons between expected breeding values obtained from best linear unbiased prediction procedures in genetic evaluations is dependent on genetic connectedness among herds. Different cattle breeding programmes have their own particular features that distinguish their database structure and can affect connectedness. Thus, the evolution of these programmes can also alter the connectedness measures. This study analysed the evolution of the genetic connectedness measures among Brazilian Nelore cattle herds from 1999 to 2008, using the French Criterion of Admission to the group of Connected Herds (CACO) method, based on coefficients of determination (CD) of contrasts. Genetic connectedness levels were analysed by using simple and multiple regression analyses on herd descriptors to understand their relationship and their temporal trends from the 1999-2003 to the 2004-2008 period. The results showed a high level of genetic connectedness, with CACO estimates higher than 0.4 for the majority of them. Evaluation of the last 5-year period showed only a small increase in average CACO measures compared with the first 5 years, from 0.77 to 0.80. The percentage of herds with CACO estimates lower than 0.7 decreased from 27.5% in the first period to 16.2% in the last one. The connectedness measures were correlated with percentage of progeny from connecting sires, and the artificial insemination spread among Brazilian herds in recent years. But changes in connectedness levels were shown to be more complex, and their complete explanation cannot consider only herd descriptors. They involve more comprehensive changes in the relationship matrix, which can be only fully expressed by the CD of contrasts.

Genetic structure of the European Charolais and Limousin cattle metapopulations using pedigree analyses

Pedigree collected by the Interbeef service allowed genetic diversity to be assessed by using pedigree analyses for the European Charolais (CHA) and Limousin (LIM) cattle populations registered in national herdbooks in Denmark (DNK), France (FRA), Ireland (IRL), Sweden (SWE), and, solely for the LIM breed, the United Kingdom (UK). The CHA data set included 2,563,189 calves with weaning performance, of which 96.1% were recorded in FRA, 3.0% in SWE, 0.5% in IRL, and 0.4% in DNK. The LIM data set included 1,652,734 calves with weaning performance, of which 91.9% were recorded in FRA, 4.9% in UK, 1.8% in DNK, 0.9% SWE, and 0.5% in IRL. Pedigree files included 3,191,132 CHA and 2,409,659 LIM animals. Gene flows were rather limited between populations, except from FRA toward other countries. Pedigree completeness was good in all subpopulations for both breeds and allowed the pedigree to be traced back to the French population. A relatively high level of genetic diversity was assessed in each CHA and LIM subpopulation by estimating either effective population sizes (N(e) >244 and N(e) >345 in the CHA and LIM subpopulations, respectively), relationship coefficients within subpopulations (<1.3% in both breeds), or probability of gene origins. However, in each subpopulation, it was shown that founders and also ancestors had unbalanced genetic contributions, leading to a moderate but continuous reduction in genetic diversity. Analyses between populations suggested that all European CHA and LIM populations were differentiated very little. The Swedish CHA population was assessed as genetically more distant from the other CHA populations because of fewer gene flows from other countries and because of the use of North American sires to introgress the polled phenotype. In each European subpopulation, most of the main ancestors, which explained 50% of gene origin, were born in FRA. However, those main ancestors were different between countries. Moreover, in both breeds, the main ancestors, which explained 50% of the gene origin in DNK, IRL, SWE, and UK for the LIM breed, were found to be infrequently used in FRA. Those results were consistent with the low relationship coefficients estimated between subpopulations (<0.6% in both the CHA and LIM breeds). Therefore, in both breeds, each subpopulation may constitute a reservoir of genetic diversity for the other ones.

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 individual's 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.1<h2≤0.4) or high (h2>0.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.

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.

Genetic and morphometric characterization of a local Vietnamese swamp buffalo population

The water buffalo plays a key role in the socio-economy of South-East Asia as it is the main draught power for paddy rice cultivation. While in the Indian subcontinent the water buffalo is the riverine type, in South-East Asia the majority of buffaloes are of the swamp type. In the poor remote northern province of Ha Giang in Vietnam, improvement of the swamp buffalo breed may be one of the best ways to increase sustainability of farming systems. Therefore, analysis of the genetic structure of the province buffalo population is a prerequisite to any conservation or improvement project. A total of 1122 animals were described for 11 body and horn measurements for morphometric characterization. From this sample set, 744 animals were genotyped for 17 microsatellite markers. Also 17 animals from southern provinces of Vietnam were genotyped as a comparative sample. The results showed that genetic diversity as well as inbreeding value in the Ha Giang was high. The F(ST) values within the province and across Vietnam were low indicating that most of the population variation is explained by individual variability. Bayesian clustering analysis did not highlight the presence of subdivided populations. These results are useful for the implementation of a conservation and improvement strategy of the swamp buffalo in order to guarantee the householders' needs for sustainability of the farming system in the Ha Giang province.

Measuring connectedness among herds in mixed linear models: from theory to practice in large-sized genetic evaluations

A procedure to measure connectedness among groups in large-sized genetic evaluations is presented. It consists of two steps: (a) computing coefficients of determination (CD) of comparisons among groups of animals; and (b) building sets of connected groups. The CD of comparisons were estimated using a sampling-based method that estimates empirical variances of true and predicted breeding values from a simulated n-sample. A clustering method that may handle a large number of comparisons and build compact clusters of connected groups was developed. An aggregation criterion (Caco) that reflects the level of connectedness of each herd was computed. This procedure was validated using a small beef data set. It was applied to the French genetic evaluation of the beef breed with most records and to the genetic evaluation of goats. Caco was more related to the type of service of sires used in the herds than to herd size. It was very sensitive to the percentage of missing sires. Disconnected herds were reliably identified by low values of Caco. In France, this procedure is the reference method for evaluating connectedness among the herds involved in on-farm genetic evaluation of beef cattle (IBOVAL) since 2002 and for genetic evaluation of goats from 2007 onwards.

Consensus genetic structuring and typological value of markers using multiple co-inertia analysis

Working with weakly congruent markers means that consensus genetic structuring of populations requires methods explicitly devoted to this purpose. The method, which is presented here, belongs to the multivariate analyses. This method consists of different steps. First, single-marker analyses were performed using a version of principal component analysis, which is designed for allelic frequencies (%PCA). Drawing confidence ellipses around the population positions enhances %PCA plots. Second, a multiple co-inertia analysis (MCOA) was performed, which reveals the common features of single-marker analyses, builds a reference structure and makes it possible to compare single-marker structures with this reference through graphical tools. Finally, a typological value is provided for each marker. The typological value measures the efficiency of a marker to structure populations in the same way as other markers. In this study, we evaluate the interest and the efficiency of this method applied to a European and African bovine microsatellite data set. The typological value differs among markers, indicating that some markers are more efficient in displaying a consensus typology than others. Moreover, efficient markers in one collection of populations do not remain efficient in others. The number of markers used in a study is not a sufficient criterion to judge its reliability. "Quantity is not quality".

Use of structured antedependence models for the genetic analysis of growth curves

Growth curve analysis is an important issue for many agricultural and laboratory species, for both phenotypic and genetic studies. The aim of this paper is to present the use of a novel statistical approach, namely the structured antedependence (SAD) models, to deal with this issue. The basic idea of these models is that an observation at time t can be explained by the previous observations. These models are especially appropriate to deal with cumulative traits such as growth, as BW at age t clearly depends on BW measures at ages (t -1), (t -2), etc. These models were applied on an INRA experimental Charolais herd data set. The data comprised BW records for 560 cows born over an 11-yr period (from 1988 to 1998) from 60 sires and 369 dams. The proposed SAD models were compared with the well-known random regression (RR) models that are already widely used in various areas of longitudinal data analysis. It was found that the SAD models fit the growth process better with far fewer parameters than the RR models (9 instead of 16 covariance parameters for the phenotypic analysis, and 14 instead of 21 for the genetic analysis). Despite this smaller number of covariance parameters, the likelihood value was found to be much higher with the SAD vs. the RR models, with a difference of 262.9 for the phenotypic analysis with a quartic polynomial for the RR and 751.5 for the genetic analysis with a cubic polynomial for both the genetic and environmental parts of the RR model. The SAD models also proved to be better able to interpolate missing values. Heritability, genetic, and environmental correlation coefficients were estimated for weights from birth to adulthood. The structured antedependence models proved, in this study, to be very appropriate to model growth data in a parsimonious and flexible way.

Evaluation models and genetic parameters for calving difficulty in beef cattle

Calving difficulty was analyzed under threshold and linear models considering either a fixed or random herd-year effect. The aim of the study was to compare models for predicting breeding values according to the size of herd-year groups. When simulating data sets with small herds, in order to obtain an unbiased evaluation under a nonrandom and negative association of sire and herd effects, the best model for a practical evaluation was the fixed linear model. Field data included 246,576 records of the largest Charolais herds in France. Models were compared using the correlations of estimated breeding values between the different models. Although the best model from a theoretical point of view was a threshold model with a fixed herd-year effect, a linear model with a fixed herd-year effect was the best choice from a practical point of view for predicting direct effects for calving difficulty in beef cattle and was a sufficient choice for predicting the associated maternal effects for data set with large herds. Correlations between direct estimated breeding values under the reference model and the fixed linear model and the random threshold model were 0.94 and 0.91, respectively. Correlations between the corresponding maternal estimated breeding values were 0.94 and 0.98. Heritabilities of direct effects were 0.27 and 0.14 under fixed threshold and fixed linear models, respectively. The corresponding heritabilities of maternal effects were 0.18 and 0.13, and the genetic correlation between direct and maternal effects were -0.36 and -0.34, respectively.

Is a Multivariate Consensus Representation of Genetic Relationships Among Populations Always Meaningful?

To determine the relationships among closely related populations or species, two methods are commonly used in the literature: phylogenetic reconstruction or multivariate analysis. The aim of this article is to assess the reliability of multivariate analysis. We describe a method that is based on principal component analysis and Mantel correlations, using a two-step process: The first step consists of a single-marker analysis and the second step tests if each marker reveals the same typology concerning population differentiation. We conclude that if single markers are not congruent, the compromise structure is not meaningful. Our model is not based on any particular mutation process and it can be applied to most of the commonly used genetic markers. This method is also useful to determine the contribution of each marker to the typology of populations. We test whether our method is efficient with two real data sets based on microsatellite markers. Our analysis suggests that for closely related populations, it is not always possible to accept the hypothesis that an increase in the number of markers will increase the reliability of the typology analysis.

A sampling method for estimating the accuracy of predicted breeding values in genetic evaluation

A sampling-based method for estimating the accuracy of estimated breeding values using an animal model is presented. Empirical variances of true and estimated breeding values were estimated from a simulated n-sample. The method was validated using a small data set from the Parthenaise breed with the estimated coefficient of determination converging to the true values. It was applied to the French Salers data file used for the 2000 on-farm evaluation (IBOVAL) of muscle development score. A drawback of the method is its computational demand. Consequently, convergence can not be achieved in a reasonable time for very large data files. Two advantages of the method are that a) it is applicable to any model (animal, sire, multivariate, maternal effects...) and b) it supplies off-diagonal coefficients of the inverse of the mixed model equations and can therefore be the basis of connectedness studies.

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 F(st) = 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.

Analysis of genetic relationships between 10 cattle breeds with 17 microsatellites

To guide genetic conservation programmes with objective criteria, general genetic variability has to be taken into account. This study was conducted to determine the genetic variation between 10 cattle breeds by using 17 microsatellite loci and 13 biochemical markers (11 blood groups, the transferrin and beta-casein loci). Microsatellite loci were amplified in 31-50 unrelated individuals from 10 cattle breeds: Charolais, Limousin, Breton Black Pied, Parthenais, Montbéliard, Vosgien, Maine-Anjou, Normande, Jersey and Holstein. Neighbor-joining trees were calculated from genetic distance estimates. The robustness of tree topology was obtained by bootstrap resampling of loci. A total of 210 alleles of the 17 microsatellites were detected in this study and average heterozygosities ranged from 0.53 in the Jersey breed to 0.66 in the Parthenais breed. In general, low bootstrap values were obtained: with the 17 microsatellites, the highest bootstrap values concerned the Holstein/Maine-Anjou grouping with an occurrence of 74%; with the biochemical markers, this node had an occurrence of 79% and the Charolais/Limousin grouping appeared with an occurrence of 74%; when microsatellites and biochemical polymorphism were analysed together, the occurrence of the Holstein/Maine-Anjou grouping was 90% and that of the Charolais/Limousin grouping was 42%. These results suggest that 30 microsatellites, a number currently considered as sufficient to distinguish closely related breeds is, in fact, probably insufficient.