Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Whole genome analysis of p38 SAPK-mediated gene expression upon stress

Background: Cells have the ability to respond and adapt to environmental changes through activation of stress-activated protein kinases (SAPKs). Although p38 SAPK signalling is known to participate in the regulation of gene expression little is known on the molecular mechanisms used by this SAPK to regulate stress-responsive genes and the overall set of genes regulated by p38 in response to different stimuli.Results: Here, we report a whole genome expression analyses on mouse embryonic fibroblasts (MEFs) treated with three different p38 SAPK activating-stimuli, namely osmostress, the cytokine TNFα and the protein synthesis inhibitor anisomycin. We have found that the activation kinetics of p38α SAPK in response to these insults is different and also leads to a complex gene pattern response specific for a given stress with a restricted set of overlapping genes. In addition, we have analysed the contribution of p38α the major p38 family member present in MEFs, to the overall stress-induced transcriptional response by using both a chemical inhibitor (SB203580) and p38α deficient (p38α-/-) MEFs. We show here that p38 SAPK dependency ranged between 60% and 88% depending on the treatments and that there is a very good overlap between the inhibitor treatment and the ko cells. Furthermore, we have found that the dependency of SAPK varies depending on the time the cells are subjected to osmostress. Conclusions: Our genome-wide transcriptional analyses shows a selective response to specific stimuli and a restricted common response of up to 20% of the stress up-regulated early genes that involves an important set of transcription factors, which might be critical for either cell adaptation or preparation for continuous extra-cellular changes. Interestingly, up to 85% of the up-regulated genes are under the transcriptional control of p38 SAPK. Thus, activation of p38 SAPK is critical to elicit the early gene expression program required for cell adaptation to stress.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Characterization of the role of AKAP-Lbc/p38 complex in the hypertrophic response of the heart

In response to chronic stress the heart undergoes an adverse remodeling process associated with cardiomyocyte hypertrophy, increased cellular apoptosis and fibrosis, which ultimately causes cardiac dysfunction and heart failure. Increasing evidence suggest the role of scaffolding and anchoring proteins in coordinating different signaling pathways that mediate the hypertrophic response of the heart. In this context, the family of Α-kinase anchoring proteins (AKAPs) emerged as important regulators of the cardiac function. During my thesis work I have conducted two independent projects, both of them aiming at elucidating the role of AKAPs in the heart. It has been shown that AKAP-Lbc, an anchoring protein that possesses an intrinsic Rho- specific exchange factor activity, organizes a signaling complex that links AKAP-Lbc- dependent activation of RhoA with the mitogen activated protein kinase (MAPK) p38. The first aim of my thesis was to study the role of this novel transduction pathway in the context of cardiac hypertrophy. Here we show that transgenic mice overexpressing in cardiomyocytes a competitor fragment of AKAP-Lbc, which specifically disrupts endogenous AKAP-Lbc / p38 complexes, developed early dilated cardiomyopathy in response to two weeks of transverse aortic constriction (TAC) as compared to controls. Interestingly, inhibition of the AKAP-Lbc / p38 transduction pathway significantly reduced the hypertrophic growth of single cardiomyocytes induced by pressure overload. Therefore, it appears that the AKAP- Lbc / p38 complex is crucially involved in the regulation of stress-induced cardiomyocyte hypertrophy and that disruption of this signaling pathway is detrimental for the heart under conditions of sustained hemodynamic stress. Secondly, in order to identify new AKAPs involved in the regulation of cardiac function, we followed a proteomic approach which allowed us to characterize AKAP2 as a major AKAP in the heart. Importantly, here we show that AKAP2 interacts with several proteins known to be involved in the control of gene transcription, such as the nuclear receptor coactivator 3 (NCoA3) or the ATP-dependent SWI/SNF chromatin remodeling complex. Thus, we propose AKAP2 as a novel mediator of cardiac gene expression through its interaction with these transcriptional regulators.

A large deletion in the adRP gene PRPF31: evidence that haploinsufficiency is the cause of disease.

