Cell cycle synchronization and BrdU incorporation as a tool to study the possible selective elimination of ErbB1 gene in the micronuclei in A549 cells

Lung cancer often exhibits molecular changes, such as the overexpression of the ErbB1 gene that encodes epidermal growth factor receptor (EGFR). ErbB1 amplification and mutation are associated with tumor aggressiveness and low response to therapy. The aim of the present study was to design a schedule to synchronize the cell cycle of A549 cell line (a non-small cell lung cancer) and to analyze the possible association between the micronuclei (MNs) and the extrusion of ErbB1 gene extra-copies. After double blocking, by the process of fetal bovine serum deprivation and vincristine treatment, MNs formation was monitored with 5-bromo-2-deoxyuridine (BrdU) incorporation, which is an S-phase marker. Statistical analyses allowed us to infer that MNs may arise both in mitosis as well as in interphase. The MNs were able to replicate their DNA and this process seemed to be non-synchronous with the main cell nuclei. The presence of ErbB1 gene in the MNs was evaluated by fluorescent in situ hybridization (FISH). ErbB1 sequences were detected in the MNs, but a relation between the MNs formation and extrusion of amplified ErbB1could not be established. The present study sought to elucidate the meaning of MNs formation and its association with the elimination of oncogenes or other amplified sequences from the tumor cells.

Partial characterization of the cloned dihydrofolate reductase gene of Saccharomyces cerevisiae

The cloned dihydrofolate reductase gene of Saccharomyces cerevisiae (DFR 1) is expressed in Escherichia coli. Bacterial strain JF1754 transformed with plasmids containing DFR 1 is at least 5X more resistant to inhibition by the folate antagonist trimethoprim. Expression of yeast DFR 1 in E. coli suggests it is likely that the gene lacks intervening sequences. The 1.8 kbp DNA fragment encoding yeast dhfr activity probably has its own promotor, as the gene is expressed in both orientations in E. coli. Expression of the yeast dhfr gene cloned into M13 viral vectors allowed positive selection of DFR 1 - M13 bacterial transfectants in medium supplemented with trimethoprim. A series of nested deletions generated by nuclease Bal 31 digestion and by restriction endonuclease cleavage of plasmids containing DFR 1 physically mapped the gene to a 930 bp region between the Pst 1 and Sal 1 cut sites. This is consistent with the 21,000 molecular weight attributed to yeast dhfr in previous reports. From preliminary DNA sequence analysis of the dhfr DNA fragment the 3' terminus of DFR 1 was assigned to a position 27 nucleotides from the Eco Rl cut site on the Bam Hi - Eco Rl DNA segment. Several putative yeast transcription termination consensus sequences were identified 3' to the opal stop codon. DFR 1 is expressed in yeast and it confers resistance to the antifolate methotrexate when the gene is present in 2 - 10 copies per cell. Plasmid-dependent resistance to methotrexate is also observed in a rad 6 background although the effect is somewhat less than that conferred to wild-type or rad 18 cells. Integration of DFR 1 into the yeast genome showed an intermediate sensitivity to folate antagonists. This may suggest a gene dosage effect. No change in petite induction in these yeast strains was observed in transformed cells containing yeast dhfr plasmids. The sensitivity of rad 6 , rad 18 and wild-type cell populations to trimethoprim were unaffected by the presence of DFR 1 in transformants. Moreover, trimethoprim did not induce petites in any strain tested, which normally results if dhfr is inhibited by other antifolates such as methotrexate. This may suggest that the dhfr enzyme is not the only possible target of trimethoprim in yeast. rad 6 mutants showed a very low level of spontaneous petite formation. Methotrexate failed to induce respiratory deficient mutants in this strain which suggested that rad 6 might be an obligate grande. However, ethidium bromide induced petites to a level approximately 50% of that exhibited by wild-type and rad 18 strains.

Starvation Induces Expression of the Plant-Adhesin Gene, Mad2, of the Entomopathogenic Fungus Metarhizium robertsii.

