Emmanuelle Rebours

Identification of large intergenic non-coding RNAs in bovine muscle using next-generation transcriptomic sequencing

BackgroundThe advent of large-scale gene expression technologies has helped to reveal in eukaryotic cells, the existence of thousands of non-coding transcripts, whose function and significance remain mostly poorly understood. Among these non-coding transcripts, long non-coding RNAs (lncRNAs) are the least well-studied but are emerging as key regulators of diverse cellular processes. In the present study, we performed a survey in bovine Longissimus thoraci of lincRNAs (long intergenic non-coding RNAs not overlapping protein-coding transcripts). To our knowledge, this represents the first such study in bovine muscle.ResultsTo identify lincRNAs, we used paired-end RNA sequencing (RNA-Seq) to explore the transcriptomes of Longissimus thoraci from nine Limousin bull calves. Approximately 14–45 million paired-end reads were obtained per library. A total of 30,548 different transcripts were identified. Using a computational pipeline, we defined a stringent set of 584 different lincRNAs with 418 lincRNAs found in all nine muscle samples. Bovine lincRNAs share characteristics seen in their mammalian counterparts: relatively short transcript and gene lengths, low exon number and significantly lower expression, compared to protein-encoding genes. As for the first time, our study identified lincRNAs from nine different samples from the same tissue, it is possible to analyse the inter-individual variability of the gene expression level of the identified lincRNAs. Interestingly, there was a significant difference when we compared the expression variation of the 418 lincRNAs with the 10,775 known selected protein-encoding genes found in all muscle samples. In addition, we found 2,083 pairs of lincRNA/protein-encoding genes showing a highly significant correlated expression. Fourteen lincRNAs were selected and 13 were validated by RT-PCR. Some of the lincRNAs expressed in muscle are located within quantitative trait loci for meat quality traits.ConclusionsOur study provides a glimpse into the lincRNA content of bovine muscle and will facilitate future experimental studies to unravel the function of these molecules. It may prove useful to elucidate their effect on mechanisms underlying the genetic variability of meat quality traits. This catalog will complement the list of lincRNAs already discovered in cattle and therefore will help to better annotate the bovine genome.

Validation of RNA isolated from milk fat globules to profile mammary epithelial cell expression during lactation and transcriptional response to a bacterial infection

Mastitis, an inflammation of the mammary gland, is the most costly infectious disease of dairy ruminants worldwide. Although it receives considerable attention, the early steps of the host response remain poorly defined. Here, we report a noninvasive method using milk fat globules (MFG) as a source of mammary RNA to follow the dynamics of the global transcriptional response of mammary epithelial cells (MEC) during the course of a bacterial infection. We first assessed that RNA isolated from MFG were representative of MEC RNA; we then evaluated whether MFG RNA could be used to monitor the MEC response to infection. Sufficiently high yields of good-quality RNA (RNA integrity numbers ranging between 6.7 and 8.7) were obtained from goat MFG for subsequent analyses. Contamination of MFG by macrophages and neutrophils, which can be trapped during creaming, was assessed and when using quantitative real-time PCR for cell-type specific markers, was shown to be weak enough (<8%) to affect MFG gene expression profiling. Using microarrays, we showed that RNA extracted from MFG and from mammary alveolar parenchyma shared approximately 90% of the highlighted probes corresponding in particular to genes encoding milk proteins (CSN, BLG, LALBA) and enzymes involved in milk fat synthesis and secretion (FASN, XDH, ADRP, SCD, and DGAT1). In addition, a gene involved in the acute-phase reaction, coding for the serum amyloid A3 (SAA3) protein, was found within the first 50 most highly expressed genes in a noninfectious context in both mammary alveolar parenchyma and MFG, strongly suggesting that SAA3 is expressed in MEC. We took advantage of this noninvasive RNA sampling to follow the early proinflammatory response of MEC during the course of a bacterial infection and showed that the levels of mRNA encoding SAA3 sharply increased at 24h postinfection. Taken together, our results demonstrate that MFG represent a unique source of MEC RNA to noninvasively sample sufficient amounts of high-quality RNA to assess the dynamics of MEC gene expression in vivo, especially during the first steps of infection, thereby paving the way for the discovery of early biomarkers for the control of intramammary infections. Furthermore, this noninvasive technique could be used to provide mammary transcriptomic data on a large scale, thus filling the gap between genomic and phenotypic data.

