Conserved gene structure and function of interleukin-10 in teleost fish

Interleukin-10 (IL-10) is an important immunoregulatory cytokine produced by various types of cells. Researchers describe here the isolation and characterization of olive flounder IL-10 (ofIL-10) cDNA and genomic organization. The ofIL-10 gene encodes a 187 amino acid protein and is composed of a five exon/four intron structure, similar to other known IL-10 genes. The ofIL-10 promoter sequence analysis shows a high level of homology in putative binding sites for transcription factors which are sufficient for transcriptional regulation ofIL-10. Important structural residues are maintained in the ofIL-10 protein including the four cysteines responsible for the two intra-chain disulfide bridges reported for human IL-10 and two extra cysteine residues that exist only in fish species. The phylogenetic analysis clustered ofIL-10 with other fish IL-10s and apart from mammalian IL-10 molecules. Quantitative real-time Polymerase Chain Reaction (PCR) analysis demonstrated ubiquitous ofIL-10 gene expression in the 13 tissues examined. Additionally, the induction of ofIL-10 gene expression was observed in the kidney tissue from olive flounder infected with bacteria (Edawardsiella tarda) or virus (Viral Hemorrhagic Septicemia Virus; VHSV). These data indicate that IL-10 is an important immune regulator that is conserved strictly genomic organization and function during the evolution of vertebrate immunity.

Computational Methods for Locating and Analyzing Conserved Gene Regulatory DNA Elements

This thesis presents methods for locating and analyzing cis-regulatory DNA elements involved with the regulation of gene expression in multicellular organisms. The regulation of gene expression is carried out by the combined effort of several transcription factor proteins collectively binding the DNA on the cis-regulatory elements. Only sparse knowledge of the 'genetic code' of these elements exists today. An automatic tool for discovery of putative cis-regulatory elements could help their experimental analysis, which would result in a more detailed view of the cis-regulatory element structure and function. We have developed a computational model for the evolutionary conservation of cis-regulatory elements. The elements are modeled as evolutionarily conserved clusters of sequence-specific transcription factor binding sites. We give an efficient dynamic programming algorithm that locates the putative cis-regulatory elements and scores them according to the conservation model. A notable proportion of the high-scoring DNA sequences show transcriptional enhancer activity in transgenic mouse embryos. The conservation model includes four parameters whose optimal values are estimated with simulated annealing. With good parameter values the model discriminates well between the DNA sequences with evolutionarily conserved cis-regulatory elements and the DNA sequences that have evolved neutrally. In further inquiry, the set of highest scoring putative cis-regulatory elements were found to be sensitive to small variations in the parameter values. The statistical significance of the putative cis-regulatory elements is estimated with the Two Component Extreme Value Distribution. The p-values grade the conservation of the cis-regulatory elements above the neutral expectation. The parameter values for the distribution are estimated by simulating the neutral DNA evolution. The conservation of the transcription factor binding sites can be used in the upstream analysis of regulatory interactions. This approach may provide mechanistic insight to the transcription level data from, e.g., microarray experiments. Here we give a method to predict shared transcriptional regulators for a set of co-expressed genes. The EEL (Enhancer Element Locator) software implements the method for locating putative cis-regulatory elements. The software facilitates both interactive use and distributed batch processing. We have used it to analyze the non-coding regions around all human genes with respect to the orthologous regions in various other species including mouse. The data from these genome-wide analyzes is stored in a relational database which is used in the publicly available web services for upstream analysis and visualization of the putative cis-regulatory elements in the human genome.

Production of bioactive diterpenoids in the Euphorbiaceae depends on evolutionarily conserved gene clusters

The Euphorbiaceae produce a diverse range of diterpenoids, many of which have pharmacological activities. These diterpenoids include ingenol mebutate, which is licensed for the treatment of a precancerous skin condition (actinic keratosis), and phorbol derivatives such as resiniferatoxin and prostratin, which are undergoing investigation for the treatment of severe pain and HIV, respectively. Despite the interest in these diterpenoids, their biosynthesis is poorly understood at present, with the only characterized step being the conversion of geranylgeranyl pyrophosphate into casbene. Here, we report a physical cluster of diterpenoid biosynthetic genes from castor (Ricinus communis), including casbene synthases and cytochrome P450s from the CYP726A subfamily. CYP726A14, CYP726A17, and CYP726A18 were able to catalyze 5-oxidation of casbene, a conserved oxidation step in the biosynthesis of this family of medicinally important diterpenoids. CYP726A16 catalyzed 7,8-epoxidation of 5-keto-casbene and CYP726A15 catalyzed 5-oxidation of neocembrene. Evidence of similar gene clustering was also found in two other Euphorbiaceae, including Euphorbia peplus, the source organism of ingenol mebutate. These results demonstrate conservation of gene clusters at the higher taxonomic level of the plant family and that this phenomenon could prove useful in further elucidating diterpenoid biosynthetic pathways.

