Comparative analysis of transcriptome remodeling in plaque-associated and plaque-distant microglia during amyloid-ß pathology progression in mice.

BACKGROUND: Research in recent years firmly established that microglial cells play an important role in the pathogenesis of Alzheimer's disease (AD). In parallel, a series of studies showed that, under both homeostatic and pathological conditions, microglia are a heterogeneous cell population. In AD, amyloid-ß (Aß) plaque-associated microglia (PAM) display a clearly distinct phenotype compared to plaque-distant microglia (PCM), suggesting that these two microglia subtypes likely differently contribute to disease progression. So far, molecular characterization of PAM was performed indirectly using single cell RNA sequencing (scRNA-seq) approaches or based on markers that are supposedly up-regulated in this microglia subpopulation. METHODS: In this study based on a well-characterized AD mouse model, we combined cell-specific laser capture microdissection and RNA-seq analysis to i) identify, without preconceived notions of the molecular and/or functional changes that would affect these cells, the genes and gene networks that are dysregulated in PAM or PCM at three critical stages of the disease, and ii) to investigate the potential contribution of both plaque-associated and plaque-distant microglia. RESULTS: First, we established that our approach allows selective isolation of microglia, while preserving spatial information and preventing transcriptome changes induced by classical purification approaches. Then, we identified, in PAM and PCM subpopulations, networks of co-deregulated genes and analyzed their potential functional roles in AD. Finally, we investigated the dynamics of microglia transcriptomic remodeling at early, intermediate and late stages of the disease and validated select findings in postmortem human AD brain. CONCLUSIONS: Our comprehensive study provides useful transcriptomic information regarding the respective contribution of PAM and PCM across the Aß pathology progression. It highlights specific pathways that would require further study to decipher their roles across disease progression. It demonstrates that the proximity of microglia to Aß-plaques dramatically alters the microglial transcriptome and reveals that these changes can have both positive and negative impacts on the surrounding cells. These opposing effects may be driven by local microglia heterogeneity also demonstrated by this study. Our approach leads to molecularly define the less well studied plaque-distant microglia. We show that plaque-distant microglia are not bystanders of the disease, although the transcriptomic changes are far less striking compared to what is observed in plaque-associated microglia. In particular, our results suggest they may be involved in Aß oligomer detection and in Aß-plaque initiation, with increased contribution as the disease progresses.

Survey of allele specific expression in bovine muscle.

Allelic imbalance is a common phenomenon in mammals that plays an important role in gene regulation. An Allele Specific Expression (ASE) approach can be used to detect variants with a cis-regulatory effect on gene expression. In cattle, this type of study has only been done once in Holstein. In our study we performed a genome-wide analysis of ASE in 19 Limousine muscle samples. We identified 5,658 ASE SNPs (Single Nucleotide Polymorphisms showing allele specific expression) in 13% of genes with detectable expression in the Longissimus thoraci muscle. Interestingly we found allelic imbalance in AOX1, PALLD and CAST genes. We also found 2,107 ASE SNPs located within genomic regions associated with meat or carcass traits. In order to identify causative cis-regulatory variants explaining ASE we searched for SNPs altering binding sites of transcription factors or microRNAs. We identified one SNP in the 3'UTR region of PRNP that could be a causal regulatory variant modifying binding sites of several miRNAs. We showed that ASE is frequent within our muscle samples. Our data could be used to elucidate the molecular mechanisms underlying gene expression imbalance.

Identification of copy number variation in French dairy and beef breeds using next-generation sequencing.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Genetic variability of the activity of bidirectional promoters: a pilot study in bovine muscle.