Metarhizium robertsii is an entomopathogenic fungus that is additionally plant rhizosphere competent. Two adhesin-encoding gens, Mad1 and Mad2, are involved in insect pathogenesis or plant root colonization, respectively. This study examined differential expression of the Mad genes for M robertsii grown on a variety of insectand plant-related substrates. Mad1 was up regulated in response to insect cuticles and up regulation of Mad2 resulted from root exudates, tomato stems and non-preferred carbohydrates. A time course analysis that compared water, minimal media, and nutrient rich broth revealed Mad2 gene expression increased as nutrient availability decreased. The regulation of Mad2 compared to known stress-related genes (Hsp30, Hsp70 and ssgA) under various stresses (nutrient, pH, osmotic, oxidative, temperature) revealed Mad2 to be generally up regulated by nutrient starvation only. Examination of the Mad2 promoter region revealed two copies of a stress-response element (S TRE) known to be regulated under the general stress response pathway.

The cumate gene-switch: a system for regulated expression in mammalian cells

Affiliation: Sophie Broussau, Amelie Pilotte & Bernard Massie : Départment de microbiologie et immunologie, Faculté de médecine, Université de Montréal

Regulation of gene expression in the dinoflagellate Lingulodinium polyedrum

Les dinoflagellés sont des eucaryotes unicellulaires que l’on retrouve autant en eau douce qu’en milieu marin. Ils sont particulièrement connus pour causer des fleurs d’algues toxiques nommées ‘marée-rouge’, ainsi que pour leur symbiose avec les coraux et pour leur importante contribution à la fixation du carbone dans les océans. Au point de vue moléculaire, ils sont aussi connus pour leur caractéristiques nucléaires uniques, car on retrouve généralement une quantité immense d’ADN dans leurs chromosomes et ceux-ci sont empaquetés et condensés sous une forme cristalline liquide au lieu de nucléosomes. Les gènes encodés par le noyau sont souvent présents en multiples copies et arrangés en tandem et aucun élément de régulation transcriptionnelle, y compris la boite TATA, n’a encore été observé. L’organisation unique de la chromatine des dinoflagellés suggère que différentes stratégies sont nécessaires pour contrôler l’expression des gènes de ces organismes. Dans cette étude, j’ai abordé ce problème en utilisant le dinoflagellé photosynthétique Lingulodinium polyedrum comme modèle. L. polyedrum est d’un intérêt particulier, car il a plusieurs rythmes circadiens (journalier). À ce jour, toutes les études sur l’expression des gènes lors des changements circadiens ont démontrées une régulation à un niveau traductionnel. Pour mes recherches, j’ai utilisé les approches transcriptomique, protéomique et phosphoprotéomique ainsi que des études biochimiques pour donner un aperçu de la mécanique de la régulation des gènes des dinoflagellés, ceci en mettant l’accent sur l’importance de la phosphorylation du système circadien de L. polyedrum. L’absence des protéines histones et des nucléosomes est une particularité des dinoflagellés. En utilisant la technologie RNA-Seq, j’ai trouvé des séquences complètes encodant des histones et des enzymes modifiant les histones. L polyedrum exprime donc des séquences conservées codantes pour les histones, mais le niveau d’expression protéique est plus faible que les limites de détection par immunodétection de type Western. Les données de séquençage RNA-Seq ont également été utilisées pour générer un transcriptome, qui est une liste des gènes exprimés par L. polyedrum. Une recherche par homologie de séquences a d’abord été effectuée pour classifier les transcrits en diverses catégories (Gene Ontology; GO). Cette analyse a révélé une faible abondance des facteurs de transcription et une surprenante prédominance, parmi ceux-ci, des séquences à domaine Cold Shock. Chez L. polyedrum, plusieurs gènes sont répétés en tandem. Un alignement des séquences obtenues par RNA-Seq avec les copies génomiques de gènes organisés en tandem a été réalisé pour examiner la présence de transcrits polycistroniques, une hypothèse formulée pour expliquer le manque d’élément promoteur dans la région intergénique de la séquence de ces gènes. Cette analyse a également démontré une très haute conservation des séquences codantes des gènes organisés en tandem. Le transcriptome a également été utilisé pour aider à l’identification de protéines après leur séquençage par spectrométrie de masse, et une fraction enrichie en phosphoprotéines a été déterminée comme particulièrement bien adapté aux approches d’analyse à haut débit. La comparaison des phosphoprotéomes provenant de deux périodes différentes de la journée a révélée qu’une grande partie des protéines pour lesquelles l’état de phosphorylation varie avec le temps est reliées aux catégories de liaison à l’ARN et de la traduction. Le transcriptome a aussi été utilisé pour définir le spectre des kinases présentes chez L. polyedrum, qui a ensuite été utilisé pour classifier les différents peptides phosphorylés qui sont potentiellement les cibles de ces kinases. Plusieurs peptides identifiés comme étant phosphorylés par la Casein Kinase 2 (CK2), une kinase connue pour être impliquée dans l’horloge circadienne des eucaryotes, proviennent de diverses protéines de liaison à l’ARN. Pour évaluer la possibilité que quelques-unes des multiples protéines à domaine Cold Shock identifiées dans le transcriptome puissent moduler l’expression des gènes de L. polyedrum, tel qu’observé chez plusieurs autres systèmes procaryotiques et eucaryotiques, la réponse des cellules à des températures froides a été examinée. Les températures froides ont permis d’induire rapidement un enkystement, condition dans laquelle ces cellules deviennent métaboliquement inactives afin de résister aux conditions environnementales défavorables. Les changements dans le profil des phosphoprotéines seraient le facteur majeur causant la formation de kystes. Les phosphosites prédits pour être phosphorylés par la CK2 sont la classe la plus fortement réduite dans les kystes, une découverte intéressante, car le rythme de la bioluminescence confirme que l’horloge a été arrêtée dans le kyste.