Spatio-temporal distribution of fatty acid-binding protein 6 (fabp6) gene transcripts in the developing and adult zebrafish (Danio rerio)

We have determined the structure of the fatty acid-binding protein 6 (fabp6) gene and the tissue-specific distribution of its transcripts in embryos, larvae and adult zebrafish (Danio rerio). Like most members of the vertebrate FABP multigene family, the zebrafish fabp6 gene contains four exons separated by three introns. The coding region of the gene and expressed sequence tags code for a polypeptide of 131 amino acids (14 kDa, pI 6.59). The putative zebrafish Fabp6 protein shared greatest sequence identity with human FABP6 (55.3%) compared to other orthologous mammalian FABPs and paralogous zebrafish Fabps. Phylogenetic analysis showed that the zebrafish Fabp6 formed a distinct clade with the mammalian FABP6s. The zebrafish fabp6 gene was assigned to linkage group (chromosome) 21 by radiation hybrid mapping. Conserved gene synteny was evident between the zebrafish fabp6 gene on chromosome 21 and the FABP6/Fabp6 genes on human chromosome 5, rat chromosome 10 and mouse chromosome 11. Zebrafish fabp6 transcripts were first detected in the distal region of the intestine of embryos at 72 h postfertilization. This spatial distribution remained constant to 7-day-old larvae, the last stage assayed during larval development. In adult zebrafish, fabp6 transcripts were detected by RT-PCR in RNA extracted from liver, heart, intestine, ovary and kidney (most likely adrenal tissue), but not in RNA from skin, brain, gill, eye or muscle. In situ hybridization of a fabp6 riboprobe to adult zebrafish sections revealed intense hybridization signals in the adrenal homolog of the kidney and the distal region of the intestine, and to a lesser extent in ovary and liver, a transcript distribution that is similar, but not identical, to that seen for the mammalian FABP6/Fabp6 gene. © 2008 The Authors.

Worldwide HLA-E nucleotide and haplotype variability reveals a conserved gene for coding and 3 ' untranslated regions

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

ECA39, a conserved gene regulated by c-Myc in mice, is involved in G1/S cell cycle regulation in yeast.

The c-myc oncogene has been shown to play a role in cell proliferation and apoptosis. The realization that myc oncogenes may control the level of expression of other genes has opened the field to search for genetic targets for Myc regulation. Recently, using a subtraction/coexpression strategy, a murine genetic target for Myc regulation, called EC439, was isolated. To further characterize the ECA39 gene, we set out to determine the evolutionary conservation of its regulatory and coding sequences. We describe the human, nematode, and budding yeast homologs of the mouse ECA39 gene. Identities between the mouse ECA39 protein and the human, nematode, or yeast proteins are 79%, 52%, and 49%, respectively. Interestingly, the recognition site for Myc binding, located 3' to the start site of transcription in the mouse gene, is also conserved in the human homolog. This regulatory element is missing in the ECA39 homologs from nematode or yeast, which also lack the regulator c-myc. To understand the function of ECA39, we deleted the gene from the yeast genome. Disruption of ECA39 which is a recessive mutation that leads to a marked alteration in the cell cycle. Mutant haploids and homozygous diploids have a faster growth rate than isogenic wild-type strains. Fluorescence-activated cell sorter analyses indicate that the mutation shortens the G1 stage in the cell cycle. Moreover, mutant strains show higher rates of UV-induced mutations. The results suggest that the product of ECA39 is involved in the regulation of G1 to S transition.

The Genome of Anopheles darlingi, the main neotropical malaria vector

