LUAD transcriptomic profile analysis of d-limonene and potential lncRNA chemopreventive target.

d-Limonene, a type of natural extract obtained from citrus oils, was reported to have anti-cancer effects and be well-tolerated by cancer patients. Despite arousing interest as a cancer chemopreventive substance, the transcriptomic profile of d-limonene in humans is poorly understood. Based on the results of the transcriptomic profiling, a lncRNA named protein disulfide isomerase family A member three pseudogene (PDIA3P1) was found to be regulated by d-limonene. PDIA3P1 is an oncogene verified by three lung adenocarcinoma (LUAD) datasets. The knockdown of PDIA3P1 with siRNA decreased the viability, invasion, migration, and proliferation of LUAD cells. Based on The Cancer Genome Atlas (TCGA) LUAD datasets, PDIA3P1 regulates functions and pathways mainly including lipid metabolism, immunity, and the change of the chromosome structure. This study comprehensively performs the transcriptomic analysis of the d-limonene regulation on LUAD, and reveals that PDIA3P1 may be the mediator in helping d-limonene to prevent and suppress LUAD via lipid metabolism, immunity pathway, and the change in the chromosome structure.

Ptr1 evolved convergently with RPS2 and Mr5 to mediate recognition of AvrRpt2 in diverse solanaceous species.

The Ptr1 (Pseudomonas tomato race 1) locus in Solanum lycopersicoides confers resistance to strains of Pseudomonas syringae pv. tomato expressing AvrRpt2 and Ralstonia pseudosolanacearum expressing RipBN. Here we describe the identification and phylogenetic analysis of the Ptr1 gene. A single recombinant among 585 F2 plants segregating for the Ptr1 locus was discovered that narrowed the Ptr1 candidates to eight nucleotide-binding leucine-rich repeat protein (NLR)-encoding genes. From analysis of the gene models in the S. lycopersicoides genome sequence and RNA-Seq data, two of the eight genes emerged as the strongest candidates for Ptr1. One of these two candidates was found to encode Ptr1 based on its ability to mediate recognition of AvrRpt2 and RipBN when it was transiently expressed with these effectors in leaves of Nicotiana glutinosa. The ortholog of Ptr1 in tomato and in Solanum pennellii is a pseudogene. However, a functional Ptr1 ortholog exists in Nicotiana benthamiana and potato, and both mediate recognition of AvrRpt2 and RipBN. In apple and Arabidopsis, recognition of AvrRpt2 is mediated by the Mr5 and RPS2 proteins, respectively. Phylogenetic analysis places Ptr1 in a distinct clade compared with Mr5 and RPS2, and it therefore appears to have arisen by convergent evolution for recognition of AvrRpt2.

Complete plastome sequencing of Allium paradoxum reveals unusual rearrangements and the loss of the ndh genes as compared to Allium ursinum and other onions.

In this work the complete chloroplast DNAs of Allium paradoxum and Allium ursinum, two edible species of Allium subg. Amerallium (the first lineage), were sequenced, assembled, annotated, and compared with complete Allium plastomes of the second and third evolutionary lines from GenBank database. The A. ursinum plastome contains 90 predicted genes (81 unique) including 5 pseudogenes, while A. paradoxum has 88 predicted genes (79 unique) including 19 pseudogenes. The comparative analysis has revealed that the A. paradoxum plastome differs markedly from those of other species. Due to many deletions, the A. paradoxum plastome is the shortest of known for Allium species, being only 145,819 bp long. The most prominent distinctions are (1) a 4825 bp long local inversion that spans from the ndhE to the rpl32 gene in the small single copy region and (2) pseudogenization, or the loss of all NADH-genes. In contrast, the plastome of A. ursinum - a species from the first evolutionary line (as well as A. paradoxum) - resembles the Allium species of the second and third evolutionary lines, showing no large rearrangements or discrepancies in gene content. It is unclear yet whether only A. paradoxum was affected by some evolutionary events or its close relatives from both sect. Briseis and other sections of Amerallium were altered as well. We speculate the sunlight-intolerant, shade-loving nature of A. paradoxum and the impairment of the ndh genes in its plastome could be interrelated phenomena.

