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.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Genetic variability of porcine parvovirus isolates revealed by analysis of partial sequences of the structural coding gene VP2

The 3'-terminal 853 nt (and the putative 283 aa) sequence of the VP2-encoding gene from 29 field strains of porcine parvovirus (PPV) were determined and compared both to each other and with other published sequences. Sequences were examined using maximum-parsimony and statistical analyses for nucleotide diversity and sequence variability. Among the nucleotide sequences of the PPV field strains, 26 polymorphic sites were encountered; 22 polymorphic sites were detected in the putative amino acid sequence. Mapping polymorphic sites of protein data onto the three-dimensional (3D) structure of PPV VP2 revealed that almost all substitutions were located on the external surface of the viral capsid. Mapping amino acid substitutions to the alignment between PPV VP2 sequences and the 3D structure of canine parvovirus (CPV) capsid, many PPV substitutions were observed to map to regions of recognized antigenicity and/or to contain phenotypically important residues for CPV and other parvoviruses. In spite of the high sequence similarity, genetic analysis has shown the existence of at least two virus lineages among the samples. In conclusion, these results highlight the need for close surveillance on PPV genetic drift, with an assessment of its potential ability to modify the antigenic make-up of the virus.

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 ∼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. ©FUNPEC-RP.

Plant protein-coding gene families: emerging bioinformatics approaches

Protein-coding gene families are sets of similar genes with a shared evolutionary origin and, generally, with similar biological functions. In plants, the size and role of gene families has been only partially addressed. However, suitable bioinformatics tools are being developed to cluster the enormous number of sequences currently available in databases. Specifically, comparative genomic databases promise to become powerful tools for gene family annotation in plant clades. In this review, I evaluate the data retrieved from various gene family databases, the ease with which they can be extracted and how useful the extracted information is.

The U5 RNA of trypanosomes deviates from the canonical U5 RNA: The Leptomonas collosoma U5 RNA and its coding gene

Fractionation of the abundant small ribonucleoproteins (RNPs) of the trypanosomatid Leptomonas collosoma revealed the existence of a group of unidentified small RNPs that were shown to fractionate differently than the well-characterized trans-spliceosomal RNPs. One of these RNAs, an 80-nt RNA, did not possess a trimethylguanosine (TMG) cap structure but did possess a 5′ phosphate terminus and an invariant consensus U5 snRNA loop 1. The gene coding for the RNA was cloned, and the coding region showed 55% sequence identity to the recently described U5 homologue of Trypanosoma brucei [Dungan, J. D., Watkins, K. P. & Agabian, N. (1996) EMBO J. 15, 4016–4029]. The L. collosoma U5 homologue exists in multiple forms of RNP complexes, a 10S monoparticle, and two subgroups of 18S particles that either contain or lack the U4 and U6 small nuclear RNAs, suggesting the existence of a U4/U6⋅U5 tri-small nuclear RNP complex. In contrast to T. brucei U5 RNA (62 nt), the L. collosoma homologue is longer (80 nt) and possesses a second stem–loop. Like the trypanosome U3, U6, and 7SL RNA genes, a tRNA gene coding for tRNACys was found 98 nt upstream to the U5 gene. A potential for base pair interaction between U5 and SL RNA in the 5′ splice site region (positions −1 and +1) and downstream from it is proposed. The presence of a U5-like RNA in trypanosomes suggests that the most essential small nuclear RNPs are ubiquitous for both cis- and trans-splicing, yet even among the trypanosomatids the U5 RNA is highly divergent.

Pseudo-messenger RNA: Phantoms of the transcriptome

The mammalian transcriptome harbours shadowy entities that resist classification and analysis. In analogy with pseudogenes, we define pseudo-messenger RNA to be RNA molecules that resemble protein- coding mRNA, but cannot encode full-length proteins owing to disruptions of the reading frame. Using a rigorous computational pipeline, which rules out sequencing errors, we identify 10,679 pseudo - messenger RNAs ( approximately half of which are transposonassociated) among the 102,801 FANTOM3 mouse cDNAs: just over 10% of the FANTOM3 transcriptome. These comprise not only transcribed pseudogenes, but also disrupted splice variants of otherwise protein- coding genes. Some may encode truncated proteins, only a minority of which appear subject to nonsense- mediated decay. The presence of an excess of transcripts whose only disruptions are opal stop codons suggests that there are more selenoproteins than currently estimated. We also describe compensatory frameshifts, where a segment of the gene has changed frame but remains translatable. In summary, we survey a large class of non- standard but potentially functional transcripts that are likely to encode genetic information and effect biological processes in novel ways. Many of these transcripts do not correspond cleanly to any identifiable object in the genome, implying fundamental limits to the goal of annotating all functional elements at the genome sequence level.

Differential diagnosis of oocysts of Hammondia-like organisms of dogs and cats by PCR-RFLP analysis of 70-kilodalton heat shock protein (HSP70) gene