PURPOSE: To report a large deletion that encompasses more than 90% of PRPF31 gene and two other neighboring genes in their entirety in an adRP pedigree that appears to show only the typical clinical features of retinitis pigmentosa. METHODS: To identify PRPF31 mutation in a dominant RP family (ADRP2) previously linked to the RP11 locus, the 14 exons of PRPF31 were screened for mutations by direct sequencing. To investigate the possibility of a large deletion, microsatellite markers near PRPF31 gene were analyzed by non-denaturing PAGE. RESULTS: Initial screening of PRPF31 gene in the ADRP2 family did not reveal an obvious mutation. A large deletion was however suspected due to lack of heterozygosity for nearly all PRPF31 intragenic single nucleotide polymorphysm (SNPs). In order to estimate the size of the deletion, SNPs and microsatellite markers spanning and flanking PRPF31 were analyzed in the entire ADRP2 family. Haplotype analysis with the above markers suggested a deletion of approximately 30 kb that included the putative promoter region of a novel gene OSCAR, the entire genomic content of genes NDUFA3, TFPT and more than 90% of PRPF31 gene. Sequence analysis of the region flanking the potential deletion showed a high presence of Alu elements implicating Alu mediated recombination as the mechanism responsible for this event. CONCLUSIONS: This mutation provides evidence that haploinsufficiency rather than aberrant function of mutated proteins is the cause of disease in these adRP patients with mutations in PRPF31 gene.

Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Dissection of the complex phenotype in cuticular mutants of Arabidopsis reveals a role of SERRATE as a mediator.

Mutations in LACERATA (LCR), FIDDLEHEAD (FDH), and BODYGUARD (BDG) cause a complex developmental syndrome that is consistent with an important role for these Arabidopsis genes in cuticle biogenesis. The genesis of their pleiotropic phenotypes is, however, poorly understood. We provide evidence that neither distorted depositions of cutin, nor deficiencies in the chemical composition of cuticular lipids, account for these features, instead suggesting that the mutants alleviate the functional disorder of the cuticle by reinforcing their defenses. To better understand how plants adapt to these mutations, we performed a genome-wide gene expression analysis. We found that apparent compensatory transcriptional responses in these mutants involve the induction of wax, cutin, cell wall, and defense genes. To gain greater insight into the mechanism by which cuticular mutations trigger this response in the plants, we performed an overlap meta-analysis, which is termed MASTA (MicroArray overlap Search Tool and Analysis), of differentially expressed genes. This suggested that different cell integrity pathways are recruited in cesA cellulose synthase and cuticular mutants. Using MASTA for an in silico suppressor/enhancer screen, we identified SERRATE (SE), which encodes a protein of RNA-processing multi-protein complexes, as a likely enhancer. In confirmation of this notion, the se lcr and se bdg double mutants eradicate severe leaf deformations as well as the organ fusions that are typical of lcr and bdg and other cuticular mutants. Also, lcr does not confer resistance to Botrytis cinerea in a se mutant background. We propose that there is a role for SERRATE-mediated RNA signaling in the cuticle integrity pathway.

Gene-environment interactions of selected pharmacogenes in arterial hypertension.

Hypertension affects approximately 1 billion people worldwide. Owing to population aging, hypertension-related cardiovascular burden is expected to rise in the near future. In addition to genetic variants influencing the blood pressure response to antihypertensive drugs, several genes encoding for drug-metabolizing or -transporting enzymes have been associated with blood pressure and/or hypertension in humans (e.g., ACE, CYP1A2, CYP3A5, ABCB1 and MTHFR) regardless of drug treatment. These genes are also involved in the metabolism and transport of endogenous substances and their effects may be modified by selected environmental factors, such as diet or lifestyle. However, little is currently known on the complex interplay between environmental factors, endogenous factors, genetic variants and drugs on blood pressure control. This review will discuss the respective role of population-based primary prevention and personalized medicine for arterial hypertension, taking a pharmacogenomics' perspective focusing on selected pharmacogenes.

PPARgamma but not PPARalpha ligands are potent repressors of major histocompatibility complex class II induction in atheroma-associated cells.

Peroxisome proliferator-activated receptors (PPARs) are essential in glucose and lipid metabolism and are implicated in metabolic disorders predisposing to atherosclerosis, such as diabetes and dyslipidemia. Conversely, antidiabetic glitazones and hypolipidemic fibrate drugs, known as PPARgamma and PPARalpha ligands, respectively, reduce the process of atherosclerotic lesion formation, which involves chronic immunoinflammatory processes. Major histocompatibility complex class II (MHC-II) molecules, expressed on the surface of specialized cells, are directly involved in the activation of T lymphocytes and in the control of the immune response. Interestingly, expression of MHC-II has recently been observed in atherosclerotic plaques, and it can be induced by the proinflammatory cytokine interferon-gamma (IFN-gamma) in vascular cells. To explore a possible role for PPAR ligands in the regulation of the immune response, we investigated whether PPAR activation affects MHC-II expression in atheroma-associated cells. In the present study, we demonstrate that PPARgamma but not PPARalpha ligands act as inhibitors of IFN-gamma-induced MHC-II expression and thus as repressors of MHC-II-mediated T-cell activation. All different types of PPARgamma ligands tested inhibit MHC-II. This effect of PPARgamma ligands is due to a specific inhibition of promoter IV of CIITA and does not concern constitutive expression of MHC-II. Thus, the beneficial effects of antidiabetic PPARgamma activators on atherosclerotic plaque development may be partly explained by their repression of MHC-II expression and subsequent inhibition of T-lymphocyte activation.

