Steady-state expression of self-regulated genes.

MOTIVATION: Regulatory gene networks contain generic modules such as feedback loops that are essential for the regulation of many biological functions. The study of the stochastic mechanisms of gene regulation is instrumental for the understanding of how cells maintain their expression at levels commensurate with their biological role, as well as to engineer gene expression switches of appropriate behavior. The lack of precise knowledge on the steady-state distribution of gene expression requires the use of Gillespie algorithms and Monte-Carlo approximations. METHODOLOGY: In this study, we provide new exact formulas and efficient numerical algorithms for computing/modeling the steady-state of a class of self-regulated genes, and we use it to model/compute the stochastic expression of a gene of interest in an engineered network introduced in mammalian cells. The behavior of the genetic network is then analyzed experimentally in living cells. RESULTS: Stochastic models often reveal counter-intuitive experimental behaviors, and we find that this genetic architecture displays a unimodal behavior in mammalian cells, which was unexpected given its known bimodal response in unicellular organisms. We provide a molecular rationale for this behavior, and we implement it in the mathematical picture to explain the experimental results obtained from this network.

The untranslated regions of genes from Trypanosoma cruzi: perspectives for functional characterization of strains and isolates

The sequencing of Trypanosoma cruzi genome has been completed and a great deal of information is now available. However, the organization of protozoa genomes is somewhat elusive and much effort must be applied to reveal all the information coded in the nucleotide sequences. Among the DNA segments that needs further investigation are the untranslated regions of genes. Many of the T. cruzi genes that were revealed by the genome sequencing lack information about the untranslated regions. In this paper, some features of these untranslated segments as well as their applications in T. cruzi populations are discussed.

Challenges of phylogenetic analyses of aDNA sequences

One of the crucial steps of authentication of aDNA sequences is phylogenetic consistency. Amplified sequences should fit into the phylogenetic framework of their supposed origin. An inherent property of aDNA sequences however, is their short sequence length. Additionally, genes for aDNA studies are often chosen by their preservation potential rather than by phylogenetically informative content. This poses potential challenges regarding their analyses, and might result in an inaccurate reflection of the supposed phylogenetic history of the sequence or organism under study. In this paper some fundamental problems of phylogenetic analysis and interpretation of aDNA datasets are discussed. Suggestions for character sampling and treatment of missing data are made. The publication is the result of a talk from the 1st PAMINSA Meeting in Rio de Janeiro, July 2005.

Re-mapping the molecular features of the human immunodeficiency virus type 1 and human T-cell lymphotropic virus type 1 Brazilian sequences using a bioinformatics unit established in Salvador, Bahia, Brazil, to give support to the viral epidemiology studies

The analysis of genetic data for human immunodeficiency virus type 1 (HIV-1) and human T-cell lymphotropic virus type 1 (HTLV-1) is essential to improve treatment and public health strategies as well as to select strains for vaccine programs. However, the analysis of large quantities of genetic data requires collaborative efforts in bioinformatics, computer biology, molecular biology, evolution, and medical science. The objective of this study was to review and improve the molecular epidemiology of HIV-1 and HTLV-1 viruses isolated in Brazil using bioinformatic tools available in the Laboratório Avançado de Sáude Pública (Lasp) bioinformatics unit. The analysis of HIV-1 isolates confirmed a heterogeneous distribution of the viral genotypes circulating in the country. The Brazilian HIV-1 epidemic is characterized by the presence of multiple subtypes (B, F1, C) and B/F1 recombinant virus while, on the other hand, most of the HTLV-1 sequences were classified as Transcontinental subgroup of the Cosmopolitan subtype. Despite the high variation among HIV-1 subtypes, protein glycosylation and phosphorylation domains were conserved in the pol, gag, and env genes of the Brazilian HIV-1 strains suggesting constraints in the HIV-1 evolution process. As expected, the functional protein sites were highly conservative in the HTLV-1 env gene sequences. Furthermore, the presence of these functional sites in HIV-1 and HTLV-1 strains could help in the development of vaccines that pre-empt the viral escape process.

Functional Variants in NOS1 and NOS2A Are Not Associated with Progressive Hearing Loss in Ménière's Disease in a European Caucasian Population

Hearing loss in Meniere's disease (MD) is associated with loss of spiral ganglion neurons and hair cells. In a guinea pig model of endolymphatic hydrops, nitric oxide synthases (NOS) and oxidative stress mediate loss of spiral ganglion neurons. To test the hypothesis that functional variants of NOS1 and NOS2A are associated with MD, wed genotyped three functional variants of NOS1 (rs41279104,rs2682826, and a cytosine-adenosine microsatellite repeat in exon 1f) and the CCTTT repeat in the promoter of NOS2A gene (rs3833912) in two independent MD sets(273 patients in total) and 550 controls. A third cohort of American patients was genotyped as replication cohort for the CCTTT repeat. Neither allele nor genotype frequencies of rs41279104 and rs2682826 were associated with MD, although longer alleles of the cytosine-adenosine microsatellite repeat were marginally significant (corrected p = 0.05) in the Mediterranean cohort but not in a second Galicia cohort. Shorter numbers of the CCTTT repeat in NOS2A were significantly more frequent in Galicia controls (OR = 0.37 [CI, 0.18-0.76], corrected p =0.04), but this finding could not be replicated in Mediterranean or American case-control populations. Meta-analysis did not support an association between CCTTT repeats and risk for MD. Severe hearing loss (>75 dB) was also not associated with any functional variants studied. Functional variants of NOS1 and and NOS2A do not confer susceptibility for MD.

