Gaining biological perspectives from schistosome genomes

Characterization of the genomic basis underlying schistosome biology is an important strategy for the development of future treatments and interventions. Genomic sequence is now available for the three major clinically relevant schistosome species, Schistosoma mansoni, S. japonicum and S. haematobium, and this information represents an invaluable resource for the future control of human schistosomiasis. The identification of a biologically important, but distinct from the host, schistosome gene product is the ultimate goal for many research groups. While the initial elucidation of the genome of an organism is critical for most biological research, continued improvement or curation of the genome construction should be an ongoing priority. In this review we will discuss prominent recent findings utilizing a systems approach to schistosome biology, as well as the increased use of interference RNA (RNAi). Both of these research strategies are aiming to place parasite genes into a more meaningful biological perspective.

Alternative splicing takes control of cytokine signaling

Thesis (Ph.D.)--University of Washington, 2016-08

CpxR activates MexAB-OprM efflux pump expression and enhances antibiotic resistance in both laboratory and clinical nalB-type isolates of Pseudomonas aeruginosa

Resistance-Nodulation-Division (RND) efflux pumps are responsible for multidrug resistance in Pseudomonas aeruginosa. In this study, we demonstrate that CpxR, previously identified as a regulator of the cell envelope stress response in Escherichia coli, is directly involved in activation of expression of RND efflux pump MexAB-OprM in P. aeruginosa. A conserved CpxR binding site was identified upstream of the mexA promoter in all genome-sequenced P. aeruginosa strains. CpxR is required to enhance mexAB-oprM expression and drug resistance, in the absence of repressor MexR, in P. aeruginosa strains PA14. As defective mexR is a genetic trait associated with the clinical emergence of nalB-type multidrug resistance in P. aeruginosa during antibiotic treatment, we investigated the involvement of CpxR in regulating multidrug resistance among resistant isolates generated in the laboratory via antibiotic treatment and collected in clinical settings. CpxR is required to activate expression of mexAB-oprM and enhances drug resistance, in the absence or presence of MexR, in ofloxacin-cefsulodin-resistant isolates generated in the laboratory. Furthermore, CpxR was also important in the mexR-defective clinical isolates. The newly identified regulatory linkage between CpxR and the MexAB-OprM efflux pump highlights the presence of a complex regulatory network modulating multidrug resistance in P. aeruginosa.

Phylum-wide transcriptome analysis of oogenesis and early embryogenesis in selected nematode species

