Searching for homozygous haplotype deficiency in Manech Tête Rousse dairy sheep revealed a nonsense variant in the MMUT gene affecting newborn lamb viability.

BACKGROUND: Recessive deleterious variants are known to segregate in livestock populations, as in humans, and some may be lethal in the homozygous state. RESULTS: We used phased 50 k single nucleotide polymorphism (SNP) genotypes and pedigree data to scan the genome of 6845 Manech Tête Rousse dairy sheep to search for deficiency in homozygous haplotypes (DHH). Five Manech Tête Rousse deficient homozygous haplotypes (MTRDHH1 to 5) were identified, with a homozygous deficiency ranging from 84 to 100%. These haplotypes are located on Ovis aries chromosome (OAR)1 (MTRDHH2 and 3), OAR10 (MTRDHH4), OAR13 (MTRDHH5), and OAR20 (MTRDHH1), and have carrier frequencies ranging from 7.8 to 16.6%. When comparing at-risk matings between DHH carriers to safe matings between non-carriers, two DHH (MTRDHH1 and 2) were linked with decreased insemination success and/or increased stillbirth incidence. We investigated the MTRDHH1 haplotype, which substantially increased stillbirth rate, and identified a single nucleotide variant (SNV) inducing a premature stop codon (p.Gln409*) in the methylmalonyl-CoA mutase (MMUT) gene by using a whole-genome sequencing approach. We generated homozygous lambs for the MMUT mutation by at-risk mating between heterozygous carriers, and most of them died within the first 24 h after birth without any obvious clinical symptoms. Reverse transcriptase-qPCR and western blotting on post-mortem liver and kidney biological samples showed a decreased expression of MMUT mRNA in the liver and absence of a full-length MMUT protein in the mutant homozygous lambs. CONCLUSIONS: We identified five homozygous deficient haplotypes that are likely to harbor five independent deleterious recessive variants in sheep. One of these was detected in the MMUT gene, which is associated with lamb lethality in the homozygous state. A specific management of these haplotypes/variants in the MTR dairy sheep selection program would help enhance the overall fertility and lamb survival.

A Nonsense Variant in CCDC65 Gene Causes Respiratory Failure Associated with Increased Lamb Mortality in French Lacaune Dairy Sheep.

We recently demonstrated that the Lacaune deficient homozygous haplotype 6 (LDHH6) potentially hosts a recessive perinatal lethal mutation in Lacaune dairy sheep mapped on OAR3. In the present study, we have analyzed the whole-genome sequences of two Lacaune ram heterozygous carriers of LDHH6. After variant calling and filtering against the variants of 86 non-carrier rams, we have identified a single nucleotide variant (SNV) in the two LDHH6 carriers whose variant allele induced a premature stop codon (p.Glu111*) in the Coiled-Coil Domain Containing 65 (CCDC65) gene. CCDC65 is involved in the assembly of the nexin-dynein regulatory complex for the formation of microtubules in ciliated cells. In order to identify the phenotype in homozygous sheep, we generated at-risk matings (n = 17) between rams and ewes heterozygous for the candidate variant in CCDC65. A total of 16 lambs were born alive with five genotyped as homozygous carriers. The homozygous lambs suffered from respiratory problems, and four of them died within the first month of life. At necropsy, we observed a broad hepatization of lung lobes possibly induced by infectious pneumonia. The management of this lethal recessive allele (frequency of 0.06) through reasoned mating in the Lacaune sheep selection schemes could reduce lamb mortality by 2%.

Genome-Wide Identification of a Regulatory Mutation in BMP15 Controlling Prolificacy in Sheep.

The search for the genetic determinism of prolificacy variability in sheep has evidenced several major mutations in genes playing a crucial role in the control of ovulation rate. In the Noire du Velay (NV) sheep population, a recent genetic study has evidenced the segregation of such a mutation named FecL L . However, based on litter size (LS) records of FecL L non-carrier ewes, the segregation of a second prolificacy major mutation was suspected in this population. In order to identify this mutation, we have combined a case/control genome-wide association study with ovine 50k SNP chip genotyping, whole genome sequencing, and functional analyses. A new single nucleotide polymorphism (OARX:50977717T > A, NC_019484) located on the X chromosome upstream of the BMP15 gene was evidenced to be highly associated with the prolificacy variability (P = 1.93E-11). The variant allele was called FecX N and shown to segregate also in the Blanche du Massif Central (BMC) sheep population. In both NV and BMC, the FecX N allele frequency was estimated close to 0.10, and its effect on LS was estimated at +0.20 lamb per lambing at the heterozygous state. Homozygous FecX N carrier ewes were fertile with increased prolificacy in contrast to numerous mutations affecting BMP15. At the molecular level, FecX N was shown to decrease BMP15 promoter activity and supposed to impact BMP15 expression in the oocyte. This regulatory action was proposed as the causal mechanism for the FecX N mutation to control ovulation rate and prolificacy in sheep.

The high prolificacy of D'man sheep is associated with the segregation of the FecLL mutation in the B4GALNT2 gene.

