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.

Spatio-Temporal Gene Expression Profiling during In Vivo Early Ovarian Folliculogenesis: Integrated Transcriptomic Study and Molecular Signature of Early Follicular Growth.

BACKGROUND: The successful achievement of early ovarian folliculogenesis is important for fertility and reproductive life span. This complex biological process requires the appropriate expression of numerous genes at each developmental stage, in each follicular compartment. Relatively little is known at present about the molecular mechanisms that drive this process, and most gene expression studies have been performed in rodents and without considering the different follicular compartments. RESULTS: We used RNA-seq technology to explore the sheep transcriptome during early ovarian follicular development in the two main compartments: oocytes and granulosa cells. We documented the differential expression of 3,015 genes during this phase and described the gene expression dynamic specific to these compartments. We showed that important steps occurred during primary/secondary transition in sheep. We also described the in vivo molecular course of a number of pathways. In oocytes, these pathways documented the chronology of the acquisition of meiotic competence, migration and cellular organization, while in granulosa cells they concerned adhesion, the formation of cytoplasmic projections and steroid synthesis. This study proposes the involvement in this process of several members of the integrin and BMP families. The expression of genes such as Kruppel-like factor 9 (KLF9) and BMP binding endothelial regulator (BMPER) was highlighted for the first time during early follicular development, and their proteins were also predicted to be involved in gene regulation. Finally, we selected a data set of 24 biomarkers that enabled the discrimination of early follicular stages and thus offer a molecular signature of early follicular growth. This set of biomarkers includes known genes such as SPO11 meiotic protein covalently bound to DSB (SPO11), bone morphogenetic protein 15 (BMP15) and WEE1 homolog 2 (S. pombe)(WEE2) which play critical roles in follicular development but other biomarkers are also likely to play significant roles in this process. CONCLUSIONS: To our knowledge, this is the first in vivo spatio-temporal exploration of transcriptomes derived from early follicles in sheep.

Design and characterization of a 52K SNP chip for goats.

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

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.

Genome-wide association studies identify two novel BMP15 mutations responsible for an atypical hyperprolificacy phenotype in sheep.

Some sheep breeds are naturally prolific, and they are very informative for the studies of reproductive genetics and physiology. Major genes increasing litter size (LS) and ovulation rate (OR) were suspected in the French Grivette and the Polish Olkuska sheep populations, respectively. To identify genetic variants responsible for the highly prolific phenotype in these two breeds, genome-wide association studies (GWAS) followed by complementary genetic and functional analyses were performed. Highly prolific ewes (cases) and normal prolific ewes (controls) from each breed were genotyped using the Illumina OvineSNP50 Genotyping Beadchip. In both populations, an X chromosome region, close to the BMP15 gene, harbored clusters of markers with suggestive evidence of association at significance levels between 1E(-05) and 1E(-07). The BMP15 candidate gene was then sequenced, and two novel non-conservative mutations called FecX(Gr) and FecX(O) were identified in the Grivette and Olkuska breeds, respectively. The two mutations were associated with the highly prolific phenotype (p FecX (Gr) = 5.98E(-06) and p FecX(O) = 2.55E(-08)). Homozygous ewes for the mutated allele showed a significantly increased prolificacy (FecX(Gr)/FecX(Gr), LS = 2.50 ± 0.65 versus FecX(+)/FecX(Gr), LS = 1.93 ± 0.42, p<1E(-03) and FecX(O)/FecX(O), OR = 3.28 ± 0.85 versus FecX(+)/FecX(O), OR = 2.02 ± 0.47, p<1E(-03)). Both mutations are located in very well conserved motifs of the protein and altered the BMP15 signaling activity in vitro using a BMP-responsive luciferase test in COV434 granulosa cells. Thus, we have identified two novel mutations in the BMP15 gene associated with increased LS and OR. Notably, homozygous FecX(Gr)/FecX(Gr) Grivette and homozygous FecX(O)/FecX(O) Olkuska ewes are hyperprolific in striking contrast with the sterility exhibited by all other known homozygous BMP15 mutations. Our results bring new insights into the key role played by the BMP15 protein in ovarian function and could contribute to a better understanding of the pathogenesis of women's fertility disorders.