Mean Squared Residue Based Biclustering Algorithms for the Analysis of Gene Expression Data

Computational Biology is the research are that contributes to the analysis of biological data through the development of algorithms which will address significant research problems.The data from molecular biology includes DNA,RNA ,Protein and Gene expression data.Gene Expression Data provides the expression level of genes under different conditions.Gene expression is the process of transcribing the DNA sequence of a gene into mRNA sequences which in turn are later translated into proteins.The number of copies of mRNA produced is called the expression level of a gene.Gene expression data is organized in the form of a matrix. Rows in the matrix represent genes and columns in the matrix represent experimental conditions.Experimental conditions can be different tissue types or time points.Entries in the gene expression matrix are real values.Through the analysis of gene expression data it is possible to determine the behavioral patterns of genes such as similarity of their behavior,nature of their interaction,their respective contribution to the same pathways and so on. Similar expression patterns are exhibited by the genes participating in the same biological process.These patterns have immense relevance and application in bioinformatics and clinical research.Theses patterns are used in the medical domain for aid in more accurate diagnosis,prognosis,treatment planning.drug discovery and protein network analysis.To identify various patterns from gene expression data,data mining techniques are essential.Clustering is an important data mining technique for the analysis of gene expression data.To overcome the problems associated with clustering,biclustering is introduced.Biclustering refers to simultaneous clustering of both rows and columns of a data matrix. Clustering is a global whereas biclustering is a local model.Discovering local expression patterns is essential for identfying many genetic pathways that are not apparent otherwise.It is therefore necessary to move beyond the clustering paradigm towards developing approaches which are capable of discovering local patterns in gene expression data.A biclusters is a submatrix of the gene expression data matrix.The rows and columns in the submatrix need not be contiguous as in the gene expression data matrix.Biclusters are not disjoint.Computation of biclusters is costly because one will have to consider all the combinations of columans and rows in order to find out all the biclusters.The search space for the biclustering problem is 2 m+n where m and n are the number of genes and conditions respectively.Usually m+n is more than 3000.The biclustering problem is NP-hard.Biclustering is a powerful analytical tool for the biologist.The research reported in this thesis addresses the problem of biclustering.Ten algorithms are developed for the identification of coherent biclusters from gene expression data.All these algorithms are making use of a measure called mean squared residue to search for biclusters.The objective here is to identify the biclusters of maximum size with the mean squared residue lower than a given threshold. All these algorithms begin the search from tightly coregulated submatrices called the seeds.These seeds are generated by K-Means clustering algorithm.The algorithms developed can be classified as constraint based,greedy and metaheuristic.Constarint based algorithms uses one or more of the various constaints namely the MSR threshold and the MSR difference threshold.The greedy approach makes a locally optimal choice at each stage with the objective of finding the global optimum.In metaheuristic approaches particle Swarm Optimization(PSO) and variants of Greedy Randomized Adaptive Search Procedure(GRASP) are used for the identification of biclusters.These algorithms are implemented on the Yeast and Lymphoma datasets.Biologically relevant and statistically significant biclusters are identified by all these algorithms which are validated by Gene Ontology database.All these algorithms are compared with some other biclustering algorithms.Algorithms developed in this work overcome some of the problems associated with the already existing algorithms.With the help of some of the algorithms which are developed in this work biclusters with very high row variance,which is higher than the row variance of any other algorithm using mean squared residue, are identified from both Yeast and Lymphoma data sets.Such biclusters which make significant change in the expression level are highly relevant biologically.