Mutation analysis of genes that control the G1/S cell cycle in melanoma: TP53, CDKN1A, CDKN2A, and CDKN2B.

BACKGROUND The role of genes involved in the control of progression from the G1 to the S phase of the cell cycle in melanoma tumors in not fully known. The aim of our study was to analyse mutations in TP53, CDKN1A, CDKN2A, and CDKN2B genes in melanoma tumors and melanoma cell lines METHODS We analysed 39 primary and metastatic melanomas and 9 melanoma cell lines by single-stranded conformational polymorphism (SSCP). RESULTS The single-stranded technique showed heterozygous defects in the TP53 gene in 8 of 39 (20.5%) melanoma tumors: three new single point mutations in intronic sequences (introns 1 and 2) and exon 10, and three new single nucleotide polymorphisms located in introns 1 and 2 (C to T transition at position 11701 in intron 1; C insertion at position 11818 in intron 2; and C insertion at position 11875 in intron 2). One melanoma tumor exhibited two heterozygous alterations in the CDKN2A exon 1 one of which was novel (stop codon, and missense mutation). No defects were found in the remaining genes. CONCLUSION These results suggest that these genes are involved in melanoma tumorigenesis, although they may be not the major targets. Other suppressor genes that may be informative of the mechanism of tumorigenesis in skin melanomas should be studied.

Modeling stochasticity and robustness in gene regulatory networks.

MOTIVATION: Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations. RESULTS: In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs. AVAILABILITY: Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/ approximately garg/genysis.html.

Analysis of membrane protein genes in a Brazilian isolate of Anaplasma marginale

The sequencing of the complete genome of Anaplasma marginale has enabled the identification of several genes that encode membrane proteins, thereby increasing the chances of identifying candidate immunogens. Little is known regarding the genetic variability of genes that encode membrane proteins in A. marginale isolates. The aim of the present study was to determine the degree of conservation of the predicted amino acid sequences of OMP1, OMP4, OMP5, OMP7, OMP8, OMP10, OMP14, OMP15, SODb, OPAG1, OPAG3, VirB3, VirB9-1, PepA, EF-Tu and AM854 proteins in a Brazilian isolate of A. marginale compared to other isolates. Hence, primers were used to amplify these genes: omp1, omp4, omp5, omp7, omp8, omp10, omp14, omp15, sodb, opag1, opag3, virb3, VirB9-1, pepA, ef-tu and am854. After polimerase chain reaction amplification, the products were cloned and sequenced using the Sanger method and the predicted amino acid sequence were multi-aligned using the CLUSTALW and MEGA 4 programs, comparing the predicted sequences between the Brazilian, Saint Maries, Florida and A. marginale centrale isolates. With the exception of outer membrane protein (OMP) 7, all proteins exhibited 92-100% homology to the other A. marginale isolates. However, only OMP1, OMP5, EF-Tu, VirB3, SODb and VirB9-1 were selected as potential immunogens capable of promoting cross-protection between isolates due to the high degree of homology (over 72%) also found with A. (centrale) marginale.

Molecular characterization of regulatory polymorphisms in the promoter region of the STAT6 gene in a Gabonese population

Parasites remain competent invaders of host immunity. Their invasion strategies have proven to impact immunorelevant genes leading to diversity among gene families. We focussed on signal transducer and activator of transcription (STAT6) factor that plays a fundamental role in signal transduction and activation of transcription. Recent studies have highlighted the role of STAT6 variants in control of infection levels. We identified and investigated regulatory single nucleotide polymorphisms (SNPs) in the promoter regions of the STAT6 gene in a group of Gabonese individuals exposed to a variety of parasitic infections. Three promoter variants were identified in 40 individual subjects. We further validated these promoter variants for their allelic gene expression using transient transfection assays. One promoter variant, rs3024944 (G/C), revealed an altered expression of the marker gene. The identification of function-altering SNPs in the promoter may facilitate studying parasite susceptibility in association studies.

Regulatory elements involved in the post-transcriptional control of stage-specific gene expression in Trypanosoma cruzi: a review

Trypanosoma cruzi, a protozoan parasite that causes Chagas disease, exhibits unique mechanisms for gene expression such as constitutive polycistronic transcription of protein-coding genes, RNA editing and trans-splicing. In the absence of mechanism controlling transcription initiation, organized subsets of T. cruzi genes must be post-transcriptionally co-regulated in response to extracellular signals. The mechanisms that regulate stage-specific gene expression in this parasite have become much clearer through sequencing its whole genome as well as performing various proteomic and microarray analyses, which have demonstrated that at least half of the T. cruzi genes are differentially regulated during its life cycle. In this review, we attempt to highlight the recent advances in characterising cis and trans-acting elements in the T. cruzi genome that are involved in its post-transcriptional regulatory machinery.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications.