Genome-wide DNA methylation changes in placenta tissues associated with small for gestational age newborns; cohort study in the Chinese population.

Aim: To investigate DNA methylation changes in placenta tissues associated with small for gestational age (SGA). Materials & methods: A prospective cohort study consisting of 1292 pregnant women from China (including 39 SGA with placenta tissues) was performed, microarray and pyrosequencing were conducted. Results: Total 2012 methylation variable positions stood out from all probes (p < 0.05; Δß > 0.2). In SGA cases, a CpG site within ANKRD20B showed lower methylation level (p = 0.032) than appropriate for gestational age in validation cohort. Five sites within FAM198A (p = 0.047, 0.050, 0.039, 0.026 and 0.043, respectively) had a reduced methylation in male newborns whose mother had preconception folic acid supplementation. Conclusion: DNA methylation changes in placenta tissues may be associated with SGA, maternal preconception folic acid supplementation status and also be fetal sex-specific.

A Carbonic Anhydrase Pseudogene Sensitizes Select Brucella Lineages to Low CO2 Tension.

Brucella spp. are intracellular pathogens that cause a disease known as brucellosis. Though the genus is highly monomorphic at the genetic level, species have animal host preferences and some defining physiologic characteristics. Of note is the requirement for CO2 supplementation to cultivate particular species, which confounded early efforts to isolate B. abortus from diseased cattle. Differences in the capacity of Brucella species to assimilate CO2 are determined by mutations in the carbonic anhydrase gene, bcaA Ancestral single-nucleotide insertions in bcaA have resulted in frameshifted pseudogenes in B. abortus and B. ovis lineages, which underlie their inability to grow under the low CO2 tension of a standard atmosphere. Incubation of wild-type B. ovis in air selects for mutations that "rescue" a functional bcaA reading frame, which enables growth under low CO2 and enhances the growth rate under high CO2 Accordingly, we show that heterologous expression of functional Escherichia coli carbonic anhydrases enables B. ovis growth in air. Growth of B. ovis is acutely sensitive to a reduction in CO2 tension, while frame-rescued B. ovis mutants are insensitive to CO2 shifts. B. ovis initiates a gene expression program upon CO2 downshift that resembles the stringent response and results in transcriptional activation of its type IV secretion system. Our study provides evidence that loss-of-function insertion mutations in bcaA sensitize the response of B. ovis and B. abortus to reduced CO2 tension relative to that of other Brucella lineages. CO2-dependent starvation and virulence gene expression programs in these species may influence persistence or transmission in natural hosts.IMPORTANCEBrucella spp. are highly related, but they exhibit differences in animal host preference that must be determined by genome sequence differences. B. ovis and the majority of B. abortus strains require high CO2 tension to be cultivated in vitro and harbor conserved insertional mutations in the carbonic anhydrase gene, bcaA, which underlie this trait. Mutants that grow in a standard atmosphere, first reported nearly a century ago, are easily selected in the laboratory. These mutants harbor varied indel polymorphisms in bcaA that restore its consensus reading frame and rescue its function. Loss of bcaA function has evolved independently in the B. ovis and B. abortus lineages and results in a dramatically increased sensitivity to CO2 limitation.

Dissimilar evolutionary histories of two resistance gene families in the genus Solanum.

Genomic analyses have shown that most genes in eukaryotic lineages belong to families. Gene families vary in terms of number of members, nucleotide similarity, gene integrity, expression, and function. Often, the members of gene families are arranged in clusters, which contribute to maintaining similarity among gene copies and also to generate duplicates through replication errors. Gene families offer us an opportunity to examine the forces involved in the evolution of the genomes and to study recombination events and genomic rearrangements. In this work, we focused on the evolution of two plant resistance gene families, Sw5 and Mi-1, and analyzed the completely sequenced nuclear genomes of potato and tomato. We first noticed that the potato genome carries larger resistance gene families than tomato, but all gene copies are pseudogenes. Second, phylogenetic analyses indicated that Sw5 and Mi-1 gene families had dissimilar evolutionary histories. In contrast to Sw5, Mi-1 homologues suffered repeated gene conversion events among the gene copies, particularly in the tomato genome.

