Frédéric Lecerf

Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Mérino ewes

Ewes from the Booroola strain of Australian Mérino sheep are characterized by high ovulation rate and litter size. This phenotype is due to the action of the FecBB allele of a major gene named FecB, as determined by statistical analysis of phenotypic data. By genetic analysis of 31 informative half-sib families from heterozygous sires, we showed that the FecB locus is situated in the region of ovine chromosome 6 corresponding to the human chromosome 4q22–23 that contains the bone morphogenetic protein receptor IB (BMPR-IB) gene encoding a member of the transforming growth factor-β (TGF-β) receptor family. A nonconservative substitution (Q249R) in the BMPR-IB coding sequence was found to be associated fully with the hyperprolificacy phenotype of Booroola ewes. In vitro, ovarian granulosa cells from FecBB/FecBB ewes were less responsive than granulosa cells from FecB+/FecB+ ewes to the inhibitory effect on steroidogenesis of GDF-5 and BMP-4, natural ligands of BMPR-IB. It is suggested that in FecBB/FecBB ewes, BMPR-IB would be inactivated partially, leading to an advanced differentiation of granulosa cells and an advanced maturation of ovulatory follicles.

DNA Tests in Prolific Sheep from Eight Countries Provide New Evidence on Origin of the Booroola (FecB) Mutation

Abstract 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 (FecXI) 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 FecXI. 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 FecXI 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.

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.

Transcriptome profiling of the feeding-to-fasting transition in chicken liver

BackgroundStarvation triggers a complex array of adaptative metabolic responses including energy-metabolic responses, a process which must imply tissue specific alterations in gene expression and in which the liver plays a central role. The present study aimed to describe the evolution of global gene expression profiles in liver of 4-week-old male chickens during a 48 h fasting period using a chicken 20 K oligoarray.ResultsA large number of genes were modulated by fasting (3532 genes with a pvalue corrected by Benjamini-Hochberg < 0.01); 2062 showed an amplitude of variation higher than +/- 40% among those, 1162 presented an human ortholog, allowing to collect functional information. Notably more genes were down-regulated than up-regulated, whatever the duration of fasting (16 h or 48 h). The number of genes differentially expressed after 48 h of fasting was 3.5-fold higher than after 16 h of fasting. Four clusters of co-expressed genes were identified by a hierarchical cluster analysis. Gene Ontology, KEGG and Ingenuity databases were then used to identify the metabolic processes associated to each cluster. After 16 h of fasting, genes involved in ketogenesis, gluconeogenesis and mitochondrial or peroxisomal fatty acid beta-oxidation, were up-regulated (cluster-1) whereas genes involved in fatty acid and cholesterol synthesis were down-regulated (cluster-2). For all genes tested, the microarray data was confirmed by quantitative RT-PCR. Most genes were altered by fasting as already reported in mammals. A notable exception was the HMG-CoA synthase 1 gene, which was up-regulated following 16 and 48 h of fasting while the other genes involved in cholesterol metabolism were down-regulated as reported in mammalian studies. We further focused on genes not represented on the microarray and candidates for the regulation of the target genes belonging to cluster-1 and -2 and involved in lipid metabolism. Data are provided concerning PPARa, SREBP1, SREBP2, NR1H3 transcription factors and two desaturases (FADS1, FADS2).ConclusionThis study evidences numerous genes altered by starvation in chickens and suggests a global repression of cellular activity in response to this stressor. The central role of lipid and acetyl-CoA metabolisms and its regulation at transcriptional level are confirmed in chicken liver in response to short-term fasting. Interesting expression modulations were observed for NR1H3, FADS1 and FADS2 genes. Further studies are needed to precise their role in the complex regulatory network controlling lipid metabolism.

QTL for resistance to Salmonella carrier state confirmed in both experimental and commercial chicken lines