Nuclear organisation and epigenetic regulation of gene expression in Dictyostelium discoideum

A series of vectors for the over-expression of tagged proteins in Dictyostelium were designed, constructed and tested. These vectors allow the addition of an N- or C-terminal tag (GFP, RFP, 3xFLAG, 3xHA, 6xMYC and TAP) with an optimized polylinker sequence and no additional amino acid residues at the N or C terminus. Different selectable markers (Blasticidin and gentamicin) are available as well as an extra chromosomal version; these allow copy number and thus expression level to be controlled, as well as allowing for more options with regard to complementation, co- and super-transformation. Finally, the vectors share standardized cloning sites, allowing a gene of interest to be easily transfered between the different versions of the vectors as experimental requirements evolve. The organisation and dynamics of the Dictyostelium nucleus during the cell cycle was investigated. The centromeric histone H3 (CenH3) variant serves to target the kinetochore to the centromeres and thus ensures correct chromosome segregation during mitosis and meiosis. A number of Dictyostelium histone H3-domain containing proteins as GFP-tagged fusions were expressed and it was found that one of them functions as CenH3 in this species. Like CenH3 from some other species, Dictyostelium CenH3 has an extended N-terminal domain with no similarity to any other known proteins. The targeting domain, comprising α-helix 2 and loop 1 of the histone fold is required for targeting CenH3 to centromeres. Compared to the targeting domain of other known and putative CenH3 species, Dictyostelium CenH3 has a shorter loop 1 region. The localisation of a variety of histone modifications and histone modifying enzymes was examined. Using fluorescence in situ hybridisation (FISH) and CenH3 chromatin-immunoprecipitation (ChIP) it was shown that the six telocentric centromeres contain all of the DIRS-1 and most of the DDT-A and skipper transposons. During interphase the centromeres remain attached to the centrosome resulting in a single CenH3 cluster which also contains the putative histone H3K9 methyltransferase SuvA, H3K9me3 and HP1 (heterochromatin protein 1). Except for the centromere cluster and a number of small foci at the nuclear periphery opposite the centromeres, the rest of the nucleus is largely devoid of transposons and heterochromatin associated histone modifications. At least some of the small foci correspond to the distal telomeres, suggesting that the chromosomes are organised in a Rabl-like manner. It was found that in contrast to metazoans, loading of CenH3 onto Dictyostelium centromeres occurs in late G2 phase. Transformation of Dictyostelium with vectors carrying the G418 resistance cassette typically results in the vector integrating into the genome in one or a few tandem arrays of approximately a hundred copies. In contrast, plasmids containing a Blasticidin resistance cassette integrate as single or a few copies. The behaviour of transgenes in the nucleus was examined by FISH, and it was found that low copy transgenes show apparently random distribution within the nucleus, while transgenes with more than approximately 10 copies cluster at or immediately adjacent to the centromeres in interphase cells regardless of the actual integration site along the chromosome. During mitosis the transgenes show centromere-like behaviour, and ChIP experiments show that transgenes contain the heterochromatin marker H3K9me2 and the centromeric histone variant H3v1. This clustering, and centromere-like behaviour was not observed on extrachromosomal transgenes, nor on a line where the transgene had integrated into the extrachromosomal rDNA palindrome. This suggests that it is the repetitive nature of the transgenes that causes the centromere-like behaviour. A Dictyostelium homolog of DET1, a protein largely restricted to multicellular eukaryotes where it has a role in developmental regulation was identified. As in other species Dictyostelium DET1 is nuclear localised. In ChIP experiments DET1 was found to bind the promoters of a number of developmentally regulated loci. In contrast to other species where it is an essential protein, loss of DET1 is not lethal in Dictyostelium, although viability is greatly reduced. Loss of DET1 results in delayed and abnormal development with enlarged aggregation territories. Mutant slugs displayed apparent cell type patterning with a bias towards pre-stalk cell types.