Heterogeneity of Human Neutrophil CD177 Expression Results from CD177P1 Pseudogene Conversion.

Most humans harbor both CD177neg and CD177pos neutrophils but 1-10% of people are CD177null, placing them at risk for formation of anti-neutrophil antibodies that can cause transfusion-related acute lung injury and neonatal alloimmune neutropenia. By deep sequencing the CD177 locus, we catalogued CD177 single nucleotide variants and identified a novel stop codon in CD177null individuals arising from a single base substitution in exon 7. This is not a mutation in CD177 itself, rather the CD177null phenotype arises when exon 7 of CD177 is supplied entirely by the CD177 pseudogene (CD177P1), which appears to have resulted from allelic gene conversion. In CD177 expressing individuals the CD177 locus contains both CD177P1 and CD177 sequences. The proportion of CD177hi neutrophils in the blood is a heritable trait. Abundance of CD177hi neutrophils correlates with homozygosity for CD177 reference allele, while heterozygosity for ectopic CD177P1 gene conversion correlates with increased CD177neg neutrophils, in which both CD177P1 partially incorporated allele and paired intact CD177 allele are transcribed. Human neutrophil heterogeneity for CD177 expression arises by ectopic allelic conversion. Resolution of the genetic basis of CD177null phenotype identifies a method for screening for individuals at risk of CD177 isoimmunisation.

Identification and Evolution of Functional Alleles of the Previously Described Pollen Specific Myrosinase Pseudogene AtTGG6 in Arabidopsis thaliana.

Myrosinases are ß-thioglucoside glucohydrolases and serve as defense mechanisms against insect pests and pathogens by producing toxic compounds. AtTGG6 in Arabidopsis thaliana was previously reported to be a myrosinase pseudogene but specifically expressed in pollen. However, we found that AlTGG6, an ortholog to AtTGG6 in A. lyrata (an outcrossing relative of A. thaliana) was functional, suggesting that functional AtTGG6 alleles may still exist in A. thaliana. AtTGG6 alleles in 29 A. thaliana ecotypes were cloned and sequenced. Results indicate that ten alleles were functional and encoded Myr II type myrosinase of 512 amino acids, and myrosinase activity was confirmed by overexpressing AtTGG6 in Pichia pastoris. However, the 19 other ecotypes had disabled alleles with highly polymorphic frame-shift mutations and diversified sequences. Thirteen frame-shift mutation types were identified, which occurred independently many times in the evolutionary history within a few thousand years. The functional allele was expressed specifically in pollen similar to the disabled alleles but at a higher expression level, suggesting its role in defense of pollen against insect pests such as pollen beetles. However, the defense function may have become less critical after A. thaliana evolved to self-fertilization, and thus resulted in loss of function in most ecotypes.

Comparative Genomics of a Plant-Parasitic Nematode Endosymbiont Suggest a Role in Nutritional Symbiosis.

