Copy number variations (CNV) are known to play a major role in genetic variability and disease pathogenesis in several species including cattle. In this study, we report the identification and characterization of CNV in eight French beef and dairy breeds using whole-genome sequence data from 200 animals. Bioinformatics analyses to search for CNV were carried out using four different but complementary tools and we validated a subset of the CNV by both in silico and experimental approaches. We report the identification and localization of 4178 putative deletion-only, duplication-only and CNV regions, which cover 6% of the bovine autosomal genome; they were validated by two in silico approaches and/or experimentally validated using array-based comparative genomic hybridization and single nucleotide polymorphism genotyping arrays. The size of these variants ranged from 334 bp to 7.7 Mb, with an average size of ~ 54 kb. Of these 4178 variants, 3940 were deletions, 67 were duplications and 171 corresponded to both deletions and duplications, which were defined as potential CNV regions. Gene content analysis revealed that, among these variants, 1100 deletions and duplications encompassed 1803 known genes, which affect a wide spectrum of molecular functions, and 1095 overlapped with known QTL regions. Our study is a large-scale survey of CNV in eight French dairy and beef breeds. These CNV will be useful to study the link between genetic variability and economically important traits, and to improve our knowledge on the genomic architecture of cattle.

Knowledge of population structure is essential to improve the management and conservation of farm animal genetic resources. Microsatellites, which have long been popular for this type of analysis, are more and more neglected in favor of whole-genome single nucleotide polymorphism (SNP) chips that are now available for the main farmed animal species. In this study, we compared genetic patterns derived from microsatellites to that inferred by SNPs, considering three pairs of datasets of sheep and cattle. Population genetic differentiation analyses (Fixation index, F_ST), as well as STRUCTURE analyses showed a very strong consistency between the two types of markers. Microsatellites gave pictures that were largely concordant with SNPs, although less accurate. The best concordance was found in the most complex dataset, which included 17 French sheep breeds (with a Pearson correlation coefficient of 0.95 considering the 136 values of pairwise F_ST, obtained with both types of markers). The use of microsatellites reduces the cost and the related analyses do not require specific computer equipment (i.e., information technology (IT) infrastructure able to provide adequate computing and storage capacity). Therefore, this tool may still be a very appropriate solution to evaluate, in a first stage, the general state of livestock at national scales. At a time when local breeds are disappearing at an alarming rate, it is urgent to improve our knowledge of them, in particular by promoting tools accessible to the greatest number.

In recent years, several bovine genome sequencing projects were carried out with the aim of developing genomic tools to improve dairy and beef production efficiency and sustainability. In this study, we describe the first French cattle genome variation dataset obtained by sequencing 274 whole genomes representing several major dairy and beef breeds. This dataset contains over 28 million single nucleotide polymorphisms (SNPs) and small insertions and deletions. Comparisons between sequencing results and SNP array genotypes revealed a very high genotype concordance rate, which indicates the good quality of our data. To our knowledge, this is the first large-scale catalog of small genomic variations in French dairy and beef cattle. This resource will contribute to the study of gene functions and population structure and also help to improve traits through genotype-guided selection.

Background: Myostatin, a member of the TGFβ superfamily, is well known as a potent and specific negative regulator of muscle growth. Targeting the myostatin signalling pathway may offer promising therapeutic strategies for the treatment of muscle-wasting disorders. In the last decade, various myostatin-binding proteins have been identified to be able to inhibit myostatin activity. One of these is GASP1 (Growth and Differentiation Factor-Associated Serum Protein-1), a protein containing a follistatin domain as well as multiple domains associated with protease inhibitors. Despite in vitro data, remarkably little is known about in vivo functions of Gasp1. To further address the role of GASP1 during mouse development and in adulthood, we generated a gain-of-function transgenic mouse model that overexpresses Gasp1 under transcriptional control of the human cytomegalovirus immediate-early promoter/enhancer. Results: Overexpression of Gasp1 led to an increase in muscle mass observed not before day 15 of postnatal life. The surGasp1 transgenic mice did not display any other gross abnormality. Histological and morphometric analysis of surGasp1 rectus femoris muscles revealed an increase in myofiber size without a corresponding increase in myofiber number. Fiber-type distribution was unaltered. Interestingly, we do not detect a change in total fat mass and lean mass. These results differ from those for myostatin knockout mice, transgenic mice overexpressing the myostatin propeptide or follistatin which exhibit both muscle hypertrophy and hyperplasia, and show minimal fat deposition. Conclusions: Altogether, our data give new insight into the in vivo functions of Gasp1. As an extracellular regulatory factor in the myostatin signalling pathway, additional studies on GASP1 and its homolog GASP2 are required to elucidate the crosstalk between the different intrinsic inhibitors of the myostatin.

Background: Linkage disequilibrium (LD) is a key parameter to study the history of populations and to identify and fine map quantitative trait loci (QTL) and it has been studied for many years in animal populations. The advent of new genotyping technologies has allowed whole-genome LD studies in most cattle populations. However, to date, long-range LD (LRLD) between distant variants on the genome has not been investigated in detail in cattle. Here, we present the first comprehensive study of LRLD in French beef cattle by analysing data on 672 Charolais (CHA), 462 Limousine (LIM) and 326 Blonde d’Aquitaine (BLA) individuals that were genotyped on the Illumina BovineHD Beadchip. Furthermore, whole-genome LD and haplotype block structure were analysed in these three breeds. Results: We computed linkage disequilibrium ( r ² ) values for 5.9, 5.6 and 6.0 billion pairs of SNPs on the 29 autosomes of CHA, LIM and BLA, respectively. Mean r ² values drop to less than 0.1 for distances between SNPs greater than 120 kb. However, for the first time, we detected the existence of LRLD in the three main French beef breeds. In total, 598, 266, and 795 LRLD events ( r ² ≥ 0.6) were detected in CHA, LIM and BLA, respectively. Each breed had predominantly population-specific LRLD interactions, although shared LRLD events occurred in a number of regions (55 LRLD events were shared between two breeds and nine between the three breeds). Examples of possible functional gene interactions and QTL co-location were observed with some of these LRLD events, which suggests epistatic selection. Conclusions: We identified long-range linkage disequilibrium for the first time in French beef cattle populations. Epistatic selection may be the main source of the observed LRLD events, but other forces may also be involved. LRLD information should be accounted for in genome-wide association studies.