Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background: The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F2 crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism. Methodology/Principal Findings: GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained. Conclusions/Significance: Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.

The evolution of alternative splicing patterns and regulatory mechanisms

Alternative splicing produces multiple isoforms from the same gene, thus increasing the number of transcripts of the species. Alternative splicing is a virtually ubiquitous mechanism in eukaryotes, for example more than 90% of protein-coding genes in human are alternatively spliced. Recent evolutionary studies showed that alternative splicing is a fast evolving and highly species- specific mechanism. The rapid evolution of alternative splicing was considered as a contribution to the phenotypic diversity between species. However, the function of many isoforms produced by alternative splicing remains unclear and they might be the result of noisy splicing. Thus, the functional relevance of alternative splicing and the evolutionary mechanisms of its rapid divergence among species are still poorly understood. During my thesis, I performed a large-scale analysis of the regulatory mechanisms that drive the rapid evolution of alternative splicing. To study the evolution of alternative splicing regulatory mechanisms, I used an extensive RNA-sequencing dataset comprising 12 tetrapod species (human, chimpanzee and bonobo, gorilla, orangutan, macaque, marmoset, mouse, opossum, platypus, chicken and frog) and 8 tissues (cerebellum, brain, heart, kidney, liver, testis, placenta and ovary). To identify the catalogue of alternative splicing eis-acting regulatory elements in the different tetrapod species, I used a previously defined computational approach. This approach is a statistical analysis of exons/introns and splice sites composition and relies on a principle of compensation between splice sites strength and the presence of additional regulators. With an evolutionary comparative analysis of the exonic eis-acting regulators, I showed that these regulatory elements are generally shared among primates and more conserved than non-regulatory elements. In addition, I showed that the usage of these regulatory elements is also more conserved than expected by chance. In addition to the identification of species- specific eis-acting regulators, these results may explain the rapid evolution of alternative splicing. I also developed a new approach based on evolutionary sequence changes and corresponding alternative splicing changes to identify potential splicing eis-acting regulators in primates. The identification of lineage-specific substitutions and corresponding lineage-specific alternative splicing changes, allowed me to annotate the genomic sequences that might have played a role in the alternative splicing pattern differences among primates. Finally, I showed that the identified splicing eis-acting regulator datasets are enriched in human disease-causing mutations, thus confirming their biological relevance.

CBS: an open platform that integrates predictive methods and epigenetics information to characterize conserved regulatory features in multiple Drosophila genomes.