Bacterial mutualists can modulate the biochemical capacity of animals. Highly coevolved nutritional mutualists do this by synthesizing nutrients missing from the host's diet. Genomics tools have advanced the study of these partnerships. Here we examined the endosymbiont Xiphinematobacter (phylum Verrucomicrobia) from the dagger nematode Xiphinema americanum, a migratory ectoparasite of numerous crops that also vectors nepovirus. Previously, this endosymbiont was identified in the gut, ovaries, and eggs, but its role was unknown. We explored the potential role of this symbiont using fluorescence in situ hybridization, genome sequencing, and comparative functional genomics. We report the first genome of an intracellular Verrucomicrobium and the first exclusively intracellular non-Wolbachia nematode symbiont. Results revealed that Xiphinematobacter had a small 0.916-Mb genome with only 817 predicted proteins, resembling genomes of other mutualist endosymbionts. Compared with free-living relatives, conserved proteins were shorter on average, and there was large-scale loss of regulatory pathways. Despite massive gene loss, more genes were retained for biosynthesis of amino acids predicted to be essential to the host. Gene ontology enrichment tests showed enrichment for biosynthesis of arginine, histidine, and aromatic amino acids, as well as thiamine and coenzyme A, diverging from the profiles of relatives Akkermansia muciniphilia (in the human colon), Methylacidiphilum infernorum, and the mutualist Wolbachia from filarial nematodes. Together, these features and the location in the gut suggest that Xiphinematobacter functions as a nutritional mutualist, supplementing essential nutrients that are depleted in the nematode diet. This pattern points to evolutionary convergence with endosymbionts found in sap-feeding insects.

Polymorphism and evolution of ribosomal DNA in tea (Camellia sinensis, Theaceae).

Ribosomal DNA (rDNA) repeats often exhibit a high level of within-species homogeneity, but intra-individual polymorphism of rDNA has been found in Camellia species. In order to reveal the evolutionary pattern of rDNA repeats in tea [Camellia sinensis (L.) O. Kuntze], we identified the 45S rDNA loci, estimated their copy number, and cloned partial regions of them from different PCR products and from digested genomic DNA. The results show that there are 3 loci (6 sites) containing 45S rDNA in the tea genome; they are located at the ends of the short arms of 6 chromosomes and consist of a maximum of ∼6500 repeat units. On one hand, 164 sequences of the 26S rDNA cloned from PCR products contained 147 haplotypes, consisting of 62% pseudogenes, 24% putative functional genes, and 14% PCR-mediated recombinants, suggesting that the 45S rDNA of tea maintains an extremely high level of polymorphism and divergence; on the other hand, rDNA fragments cloned directly from genomic DNA exhibited a very high level of homogeneity: only one of 33 rDNA fragments was from a pseudogene. These results demonstrate that although over 60% of the 26S rDNA sequences identified in PCR products belong to pseudogenes, most 45S rDNA repeats are functional genes and have undergone concerted evolution. In this study, strong PCR bias and PCR-mediated recombination greatly increased the apparent proportion of pseudogenes in PCR products. Phylogenetic analysis and genetic divergence values for 26S rDNA sequences obtained in this study show that many pseudogenes have originated independently from functional genes at different times, and despite thus escaping from concerted evolution, they have failed to be eliminated from the tea genome over a long period, some of them having even produced addition copies by rapid expansion. Importantly, our study suggests that in order to determine the true pattern of evolution of rDNA it is necessary to combine data from more than one method rather than relying only on sequences from PCR products.

The complete chloroplast genome of colchicine plants (Colchicum autumnale L. and Gloriosa superba L.) and its application for identifying the genus.

MAIN CONCLUSION: The complete chloroplast genome of two colchicine medicinal plants is reported for the first time. Deletion of ycf 15 gene occurred only in Colchicum but not in Gloriosa and suggests this as a potential marker for delineating the two species. Colchicum autumnale L. and Gloriosa superba L. are well-known sources of colchicine, a type of alkaloid and an ancient anti-inflammatory drug used to prevent gout. Accordingly, this alkaloid has been used as a chemical marker for identifying the expanded Colchicaceae family. In the present study, we report the complete chloroplast genome (cpDNA) sequence of two colchicine medicinal plants (G. superba and C. autumnale) that belong to the tribe Colchiceae of the Colchicaceae family. In C. autumnale, the circular double-stranded cpDNA sequence of 156,462 bp consists of two inverted repeat (IR) regions of 27,741 bp each, a large single-copy region (LSC) of 84,246 bp, and a small single-copy region (SSC) of 16,734 bp. The cpDNA sequence of G. superba is longer than that of C. autumnale (157,924 bp), which consists of two IRs (28,063 bp), an SSC (16,786 bp), and an LSC (85,012 bp). Significant structural differences between them were observed in the ycf15 gene. ycf15 gene was absent from C. autumnale cpDNA and affected the length of the chloroplast genome between the species. Furthermore, this gene loss event was specific to the expanded genus of Colchicum sensu Vinnersten and Manning. Therefore, this gene may be an effective and powerful molecular marker for identifying the Colchicum genus within the family.