Mutations in the FecL locus are associated with large variation in ovulation rate and litter size in the French Lacaune sheep breed. It has been shown that the B4GALNT2 gene within the FecL locus is most likely responsible for the high fecundity in the French breed. In this study, we have highlighted the segregation of the FecLL mutation within the B4GALNT2 gene in North African sheep breeds and notably in the highly prolific D'man breed. Genotyping of a sample of 183 Tunisian D'man individuals revealed a high frequency (0.65) of the prolific allele FecLL which was attributed to the adoption of a decades-old breeding strategy based on the selection of ewe lambs born from large litter size. Homozygous LL ewes showed a significantly increased litter size compared to heterozygous and non-carrier ewes (FecLL /FecLL  = 2.47 ± 0.09 vs. FecLL /FecL+  = 2.23 ± 0.09, p < 0.05 and FecL+ /FecL+  = 1.93 ± 0.18, p < 0.01). The presence of the FecLL polymorphism in both D'man and Lacaune breeds argues for an ancestral origin of this mutation and brings an answer to the old question of the genetic determinism of the extreme prolificacy of the D'man ewes. The results of this study can help to establish planned genotype-based mating allowing both higher profit for the breeders and an optimal management of the FecLL mutation in D'man sheep populations.

Author Correction: A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content.

The quantity of milk and milk fat and proteins are particularly important traits in dairy livestock. However, little is known about the regions of the genome that influence these traits in goats. We conducted a genome wide association study in French goats and identified 109 regions associated with dairy traits. For a major region on chromosome 14 closely associated with fat content, the Diacylglycerol O-Acyltransferase 1 (DGAT1) gene turned out to be a functional and positional candidate gene. The caprine reference sequence of this gene was completed and 29 polymorphisms were found in the gene sequence, including two novel exonic mutations: R251L and R396W, leading to substitutions in the protein sequence. The R251L mutation was found in the Saanen breed at a frequency of 3.5% and the R396W mutation both in the Saanen and Alpine breeds at a frequencies of 13% and 7% respectively. The R396W mutation explained 46% of the genetic variance of the trait, and the R251L mutation 6%. Both mutations were associated with a notable decrease in milk fat content. Their causality was then demonstrated by a functional test. These results provide new knowledge on the genetic basis of milk synthesis and will help improve the management of the French dairy goat breeding program.

FecX Bar a Novel BMP15 mutation responsible for prolificacy and female sterility in Tunisian Barbarine Sheep.

BACKGROUND: Naturally occurring mutations in growth and differentiation factor 9 (GDF9) or bone morphogenetic protein 15 (BMP15) genes are associated with increased ovulation rate (OR) and litter size (LS) but also sterility. Observing the Tunisian Barbarine ewes of the "W" flock selected for improved prolificacy, we found prolific and infertile ewes with streaky ovaries. Blood genomic DNA was extracted from a subset of low-ovulating, prolific and infertile ewes of the "W" flock, and the entire coding sequences of GDF9 and BMP15 were sequenced. RESULTS: We evidenced a novel polymorphism in the exon 1 of the BMP15 gene associated with increased prolificacy and sterility. This novel mutation called FecX Bar is a composite polymorphism associating a single nucleotide substitution (c.301G > T), a 3 bp deletion (c.302_304delCTA) and a C insertion (c.310insC) in the ovine BMP15 cDNA leading to a frame shift at protein position 101. Calculated in the "W" flock, the FecX Bar allele increased OR by 0.7 ova and LS by 0.3 lambs (p = 0.08). As for already identified mutations, homozygous females carrying FecX Bar exhibited streaky ovaries with a blockade at the primary stage of folliculogenesis as shown by histochemistry. CONCLUSIONS: Our investigation demonstrates a new mutation in the BMP15 gene providing a valuable genetic tool to control fecundity in Tunisian Barbarine, usable for diffusion program into conventional flocks looking for prolificacy improvement.

Genome-Wide Identification of the Mutation Underlying Fleece Variation and Discriminating Ancestral Hairy Species from Modern Woolly Sheep.

The composition and structure of fleece variation observed in mammals is a consequence of a strong selective pressure for fiber production after domestication. In sheep, fleece variation discriminates ancestral species carrying a long and hairy fleece from modern domestic sheep (Ovis aries) owning a short and woolly fleece. Here, we report that the "woolly" allele results from the insertion of an antisense EIF2S2 retrogene (called asEIF2S2) into the 3' UTR of the IRF2BP2 gene leading to an abnormal IRF2BP2 transcript. We provide evidence that this chimeric IRF2BP2/asEIF2S2 messenger 1) targets the genuine sense EIF2S2 RNA and 2) creates a long endogenous double-stranded RNA which alters the expression of both EIF2S2 and IRF2BP2 mRNA. This represents a unique example of a phenotype arising via a RNA-RNA hybrid, itself generated through a retroposition mechanism. Our results bring new insights on the sheep population history thanks to the identification of the molecular origin of an evolutionary phenotypic variation.

