Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is similar to 85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

The role of exon shuffling in shaping protein-protein interaction networks

Background: Physical protein-protein interaction (PPI) is a critical phenomenon for the function of most proteins in living organisms and a significant fraction of PPIs are the result of domain-domain interactions. Exon shuffling, intron-mediated recombination of exons from existing genes, is known to have been a major mechanism of domain shuffling in metazoans. Thus, we hypothesized that exon shuffling could have a significant influence in shaping the topology of PPI networks. Results: We tested our hypothesis by compiling exon shuffling and PPI data from six eukaryotic species: Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Cryptococcus neoformans and Arabidopsis thaliana. For all four metazoan species, genes enriched in exon shuffling events presented on average higher vertex degree (number of interacting partners) in PPI networks. Furthermore, we verified that a set of protein domains that are simultaneously promiscuous (known to interact to multiple types of other domains), self-interacting (able to interact with another copy of themselves) and abundant in the genomes presents a stronger signal for exon shuffling. Conclusions: Exon shuffling appears to have been a recurrent mechanism for the emergence of new PPIs along metazoan evolution. In metazoan genomes, exon shuffling also promoted the expansion of some protein domains. We speculate that their promiscuous and self-interacting properties may have been decisive for that expansion.

Conservation of dual-targeted proteins in Arabidopsis and rice points to a similar pattern of gene-family evolution

Gene duplication followed by acquisition of specific targeting information and dual targeting were evolutionary strategies enabling organelles to cope with overlapping functions. We examined the evolutionary trend of dual-targeted single-gene products in Arabidopsis and rice genomes. The number of paralogous proteins encoded by gene families and the dual-targeted orthologous proteins were analysed. The number of dual-targeted proteins and the corresponding gene-family sizes were similar in Arabidopsis and rice irrespective of genome sizes. We show that dual targeting of methionine aminopeptidase, monodehydroascorbate reductase, glutamyl-tRNA synthetase, and tyrosyl-tRNA synthetase was maintained despite occurrence of whole-genome duplications in Arabidopsis and rice as well as a polyploidization followed by a diploidization event (gene loss) in the latter.

Dissemination of bla(IMP-1)-carrying integron In86 among Klebsiella pneumoniae isolates harboring a new trimethoprim resistance gene dfr23