Genetic control of multiple births in low ovulating mammalian species.

In mammals, litter size is a highly variable trait. Some species such as humans or cattle are monotocous, with one or sometimes two newborns per birth, whereas others, the polytocous species such as mice or pigs, are highly prolific and often produce a dozen newborns at each farrowing. In monotocous species, however, two or three newborns per birth may sometime be unwanted. In more polytocous species such as sheep or pigs, litter size is studied in order to increase livestock prolificacy. By contrast, twinning rates in humans or cattle may increase birth difficulties and health problems in the newborns. In this context, the aim of our review was to provide a clearer understanding of the genetic and physiological factors that control multiple births in low-ovulating mammalian species, with particular focus on three species: sheep, cattle, and humans, where knowledge of the ovulation rate in one may enlighten findings in the others. This article therefore reviews the phenotypic and genetic variability observed with respect to ovulation and twinning rates. It then presents the QTL and major genes that have been identified in each species. Finally, we draw a picture of the diversity of the physiological mechanisms underlying multiple ovulation. Although several major genes have been discovered in sheep, QTL detection methods in humans or cattle have suggested that the determinism of litter size is complex and probably involves several genes in order to explain variations in the number of ovulations.

Endocrine characterization of the reproductive axis in highly prolific lacaune sheep homozygous for the FecLL mutation.

A prolific allele named FecL(L) is known to segregate in the meat breed of the French Lacaune sheep on the basis of ovulation rate record. To gain more knowledge about the physiological effects of FecL(L), ewes homozygous for FecL(L) (L/L) were compared to wild-type ewes (+/+) for ovarian phenotype and reproductive endocrine profiles. At the ovarian level, the increased ovulation rate in L/L ewes was associated with an increased number of antral follicles of greater than 3 mm and with preovulatory follicles being, on average, 1 mm smaller. Intrafollicular estradiol and testosterone concentrations were not significantly different between the two genotypes. In contrast, L/L large follicles (>or=6 mm) had lower intrafollicular progesterone concentration. At the molecular level, expressions of ovarian markers, such as CYP19A1, CYP11A1, CYP17A1, LHR, and INHA, were not significantly different between the two genotypes. In contrast, FSHR and STAR mRNA levels increased in granulosa cells from L/L ewes. Plasma concentrations of estradiol, luteinizing hormone (LH), follicle-stimulating hormone (FSH), and progesterone measured across a synchronized estrous cycle revealed a significant increase in estradiol levels during the follicular phase, a precocious LH surge, and an increase in progesterone level during the luteal phase of L/L ewes compared to +/+ ewes. Circulating concentrations of FSH were not different between the two genotypes. The precocious LH surge was associated with an increase in frequency of LH pulsatility during the follicular phase. At the pituitary level, mRNA levels for LHB, FSHB, GNRHR, and ESR1 were not significantly different between the two genotypes. In contrast, ESR2 mRNA expression was decreased in L/L ewes compared to +/+ ewes. Based on ovarian phenotype and endocrine profiles, these findings suggest that the mutation in the FecL gene affects ovarian function in a different way compared to other known prolific mutations affecting the bone morphogenetic protein signaling system in the ovine species.

A novel mutation in the bone morphogenetic protein 15 gene causing defective protein secretion is associated with both increased ovulation rate and sterility in Lacaune sheep.