Anopheles darlingi is the principal neotropical malaria vector, responsible for more than a million cases of malaria per year on the American continent. Anopheles darlingi diverged from the African and Asian malaria vectors ∼100 million years ago (mya) and successfully adapted to the New World environment. Here we present an annotated reference A. darlingi genome, sequenced from a wild population of males and females collected in the Brazilian Amazon. A total of 10 481 predicted protein-coding genes were annotated, 72% of which have their closest counterpart in Anopheles gambiae and 21% have highest similarity with other mosquito species. In spite of a long period of divergent evolution, conserved gene synteny was observed between A. darlingi and A. gambiae. More than 10 million single nucleotide polymorphisms and short indels with potential use as genetic markers were identified. Transposable elements correspond to 2.3% of the A. darlingi genome. Genes associated with hematophagy, immunity and insecticide resistance, directly involved in vectorhuman and vectorparasite interactions, were identified and discussed. This study represents the first effort to sequence the genome of a neotropical malaria vector, and opens a new window through which we can contemplate the evolutionary history of anopheline mosquitoes. It also provides valuable information that may lead to novel strategies to reduce malaria transmission on the South American continent. The A. darlingi genome is accessible at www.labinfo.lncc.br/index.php/anopheles- darlingi. © 2013 The Author(s).

Origin and evolution of the RIG-I like RNA helicase gene family

Background: The DExD/H domain containing RNA helicases such as retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5) are key cytosolic pattern recognition receptors (PRRs) for detecting nucleotide pathogen associated molecular patterns (PAMPs) of invading viruses. The RIG-I and MDA5 proteins differentially recognise conserved PAMPs in double stranded or single stranded viral RNA molecules, leading to activation of the interferon system in vertebrates. They share three core protein domains including a RNA helicase domain near the C terminus (HELICc), one or more caspase activation and recruitment domains (CARDs) and an ATP dependent DExD/H domain. The RIG-I/MDA5 directed interferon response is negatively regulated by laboratory of genetics and physiology 2 (LGP2) and is believed to be controlled by the mitochondria antiviral signalling protein (MAVS), a CARD containing protein associated with mitochondria. Results: The DExD/H containing RNA helicases including RIG-I, MDA5 and LGP2 were analysed in silico in a wide spectrum of invertebrate and vertebrate genomes. The gene synteny of MDA5 and LGP2 is well conserved among vertebrates whilst conservation of the gene synteny of RIG-I is less apparent. Invertebrate homologues had a closer phylogenetic relationship with the vertebrate RIG-Is than the MDA5/LGP2 molecules, suggesting the RIG-I homologues may have emerged earlier in evolution, possibly prior to the appearance of vertebrates. Our data suggest that the RIG-I like helicases possibly originated from three distinct genes coding for the core domains including the HELICc, CARD and ATP dependent DExD/H domains through gene fusion and gene/domain duplication. Furthermore, presence of domains similar to a prokaryotic DNA restriction enzyme III domain (Res III), and a zinc finger domain of transcription factor (TF) IIS have been detected by bioinformatic analysis. Conclusion: The RIG-I/MDA5 viral surveillance system is conserved in vertebrates. The RIG-I like helicase family appears to have evolved from a common ancestor that originated from genes encoding different core functional domains. Diversification of core functional domains might be fundamental to their functional divergence in terms of recognition of different viral PAMPs.

Morphological and genomic characterization of Filobasidiella depauperata: a homothallic sibling species of the pathogenic cryptococcus species complex.

The fungal species Cryptococcus neoformans and Cryptococcus gattii cause respiratory and neurological disease in animals and humans following inhalation of basidiospores or desiccated yeast cells from the environment. Sexual reproduction in C. neoformans and C. gattii is controlled by a bipolar system in which a single mating type locus (MAT) specifies compatibility. These two species are dimorphic, growing as yeast in the asexual stage, and producing hyphae, basidia, and basidiospores during the sexual stage. In contrast, Filobasidiella depauperata, one of the closest related species, grows exclusively as hyphae and it is found in association with decaying insects. Examination of two available strains of F. depauperata showed that the life cycle of this fungal species shares features associated with the unisexual or same-sex mating cycle in C. neoformans. Therefore, F. depauperata may represent a homothallic and possibly an obligately sexual fungal species. RAPD genotyping of 39 randomly isolated progeny from isolate CBS7855 revealed a new genotype pattern in one of the isolated basidiospores progeny, therefore suggesting that the homothallic cycle in F. depauperata could lead to the emergence of new genotypes. Phylogenetic analyses of genes linked to MAT in C. neoformans indicated that two of these genes in F. depauperata, MYO2 and STE20, appear to form a monophyletic clade with the MATa alleles of C. neoformans and C. gattii, and thus these genes may have been recruited to the MAT locus before F. depauperata diverged. Furthermore, the ancestral MATa locus may have undergone accelerated evolution prior to the divergence of the pathogenic Cryptococcus species since several of the genes linked to the MATa locus appear to have a higher number of changes and substitutions than their MATalpha counterparts. Synteny analyses between C. neoformans and F. depauperata showed that genomic regions on other chromosomes displayed conserved gene order. In contrast, the genes linked to the MAT locus of C. neoformans showed a higher number of chromosomal translocations in the genome of F. depauperata. We therefore propose that chromosomal rearrangements appear to be a major force driving speciation and sexual divergence in these closely related pathogenic and saprobic species.