Oogenesis is a prerequisite for embryogenesis in Metazoa. During both biological processes important decisions must be made to form the embryo and hence ensure the next generation: (1) Maternal gene products (mRNAs, proteins and nutrients) must be supplied to the embryo. (2) Polarity must be established and axes must be specified. While incorporation of maternal gene products occurs during oogenesis, the time point of polarity establishment and axis specification varies among species, as it is accomplished either prior, during, or after fertilisation. But not only the time point when these events take place varies among species but also the underlying mechanisms by which they are triggered. For the nematode model Caenorhabditis elegans the underlying pathways and gene regulatory networks (GRNs) are well understood. It is known that there the sperm entry point initiates a primary polarity in the 1-celled egg and with it the establishment of the anteroposterior axis. However, studies of other nematodes demonstrated that polarity establishment can be independent of sperm entry (Goldstein et al., 1998; Lahl et al., 2006) and that cleavage patterns, symmetry formation and cell specification also differ from C. elegans. In contrast to the studied Chromadorea (more derived nematodes including C. elegans), embryos of some marine Enoplea (more basal representatives) even show no discernible early polarity and blastomeres can adopt variable cell fates (Voronov and Panchin 1998). The underlying pathways controlling the obviously variant embryonic processes in non-Caenorhabditis nematodes are essentially unknown. In this thesis I addressed this issue by performing a detailed unbiased comparative transcriptome analysis based on microarrays and RNA sequencing of selected developmental stages in a variety of nematodes from different phylogenetic branches with C. elegans as a reference system and a nematomorph as an outgroup representative. In addition, I made use of available genomic data to determine the presence or absence of genes for which no expression had been detected. In particular, I focussed on components of selected pathways or GRNs which are known to play essential roles during C. elegans development and/or other invertebrate or vertebrate model systems. Oogenesis must be regulated differently in non-Caenorhabditis nematodes, as crucial controlling components of Wnt and sex determination signaling are absent in these species. In this respect, I identified female-specific expression of potential polarity associated genes during gonad development and oogenesis in the Enoplean nematode Romanomermis culicivorax. I could show that known downstream components of the polarity complexes PAR-3/-6/PKC-3 and PAR-1/-2 are absent in non-Caenorhabditis species. Even PAR-2 as part of the polarity complex does not exist in these nematodes. Instead, transcriptomes of nematodes (including C. elegans), show expression of other polarity-associated complexes such as the Lgl (Lethal giant larvae) complex. This result could pose an alternative route for nematodes and nematomorphs to initiate polarity during early embryogenesis. I could show that crucial pathways of axis specification, such as Wnt and BMP are very different in C. elegans compared to other nematodes. In the former, Wnt signaling, for instance, is mediated by four paralogous beta-catenins, while other Chromadorea have fewer and Enoplea only one beta-catenin. The transcriptomes of R. culicivorax and the nematomorph show that regulators of BMP (e.g. Chordin), are specifically expressed during early embryogenesis only in Enoplea and the close outgroup of nematomorphs. In conclusion, my results demonstrate that the molecular machinery controlling oogenesis and embryogenesis in nematodes is unexpectedly variable and C. elegans cannot be taken as a general model for nematode development. Under this perspective, Enoplean nematodes show more similarities with outgroups than with C. elegans. It appears that certain pathway components were lost or gained during evolution and others adopted new functions. Based on my findings I can conjecture, which pathway components may be ancestral and which were newly acquired in the course of nematode evolution.

Mutations in SNX14 cause a distinctive autosomal-recessive cerebellar ataxia and intellectual disability syndrome

Intellectual disability and cerebellar atrophy occur together in a large number of genetic conditions and are frequently associated with microcephaly and/or epilepsy. Here we report the identification of causal mutations in Sorting Nexin 14 (SNX14) found in seven affected individuals from three unrelated consanguineous families who presented with recessively inherited moderate-severe intellectual disability, cerebellar ataxia, early-onset cerebellar atrophy, sensorineural hearing loss, and the distinctive association of progressively coarsening facial features, relative macrocephaly, and the absence of seizures. We used homozygosity mapping and whole-exome sequencing to identify a homozygous nonsense mutation and an in-frame multiexon deletion in two families. A homozygous splice site mutation was identified by Sanger sequencing of SNX14 in a third family, selected purely by phenotypic similarity. This discovery confirms that these characteristic features represent a distinct and recognizable syndrome. SNX14 encodes a cellular protein containing Phox (PX) and regulator of G protein signaling (RGS) domains. Weighted gene coexpression network analysis predicts that SNX14 is highly coexpressed with genes involved in cellular protein metabolism and vesicle-mediated transport. All three mutations either directly affected the PX domain or diminished SNX14 levels, implicating a loss of normal cellular function. This manifested as increased cytoplasmic vacuolation as observed in cultured fibroblasts. Our findings indicate an essential role for SNX14 in neural development and function, particularly in development and maturation of the cerebellum.

Transcriptional regulation os phenylalaline biosynthesis and utilization