The ability of chickens to carry Salmonella without displaying disease symptoms is responsible for Salmonella propagation in poultry stocks and for subsequent human contamination through the consumption of contaminated eggs or meat. The selection of animals more resistant to carrier state might be a way to decrease the propagation of Salmonella in poultry stocks and its transmission to humans. Five QTL controlling variation for resistance to carrier state in a chicken F(2) progeny derived from the White Leghorn inbred lines N and 6(1) had been previously identified using a selective genotyping approach. Here, a second analysis on the whole progeny was performed, which led to the confirmation of two QTL on chromosomes 2 and 16. To assess the utility of these genomic regions for selection in commercial lines, we tested them together with other QTL identified in an [Nx6(1)] x N backcross progeny and with the candidate genes SLC11A1 and TLR4. We used a commercial line divergently selected for either low or high carrier-state resistance both in young chicks and in adult hens. In divergent chick lines, one QTL on chromosome 1 and one in the SLC11A1 region were significantly associated with carrier-state resistance variations; in divergent adult lines, one QTL located in the major histocompatibility complex on chromosome 16 and one in the SLC11A1 region were involved in these variations. Genetic studies conducted on experimental lines can therefore be of potential interest for marker-assisted selection in commercial lines.

The duplicated genes database: identification and functional annotation of co-localised duplicated genes across genomes.

Background There has been a surge in studies linking genome structure and gene expression, with special focus on duplicated genes. Although initially duplicated from the same sequence, duplicated genes can diverge strongly over evolution and take on different functions or regulated expression. However, information on the function and expression of duplicated genes remains sparse. Identifying groups of duplicated genes in different genomes and characterizing their expression and function would therefore be of great interest to the research community. The ‘Duplicated Genes Database’ (DGD) was developed for this purpose. Methodology Nine species were included in the DGD. For each species, BLAST analyses were conducted on peptide sequences corresponding to the genes mapped on a same chromosome. Groups of duplicated genes were defined based on these pairwise BLAST comparisons and the genomic location of the genes. For each group, Pearson correlations between gene expression data and semantic similarities between functional GO annotations were also computed when the relevant information was available. Conclusions The Duplicated Gene Database provides a list of co-localised and duplicated genes for several species with the available gene co-expression level and semantic similarity value of functional annotation. Adding these data to the groups of duplicated genes provides biological information that can prove useful to gene expression analyses. The Duplicated Gene Database can be freely accessed through the DGD website at http://dgd.genouest.org.

Limited RNA Editing in Exons of Mouse Liver and Adipose

Several studies have investigated RNA–DNA differences (RDD), presumably due to RNA editing, with conflicting results. We report a rigorous analysis of RDD in exonic regions in mice, taking into account critical biases in RNA-Seq analysis. Using deep-sequenced F1 reciprocal inbred mice, we mapped 40 million RNA-Seq reads per liver sample and 180 million reads per adipose sample. We found 7300 apparent hepatic RDDs using a multiple-site mapping procedure, compared with 293 RDD found using a unique-site mapping procedure. After filtering for repeat sequence, splice junction proximity, undirectional strand, and extremity read bias, 63 RDD remained. In adipose tissue unique-site mapping identified 1667 RDD, and after applying the same four filters, 188 RDDs remained. In both tissues, the filtering procedure increased the proportion of canonical (A-to-I and C-to-U) editing events. The genomic DNA of 12 RDD sites among the potential 63 hepatic RDD was tested by Sanger sequencing, three of which proved to be due to unreferenced SNPs. We validated seven liver RDD with Sequenom technology, including two noncanonical, Gm5424 C-to-I(G) and Pisd I(G)-to-A RDD. Differences in diet, sex, or genetic background had very modest effects on RDD occurrence. Only a small number of apparent RDD sites overlapped between liver and adipose, indicating a high degree of tissue specificity. Our findings underscore the importance of properly filtering for bias in RNA-Seq investigations, including the necessity of confirming the DNA sequence to eliminate unreferenced SNPs. Based on our results, we conclude that RNA editing is likely limited to hundreds of events in exonic RNA in liver and adipose.