Bidirectional promoters are regulatory regions co-regulating the expression of two neighbouring genes organized in a head-to-head orientation. In recent years, these regulatory regions have been studied in many organisms; however, no investigation to date has been done to analyse the genetic variation of the activity of this type of promoter regions. In our study, we conducted an investigation to first identify bidirectional promoters sharing genes expressed in bovine Longissimus thoracis and then to find genetic variants affecting the activity of some of these bidirectional promoters. Combining bovine gene information and expression data obtained using RNA-Seq, we identified 120 putative bidirectional promoters active in bovine muscle. We experimentally validated in vitro 16 of these bidirectional promoters. Finally, using gene expression and whole-genome genotyping data, we explored the variability of the activity in muscle of the identified bidirectional promoters and discovered genetic variants affecting their activity. We found that the expression level of 77 genes is correlated with the activity of 12 bidirectional promoters. We also identified 57 single nucleotide polymorphisms associated with the activity of 5 bidirectional promoters. To our knowledge, our study is the first analysis in any species of the genetic variability of the activity of bidirectional promoters.

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

BACKGROUND: The 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. RESULTS: To 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. CONCLUSIONS: Our 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.

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.

Bovine TWINKLE and mitochondrial ribosomal protein L43 genes are regulated by an evolutionary conserved bidirectional promoter.

TWINKLE is a mitochondrial DNA helicase playing an important role in mitochondrial DNA replication. In human, mutations in this gene cause progressive external ophtalmoplegia and mitochondrial DNA depletion syndrome-7. TWINKLE is well conserved among multicellular eukaryotes and is believed to be a key regulator of mitochondrial DNA copy number in mammals. Despite its involvement in several diseases and its important function in mitochondrial DNA metabolism, nothing is known about the regulation of the expression of TWINKLE. We have analysed the 5'-flanking genomic region of the bovine TWINKLE gene and found it was localised adjacent to the MRPL43 gene in a head-to-head orientation, suggesting that both genes are regulated by a shared bidirectional promoter. The bovine 75-bp long intergenic region shows substantial homology across different species and contains several conserved putative transcription factor binding sites. A TATA box, however, was lacking. Using a dual fluorescent reporter system and transient transfection assays, we have analysed the bovine intergenic region between TWINKLE and MRPL43. This small genomic fragment showed a bidirectional promoter activity. As the TWINKLE/MRPL43 bidirectional promoter tested was highly conserved, it is likely that the results we obtained here in cattle may be extended to the other species.

Gene-based single nucleotide polymorphism discovery in bovine muscle using next-generation transcriptomic sequencing.

BACKGROUND: Genetic information based on molecular markers has increasingly being used in cattle breeding improvement programmes, as a mean to improve conventionally phenotypic selection. Advances in molecular genetics have led to the identification of several genetic markers associated with genes affecting economic traits. Until recently, the identification of the causative genetic variants involved in the phenotypes of interest has remained a difficult task. The advent of novel sequencing technologies now offers a new opportunity for the identification of such variants. Despite sequencing costs plummeting, sequencing whole-genomes or large targeted regions is still too expensive for most laboratories. A transcriptomic-based sequencing approach offers a cheaper alternative to identify a large number of polymorphisms and possibly to discover causative variants. In the present study, we performed a gene-based single nucleotide polymorphism (SNP) discovery analysis in bovine Longissimus thoraci, using RNA-Seq. To our knowledge, this represents the first study done in bovine muscle. RESULTS: Messenger RNAs from Longissimus thoraci from three Limousin bull calves were subjected to high-throughput sequencing. Approximately 36-46 million paired-end reads were obtained per library. A total of 19,752 transcripts were identified and 34,376 different SNPs were detected. Fifty-five percent of the SNPs were found in coding regions and ~22% resulted in an amino acid change. Applying a very stringent SNP quality threshold, we detected 8,407 different high-confidence SNPs, 18% of which are non synonymous coding SNPs. To analyse the accuracy of RNA-Seq technology for SNP detection, 48 SNPs were selected for validation by genotyping. No discrepancies were observed when using the highest SNP probability threshold. To test the usefulness of the identified SNPs, the 48 selected SNPs were assessed by genotyping 93 bovine samples, representing mostly the nine major breeds used in France. Principal component analysis indicates a clear separation between the nine populations. CONCLUSIONS: The RNA-Seq data and the collection of newly discovered coding SNPs improve the genomic resources available for cattle, especially for beef breeds. The large amount of variation present in genes expressed in Limousin Longissimus thoracis, especially the large number of non synonymous coding SNPs, may prove useful to study the mechanisms underlying the genetic variability of meat quality traits.