Genetic mutations with major effects on ovulation rate and litter size in sheep were recently identified in three genes belonging to the TGFbeta superfamily pathway: the bone morphogenetic protein 15 (BMP15, also known as GDF9b), growth differentiation factor 9 (GDF9), and BMP receptor type IB (also known as activin-like kinase 6). Homozygous BMP15 or GDF9 mutations raise female sterility due to a failure of normal ovarian follicle development, whereas heterozygous animals for BMP15 or GDF9 as well as heterozygous and homozygous animals for BMP receptor type IB show increased ovulation rates. In the present work, a new naturally occurring mutation in the BMP15 gene in the high prolific Lacaune sheep breed is described. The identified variant is a C53Y missense nonconservative substitution leading to the aminoacidic change of a cysteine with a tyrosine in the mature peptide of the protein. As for other mutations found in the same gene, this is associated with an increased ovulation rate and sterility in heterozygous and homozygous animals, respectively. Further in vitro studies showed that the C53Y mutation was responsible for the impairment of the maturation process of the BMP15 protein, resulting in a defective secretion of both the precursor and mature peptide. Overall, our findings confirm the essential role of the BMP15 factor in the ovarian folliculogenesis and control of ovulation rate in sheep.

Regulation of ovulation rate in mammals: contribution of sheep genetic models.

Ovarian folliculogenesis in mammals from the constitution of primordial follicles up to ovulation is a reasonably well understood mechanism. Nevertheless, underlying mechanisms that determine the number of ovulating follicles were enigmatic until the identification of the fecundity genes affecting ovulation rate in sheep, bone morphogenetic protein-15 (BMP-15), growth and differentiation factor-9 (GDF-9) and BMP receptor-1B (BMPR-1B). In this review, we focus on the use of these sheep genetic models for understanding the role of the BMP system as an intra-ovarian regulator of follicular growth and maturation, and finally, ovulation rate.

Mutations in the genes for oocyte-derived growth factors GDF9 and BMP15 are associated with both increased ovulation rate and sterility in Cambridge and Belclare sheep (Ovis aries).

Belclare and Cambridge are prolific sheep breeds, the origins of which involved selecting ewes with exceptionally high litter size records from commercial flocks. The variation in ovulation rate in both breeds is consistent with segregation of a gene (or genes) with a large effect on this trait. Sterile ewes, due to a failure of normal ovarian follicle development, occur in both breeds. New naturally occurring mutations in genes for the oocyte-derived growth factors growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP15) are described. These mutations are associated with increased ovulation rate in heterozygous carriers and sterility in homozygous carriers in both breeds. This is the first time that a mutation in the gene for GDF9 has been found that causes increased ovulation rate and infertility in a manner similar to inactivating mutations in BMP15, and shows that GDF9 is essential for normal folliculogenesis in sheep. Furthermore, it is shown, for the first time in any species, that individuals with mutations in both GDF9 and BMP15 have a greater ovulation rate than sheep with either of the mutations separately.

Segregation of a major gene influencing ovulation in progeny of Lacaune meat sheep.

Inheritance of the ovulation rate (OR) in the Lacaune meat breed was studied through records from a small nucleus of 36 hyper-prolific ewes screened on farms on the basis of their natural litter size, and from progeny data of three selected Lacaune sires. These sires were chosen at the AI centre according to their breeding values estimated for the mean and the variability of their daughters' litter size. Non-carrier Lacaune dairy ewes were inseminated to produce 121 F1 daughters and 27 F1 sons. Twelve sons (four from each sire) were used in turn to inseminate non-carrier Lacaune dairy ewes providing 260 BC progeny ewes. F1 and BC progeny were brought from private farms and gathered after weaning on an experimental farm where ovulation rates were recorded in the first and second breeding seasons. With an average of 6.5 records each, the mean OR of hyper-prolific ewes was very high (5.34), and 38.4% of records showed a rate of 6 or more. F1 data showed high repeatability of OR (r = 0.54) within ewe, with significant variability among ewes. High OR (> or = 4) were observed in each family. A segregation analysis provided a significant likelihood ratio and classified the three founders as heterozygous. BC ewes also displayed high repeatability of OR (r = 0.47) and the mean OR varied considerably between families (from 1.24 to 1.78). Seven of the 12 BC families presented high-ovulating ewes (at least one record > or = 4) and segregation analysis yielded a highly significant likelihood ratio as compared to an empirical test distribution. The high variability of the mean ovulation rate shown by a small group of daughters of BC ewes inseminated by putative carrier F1 rams, and the very high ovulation rate observed for some of these ewe lambs, confirmed the segregation of a major gene with two co-dominant alleles borne by an autosome. The difference between homozygous non-carriers and heterozygous ewes was about one ovulation on the observed scale and 2.2 standard deviations on the underlying scale.