Conifer trees divert large quantities of carbon into the biosynthesis of phenylpropanoids, particularly to generate lignin, an important constituent of wood. Since phenylalanine is the precursor for phenylpropanoid biosynthesis, the precise regulation of phenylalanine synthesis and utilization should occur simultaneously. This crucial pathway is finely regulated primarily at the transcriptional level. Transcriptome analyses indicate that the transcription factors (TFs) preferentially expressed during wood formation in plants belong to the MYB and NAC families. Craven-Bartle et al. (2013) have shown in conifers that Myb8 is a candidate regulator of key genes in phenylalanine biosynthesis involved in the supply of the phenylpropane carbon skeleton necessary for lignin biosynthesis. This TF is able to bind AC elements present in the promoter regions of these genes to activate transcription. Constitutive overexpression of Myb8 in white spruce increased secondary-wall thickening and led to ectopic lignin deposition (Bomal et al. 2008). In Arabidopsis, the transcriptional network controlling secondary cell wall involves NAC-domain regulators operating upstream Myb transcription factors. Functional orthologues of members of this network described have been identified in poplar and eucalyptus, but in conifers functional evidence had only been obtained for MYBs. We have identified in the P. pinaster genome 37 genes encoding NAC proteins, which 3 NAC proteins could be potential candidates to be involved in vascular development (Pascual et al. 2015). The understanding of the transcriptional regulatory network associated to phenylpropanoids and lignin biosynthesis in conifers is crucial for future applications in tree improvement and sustainable forest management. This work is supported by the projects BIO2012-33797, BIO2015-69285-R and BIO-474 References: Bomal C, et al. (2008) Involvement of Pinus taeda MYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis. J Exp Bot. 59: 3925-3939. Craven-Bartle B, et al. (2013) A Myb transcription factor regulates genes of the phenylalanine pathway in maritime pine. Plant J, 74: 755-766. Pascual MB, et al. (2015) The NAC transcription factor family in maritime pine (Pinus pinaster): molecular regulation of two genes involved in stress responses. BMC Plant Biol, 15: 254.

NAC-MYB basedtranscriptional networkinvolved in the regulation of phenylalanine biosynthesis in P. pinaster

P2-2 NAC-MYB-BASED TRANSCRIPCIONAL NETWORK INVOLVED IN THE REGULATION OF PHENYLALANINE BIOSYNTHESIS IN P. PINASTER Mª Belén Pascual, Rafael A. Cañas, Blanca Craven-Bartle, Francisco M. Cánovas and Concepción Ávila Departamento de Biología Molecular y Bioquímica. Facultad de Ciencias. Universidad de Málaga. Campus de teatinos s/n, Málaga, Spain Email: cavila@uma.es Conifer trees divert large quantities of carbon into the biosynthesis of phenylpropanoids, particularly to generate lignin, an important constituent of wood. Since phenylalanine is the precursor for phenylpropanoid biosynthesis, the precise regulation of phenylalanine synthesis and use should occur simultaneously. This crucial pathway is finely regulated primarily at the transcriptional level. Transcriptome analyses indicate that the transcription factors (TFs) preferentially expressed during wood formation in plants belong to the MYB and NAC families. Craven-Bartle et al. (2013) have shown that Myb8 is a candidate regulator of key genes in phenylalanine biosynthesis involved in the supply of the phenylpropane carbon skeleton necessary for lignin biosynthesis. This TF is able to bind AC elements present in the promoter regions of these genes to activate transcription. In Arabidopsis, the transcriptional network controlling secondary cell wall involves NAC-domain regulators operating upstream Myb transcription factors. We have identified in the P. pinaster genome three NAC proteins as potential candidates to be involved in vascular development. One of them, PpNAC1 is expressed both in xylem and compression wood from adult trees and has been thoroughly characterized. Its role upstream the transcriptional network involving Myb8 will be discussed. The understanding of the transcriptional regulatory network associated to phenylpropanoids and lignin biosynthesis in conifers is crucial for future applications in tree improvement and sustainable forest management.

Revealing differences in gene network inference algorithms on the network level by ensemble methods

Motivation: The inference of regulatory networks from large-scale expression data holds great promise because of the potentially causal interpretation of these networks. However, due to the difficulty to establish reliable methods based on observational data there is so far only incomplete knowledge about possibilities and limitations of such inference methods in this context.

Results: In this article, we conduct a statistical analysis investigating differences and similarities of four network inference algorithms, ARACNE, CLR, MRNET and RN, with respect to local network-based measures. We employ ensemble methods allowing to assess the inferability down to the level of individual edges. Our analysis reveals the bias of these inference methods with respect to the inference of various network components and, hence, provides guidance in the interpretation of inferred regulatory networks from expression data. Further, as application we predict the total number of regulatory interactions in human B cells and hypothesize about the role of Myc and its targets regarding molecular information processing.

Structure and dynamics of a gene network model incorporating small RNAs