Functional redundancy and/or ongoing pseudogenization among F-box protein genes expressed in Arabidopsis male gametophyte.

F-box protein genes family is one of the largest gene families in plants, with almost 700 predicted genes in the model plant Arabidopsis. F-box proteins are key components of the ubiquitin proteasome system that allows targeted protein degradation. Transcriptome analyses indicate that half of these F-box protein genes are found expressed in microspore and/or pollen, i.e., during male gametogenesis. To assess the role of F-box protein genes during this crucial developmental step, we selected 34 F-box protein genes recorded as highly and specifically expressed in pollen and isolated corresponding insertion mutants. We checked the expression level of each selected gene by RT-PCR and confirmed pollen expression for 25 genes, but specific expression for only 10 of the 34 F-box protein genes. In addition, we tested the expression level of selected F-box protein genes in 24 mutant lines and showed that 11 of them were null mutants. Transmission analysis of the mutations to the progeny showed that none of the single mutations was gametophytic lethal. These unaffected transmission efficiencies suggested leaky mutations or functional redundancy among F-box protein genes. Cytological observation of the gametophytes in the mutants confirmed these results. Combinations of mutations in F-box protein genes from the same subfamily did not lead to transmission defect either, further highlighting functional redundancy and/or a high proportion of pseudogenes among these F-box protein genes.

Genome-wide identification and mapping of NBS-encoding resistance genes in Solanum tuberosum group phureja.

The majority of disease resistance (R) genes identified to date in plants encode a nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domain containing protein. Additional domains such as coiled-coil (CC) and TOLL/interleukin-1 receptor (TIR) domains can also be present. In the recently sequenced Solanum tuberosum group phureja genome we used HMM models and manual curation to annotate 435 NBS-encoding R gene homologs and 142 NBS-derived genes that lack the NBS domain. Highly similar homologs for most previously documented Solanaceae R genes were identified. A surprising ∼41% (179) of the 435 NBS-encoding genes are pseudogenes primarily caused by premature stop codons or frameshift mutations. Alignment of 81.80% of the 577 homologs to S. tuberosum group phureja pseudomolecules revealed non-random distribution of the R-genes; 362 of 470 genes were found in high density clusters on 11 chromosomes.

Composition and evolution of the vertebrate and mammalian selenoproteomes.

BACKGROUND: Selenium is an essential trace element in mammals due to its presence in proteins in the form of selenocysteine (Sec). Human genome codes for 25 Sec-containing protein genes, and mouse and rat genomes for 24. METHODOLOGY/PRINCIPAL FINDINGS: We characterized the selenoproteomes of 44 sequenced vertebrates by applying gene prediction and phylogenetic reconstruction methods, supplemented with the analyses of gene structures, alternative splicing isoforms, untranslated regions, SECIS elements, and pseudogenes. In total, we detected 45 selenoprotein subfamilies. 28 of them were found in mammals, and 41 in bony fishes. We define the ancestral vertebrate (28 proteins) and mammalian (25 proteins) selenoproteomes, and describe how they evolved along lineages through gene duplication (20 events), gene loss (10 events) and replacement of Sec with cysteine (12 events). We show that an intronless selenophosphate synthetase 2 gene evolved in early mammals and replaced functionally the original multiexon gene in placental mammals, whereas both genes remain in marsupials. Mammalian thioredoxin reductase 1 and thioredoxin-glutathione reductase evolved from an ancestral glutaredoxin-domain containing enzyme, still present in fish. Selenoprotein V and GPx6 evolved specifically in placental mammals from duplications of SelW and GPx3, respectively, and GPx6 lost Sec several times independently. Bony fishes were characterized by duplications of several selenoprotein families (GPx1, GPx3, GPx4, Dio3, MsrB1, SelJ, SelO, SelT, SelU1, and SelW2). Finally, we report identification of new isoforms for several selenoproteins and describe unusually conserved selenoprotein pseudogenes. CONCLUSIONS/SIGNIFICANCE: This analysis represents the first comprehensive survey of the vertebrate and mammal selenoproteomes, and depicts their evolution along lineages. It also provides a wealth of information on these selenoproteins and their forms.