A Point Mutation in Suppressor of Cytokine Signalling 2 (Socs2) Increases the Susceptibility to Inflammation of the Mammary Gland while Associated with Higher Body Weight and Size and Higher Milk Production in a Sheep Model.

Mastitis is an infectious disease mainly caused by bacteria invading the mammary gland. Genetic control of susceptibility to mastitis has been widely evidenced in dairy ruminants, but the genetic basis and underlying mechanisms are still largely unknown. We describe the discovery, fine mapping and functional characterization of a genetic variant associated with elevated milk leukocytes count, or SCC, as a proxy for mastitis. After implementing genome-wide association studies, we identified a major QTL associated with SCC on ovine chromosome 3. Fine mapping of the region, using full sequencing with 12X coverage in three animals, provided one strong candidate SNP that mapped to the coding sequence of a highly conserved gene, suppressor of cytokine signalling 2 (Socs2). The frequency of the SNP associated with increased SCC was 21.7% and the Socs2 genotype explained 12% of the variance of the trait. The point mutation induces the p.R96C substitution in the SH2 functional domain of SOCS2 i.e. the binding site of the protein to various ligands, as well-established for the growth hormone receptor GHR. Using surface plasmon resonance we showed that the p.R96C point mutation completely abrogates SOCS2 binding affinity for the phosphopeptide of GHR. Additionally, the size, weight and milk production in p.R96C homozygote sheep, were significantly increased by 24%, 18%, and 4.4%, respectively, when compared to wild type sheep, supporting the view that the point mutation causes a loss of SOCS2 functional activity. Altogether these results provide strong evidence for a causal mutation controlling SCC in sheep and highlight the major role of SOCS2 as a tradeoff between the host's inflammatory response to mammary infections, and body growth and milk production, which are all mediated by the JAK/STAT signaling pathway.

An overview of gene expression dynamics during early ovarian folliculogenesis: specificity of follicular compartments and bi-directional dialog.

BACKGROUND: Successful early folliculogenesis is crucial for female reproductive function. It requires appropriate gene specific expression of the different types of ovarian cells at different developmental stages. To date, most gene expression studies on the ovary were conducted in rodents and did not distinguish the type of cell. In mono-ovulating species, few studies have addressed gene expression profiles and mainly concerned human oocytes. RESULTS: We used a laser capture microdissection method combined with RNA-seq technology to explore the transcriptome in oocytes and granulosa cells (GCs) during development of the sheep ovarian follicle. We first documented the expression profile of 15 349 genes, then focused on the 5 129 genes showing differential expression between oocytes and GCs. Enriched functional categories such as oocyte meiotic arrest and GC steroid synthesis reflect two distinct cell fates. We identified the implication of GC signal transduction pathways such as SHH, WNT and RHO GTPase. In addition, signaling pathways (VEGF, NOTCH, IGF1, etc.) and GC transzonal projections suggest the existence of complex cell-cell interactions. Finally, we highlighted several transcription regulators and specifically expressed genes that likely play an important role in early folliculogenesis. CONCLUSIONS: To our knowledge, this is the first comprehensive exploration of transcriptomes derived from in vivo oocytes and GCs at key stages in early follicular development in sheep. Collectively, our data advance our understanding of early folliculogenesis in mono-ovulating species and will be a valuable resource for unraveling human ovarian dysfunction such as premature ovarian failure (POF).

The highly prolific phenotype of Lacaune sheep is associated with an ectopic expression of the B4GALNT2 gene within the ovary.

Prolific sheep have proven to be a valuable model to identify genes and mutations implicated in female fertility. In the Lacaune sheep breed, large variation in litter size is genetically determined by the segregation of a fecundity major gene influencing ovulation rate, named FecL and its prolific allele FecL(L) . Our previous work localized FecL on sheep chromosome 11 within a locus of 1.1 Mb encompassing 20 genes. With the aim to identify the FecL gene, we developed a high throughput sequencing strategy of long-range PCR fragments spanning the locus of FecL(L) carrier and non-carrier ewes. Resulting informative markers defined a new 194.6 kb minimal interval. The reduced FecL locus contained only two genes, insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) and beta-1,4-N-acetyl-galactosaminyl transferase 2 (B4GALNT2), and we identified two SNP in complete linkage disequilibrium with FecL(L) . B4GALNT2 appeared as the best positional and expressional candidate for FecL, since it showed an ectopic expression in the ovarian follicles of FecL(L) /FecL(L) ewes at mRNA and protein levels. In FecL(L) carrier ewes only, B4GALNT2 transferase activity was localized in granulosa cells and specifically glycosylated proteins were detected in granulosa cell extracts and follicular fluids. The identification of these glycoproteins by mass spectrometry revealed at least 10 proteins, including inhibin alpha and betaA subunits, as potential targets of B4GALNT2 activity. Specific ovarian protein glycosylation by B4GALNT2 is proposed as a new mechanism of ovulation rate regulation in sheep, and could contribute to open new fields of investigation to understand female infertility pathogenesis.