BACKGROUND: The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes. RESULTS: We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons. CONCLUSION: There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular - but possibly clusters of genes more generally - might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Kallikrein genes are associated with lupus and glomerular basement membrane-specific antibody-induced nephritis in mice and humans.

Immune-mediated nephritis contributes to disease in systemic lupus erythematosus, Goodpasture syndrome (caused by antibodies specific for glomerular basement membrane [anti-GBM antibodies]), and spontaneous lupus nephritis. Inbred mouse strains differ in susceptibility to anti-GBM antibody-induced and spontaneous lupus nephritis. This study sought to clarify the genetic and molecular factors that maybe responsible for enhanced immune-mediated renal disease in these models. When the kidneys of 3 mouse strains sensitive to anti-GBM antibody-induced nephritis were compared with those of 2 control strains using microarray analysis, one-fifth of the underexpressed genes belonged to the kallikrein gene family,which encodes serine esterases. Mouse strains that upregulated renal and urinary kallikreins exhibited less evidence of disease. Antagonizing the kallikrein pathway augmented disease, while agonists dampened the severity of anti-GBM antibody-induced nephritis. In addition, nephritis-sensitive mouse strains had kallikrein haplotypes that were distinct from those of control strains, including several regulatory polymorphisms,some of which were associated with functional consequences. Indeed, increased susceptibility to anti-GBM antibody-induced nephritis and spontaneous lupus nephritis was achieved by breeding mice with a genetic interval harboring the kallikrein genes onto a disease-resistant background. Finally, both human SLE and spontaneous lupus nephritis were found to be associated with kallikrein genes, particularly KLK1 and the KLK3 promoter, when DNA SNPs from independent cohorts of SLE patients and controls were compared. Collectively, these studies suggest that kallikreins are protective disease-associated genes in anti-GBM antibody-induced nephritis and lupus.

Molecular divergence in the timeless and cpr genes among three sympatric cryptic species of the Anopheles triannulatus complex

Anopheles triannulatus s.l. is a malaria vector with a wide geographic distribution, ranging from Argentina-Nicaragua and Trinidad. Here we analysed sequences of two genes, timeless and cpr, to assess the genetic variability and divergence among three sympatric cryptic species of this complex from Salobra, central-western Brazil. The timeless gene sequences did not conclusively differentiate Anopheles halophylus and An. triannulatus species "C". However, a partial separation has been observed between these species and An. triannulatus s.s. Importantly, the analysis of the cpr gene sequences revealed fixed differences, no shared polymorphisms and considerable genetic differentiation among the three species of the An. triannulatus complex. The results confirm that An. triannulatus s.s., An. halophylus and An. triannulatus species C are distinct taxa, with the latter two likely representing a more recent speciation event.

Control of hair follicle cell fate by underlying mesenchyme through a CSL-Wnt5a-FoxN1 regulatory axis.

Epithelial-mesenchymal interactions are key to skin morphogenesis and homeostasis. We report that maintenance of the hair follicle keratinocyte cell fate is defective in mice with mesenchymal deletion of the CSL/RBP-Jkappa gene, the effector of "canonical" Notch signaling. Hair follicle reconstitution assays demonstrate that this can be attributed to an intrinsic defect of dermal papilla cells. Similar consequences on hair follicle differentiation result from deletion of Wnt5a, a specific dermal papilla signature gene that we found to be under direct Notch/CSL control in these cells. Functional rescue experiments establish Wnt5a as an essential downstream mediator of Notch-CSL signaling, impinging on expression in the keratinocyte compartment of FoxN1, a gene with a key hair follicle regulatory function. Thus, Notch/CSL signaling plays a unique function in control of hair follicle differentiation by the underlying mesenchyme, with Wnt5a signaling and FoxN1 as mediators.

Chromosomal divergence and evolutionary inferences in Rhodniini based on the chromosomal location of ribosomal genes

In this study, we used fluorescence in situ hybridisation to determine the chromosomal location of 45S rDNA clusters in 10 species of the tribe Rhodniini (Hemiptera: Reduviidae: Triatominae). The results showed striking inter and intraspecific variability, with the location of the rDNA clusters restricted to sex chromosomes with two patterns: either on one (X chromosome) or both sex chromosomes (X and Y chromosomes). This variation occurs within a genus that has an unchanging diploid chromosome number (2n = 22, including 20 autosomes and 2 sex chromosomes) and a similar chromosome size and genomic DNA content, reflecting a genome dynamic not revealed by these chromosome traits. The rDNA variation in closely related species and the intraspecific polymorphism in Rhodnius ecuadoriensis suggested that the chromosomal position of rDNA clusters might be a useful marker to identify recently diverged species or populations. We discuss the ancestral position of ribosomal genes in the tribe Rhodniini and the possible mechanisms involved in the variation of the rDNA clusters, including the loss of rDNA loci on the Y chromosome, transposition and ectopic pairing. The last two processes involve chromosomal exchanges between both sex chromosomes, in contrast to the widely accepted idea that the achiasmatic sex chromosomes of Heteroptera do not interchange sequences.