Background: Information about the composition of regulatory regions is of great value for designing experiments to functionally characterize gene expression. The multiplicity of available applications to predict transcription factor binding sites in a particular locus contrasts with the substantial computational expertise that is demanded to manipulate them, which may constitute a potential barrier for the experimental community. Results: CBS (Conserved regulatory Binding Sites, http://compfly.bio.ub.es/CBS) is a public platform of evolutionarily conserved binding sites and enhancers predicted in multiple Drosophila genomes that is furnished with published chromatin signatures associated to transcriptionally active regions and other experimental sources of information. The rapid access to this novel body of knowledge through a user-friendly web interface enables non-expert users to identify the binding sequences available for any particular gene, transcription factor, or genome region. Conclusions: The CBS platform is a powerful resource that provides tools for data mining individual sequences and groups of co-expressed genes with epigenomics information to conduct regulatory screenings in Drosophila.

Epigenetic regulatory elements: Recent advances in understanding their mode of action and use for recombinant protein production in mammalian cells.

Successful generation of high producing cell lines requires the generation of cell clones expressing the recombinant protein at high levels and the characterization of the clones' ability to maintain stable expression levels. The use of cis-acting epigenetic regulatory elements that improve this otherwise long and uncertain process has revolutionized recombinant protein production. Here we review and discuss new insights into the molecular mode of action of the matrix attachment regions (MARs) and ubiquitously-acting chromatin opening elements (UCOEs), i.e. cis-acting elements, and how these elements are being used to improve recombinant protein production. These elements can help maintain the chromatin environment of the transgene genomic integration locus in a transcriptionally favorable state, which increases the numbers of positive clones and the transgene expression levels. Moreover, the high producing clones tend to be more stable in long-term cultures even in the absence of selection pressure. Therefore, by increasing the probability of isolating a high producing clone, as well as by increasing transcription efficiency and stability, these elements can significantly reduce the time and cost required for producing large quantities of recombinant proteins.

Identification and characterization of the two-component NtrY/NtrX regulatory system in Azospirillum brasilense

Two Azospirillum brasilense open reading frames (ORFs) exhibited homology with the two-component NtrY/NtrX regulatory system from Azorhizobium caulinodans. These A. brasilense ORFs, located downstream to the nifR3ntrBC operon, were isolated, sequenced and characterized. The present study suggests that ORF1 and ORF2 correspond to the A. brasilense ntrY and ntrX genes, respectively. The amino acid sequences of A. brasilense NtrY and NtrX proteins showed high similarity to sensor/kinase and regulatory proteins, respectively. Analysis of lacZ transcriptional fusions by the ß-galactosidase assay in Escherichia coli ntrC mutants showed that the NtrY/NtrX proteins failed to activate transcription of the nifA promoter of A. brasilense. The ntrYX operon complemented a nifR3ntrBC deletion mutant of A. brasilense for nitrate-dependent growth, suggesting a possible cross-talk between the NtrY/X and NtrB/C sensor/regulator pairs. Our data support the existence of another two-component regulatory system in A. brasilense, the NtrY/NtrX system, probably involved in the regulation of nitrate assimilation.

Influence of the polymorphisms of the α-major regulatory element HS-40 on in vitro gene expression

The α-MRE is the major regulatory element responsible for the expression of human α-like globin genes. It is genetically polymorphic, and six different haplotypes, named A to F, have been identified in some population groups from Europe, Africa and Asia and in native Indians from two Brazilian Indian tribes. Most of the mutations that constitute the α-MRE haplotypes are located in flanking sequences of binding sites for nuclear factors. To our knowledge, there are no experimental studies evaluating whether such variability may influence the α-MRE enhancer activity. We analyzed and compared the expression of luciferase of nine constructs containing different α-MRE elements as enhancers. Genomic DNA samples from controls with A (wild-type α-MRE) and B haplotypes were used to generate C-F haplotypes by site-directed mutagenesis. In addition, three other elements containing only the G→A polymorphism at positions +130, +199, and +209, separately, were also tested. The different α-MRE elements were amplified and cloned into a plasmid containing the luciferase reporter gene and the SV40 promoter and used to transiently transfect K562 cells. A noticeable reduction in luciferase expression was observed with all constructs compared with the A haplotype. The greatest reductions occurred with the F haplotype (+96, C→A) and the isolated polymorphism +209, both located near the SP1 protein-binding sites believed not to be active in vivo. These are the first analyses of α-MRE polymorphisms on gene expression and demonstrate that these single nucleotide polymorphisms, although outside the binding sites for nuclear factors, are able to influence in vitro gene expression.

Essential regulatory elements for NHE3 gene transcription in renal proximal tubule cells

The objectives of the present study were to identify the cis-elements of the promoter absolutely required for the efficient rat NHE3 gene transcription and to locate positive and negative regulatory elements in the 5’-flanking sequence (5’FS), which might modulate the gene expression in proximal tubules, and to compare this result to those reported for intestinal cell lines. We analyzed the promoter activity of different 5’FS segments of the rat NHE3 gene, in the OKP renal proximal tubule cell line by measuring the activity of the reporter gene luciferase. Because the segment spanning the first 157 bp of 5’FS was the most active it was studied in more detail by sequential deletions, point mutations, and gel shift assays. The essential elements for gene transcription are in the region -85 to -33, where we can identify consensual binding sites for Sp1 and EGR-1, which are relevant to NHE3 gene basal transcription. Although a low level of transcription is still possible when the first 25 bp of the 5’FS are used as promoter, efficient transcription only occurs with 44 bp of 5’FS. There are negative regulatory elements in the segments spanning -1196 to -889 and -467 to -152, and positive enhancers between -889 and -479 bp of 5’FS. Transcription factors in the OKP cell nuclear extract efficiently bound to DNA elements of rat NHE3 promoter as demonstrated by gel shift assays, suggesting a high level of similarity between transcription factors of both species, including Sp1 and EGR-1.

Melanosomal targeting sequences from gp100 are essential for MHC class II-restricted endogenous epitope presentation and mobilization to endosomal compartments

CD4+ T lymphocytes play an important role in CD8+ T cell-mediated responses against tumors. Considering that about 20% of melanomas express major histocompatibility complex (MHC) class II, it is plausible that concomitant antigenic presentation by MHC class I and class II complexes shapes positive (helper T cells) or negative (regulatory T cells) anti-tumor responses. Interestingly, gp100, a melanoma antigen, can be presented by both MHC class I and class II when expressed endogenously, suggesting that it can reach endosomal/MHC class II compartments (MIIC). Here, we demonstrated that the gp100 putative amino-terminal signal sequence and the last 70 residues in carboxy-terminus, are essential for MIIC localization and MHC class II presentation. Confocal microscopy analyses confirmed that gp100 was localized in LAMP-1+ endosomal/MIIC. Gp100-targeting sequences were characterized by deleting different sections in the carboxy-terminus (residues 590 to 661). Transfection in 293T cells, expressing MHC class I and class II molecules, revealed that specific deletions in carboxy-terminus resulted in decreased MHC class II presentation, without effects on MHC class I presentation, suggesting a role in MIIC trafficking for these deleted sections. Then, we used these gp100-targeting sequences to mobilize the green fluorescent protein (GFP) to endosomal compartments, and to allow MHC class II and class I presentation of minimal endogenous epitopes. Thus, we concluded that these specific sequences are MIIC targeting motifs. Consequently, these sequences could be included in expression cassettes for endogenously expressed tumor or viral antigens to promote MHC class II and class I presentation and optimize in vivo T cell responses, or as an in vitro tool for characterization of new MHC class II epitopes.

DNA interaction and phosphotransfer of the C-4-dicarboxylate-responsive DcuS-DcuR two-component regulatory system from Escherichia coli

The DcuS-DcuR system of Escherichia coli is a two-component sensor-regulator that controls gene expression in response to external C-4-dicarboxylates and citrate. The DcuS protein is particularly interesting since it contains two PAS domains, namely a periplasmic C-4-dicarboxylate-sensing PAS domain (PASp) and a cytosolic PAS domain (PASc) of uncertain function. For a study of the role of the PASc domain, three different fragments of DcuS were overproduced and examined: they were PASc-kinase, PASc, and kinase. The two kinase-domain-containing fragments were autophosphorylated by [gamma-P-32]ATP. The rate was not affected by fumarate or succinate, supporting the role of the PASp domain in C-4-dicarboxylate sensing. Both of the phosphorylated DcuS constructs were able to rapidly pass their phosphoryl groups to DcuR, and after phosphorylation, DcuR dephosphorylated rapidly. No prosthetic group or significant quantity of metal was found associated with either of the PASc-containing proteins. The DNA-binding specificity of DcuR was studied by use of the pure protein. It was found to be converted from a monomer to a dimer upon acetylphosphate treatment, and native polyacrylamide gel electrophoresis suggested that it can oligomerize. DcuR specifically bound to the promoters of the three known DcuSR-regulated genes (dctA, dcuB, and frdA), with apparent K(D)s of 6 to 32 muM for untreated DcuR and less than or equal to1 to 2 muM for the acetylphosphate-treated form. The binding sites were located by DNase I footprinting, allowing a putative DcuR-binding motif [tandemly repeated (T/A)(A/T)(T/C)(A/T)AA sequences] to be identified. The DcuR-binding sites of the dcuB, dctA, and frdA genes were located 27, 94, and 86 bp, respectively, upstream of the corresponding +1 sites, and a new promoter was identified for dcuB that responds to DcuR.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

14-3-3: Defense and regulatory proteins coding by Eucalyptus genome

The data mining of Eucalyptus ESTs genome finds four clusters (EGCEST2257E11.g, EGBGRT3213F11.g, and EGCCFB1223H11.g) from highly conservative 14-3-3 protein family which modulates a wide variety of cellular processes. Multiple alignments were built from twenty four sequences of 14-3-3 proteins searched into the GenBank databases and into the four pools of Eucalyptus genome programs. The alignment has shown two regions highly conservative on the sequences corresponding to the motifs of protein phosphorylation and nine highly conservative regions on the sequence corresponding to the linkage regions of alpha helices structure based on three dimensional of dimer functional structure. The differences of amino acid into the structural and functional domains of 14-3-3 plant protein were identified and can explain the functional diversity of different isoforms. The phylogenic protein trees were built by the maximum parsimony and neighborjoining procedures of Clustal X alignments and PAUP software for phylogenic analysis.

Identificação de Single Nucleotilde Polymorphisms (SMPs) no gene Nove-cis-epoxicaroteníde dioxigenase (NCED) em Eucalyptus

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Molecular characterization of interferon regulatory factor 1 in Bubalus bubalis

Interferon regulatory factor 1 (IRF1) is functionally diverse in the regulation of immune response and is considered to be an important candidate gene for studying disease susceptibility in mammals. In this paper, we characterized the whole sequence of the IRF1 gene in river buffalo (Bubalus bubalis) and compared genomic and the amino acid sequences between different species. The buffalo IRF1 gene was 7099 bp long and organized into 10 exons and nine introns. Its molecular structure showed exactly the same number of exons (10) and introns (nine) in bovids, mice, horses, humans, and chickens. However, rats did not have exon 5, but had the largest exon 4, which suggests that exon 5 was incorporated into exon 4. The coding and the amino acid sequences of the gene showed that identity varied from 73 to 99% at the coding sequence level and from 61 to 100% at the amino acid level when compared with other mammals and chickens. Comparative analysis of the gene sequence between two different buffalo breeds, Murrah and Mediterranean, revealed six potential SNPs that are primarily located in the 5' and 3'UTRs.

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations

Abstract Background Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC Methods Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences. Results A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) <5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR <5%, p-value ≤0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p <0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes’ as the main enriched category. A module map analysis of the protein-coding genes significantly (p <0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p <0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p <0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation. Conclusion Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.

Identifizierung und Charakterisierung cis regulatorischer Elemente des humanen Gens für Keratin 20 in vitro und im transgenen Mausmodell

ZusammenfassungKeratin 20 (K20) ist ein Intermediärfilament, das als Strukturelement in den Epithelien des Intestinaltrakts und den Merkelzellen der Haut exprimiert wird. Im Rahmen dieser Arbeit wurden verschiedene K20 Expressionsvektoren generiert, mit denen regulatorische Elemente des humanen Gens für K20 charakterisiert wurden. Analysiert wurde der Bereich von –4,8 kb bzw. –21,5 kb bis 12,9 kb vom Transkriptionsstart. Die Vektoren, welche die Sequenz von –21,5 kb bis 12,9 kb umfassten, konnten die Expression des EGFP Reportergens in HT-29 Zellen signifikant steigern. Zwei Enhancer-Elemente zwischen –21,5 und –18,4 kb bzw. –18,4 und –14,9 kb verstärkten die Expression der Reporterkonstrukte in vitro signifikant. Die Analyse der Vektoren in transgenen Mäusen zeigte, dass diese das Transgen gewebespezifisch exprimieren. Mit dem Bereich von –4,8 kb bis 12,9 kb ließ sich transgenspezifische mRNA Expression in intestinalen Geweben mittels RT-PCR nachweisen. Die Verwendung der Sequenz von –21,5 kb bis 21,8 kb steigerte die Expression in vivo nicht weiter.Für die gewebespezifische Expression des K20 Vektors reicht die 4,8 kb 5´upstream Sequenz aus, der Bereich bis –21,5 kb verstärkt die Expression in vitro, allerdings fehlen für die starke gewebespezifische Expression von Transgenen in vivo noch weitere Kontrollelemente.