Regulation of ovarian folliculogenesis by IGF and BMP system in domestic animals.

Involvement of insulin-like growth factors (IGFs) and IGF binding proteins (IGFBPs) in ovarian folliculogenesis has been extensively studied during the last decade. In all mammalian species, IGF-I stimulates granulosa cell proliferation and steroidogenesis. The concentrations of IGF-I and -II do not vary during terminal follicular growth and atresia. In contrast, the levels of IGFBP-2 and -4, as well as IGFBP-5 in ruminants, dramatically decrease and increase during terminal follicular growth and atresia, respectively. These changes are responsible for an increase and a decrease in IGF bioavailability during follicular growth and atresia, respectively. They are partly explained by changes in ovarian expression. In particular, expression of IGFBP-2 mRNA decreases during follicular growth in ovine, bovine and porcine ovaries, and expression of IGFBP-5 mRNA dramatically increases in granulosa cells of bovine and ovine atretic follicles. Changes in IGFBP-2 and -4 levels are also due to changes in intrafollicular levels of specific proteases. Recently, we have shown that the pregnancy-associated plasma protein-A (PAPP-A) is responsible for the degradation of IGFBP-4 in preovulatory follicles of domestic animals. Expression of PAPP-A mRNA is restricted to the granulosa cell compartment, and is positively correlated to expression of aromatase and LH receptor. From recent evidence, the bone morphogenetic protein (BMP) family would also play a key role in ovarian physiology of domestic animals. In particular, we and others have recently shown that a non-conservative substitution (Q249R) in the bone morphogenetic protein-receptor type IB (BMPR-IB) coding sequence is fully associated with the hyperprolific phenotype of FecB(B)/FecB(B) Booroola ewes. BMP-4 and GDF-5, natural ligands of BMPR-IB, strongly inhibit secretion of progesterone by ovine granulosa cells in vitro, but granulosa cells from FecB(B)/FecB(B) ewes are less responsive than those from FecB(+)/FecB(+) to the action of these peptides. It is suggested that in FecB(B)/FecB(B) ewes, Q249R substitution would impair the function of BMPR-IB, leading to a precocious differentiation of granulosa cells and of follicular maturation. Interestingly, recent findings have described mutations in BMP-15 gene associated with hyperprolific phenotypes in Inverdale and Hanna ewes, suggesting that the BMP pathway plays a crucial role in the control of ovulation rate.

DNA tests in prolific sheep from eight countries provide new evidence on origin of the Booroola (FecB) mutation.

Recent discoveries that high prolificacy in sheep carrying the Booroola gene (FecB) is the result of a mutation in the BMPIB receptor and high prolificacy in Inverdale sheep (FecX(I)) is the result of a mutation in the BMP15 oocyte-derived growth factor gene have allowed direct marker tests to be developed for FecB and FecX(I). These tests were carried out in seven strains of sheep (Javanese, Thoka, Woodlands, Olkuska, Lacaune, Belclare, and Cambridge) in which inheritance patterns have suggested the presence of major genes affecting prolificacy and in the prolific Garole sheep of India, which have been proposed as the ancestor of Australian Booroola Merinos. The FecB mutation was found in the Garole and Javanese sheep but not in Thoka, Woodlands, Olkuska, Lacaune, Belclare, and Cambridge sheep. None of the sheep tested had the FecX(I) mutation. These findings present strong evidence to support historical records that the Booroola gene was introduced into Australian flocks from Garole (Bengal) sheep in the late 18th century. It is unknown whether Javanese Thin-tailed sheep acquired the Booroola gene directly from Garole sheep from India or via Merinos from Australia. The DNA mutation test for FecB will enable breeding plans to be developed that allow the most effective use of this gene in Garole and Javanese Thin-tailed sheep and their crosses.