Molecular characterization of Gla-rich protein (GRP) gene expression and function

Gla-rich protein (GRP) is a vitamin K-dependent protein related to bone and cartilage recently described. This protein is characterized by a large number of Gla (γ-carboxyglutamic acid) residues being the protein with the highest Gla content of any known protein. It was found in a widely variety of tissues but highest levels was found in skeletal and cartilaginous tissues. This small secreted protein was also expressed and accumulated in soft tissues and it was clearly associated with calcification pathologies in the same tissues. Although the biological importance of GRP remains to be elucidated, it was suggested a physiological role in cartilage development and calcification process during vertebrate skeleton formation. Using zebrafish, an accepted model to study skeletal development, we have described two grp paralog genes, grp1 and grp2, which exhibited distinct patterns of expression, suggesting different regulatory pathways for each gene. Gene synteny analysis showed that grp2 gene is more closely related to tetrapod grp, although grp1 gene was proposed to be the vertebrate ortholog by sequence comparison. In addition, we identified a functional promoter of grp2 gene and using a functional approach we confirmed the involvement of transcription factors from Sox family (Sox9b and Sox10) in the regulation of grp2 expression. In an effort to provide more information about the function of grp isoforms, we generated two zebrafish transgenic lines capable to overexpress conditionally grp genes and possible roles in the skeleton development were studied. To better understand GRP function a mammalian system was used and the analysis of knockout mice showed that GRP is involved in chondrocyte maturation and the absence of GRP is associated to proteoglycans loss in calcified articular cartilage. In addition, we detected differences in chondrogenesis markers in articular chondrocyte primary culture. Overall, our data suggest a main role for GRP on chondrocyte differentiation.

Annotation and characterization of the Plasmodium vivax rhoptry neck protein 4 (PvRON4)

Background: The tight junction (TJ) is one of the most important structures established during merozoite invasion of host cells and a large amount of proteins stored in Toxoplasma and Plasmodium parasites’ apical organelles are involved in forming the TJ. Plasmodium falciparum and Toxoplasma gondii apical membrane antigen 1 (AMA-1) and rhoptry neck proteins (RONs) are the two main TJ components. It has been shown that RON4 plays an essential role during merozoite and sporozoite invasion to target cells. This study has focused on characterizing a novel Plasmodium vivax rhoptry protein, RON4, which is homologous to PfRON4 and PkRON4. Methods: The ron4 gene was re-annotated in the P. vivax genome using various bioinformatics tools and taking PfRON4 and PkRON4 amino acid sequences as templates. Gene synteny, as well as identity and similarity values between open reading frames (ORFs) belonging to the three species were assessed. The gene transcription of pvron4, and the expression and localization of the encoded protein were also determined in the VCG-1 strain by molecular and immunological studies. Nucleotide and amino acid sequences obtained for pvron4 in VCG-1 were compared to those from strains coming from different geographical areas. Results: PvRON4 is a 733 amino acid long protein, which is encoded by three exons, having similar transcription and translation patterns to those reported for its homologue, PfRON4. Sequencing PvRON4 from the VCG-1 strain and comparing it to P. vivax strains from different geographical locations has shown two conserved regions separated by a low complexity variable region, possibly acting as a “smokescreen”. PvRON4 contains a predicted signal sequence, a coiled-coil α-helical motif, two tandem repeats and six conserved cysteines towards the carboxyterminus and is a soluble protein lacking predicted transmembranal domains or a GPI anchor. Indirect immunofluorescence assays have shown that PvRON4 is expressed at the apical end of schizonts and co-localizes at the rhoptry neck with PvRON2.

Estudo de genômica comparativa de corynebacterium pseudotuberculosis linhagem 226 (biovar ovis)