Identification of major milk fat globule membrane proteins from pony mare milk highlights the molecular diversity of lactadherin across species

Although several studies have been devoted to the colloidal and soluble protein fractions of mare milk (caseins and whey proteins), to date little is known about the milk fat globule membrane (MFGM) protein fraction from mare milk. The objective of this study was thus to describe MFGM proteins from Equidae milk and to compare those proteins to already described MFGM proteins from cow and goat milk. Major MFGM proteins (namely, xanthine oxidase, butyrophilin, lactadherin, and adipophilin) already described in cow or goat milk were identified in mare milk using mass spectrometry. However, species-specific peculiarities were observed for 2 MFGM proteins: butyrophilin and lactadherin. A highly glycosylated 70-kDa protein was characterized for equine butyrophilin, whereas proteins of 64 and 67 kDa were characterized for cow and goat butyrophilin, respectively. Prominent differences across species were highlighted for lactadherin. Indeed, whereas 1 or 2 polypeptide chains were identified, respectively, by peptide mass fingerprinting matrix-assisted laser desorption/ionization-time of flight analysis for caprine and bovine lactadherin, 4 isoforms (60, 57, 48, and 45 kDa) for lactadherin from mare milk were identified by 10% sodium dodecyl sulfate-PAGE. Polymerase chain reaction experiments on lactadherin transcripts isolated from milk fat globules revealed the existence of 2 distinct lactadherin transcripts in the horse mammary gland. Cloning and sequencing of both transcripts encoding lactadherin showed an alternative use of a cryptic splice site located at the end of intron 5 of the equine lactadherin-encoding gene. This event results in the occurrence of an additional alanine (A) residue in the protein that disrupts a putative atypical N-glycosylation site (VNGC/VNAGC) described in human lactadherin. Liquid chromatography coupled with tandem mass spectrometry analyses confirmed the existence of both lactadherin variants in mare MFGM. We show here that lactadherin from Equidae milk is much more complex than that from Bovidae milk (i.e., cow and goat milk), therefore raising questions regarding the precise function of these different isoforms, if any, in the equine mammary gland.

Insights into gene expression profiling of natural resistance to coccidiosis in contrasting chicken lines

Coccidiosis is a parasitic disease with major economic impact, one of whose main causative agents is Eimeria tenella. Chicken breeds display variable natural resistance to this disease. Unravelling the genetic bases of such variations could provide new clues for protection strategies. Transcriptomic experiments were conducted comparing resistant (Fayoumi) and susceptible (Leghorn) lines. Caecum and caecal tonsils were analysed. A global increase in differential gene expression following infection was observed for caecum comparisons, whereas a global decrease following infection was observed for caecal tonsils.Gene lists for infected tissues display 40 genes in common across breeds, 20 of which were specific to infected tissues. Among these specific genes, 9 belong to the 100 more differentially expressed genes of the infected caecum comparison. Gene expression networks were constructed in parallel, identifying highly connected genes. Comparing information from differential gene lists and gene network analysis allows one to highlight potential pivotal genes in the infection process, one of which was located in a putative significant QTL region for infection associated lesions.

Construction and validation of a novel dual reporter vector for studying mammalian bidirectional promoters.

Regulation of gene expression plays important role in cellular functions. With the development of sequencing techniques, more and more genomes are available and genome-wide analyses of genomic structures that may affect gene expression regulation are now possible. Analyses of several genomes have found a class of regulatory regions that contain elements that initiate transcription of two different genes positioned with a head-to-head arrangement in two opposite directions. These regulatory regions are known as bidirectional promoters. Although bidirectional promoters have been known for years, recent genome-scale studies have shown that the regulation of the expression of up to 10% of the genes are controlled by bidirectional promoters. These findings are based mostly on computational work and only a limited number of putative bidirectional promoters have been experimentally validated. Developing methods to study bidirectional promoters will allow researchers to understand how these regions are regulated and the roles that divergent transcription plays in the expression of genes. Here, we have developed a novel dual-fluorescence reporter gene vector to study the transcriptional output of mammalian bidirectional promoters. We demonstrate that this vector is capable of expressing reporter genes under the control of bidirectional promoters, using the known human OSGEP/APEX bidirectional promoter.