AnnotQTL: a new tool to gather functional and comparative information on a genomic region

AnnotQTL is a web tool designed to aggregate functional annotations from different prominent web sites by minimizing the redundancy of information. Although thousands of QTL regions have been identified in livestock species, most of them are large and contain many genes. This tool was therefore designed to assist the characterization of genes in a QTL interval region as a step towards selecting the best candidate genes. It localizes the gene to a specific region (using NCBI and Ensembl data) and adds the functional annotations available from other databases (Gene Ontology, Mammalian Phenotype, HGNC and Pubmed). Both human genome and mouse genome can be aligned with the studied region to detect synteny and segment conservation, which is useful for running inter-species comparisons of QTL locations. Finally, custom marker lists can be included in the results display to select the genes that are closest to your most significant markers. We use examples to demonstrate that in just a couple of hours, AnnotQTL is able to identify all the genes located in regions identified by a full genome scan, with some highlighted based on both location and function, thus considerably increasing the chances of finding good candidate genes. AnnotQTL is available at http://annotqtl.genouest.org.

MarkerSet: a marker selection tool based on markers location and informativity in experimental designs

BackgroundThe recent sequencing of full genomes has led to the availability of many SNP markers which are very useful for the mapping of complex traits. In livestock production, there are still no commercial arrays and many studies use home-made sets of SNPs. Thus, the current methodologies for SNP genotyping are still expensive and it is a crucial step to select the SNPs to use. Indeed, the main factors affecting the power of the linkage analyses are the density of the genetic map and the heterozygosity of markers in tested animal parents.FindingsThis is why we have developed a PERL program selecting a defined number of markers based on their locations on the genome and their informativity in specific experimental designs. As an option, different experimental designs can be combined in order to select the best possible common marker set. The program has been tested using different conditions of marker informativity and density with both real and simulated datasets. The results show the efficiency of our program to select the most informative markers even if there is a wide range of informativity for whole genome scan mapping analyses. In case of combination of different experimental crosses, the multidesign mode can optimize the SNP markers selection.ConclusionWritten in PERL, it assures a maximum portability to other operating systems (OS) and the source code availability for user modifications. Except for the simulation mode which could be time consuming, MarkerSet can compute results in a very short time.

Combined QTL and Selective Sweep Mappings with Coding SNP Annotation and cis-eQTL Analysis Revealed PARK2 and JAG2 as New Candidate Genes for Adiposity Regulation

Very few causal genes have been identified by quantitative trait loci (QTL) mapping because of the large size of QTL, and most of them were identified thanks to functional links already known with the targeted phenotype. Here, we propose to combine selection signature detection, coding SNP annotation, and cis-expression QTL analyses to identify potential causal genes underlying QTL identified in divergent line designs. As a model, we chose experimental chicken lines divergently selected for only one trait, the abdominal fat weight, in which several QTL were previously mapped. Using new haplotype-based statistics exploiting the very high SNP density generated through whole-genome resequencing, we found 129 significant selective sweeps. Most of the QTL colocalized with at least one sweep, which markedly narrowed candidate region size. Some of those sweeps contained only one gene, therefore making them strong positional causal candidates with no presupposed function. We then focused on two of these QTL/sweeps. The absence of nonsynonymous SNPs in their coding regions strongly suggests the existence of causal mutations acting in cis on their expression, confirmed by cis-eQTL identification using either allele-specific expression or genetic mapping analyses. Additional expression analyses of those two genes in the chicken and mice contrasted for adiposity reinforces their link with this phenotype. This study shows for the first time the interest of combining selective sweeps mapping, coding SNP annotation and cis-eQTL analyses for identifying causative genes for a complex trait, in the context of divergent lines selected for this specific trait. Moreover, it highlights two genes, JAG2 and PARK2, as new potential negative and positive key regulators of adiposity in chicken and mice.