Functional test of Brassica self-incompatibility modifiers in Arabidopsis thaliana.

The self-incompatibility (SI) system of the Brassicaceae is based on allele-specific interactions among haplotypes of the S locus. In all tested self-incompatible Brassicaceae, the S haplotype encompasses two linked genes, one encoding the S-locus receptor kinase (SRK), a transmembrane kinase displayed at the surface of stigma epidermal cells, and the other encoding its ligand, the S-locus cysteine-rich (SCR) protein, which is localized in the pollen coat. Transfer of the two genes to self-fertile Arabidopsis thaliana allowed the establishment of robust SI in several accessions, indicating that the signaling cascade triggered by this receptor-ligand interaction and the resulting inhibition of "self" pollen by the stigma have been maintained in extant A. thaliana. Based on studies in Brassica species, the membrane-tethered kinase MLPK, the ARM repeat-containing U-box protein ARC1, and the exocyst subunit Exo70A1 have been proposed to function as components of an SI signaling cascade. Here we tested the role of these molecules in the SI response of A. thaliana SRK-SCR plants. We show that the A. thaliana ARC1 ortholog is a highly decayed pseudogene. We also show that, unlike reports in Brassica, inactivation of the MLPK ortholog AtAPK1b and overexpression of Exo70A1 neither abolish nor weaken SI in A. thaliana SRK-SCR plants. These results do not support a role for these molecules in the SI response of A. thaliana.

Identification and characterization of plastid trnF(GAA) pseudogenes in four species of Solanum (Solanaceae).

Non-functional trnF pseudogenes that rarely occur in embryophytes have been found in Solanaceae. We have sequenced the trnL-F intergenic spacer of four species of Solanum, and found duplicated regions of the original trnF gene. These repeats were 94-260 bp long causing large length variation in the trnL-F intergenic spacer resulting from differences in pseudogene copy number (2-4). The duplicated trnF regions are comprised of several highly structured motifs, which were partial residues, or entire parts of the Anticodon, T- and D-domains of the original gene, but all lacked the acceptor stems at the 5'- or 3'-end. Pseudogenes included several transitions and transversions in their sequences compared to the original trnF gene. Among pseudogene copies, T-domains were more frequent and fragmented than D-domain elements. Our results demonstrate that although chloroplast evolution is uniform such structural duplications in the sequences used for phylogenetic reconstructions should be treated with great caution.

Gene expansion and retention leads to a diverse tyrosine kinase superfamily in amphioxus.

Tyrosine kinase (TK) proteins play a central role in cellular behavior and development of animals. The expansion of this superfamily is regarded as a key event in the evolution of the complex signaling pathways and gene networks of metazoans and is a prominent example of how shuffling of protein modules may generate molecular novelties. Using the intron/exon structure within the TK domain (TK intron code) as a complementary tool for the assignment of orthology and paralogy, we identified and studied the 118 TK proteins of the amphioxus Branchiostoma floridae genome to elucidate TK gene family evolution in metazoans and chordates in particular. Unlike all characterized metazoans to date, amphioxus has members of all known widespread TK families, with not a single loss. Putting amphioxus TKs in an evolutionary context, including new data from the cnidarian Nematostella vectensis, the echinoderm Strongylocentrotus purpuratus, and the ascidian Ciona intestinalis, we suggest new evolutionary histories for different TK families and draw a new global picture of gene loss/gain in the different phyla. Surprisingly, our survey also detected an unprecedented expansion of a group of closely related TK families, including TIE, FGFR, PDGFR, and RET, due most probably to massive gene duplication and exon shuffling. Based on their highly similar intron/exon structure at the TK domain, we suggest that this group of TK families constitute a superfamily of TK proteins, which we termed EXpanding TK, after their seemingly unique propensity to gene duplication and exon shuffling, not only in amphioxus but also across all metazoan groups. Due to this extreme tendency to both retention and expansion of TK genes, amphioxus harbors the richest and most diverse TK repertoire among all metazoans studied so far, retaining most of the gene complement of its ancestors, but having evolved its own repertoire of genetic novelties.