As advances in molecular biology continue to reveal additional layers of complexity in gene regulation, computational models need to incorporate additional features to explore the implications of new theories and hypotheses. It has recently been suggested that eukaryotic organisms owe their phenotypic complexity and diversity to the exploitation of small RNAs as signalling molecules. Previous models of genetic systems are, for several reasons, inadequate to investigate this theory. In this study, we present an artificial genome model of genetic regulatory networks based upon previous work by Torsten Reil, and demonstrate how this model generates networks with biologically plausible structural and dynamic properties. We also extend the model to explore the implications of incorporating regulation by small RNA molecules in a gene network. We demonstrate how, using these signals, highly connected networks can display dynamics that are more stable than expected given their level of connectivity.

Stem cell regulatory gene expression in human adult dental pulp and periodontal ligament cells undergoing odontogenic/osteogenic differentiation

Introduction During development and regeneration, odontogenesis and osteogenesis are initiated by a cascade of signals driven by several master regulatory genes. Methods In this study, we investigated the differential expression of 84 stem cell–related genes in dental pulp cells (DPCs) and periodontal ligament cells (PDLCs) undergoing odontogenic/osteogenic differentiation. Results Our results showed that, although there was considerable overlap, certain genes had more differential expression in PDLCs than in DPCs. CCND2, DLL1, and MME were the major upregulated genes in both PDLCs and DPCs, whereas KRT15 was the only gene significantly downregulated in PDLCs and DPCs in both odontogenic and osteogenic differentiation. Interestingly, a large number of regulatory genes in odontogenic and osteogenic differentiation interact or crosstalk via Notch, Wnt, transforming growth factor β (TGF-β)/bone morphogenic protein (BMP), and cadherin signaling pathways, such as the regulation of APC, DLL1, CCND2, BMP2, and CDH1. Using a rat dental pulp and periodontal defect model, the expression and distribution of both BMP2 and CDH1 have been verified for their spatial localization in dental pulp and periodontal tissue regeneration. Conclusions This study has generated an overview of stem cell–related gene expression in DPCs and PDLCs during odontogenic/osteogenic differentiation and revealed that these genes may interact through the Notch, Wnt, TGF-β/BMP, and cadherin signalling pathways to play a crucial role in determining the fate of dental derived cell and dental tissue regeneration. These findings provided a new insight into the molecular mechanisms of the dental tissue mineralization and regeneration

Cluster-based network model for time-course gene expression data

We propose a model-based approach to unify clustering and network modeling using time-course gene expression data. Specifically, our approach uses a mixture model to cluster genes. Genes within the same cluster share a similar expression profile. The network is built over cluster-specific expression profiles using state-space models. We discuss the application of our model to simulated data as well as to time-course gene expression data arising from animal models on prostate cancer progression. The latter application shows that with a combined statistical/bioinformatics analyses, we are able to extract gene-to-gene relationships supported by the literature as well as new plausible relationships.

A transcription factor map as revealed by a genome-wide gene expression analysis of whole-blood mRNA transcriptome in multiple sclerosis

Background Several lines of evidence suggests that transcription factors are involved in the pathogenesis of Multiple Sclerosis (MS) but a complete mapping the whole network has been elusive. One of the reasons is that there are several clinical subtypes of MS and transcription factors which may be involved in one subtype may not be in others. We investigated the possibility that this network could be mapped using microarray technologies and modern bioinformatics methods on a dataset from whole blood in 99 untreated MS patients (36 Relapse Remitting MS, 43 Primary Progressive MS, and 20 Secondary Progressive MS) and 45 age-matched healthy controls, Methodology/Principal Findings We have used two different analytical methodologies: a differential expression analysis and a differential co-expression analysis, which have converged on a significant number of regulatory motifs that seem to be statistically overrepresented in genes which are either differentially expressed (or differentially co-expressed) in cases and controls (e.g. V$KROX_Q6, p-value < 3.31E-6; V$CREBP1_Q2, p-value < 9.93E-6, V$YY1_02, p-value < 1.65E-5). Conclusions/significance: Our analysis uncovered a network of transcription factors that potentially dysregulate several genes in MS or one or more of its disease subtypes. Analysing the published literature we have found that these transcription factors are involved in the early T-lymphocyte specification and commitment as well as in oligodendrocytes dedifferentiation and development. The most significant transcription factors motifs were for the Early Growth response EGR/KROX family, ATF2, YY1 (Yin and Yang 1), E2F-1/DP-1 and E2F-4/DP-2 heterodimers, SOX5, and CREB and ATF families.