Reduced development of CD4-8-B220+ T cells but normal autoantibody production in lpr/lpr mice lacking major histocompatibility complex class I molecules.

The lpr gene has recently been shown to encode a functional mutation in the Fas receptor, a molecule involved in transducing apoptotic signals. Mice homozygous for the lpr gene develop an autoimmune syndrome accompanied by massive accumulation of double-negative (DN) CD4-8-B220+ T cell receptor-alpha/beta+ cells. In order to investigate the origin of these DN T cells, we derived lpr/lpr mice lacking major histocompatibility complex (MHC) class I molecules by intercrossing them with beta 2-microglobulin (beta 2m)-deficient mice. Interestingly, these lpr beta 2m-/- mice develop 13-fold fewer DNT cells in lymph nodes as compared to lpr/lpr wild-type (lprWT) mice. Analysis of anti-DNA antibodies and rheumatoid factor in serum demonstrates that lpr beta 2m-/- mice produce comparable levels of autoantibodies to lprWT mice. Collectively our data indicate that MHC class I molecules control the development of DN T cells but not autoantibody production in lpr/lpr mice and support the hypothesis that the majority of DN T cells may be derived from cells of the CD8 lineage.

Death and resurrection of the human IRGM gene

Immunity-related GTPases (IRG) play an important role in defense against intracellular pathogens. One member of this gene family in humans, IRGM, has been recently implicated as a risk factor for Crohn's disease. We analyzed the detailed structure of this gene family among primates and showed that most of the IRG gene cluster was deleted early in primate evolution, after the divergence of the anthropoids from prosimians ( about 50 million years ago). Comparative sequence analysis of New World and Old World monkey species shows that the single-copy IRGM gene became pseudogenized as a result of an Alu retrotransposition event in the anthropoid common ancestor that disrupted the open reading frame (ORF). We find that the ORF was reestablished as a part of a polymorphic stop codon in the common ancestor of humans and great apes. Expression analysis suggests that this change occurred in conjunction with the insertion of an endogenous retrovirus, which altered the transcription initiation, splicing, and expression profile of IRGM. These data argue that the gene became pseudogenized and was then resurrected through a series of complex structural events and suggest remarkable functional plasticity where alleles experience diverse evolutionary pressures over time. Such dynamism in structure and evolution may be critical for a gene family locked in an arms race with an ever-changing repertoire of intracellular parasites.

Transcriptional activation of the GLUT2 gene by the IPF-1/STF-1/IDX-1 homeobox factor.

The homeodomain protein PDX-1, referred as IPF-1/STF-1/IDX-1, is a transcriptional factor that plays a critical role in the control of several genes expressed in the pancreatic islet. PDX-1 gene expression has been previously shown to be reduced in cultured beta-cell lines chronically exposed to high glucose concentrations. As the glucose transporter type 2 (GLUT2) gene expression is selectively decreased in the beta-pancreatic cells of experimental models of diabetes, we postulated that the loss of GLUT2 gene expression in the pancreatic islets of diabetic animals may be due to the loss of PDX-1 transacting function on the GLUT2 gene. We, therefore, investigated the potential role of PDX-1 in the transcriptional control of GLUT2. We have identified a repeat of a TAAT motif (5'-TAATA-ATAACA-3') conserved in the sequence of the human and murine GLUT2 promoters. Recombinant PDX-1 binds to this GLUT2TAAT motif in electrophoretic mobility shift experiments. PDX-1 antiserum detects the formation of the complex of PDX-1 with the GLUT2TAAT motif in nuclear extracts from the pancreatic insulin-secreting cell line, beta TC3. The GLUT2TAAT motif was mutated in the murine GLUT2 promoter (-1308/+49 bp) linked to a luciferase reporter gene and transfected into beta TC3 cells. Compared with the transcriptional activity of the wild type promoter, that of the mutated promoter decreases by 41%. Multiple copies of the GLUT2TAAT motif were ligated 5' to a heterologous promoter and transfected into a PDX-1-expressing cell line (beta TC3) and into cell lines lacking the homeobox factor (InR1-G9 and JEG-3). The GLUT2TAAT motif mediates the activation of the heterologous promoter in the PDX-1-expressing cell line but not in InR1-G9 or JEG-3 cell lines. Furthermore, cotransfection in a PDX-1-deficient cell line with the expression vector encoding PDX-1 transactivates specifically the heterologous promoter containing the multimerized GLUT2TAAT motif. These data demonstrate that the murine GLUT2 promoter is controlled by the PDX-1 homeobox factor through the identified GLUT2TAAT motif.

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.