The genetic context of the bla(IMP-1) gene was evaluated in 9 Klebsiella pneumoniae isolates recovered from 2 hospitals in Sao Paulo, Brazil. All isolates harbored a copy of In86 carrying bla(Imp-1), aac(6`)-31, and aadAl. Eight strains from the same hospital also carried another class I integron harboring a new trimethoprim resistance gene (dfr23) that was chromosomally embedded. In86 was likely to be in a 30-kb nontransferable plasmid and was flanked upstream by a sequence identical to one identified in an IMP-1-producing Pseudomonas putida isolate. The bla(IMP-1)-carrying integron In86 was recently reported from nonfermentative bacilli isolated in Sao Paulo. These isolates appear to be the Source of this integron now acquired by K. pneumoniae strains from different hospitals in the same city. Metallo-beta-lactamase production is still rare among Enterobacteriaceae isolates in Brazil, but the acquisition of genetic structures carrying these mobile resistance determinants is worrisome and could lead to an increase in the prevalence of these phenotypes of resistance. (C) 2008 Elsevier Inc. All rights reserved.

Crotalus durissus collilineatus venom gland transcriptome: Analysis of gene expression profile

Crotalus durissus rattlesnakes are responsible for the most lethal cases of snakebites in Brazil. Crotalus durissus collilineatus subspecies is related to a great number of accidents in Southeast and Central West regions, but few studies on its venom composition have been carried out to date. In an attempt to describe the transcriptional profile of the C. durissus collilineatus venom gland, we generated a cDNA library and the sequences obtained could be identified by similarity searches on existing databases. Out of 673 expressed sequence tags (ESTs) 489 produced readable sequences comprising 201 singletons and 47 clusters of two or more ESTs. One hundred and fifty reads (60.5%) produced significant hits to known sequences. The results showed a predominance of toxin-coding ESTs instead of transcripts coding for proteins involved in all cellular functions. The most frequent toxin was crotoxin, comprising 88% of toxin-coding sequences. Crotoxin B, a basic phospholipase A(2) (PLA(2)) subunit of crotoxin, was represented in more variable forms comparing to the non-enzymatic subunit (crotoxin A), and most sequences coding this molecule were identified as CB1 isoform from Crotalus durissus terrificus venom. Four percent of toxin-related sequences in this study were identified as growth factors, comprising five sequences for vascular endothelial growth factor (VEGF) and one for nerve growth factor (NGF) that showed 100% of identity with C. durissus terrificus NGF. We also identified two clusters for metalloprotease from PII class comprising 3% of the toxins, and two for serine proteases, including gyroxin (2.5%). The remaining 2.5% of toxin-coding ESTs represent singletons identified as homologue sequences to cardiotoxin, convulxin, angiotensin-converting enzyme inhibitor and C-type natriuretic peptide, Ohanin, crotamin and PLA(2) inhibitor. These results allowed the identification of the most common classes of toxins in C. durissus collilineatus snake venom, also showing some unknown classes for this subspecies and even for C. durissus species, such as cardiotoxins and VEGF. (C) 2009 Published by Elsevier Masson SAS.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Babesia bovis: Genome size, number of chromosomes and telomeric probe hybridisation

Pulsed field gel electrophoresis of intact chromosomes of Babesia bovis revealed four chromosomes in the haploid genome. A telomere probe, derived from Plasmodium berghei, hybridised to eight SfiI restriction fragments of genomic B. bovis DNA digests indicating the presence of four chromosomes. A small subunit (18S) ribosomal RNA gene probe hybridised to the third chromosome only. The genome size of B. bovis is estimated to be 9.4 million base pairs. The sizes of chromosomes 1, 2, 3 and 4 are estimated to be 1.4, 2.0, 2.8 and 3.2 million base pairs, respectively. (C) 1997 Australian Society for Parasitology. Published by Elsevier Science Ltd.

Persistence and expression of Cotesia congregata polydnavirus in host larvae of the tobacco hornworm, Manduca sexta

The gregarious braconid wasp Cotesia congregata parasitizes host larvae of Manduca sexta, and several other sphingid species. Parasitism induces host immunosuppression due to the disruptive action of the wasp's polydnavirus (PDV) on host blood cells. During the initial stages of parasitism, these cells undergo apoptosis followed by cell clumping, which clears the hemolymph of a large number of cells. In this study, the persistence and expression of Cotesia congregata PDV (CcPDV) were examined using Southern and Nor-them blots, respectively. Digoxygenin-labelled total polydnaviral DNA was used to probe genomic DNA isolated from fat body and brains of hosts with emerged wasps taken 6 days following egress of the parasitoids, and significant cross-hybridization between the host fat body genomic DNA with viral DNA was seen. Thus, the virus persists in the host for the duration of parasitism. even during the post-emergence period, and may even be integrated in the host caterpillar DNA. Viral gene expression was examined using Northern blots and probes to the Cotesia rubecula CrV1 homolog, and the CrV1-like mRNAs were expressed as early as 4 h post-parasitization for at least 72 h and faint hybrization is even seen at the time the wasps eclose. In contrast, in Pieris rapae larvae the CrV1 transcript is expressed only for a brief time, during which time hemocyte function is disrupted. The effect is transitory, and hemocytes regain their normal functions after the parasites emerge as first instars. The genome of CcPDV contains one copy of the CrV1-like homolog as shown on Southern blots of viral genomic DNA. In conjunction with our earlier studies of the PDV-encoded early protein 1, the current work suggests multiple viral transcripts are produced following parasitization of the host. and likely target host hemocytes to induce their apoptosis, thereby preventing encapsulation of the parasitoid's eggs. Whether viral DNAs are integrated in the host's genomic DNA remains to be proven, but our results provide preliminary evidence that viral DNAs are detected in the host's fat body cells examined at the time of wasp ernergence and several days later. (C) 2003 Elsevier Science Ltd. All rights reserved.

The Arabidopsis mutant iop1 exhibits induced over-expression of the plant defensin gene PDF1.2 and enhanced pathogen resistance

Jasmonate and ethylene are concomitantly involved in the induction of the Arabidopsis plant defensin gene PDF1.2. To define genes in the signal transduction pathway leading to the induction of PDF1.2, we screened for-mutants with induced over-expression of a beta-glucuronidase reporter, under the control of the PDF1.2 promoter. One mutant, iop1 (induced over-expressor of PDF1.2) produced small plants that showed induced over-expression of the pathogenesis-related genes PR-3, PR-4 and PR-1,2 (PDF1.2), combined with a down-regulated induction of PR-1 upon pathogen inoculation. The iop1 mutant showed enhanced resistance to a number of necrotrophic pathogens.

Folate-sensitive fragile site FRA10A is due to an expansion of a CGG repeat in a novel gene, FRA10AC1, encoding a nuclear protein

Fragile sites appear visually as nonstaining gaps on chromosomes that are inducible by specific cell culture conditions. Expansion of CGG/ CCG repeats has been shown to be the molecular basis of all five folate-sensitive fragile sites characterized molecularly so far, i.e., FRAXA, FRAXE, FRAXF, FRA11B, and FRA16A. In the present study we have refined the localization of the FRA10A folate-sensitive fragile site by fluorescence in situ hybridization. Sequence analysis of a BAC clone spanning FRA10A identified a single, imperfect, but polymorphic CGG repeat that is part of a CpG island in the 5'UTR of a novel gene named FRA10ACl. The number of CGG repeats varied in the population from 8 to 13. Expansions exceeding 200 repeat units were methylated in all FRA10A fragile site carriers tested. The FRA10ACl gene consists of 19 exons and is transcribed in the centromeric direction from the FRA10A repeat. The major transcript of similar to 1450 nt is ubiquitously expressed and codes for a highly conserved protein, FRA10ACl, of unknown function. Several splice variants leading to alternative 3' ends were identified (particularly in testis). These give rise to FRA10ACl proteins with altered COOH-termini. Immunofluorescence analysis of full-length, recombinant EGFP-tagged FRA10ACl protein showed that it was present exclusively in the nucleoplasm. We show that the expression of FRA10A, in parallel to the other cloned folate-sensitive fragile sites, is caused by an expansion and subsequent methylation of an unstable CGG trinucleotide repeat. Taking advantage of three cSNPs within the FRA10ACl gene we demonstrate that one allele of the gene is not transcribed in a FRA10A carrier. Our data also suggest that in the heterozygous state FRA10A is likely a benign folate-sensitive fragile site. (C) 2004 Elsevier Inc. All rights reserved.

A non-invasive method for silencing gene transcription in honeybees maintained under natural conditions

in the Apis mellifera post-genomic era, RNAi protocols have been used in functional approaches. However, sample manipulation and invasive methods such as injection of double-stranded RNA (dsRNA) can compromise physiology and survival. To circumvent these problems, we developed a non-invasive method for honeybee gene knockdown, using a well-established vitellogenin RNAi system as a model. Second instar larvae received dsRNA for vitellogenin (dsVg-RNA) in their natural diet. For exogenous control, larvae received dsRNA for GFP (dsGFP-RNA). Untreated larvae formed another control group. Around 60% of the treated larvae naturally developed until adult emergence when 0.5 mu g of dsVg-RNA or dsGFP-RNA was offered while no larvae that received 3.0 mu g of dsRNA reached pupal stages. Diet dilution did not affect the removal rates. Viability depends not only on the delivered doses but also on the internal conditions of colonies. The weight of treated and untreated groups showed no statistical differences. This showed that RNAi ingestion did not elicit drastic collateral effects. Approximately 90% of vitellogenin transcripts from 7-day-old workers were silenced compared to controls. A large number of samples are handled in a relatively short time and smaller quantities of RNAi molecules are used compared to invasive methods. These advantages culminate in a versatile and a cost-effective approach. (c) 2008 Elsevier Ltd. All rights reserved.

Mannose-binding lectin gene polymorphism and resistance to therapy in women with recurrent vulvovaginal candidiasis

Precis Women with recurrent vulvovaginal candidiasis (RVC) due to a polymorphism in codon 54 of the MBL2 gene respond better to fluconazole maintenance therapy than do women with other underlying causes. Objective To explain differences in response rates to maintenance therapy with fluconazole in women suffering from RVC by evaluating associations with a polymorphism in the gene coding for mannose-binding lectin (MBL). Design Follow-up study, neted case-control group. Setting Women attending vulvoginitis clinic for RVC. Population Women participating in a multicentric study in Belgium with a degressive dose of fluconazole for RVC (the ReCiDiF trial) were divided into good responders, intermediate responders and nonresponders according to the number of relapses they experienced during therapy. From 109 of these women with adequate follow-up data, vaginal lavage with 2 ml of saline were performed at the moment of a proven acute attack at inclusion in the study, before maintenance treatment was started. A buccal swab was obtained from 55 age-matched women without a history of Candida infections, serving as a control group. Methods Extracted DNA from buccal or vaginal cells was tested for codon 54 MBL2 gene polymorphism by polymerase chain reaction and endonuclease digestion. Main outcome measures Frequency of MBL2 condon 54 allele B in women with optimal or poor response to maintenance therapy in composition with controls. Results Women (n = 109) suffering from RVC were more likely to carry the variant MBL2 codon 54 allele B than control women (20 versus 6.6%, OR 3.4 [95% CI 1.3-8.2], P = 0.01). B alleles were present in 25% of the 36 women not suffering from any recurrence during the maintenance therapy with decreasing doses of fluconazole (OR 4.9 [95% CI 1.9-12.5], P = 0.0007 versus controls), in 20% of the 43 women with sporadic recurrences (OR 3.6 [95% CI 1.4-9.2], P = 0.007 versus controls) and in 15% of the 30 women who had to interrupt the treatment regimen due to frequent relapses (P = 0.097 versus controls). Conclusions The MBL2 codon 54 gene polymorphism is more frequent in Belgian women suffering from RVC than in controls. The presence of the B allele is associated with a superior response to fluconazole maintenance therapy as compared with RVC patients without this polymorphism. We conclude that RVC due to deficient MBL production is more easily helped with antifungal medication than is RVC due to some other mechanism.

DNA binding-independent transcriptional activation of the vascular endothelial growth factor gene (VEGF) by the Myb oncoprotein

Myb is a key transcription factor that can regulate proliferation, differentiation, and apoptosis, predominantly in the haemopoietic system. Abnormal expression of Myb is associated with a number of cancers, both haemopoietic and non-haemopoietic. In order to better understand the role of Myb in normal and tumorigenic processes, we undertook a cDNA array screen to identify genes that are regulated by this factor. In this way, we identified the gene encoding vascular endothelial growth factor (VEGF) as being potentially regulated by the Myb oncoprotein in myeloid cells. To determine whether this was a direct effect on VEGF gene transcription, we examined the activity of the murine VEGF promoter in the presence of either wild-type (WT) or mutant forms of Myb. It was found that WT Myb was able to activate the VEGF promoter and that a minimal promoter region of 120 bp was sufficient to confer Myb responsiveness. Surprisingly, activation of the VEGF promoter was independent of DNA binding by Myb. This was shown by the use of DNA binding-defective Myb mutants and by mutagenesis of a potential Myb-binding site in the minimal promoter. Mutation of Sp1 sites within this region abolished Myb-mediated regulation of a reporter construct, suggesting that Myb DNA binding-independent activation of VEGF expression occurs via these Sp1 binding elements. Regulation of VEGF production by Myb has implications for the potential role of Myb in myeloid leukaemias and in solid tumours where VEGF may be functioning as an autocrine growth factor. (c) 2006 Elsevier Inc. All rights reserved.