Hierarchical clustering of the genetic connectivity matrix reveals the network topology of gene action on brain microstructure: An N=531 twin study

Genetic correlation (rg) analysis determines how much of the correlation between two measures is due to common genetic influences. In an analysis of 4 Tesla diffusion tensor images (DTI) from 531 healthy young adult twins and their siblings, we generalized the concept of genetic correlation to determine common genetic influences on white matter integrity, measured by fractional anisotropy (FA), at all points of the brain, yielding an NxN genetic correlation matrix rg(x,y) between FA values at all pairs of voxels in the brain. With hierarchical clustering, we identified brain regions with relatively homogeneous genetic determinants, to boost the power to identify causal single nucleotide polymorphisms (SNP). We applied genome-wide association (GWA) to assess associations between 529,497 SNPs and FA in clusters defined by hubs of the clustered genetic correlation matrix. We identified a network of genes, with a scale-free topology, that influences white matter integrity over multiple brain regions.

Gene network effects on brain microstructure and intellectual performance identified in 472 twins

A major challenge in neuroscience is finding which genes affect brain integrity, connectivity, and intellectual function. Discovering influential genes holds vast promise for neuroscience, but typical genome-wide searches assess approximately one million genetic variants one-by-one, leading to intractable false positive rates, even with vast samples of subjects. Even more intractable is the question of which genes interact and how they work together to affect brain connectivity. Here, we report a novel approach that discovers which genes contribute to brain wiring and fiber integrity at all pairs of points in a brain scan. We studied genetic correlations between thousands of points in human brain images from 472 twins and their nontwin siblings (mean age: 23.7 2.1 SD years; 193 male/279 female).Wecombined clustering with genome-wide scanning to find brain systems withcommongenetic determination.Wethen filtered the image in a new way to boost power to find causal genes. Using network analysis, we found a network of genes that affect brain wiring in healthy young adults. Our new strategy makes it computationally more tractable to discover genes that affect brain integrity. The gene network showed small-world and scale-free topologies, suggesting efficiency in genetic interactions and resilience to network disruption. Genetic variants at hubs of the network influence intellectual performance by modulating associations between performance intelligence quotient and the integrity of major white matter tracts, such as the callosal genu and splenium, cingulum, optic radiations, and the superior longitudinal fasciculus.

A new regulatory variant in the interleukin-6 receptor gene associates with asthma risk

The main genetic determinant of soluble interleukin 6 receptor (sIL-6R) levels is the missense variant rs2228145 that maps to the cleavage site of IL-6R. For each Ala allele, sIL-6R serum levels increase by ∼20 ng ml -1 and asthma risk by 1.09-fold. However, this variant does not explain the total heritability for sIL-6R levels. Additional independent variants in IL6R may therefore contribute to variation in sIL-6R levels and influence asthma risk. We imputed 471 variants in IL6R and tested these for association with sIL-6R serum levels in 360 individuals. An intronic variant (rs12083537) was associated with sIL-6R levels independently of rs4129267 (P=0.0005), a proxy single-nucleotide polymorphism for rs2228145. A significant and consistent association for rs12083537 was observed in a replication panel of 354 individuals (P=0.033). Each rs12083537:A allele increased sIL-6R serum levels by 2.4 ng ml -1. Analysis of mRNA levels in two cohorts did not identify significant associations between rs12083537 and IL6R transcription levels. On the other hand, results from 16 705 asthmatics and 30 809 controls showed that the rs12083537:A allele increased asthma risk by 1.04-fold (P=0.0419). Genetic risk scores based on IL6R regulatory variants may prove useful in explaining variation in clinical response to tocilizumab, an anti-IL-6R monoclonal antibody.