Organisms of the genera Toxoplasma, Hammondia and Neospora, the Hammondia-like organisms, are closely related coccidian with similarly sized oocysts. Therefore, a diagnosis based on microscopy of oocysts in feces is not a method of choice for species identification of these important parasites. In this paper, we present a polymerase chain reaction coupled with restriction fragment length polymorphism (PCR-RFLP) method to differentially diagnose oocysts of Toxoplasma gondii from oocyst of Hammondia hammondi. Another PCR-RFLP was designed to differentiate oocysts of Hammondia heydorni from oocysts of Neospora spp. Both PCR-RFLP are based on nucleotide sequences of the Hsp70 coding gene. In conclusion, we presented two alternative molecular diagnostic assays that can be successfully applied for the differentiation of oocysts of Hammondia-like organisms shed by felids and canids.

Crab-eating fox (Cerdocyon thous), a South American canid, as a definitive host for Hammondia heydorni

Hammondia heydorni is a cyst forming coccidia closely related to other apicomplexans, such as Toxoplasma gondii, Neospora caninum and Hammondia hammondi with a two-host life cycle. Dogs and other canids as red foxes (Vulpes vulpes) and coyotes (Canis latrans) may serve as definitive hosts for H. heydorni. Sporulated oocysts are infective for cattle, sheep and goats, which may serve as intermediate hosts. Herein, we describe the ability of crab-eating fox (Cerdocyon thous), a wild carnivore that is commonly found from northern Argentina to northern South America, to serve as definitive host of H. heydorni. The whole masseter muscle and brain from two 2-year-old bovines were collected, minced and pooled together for the fox infection. The bovine pooled tissues were equally administered to four foxes, in two consecutive days. Two foxes shed subspherical unsporulated oocysts measuring 10-15 mu m, after 8 and 9 days post-infection, respectively. One of the foxes eliminated oocysts for 5 days, while the other fox shed oocysts for 9 days. A DNA sample of oocysts detected at each day of oocyst elimination was tested by two PCRs, one of them carried out employing primers directed to the common toxoplasmatiid 18S and 5.8S ribosomal RNA coding genes (PCR-ITS1) and the other based on heat-shock protein 70 kDa coding gene (PCR-HSP70). These samples were also submitted to a N. caninum specific nested-PCR protocol based on a N. caninum specific gene (Nc5-nPCR). All of them were positive by PCR-ITS1 and PCR-HSP70 but negative by Nc5-nPCR. The PCR-ITS1 and PCR-HSP70 nucleotide sequences amplified from the oocysts shed by the foxes revealed 100% identity with homologous sequences of H. heydorni. In conclusion, it is clear that H. heydorni also uses the crab-eating fox as a definitive host. The crab-eating fox is usually reported to live in close contact with livestock in several regions of Brazil. Therefore, it is reasonable to infer that such carnivores may play an important role in the sylvatic and domestic cycles of H. heydorni infection. (C) 2009 Elsevier B.V. All rights reserved.

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Identifying protein-coding genes in genomic sequences

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Identifying protein-coding genes in genomic sequences.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

A universal primer for isolation of fragments of a gene encoding phytoene desaturase for use in virus-induced gene silencing (VIGS) studies

We have been using Virus-Induced Gene Silencing (VIGS) to test the function of genes that are candidates for involvement in floral senescence. Although VIGS is a powerful tool for assaying the effects of gene silencing in plants, relatively few taxa have been studied using this approach, and most that have are in the Solanaceae. We typically use silencing of phytoene desaturase (PDS) in preliminary tests of the feasibility of using VIGS. Silencing this gene, whose product is involved in carotene biosynthesis, results in a characteristic photobleaching phenotype in the leaves. We have found that efficient silencing requires the use of fragments that are more than 90% homologous to the target gene. To simplify testing the effectiveness of VIGS in a range of species, we designed a set of universal primers to a region of the PDS gene that is highly conserved among species, and that therefore allows an investigator to isolate a fragment of the homologous PDS gene from the species of interest. We report the sequences of these primers and the results of VIGS experiments in horticultural species from the Asteraceae, Leguminosae, Balsaminaceae and Solanaceae.

Dysregulation of REST-regulated coding and non-coding RNAs in a cellular model of Huntington's disease

Huntingtin (Htt) protein interacts with many transcriptional regulators, with widespread disruption to the transcriptome in Huntington's disease (HD) brought about by altered interactions with the mutant Htt (muHtt) protein. Repressor Element-1 Silencing Transcription Factor (REST) is a repressor whose association with Htt in the cytoplasm is disrupted in HD, leading to increased nuclear REST and concomitant repression of several neuronal-specific genes, including brain-derived neurotrophic factor (Bdnf). Here, we explored a wide set of HD dysregulated genes to identify direct REST targets whose expression is altered in a cellular model of HD but that can be rescued by knock-down of REST activity. We found many direct REST target genes encoding proteins important for nervous system development, including a cohort involved in synaptic transmission, at least two of which can be rescued at the protein level by REST knock-down. We also identified several microRNAs (miRNAs) whose aberrant repression is directly mediated by REST, including miR-137, which has not previously been shown to be a direct REST target in mouse. These data provide evidence of the contribution of inappropriate REST-mediated transcriptional repression to the widespread changes in coding and non-coding gene expression in a cellular model of HD that may affect normal neuronal function and survival.