Gene conversion and evolution of Xq28 duplicons involved in recurring inversions causing severe hemophilia A

Inversions breaking the 1041 bp int1h-1 or the 9.5-kb int22h-1 sequence of the F8 gene cause hemophilia A in 1/30,000 males. These inversions are due to homologous recombination between the above sequences and their inverted copies on the same DNA molecule, respectively, int1h-2 and int22h-2 or int22h-3. We find that (1) int1h and int22h duplicated more than 25 million years ago; (2) the identity of the copies (>99%) of these sequences in humans and other primates is due to gene conversion; (3) gene conversion is most frequent in the internal regions of int22h; (4) breakpoints of int22h-related inversions also tend to involve the internal regions of int22h; (5) sequence variations in a sample of human X chromosomes defined eight haplotypes of int22h-1 and 27 of int22h-2 plus int22h-3; (6) the latter two sequences, which lie, respectively, 500 and 600 kb telomeric to int22h-1 are five-fold more identical when in cis than when in trans, thus suggesting that gene conversion may be predominantly intrachromosomal; (7) int1h, int22h, and flanking sequences evolved at a rate of about 0.1% substitutions per million years during the divergence between humans and other primates, except for int1h during the human-chimpanzee divergence, when its rate of evolution was significantly lower. This is reminiscent of the slower evolution of palindrome arms in the male specific regions of the Y chromosome and we propose, as an explanation, that intrachromosomal gene conversion and cosegregation of the duplicated regions favors retention of the ancestral sequence and thus reduces the evolution rate.

Neurotrophic gene polymorphisms and response to psychological therapy

Therapygenetics, the study of genetic determinants of response to psychological therapies, is in its infancy. Here, we investigate whether single-nucleotide polymorphisms in nerve growth factor (NGF) (rs6330) and brain-derived neutrotrophic factor (BDNF) (rs6265) genes predict the response to cognitive behaviour therapy (CBT). Neurotrophic genes represent plausible candidate genes: they are implicated in synaptic plasticity, response to stress, and are widely expressed in brain areas involved in mood and cognition. Allelic variation at both loci has shown associations with anxiety-related phenotypes. A sample of 374 anxiety-disordered children with white European ancestry was recruited from clinics in Reading, UK, and in Sydney, Australia. Participants received manualised CBT treatment and DNA was collected from buccal cells using cheek swabs. Treatment response was assessed at post-treatment and follow-up time points. We report first evidence that children with one or more copies of the T allele of NGF rs6330 were significantly more likely to be free of their primary anxiety diagnosis at follow-up (OR=0.60 (0.42–0.85), P=0.005). These effects remained even when other clinically relevant covariates were accounted for (OR=0.62 (0.41–0.92), P=0.019). No significant associations were observed between BDNF rs6265 and response to psychological therapy. These findings demonstrate that knowledge of genetic markers has the potential to inform clinical treatment decisions for psychotherapeutic interventions.

Ace gene dosage influences the development of renovascular hypertension