Patterns of partial RNA editing in mitochondrial genes of Beta vulgaris.

RNA editing is a process that modifies the information in transcripts of almost all angiosperm mitochondrial protein-coding genes. In order to determine the frequency and distribution of mitochondrial RNA editing in Beta vulgaris, cDNAs were sequenced and compared to the published genome sequence. 357 C to U conversions were identified across the 31 known protein genes and pseudogenes in Beta, the fewest so far for a plant mitochondrial genome. Editing patterns in the putative gene orf518 indicate that it is most likely a functional ccmC homolog, indicating that patterns of editing can be a useful determinant of gene functionality. orf518 also contains a triplicated repeat region whose members are nearly identical yet differentially edited, most likely due to differences in the sequence context of the editing sites. In addition, we show that partial editing in Beta is common at silent editing sites but rare at nonsilent editing sites, extending previous observations to a complete plant mitochondrial genome. Finally, the degree of partial editing observed for certain genes was dependent on the choice of primers used, demonstrating that care must be taken when designing primers for use in editing studies.

On the origin of leprosy.

Leprosy, a chronic human disease with potentially debilitating neurological consequences, results from infection with Mycobacterium leprae. This unculturable pathogen has undergone extensive reductive evolution, with half of its genome now occupied by pseudogenes. Using comparative genomics, we demonstrated that all extant cases of leprosy are attributable to a single clone whose dissemination worldwide can be retraced from analysis of very rare single-nucleotide polymorphisms. The disease seems to have originated in Eastern Africa or the Near East and spread with successive human migrations. Europeans or North Africans introduced leprosy into West Africa and the Americas within the past 500 years.

Mice have a transcribed L-threonine aldolase/GLY1 gene, but the human GLY1 gene is a non-processed pseudogene.

BACKGROUND: There are three pathways of L-threonine catabolism. The enzyme L-threonine aldolase (TA) has been shown to catalyse the conversion of L-threonine to yield glycine and acetaldehyde in bacteria, fungi and plants. Low levels of TA enzymatic activity have been found in vertebrates. It has been suggested that any detectable activity is due to serine hydroxymethyltransferase and that mammals lack a genuine threonine aldolase. RESULTS: The 7-exon murine L-threonine aldolase gene (GLY1) is located on chromosome 11, spanning 5.6 kb. The cDNA encodes a 400-residue protein. The protein has 81% similarity with the bacterium Thermotoga maritima TA. Almost all known functional residues are conserved between the two proteins including Lys242 that forms a Schiff-base with the cofactor, pyridoxal-5'-phosphate. The human TA gene is located at 17q25. It contains two single nucleotide deletions, in exons 4 and 7, which cause frame-shifts and a premature in-frame stop codon towards the carboxy-terminal. Expression of human TA mRNA was undetectable by RT-PCR. In mice, TA mRNA was found at low levels in a range of adult tissues, being highest in prostate, heart and liver. In contrast, serine/threonine dehydratase, another enzyme that catabolises L-threonine, is expressed very highly only in the liver. Serine dehydratase-like 1, also was most abundant in the liver. In whole mouse embryos TA mRNA expression was low prior to E-15 increasing more than four-fold by E-17. CONCLUSION: Mice, the western-clawed frog and the zebrafish have transcribed threonine aldolase/GLY1 genes, but the human homolog is a non-transcribed pseudogene. Serine dehydratase-like 1 is a putative L-threonine catabolising enzyme.