Corynebacterium pseudotuberculosis é uma bactéria Gram-positiva, intracelular facultativa, nãoesporulante, não-capsulada e sem mobilidade, contudo possui fímbria, e pode assumir formas cocóides e filamentosas (pleomórfica), além disto, apresenta crescimento ótimo à 37°C. Este patógeno apresenta dois biovares: ovis que geralmente acomete pequenos ruminantes, e causa a doença linfadenite caseosa, e biovar equi, mais comum em equinos, bovinos, camelídeos, e bubalinos causando a Linfangite ulcerativa. A infecção por esta bactéria pode levar a condenação das carcaças e redução de lã (em ovinos e caprinos), leite e carne destes animais, e consequentemente a perdas econômicas para a indústria agropecuária mundial. Atualmente, ainda não existe uma vacina eficaz para estas doenças. A fim de obter um maior entendimento biológico entre as espécies o presente trabalho tem como objetivo principal analisar, por meio da genômica comparativa a linhagem C. pseudotuberculosis 226 biotipo ovis isolada de um caprino na Califórnia com outras linhagens do biovarar ovis e equi. Na análise de sintenia entre as linhagens foi possível identificar que a linhagem 226 apresenta alta conservação da ordem gênica entre as linhagens do biótipo ovis. Através de análises filogenômicas foi possível identificar que as linhagens I19 e 267 apresentaram maior e menor proximidade filogenética com a linhagem 226. A linhagem 1/06-A foi a que apresentou maior proximidade filogenômica entre as linhagens do biovar equi, quando comparadas a linhagem 226. Foram preditas 8 ilhas de patogenicidade, estando presente na ilha 1 os genes relacionados a virulência de C. pseudotuberculosis mais bem descritos na literatura. Não houveram regiões novas relacionadas a genes de virulência entre nenhuma das linhagens. Foram identificados 248 genes ortólogos entre as linhagens I19, 267 e 226 e 282 genes ortólogos entre as linhagens 258,1/06-A e 226. Com base nesse estudo é possível inferir que as linhagens do biovar ovis possuem um repertório gênico pouco variado e que as linhagens do biovar equi apresentam uma quantidade menor de genes compartilhados com a linhagem 226, corroborando com a diversidade gênica entre os biovares.

Amplification of the flgE gene provides evidence for the existence of a Brazilian borreliosis

INTRODUCTION: The symptoms of Brazilian borreliosis resemble the clinical manifestations of Lyme disease (LD). However, there are differences between the two in terms of epidemiological and laboratory findings. Primers usually employed to diagnose LD have failed to detect Borrelia strains in Brazil. OBJECTIVE: We aimed to identify the Brazilian Borrelia using a conserved gene that synthesizes the flagellar hook (flgE) of Borrelia burgdorferi sensu lato. METHOD: Three patients presenting with erythema migrans and positive epidemiological histories were recruited for the study. Blood samples were collected, and the DNA was extracted by commercial kits. RESULTS: The gene flgE was amplified from DNA of all selected patients. Upon sequencing, these positive samples revealed 99% homology to B. burgdorferi flgE. CONCLUSION: These results support the existence of borreliosis in Brazil. However, it is unclear whether this borreliosis is caused by a genetically modified B. burgdorferi sensu stricto or by a new species of Borrelia spp.

IL26 gene inactivation in Equidae

Interleukin-26 (IL26) is a member of the IL10 cytokine family. The IL26 gene is located between two other well-known cytokines genes of this family encoding interferon-gamma (IFNG) and IL22 in an evolutionary conserved gene cluster. In contrast to humans and most other mammals, mice lack a functional Il26 gene. We analyzed the genome sequences of other vertebrates for the presence or absence of functional IL26 orthologs and found that the IL26 gene has also become inactivated in several equid species. We detected a one-base pair frameshift deletion in exon 2 of the IL26 gene in the domestic horse (Equus caballus), Przewalski horse (Equus przewalskii) and donkey (Equus asinus). The remnant IL26 gene in the horse is still transcribed and gives rise to at least five alternative transcripts. None of these transcripts share a conserved open reading frame with the human IL26 gene. A comparative analysis across diverse vertebrates revealed that the IL26 gene has also independently been inactivated in a few other mammals, including the African elephant and the European hedgehog. The IL26 gene thus appears to be highly variable, and the conserved open reading frame has been lost several times during mammalian evolution.