P>1. Clinical and experimental evidence highlights the importance of the renin-angiotensin system in renovascular hypertension. Furthermore, genetic factors affecting angiotensin-converting enzyme (ACE) could influence the development of renovascular hypertension. 2. To test the effect of small gene perturbations on the development of renovascular hypertension, mice harbouring two or three copies of the Ace gene were submitted to 4 weeks of two-kidney, one-clip (2K1C) hypertension. Blood pressure (BP), cardiac hypertrophy, baroreflex sensitivity and blood pressure and heart rate variability were assessed and compared between the different groups. 3. The increase in BP induced by 2K1C was higher in mice with three copies of the Ace gene compared with mice with only two copies (46 vs 23 mmHg, respectively). Moreover, there was a 3.8-fold increase in the slope of the left ventricle mass/BP relationship in mice with three copies of the Ace gene. Micewith three copies of the Ace gene exhibited greater increases in cardiac and serum ACE activity than mice with only two copies of the gene. Both baroreflex bradycardia and tachycardia were significantly depressed in mice with three copies of the Ace gene after induction of 2K1C hypertension. The variance in basal systolic BP was greater in mice with three copies of the Ace gene after 2K1C hypertension compared with those with only two copies of the gene (106 vs 54%, respectively). In addition, the low-frequency component of the pulse interval was higher mice with three copies of the Ace gene after 2K1C hypertension compared with those with only two (168 vs 86%, respectively). Finally, in mice with three copies of the Ace gene, renovascular hypertension induced a 6.1-fold increase in the sympathovagal balance compared with a 3.2-fold increase in mice with only two copies of the gene. 4. Collectively, these data provide direct evidence that small genetic disturbances in ACE levels per se have an influence on haemodynamic, cardiac mass and autonomic nervous system responses in mice under pathological perturbation.

Protein variation in blood-dwelling schistosome worms generated by differential splicing of micro-exon gene transcripts

Schistosoma mansoni is a well-adapted blood-dwelling parasitic helminth, persisting for decades in its human host despite being continually exposed to potential immune attack. Here, we describe in detail micro-exon genes (MEG) in S. mansoni, some present in multiple copies, which represent a novel molecular system for creating protein variation through the alternate splicing of short (<= 36 bp) symmetric exons organized in tandem. Analysis of three closely related copies of one MEG family allowed us to trace several evolutionary events and propose a mechanism for micro-exon generation and diversification. Microarray experiments show that the majority of MEGs are up-regulated in life cycle stages associated with establishment in the mammalian host after skin penetration. Sequencing of RT-PCR products allowed the description of several alternate splice forms of micro-exon genes, highlighting the potential use of these transcripts to generate a complex pool of protein variants. We obtained direct evidence for the existence of such pools by proteomic analysis of secretions from migrating schistosomula and mature eggs. Whole-mount in situ hybridization and immunolocalization showed that MEG transcripts and proteins were restricted to glands or epithelia exposed to the external environment. The ability of schistosomes to produce a complex pool of variant proteins aligns them with the other major groups of blood parasites, but using a completely different mechanism. We believe that our data open a new chapter in the study of immune evasion by schistosomes, and their ability to generate variant proteins could represent a significant obstacle to vaccine development.

A novel system for large-scale gene expression analysis: bacterial colonies array

In the present work, we report the use of bacterial colonies to optimize macroarray technique. The devised system is significantly cheaper than other methods available to detect large-scale differential gene expression. Recombinant Escherichia coli clones containing plasmid-encoded copies of 4,608 individual expressed sequence tag (ESTs) were robotically spotted onto nylon membranes that were incubated for 6 and 12 h to allow the bacteria to grow and, consequently, amplify the cloned ESTs. The membranes were then hybridized with a beta-lactamase gene specific probe from the recombinant plasmid and, subsequently, phosphorimaged to quantify the microbial cells. Variance analysis demonstrated that the spot hybridization signal intensity was similar for 3,954 ESTs (85.8%) after 6 h of bacterial growth. Membranes spotted with bacteria colonies grown for 12 h had 4,017 ESTs (87.2%) with comparable signal intensity but the signal to noise ratio was fivefold higher. Taken together, the results of this study indicate that it is possible to investigate large-scale gene expression using macroarrays based on bacterial colonies grown for 6 h onto membranes.

Associations between LEP, DGAT1 and FABP4 gene polymorphisms and carcass and meat traits in Nelore and crossbred beef cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evolutionary dynamics of rRNA gene clusters in cichlid fish

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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.