Search

Estrogen drives key transcriptional changes in breast cancer and stimulates breast cancer cells' growth with multiple mechanisms to coordinate transcription and translation. In addition to protein-coding transcripts, estrogen can regulate long non-coding RNA (lncRNA) transcripts, plus diverse non-coding RNAs including antisense, enhancer, and intergenic. LncRNA genes comprise the majority of human genes. The accidental, or regulated, translation of their short open reading frames by ribosomes remains a controversial topic. Here we report for the first time an integrated analysis of RNA abundance and ribosome occupancy level, using Ribo-seq combined with RNA-Seq, in the estrogen-responsive, estrogen receptor α positive, human breast cancer cell model MCF7, before and after hormone treatment. Translational profiling can determine, in an unbiased manner, which fraction of the genome is actually translated into proteins, as well as resolving whether transcription and translation respond concurrently, or differentially, to estrogen treatment. Our data showed specific transcripts more robustly detected in RNA-Seq than in the ribosome-profiling data, and vice versa, suggesting distinct gene-specific estrogen responses at the transcriptional and the translational level, respectively. Here, we showed that estrogen stimulation affects the expression levels of numerous lncRNAs, but not their association with ribosomes, and that most lncRNAs are not ribosome-bound. For the first time, we also demonstrated the transcriptional and translational response of expressed pseudogenes to estrogen, pointing to new perspectives for drug-target development in breast cancer in the future., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

The loss of teeth and evolution of baleen racks in Mysticeti was a profound transformation that permitted baleen whales to radiate and diversify into a previously underutilized ecological niche of bulk filter-feeding on zooplankton and other small prey. Ancestral state reconstructions suggest that postnatal teeth were lost in the common ancestor of crown Mysticeti. Genomic studies provide some support for this hypothesis and suggest that the genetic toolkit for enamel production was inactivated in the common ancestor of living baleen whales. However, molecular studies to date have not provided direct evidence for the complete loss of teeth, including their dentin component, on the stem mysticete branch. Given these results, several questions remain unanswered: (1) Were teeth lost in a single step or did enamel loss precede dentin loss? (2) Was enamel lost early or late on the stem mysticete branch? (3) If enamel and dentin/tooth loss were decoupled in the ancestry of baleen whales, did dentin loss occur on the stem mysticete branch or independently in different crown mysticete lineages? To address these outstanding questions, we compiled and analyzed complete protein-coding sequences for nine tooth-related genes from cetaceans with available genome data. Seven of these genes are associated with enamel formation (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) whereas two other genes are either dentin-specific (DSPP) or tooth-specific (ODAPH) but not enamel-specific. Molecular evolutionary analyses indicate that all seven enamel-specific genes have inactivating mutations that are scattered across branches of the mysticete tree. Three of the enamel genes (ACP4, KLK4, MMP20) have inactivating mutations that are shared by all mysticetes. The two genes that are dentin-specific (DSPP) or tooth-specific (ODAPH) do not have any inactivating mutations that are shared by all mysticetes, but there are shared mutations in Balaenidae as well as in Plicogulae (Neobalaenidae + Balaenopteroidea). These shared mutations suggest that teeth were lost at most two times. Shared inactivating mutations and dN/dS analyses, in combination with cetacean divergence times, were used to estimate inactivation times of genes and by proxy enamel and tooth phenotypes at ancestral nodes. The results of these analyses are most compatible with a two-step model for the loss of teeth in the ancestry of living baleen whales: enamel was lost very early on the stem Mysticeti branch followed by the independent loss of dentin (and teeth) in the common ancestors of Balaenidae and Plicogulae, respectively. These results imply that some stem mysticetes, and even early crown mysticetes, may have had vestigial teeth comprised of dentin with no enamel. Our results also demonstrate that all odontocete species (in our study) with absent or degenerative enamel have inactivating mutations in one or more of their enamel genes., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Pancreatic cancer, the most lethal malignant tumor, is notorious for its poor prognosis and metastatic potential. Non-coding RNAs (ncRNAs) are reported to play key roles in cancer metastasis. In this study, miRNA and gene expression profiles between metastatic pancreatic cancer cell M8 and its parental cell BxPC.3 were determined. Using differential expression analysis, survival analysis, target gene prediction, pathway enrichment analysis, intersection analysis and correlation analysis, hsa-miR-30d-5p/GJA1 axis was identified as the most potential pathway involved in metastasis of pancreatic cancer. Subsequently, two upstream lncRNAs (HELLPAR and OIP-AS1) and four upstream pseudogenes (AC093616.1, AC009951.1, TMEM183B and PABPC1P4) of hsa-miR-30d-5p/GJA1 axis were predicted and were then identified via assessment of RNA-RNA expression relationship. Furthermore, CTNNA1, CTNNB1 and CTNND1 were regarded as three crucial molecules to be participated in hsa-miR-30d-5p/GJA1-mediated metastatic potential in pancreatic cancer. In conclusion, we established a novel lncRNA/pseudogene-hsa-miR-30d-5p-GJA1 regulatory network linked to metastasis of pancreatic cancer., (Copyright © 2021 Elsevier Inc. All rights reserved.)

The ancestors of Cetacea underwent profound morpho-physiological alterations. By displaying an exclusive aquatic existence, cetaceans evolved unique patterns of locomotor activity, vigilant behaviour, thermoregulation and circadian rhythmicity. Deciphering the molecular landscape governing many of these adaptations is key to understand the evolution of phenotypes. Here, we investigate Cortistatin (CORT), a neuropeptide displaying an important role in mammalian biorhythm regulation. This neuropeptide is a known neuroendocrine factor, stimulating slow-wave sleep, but also involved in the regulation of energy metabolism and hypomotility inducement. We assessed the functional status of CORT in 359 mammalian genomes (25 orders), including 30 species of Cetacea. Our findings indicate that cetaceans and other mammals with atypical biorhythms, thermal constraints and/or energy metabolism, have accumulated deleterious mutations in CORT. In light of the pleiotropic action of this neuropeptide, we suggest that this inactivation contributed to a plethora of phenotypic adjustments to accommodate adaptive solutions to specific ecological niches., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The protein-coding genes and pseudogenes of Cuscuta australis had the diverse contribution to the formation and evolution of parasitism. The codon usage pattern analysis of these two type genes could be used to understand the gene transcription and translation. In this study, we systematically analyzed the codon usage patterns of protein-coding sequences and pseudogenes sequences in C. australis. The results showed that the high frequency codons of protein coding sequences and pseudogenes had the same A/U bias in the third position. However, these two sequences had converse bias at the third base in optimal codons: the protein coding sequences preferred G/C-ending codons while pseudogene sequences preferred A/U-ending codons. Neutrality plot and effective number of codons plot revealed that natural selection played a more important role than mutation pressure in two sequences codon usage bias. Furthermore, the gene expression level had a significant positive correlation with codon usage bias in C. australis. Highly-expressed protein coding genes exhibited a higher codon bias than lowly-expressed genes. Meanwhile, the high-expression genes tended to use G/C-ending synonymous codons. This result further verified the optimal codons usage bias and its correlation with the gene expression in C. australis., Competing Interests: Declaration of Competing Interest There were no competing interests., (Copyright © 2020. Published by Elsevier Inc.)

Described as "junk" DNA, pseudogenes are dead structures of previously active genes present in genomes. Pseudogenes are categorized into two main classes: processed pseudogenes, formed through retrotransposition, and non-processed pseudogenes, typically originated from gene decay following duplication events. The term "processed pseudogene" has changed to "retrocopy" since they are likely to evolve new functional roles and became a retrogene. Here, we surveyed 38,080 retrocopies from chimpanzee, dog, human, mouse, and rat genomes to assess their potential adaptive value. The retrocopies inserted in the same chromosome of the parental gene have higher chances of remain potentially "active" (absence of premature stop codons and frameshifts) (~26.1%), while those placed into a different chromosome have a twofold decrease chance of continuing potentially "active" (~7.52%). The genomic context of their placement seems associated with their expression. Retrocopies placed in intragenic regions and the same sense of the "host" gene have higher chances of being expressed relative to other genomic contexts. The proximity of retrocopies to their parental gene is associated with a lower decay rate, and their location likely influence their expression. Thus, despite their unclear role, retrocopies are probably involved in adaptive processes. Our results evidence natural selection acting in retrocopies., (Copyright © 2020. Published by Elsevier Inc.)

CACX is one of the most common cancer affecting women world-wide. Here, expression microarray analysis revealed 8 over-expressed transcribed pseudogenes (GBP1P1, HLA-DRB6, HLA-H, SLC6A10P, NAPSB, KRT16P2, PTTG3P and RNF126P1), down-regulated 7 lincRNAs (H19, MIR100HG, MEG3, DIO3OS, HOXA11-AS, CD27-AS1 and EPB41L4A-AS) and 6 snoRNAs (SNORD97, SNORD3A, SNORD3C, SNORD3D, SNORA12 and SCARNA9) as DEncGs (log2 fold-change ≥ ±1.0) in CACX. Consequently, down-regulation of lincRNA MEG3 and over-expression of pseudogenes, GBP1P1 and PTTG3P in the microarray analysis were found concordant with the real-time quantitative PCR results upon validation. Then, Ingenuity® Pathway analysis (IPA®) analysis with deregulated DEncGs identified functionally important gene, H19. Further, validation (n = 52) of expression confirmed frequent downregulation of H19 with significant association with its deletion (LOH) and promoter methylation (n = 128) in CACX. Moreover, clinicopathological analysis found Indian CACX patients (n = 26) with alterations of H19 by deletion or, promoter methylation with concomitant low expression have poor prognosis., (Copyright © 2019. Published by Elsevier Inc.)

Actinobacillus spp. are Gram-negative bacteria associated with mucosal membranes. While some are commensals, others can cause important human and animal diseases. A. pleuropneumoniae causes severe fibrinous hemorrhagic pneumonia in swine but not systemic disease whereas other species invade resulting in septicemia and death. To understand the invasive phenotype of Actinobacillus spp., complete genomes of eight isolates were obtained and pseudogenomes of five isolates were assembled and annotated. Phylogenetically, A. suis isolates clustered by surface antigen type and were more closely related to the invasive A. ureae, A. equuli equuli, and A. capsulatus than to the other swine pathogen, A. pleuropneumoniae. Using the LS-BSR pipeline, 251 putative virulence genes associated with serum resistance and invasion were detected. To our knowledge, this is the first genome-wide study of the genus Actinobacillus and should contribute to a better understanding of host tropism and mechanisms of invasion of pathogenic Actinobacillus and related genera., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Mycobacterium leprae is an unculturable obligate pathogen and causative agent for debilitating human disease leprosy. Due to reductive genome evolution M leprae genome harbours large number of pseudogenes and small number of genes (∼1600 genes and ∼1300 pseudogenes). How M leprae remained a successful human parasite with small set of genes remains poorly understood and provided us the impetus to investigate the intergenic regions of M leprae genome for the presence of possible open reading frames (ORFs). In this work, we have manually scanned all the intergenic regions of M leprae genome and identified 106 potential ORFs. Among these, 12 are large ORFs: encoding hypothetical proteins (HP) of more than 100 amino acids. We have also found 67 ORFs encoding 50-100 amino acids proteins and another 27 ORFs for 30-50 amino acids peptides. We have validated the presence of transcripts for large HPs by quantitative reverse transcriptase PCR (qRT-PCR). Our results suggest that some of the M leprae large HPs are indeed expressed at low level in leprosy patients. The present results will shed light on the intergenic ORFs of M leprae and further our understanding of the pathogenesis of leprosy., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

MC5R is one of five melanocortin receptor genes found in placental mammals. MC5R plays an important role in energy homeostasis and is also expressed in the terminal differentiation of sebaceous glands. Among placental mammals there are multiple lineages that either lack or have degenerative sebaceous glands including Cetacea (whales, dolphins, and porpoises), Hippopotamidae (hippopotamuses), Sirenia (manatees and dugongs), Proboscidea (elephants), Rhinocerotidae (rhinos), and Heterocephalus glaber (naked mole rat). Given the loss or diminution of sebaceous glands in these taxa, we procured MC5R sequences from publicly available genomes and transcriptomes, supplemented by a newly generated sequence for Choeropsis liberiensis (pygmy hippopotamus), to determine if this gene remains intact or is inactivated in association with loss/reduction of sebaceous glands. Our data set includes complete MC5R sequences for 114 placental mammal species including two individuals of Mammuthus primigenius (woolly mammoth) from Oimyakon and Wrangel Island. Complete loss or inactivation of the MC5R gene occurs in multiple placental lineages that have lost sebaceous glands (Cetacea, West Indian manatee, African elephant, white rhinoceros) or are characterized by unusual skin (pangolins, aardvarks). Both M. primigenius individuals share inactivating mutations with the African elephant even though sebaceous glands have been reported in the former. MC5R remains intact in hippopotamuses and the naked mole rat, although slightly elevated dN/dS ratios in these lineages allow for the possibility that the accumulation of inactivating mutations in MC5R may lag behind the relaxation of purifying selection. For Cetacea and Hippopotamidae, the absence of shared inactivating mutations in two different skin genes (MC5R, PSORS1C2) is consistent with the hypothesis that semi-aquatic lifestyles were acquired independently in these clades following divergence from a common ancestor., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Nuclear copies of mitochondrial genes (numts) are a well-known feature of eukaryotic genomes and a concern in systematics, as they can mislead phylogenetic inferences when inadvertently used. Studies on avian numts initially based on the chicken genome suggest that numts may be uncommon and relatively short among birds. Here we ask how common numts are in falcons, based on recently sequenced genomes of the Saker falcon (Falco cherrug) and Peregrine falcon (F. peregrinus). We identified numts by BLASTN searches and then extracted CYTB, ND2 and COI sequences from them, which were then used for phylogeny inference along with several sequences from other species in Falconiformes. Our results indicate that avian numts may be much more frequent and longer than previously thought. Phylogenetic inferences revealed multiple independent nuclear insertions throughout the history of the Falconiformes, including cases of sequences available in public databases and wrongly identified as authentic mtDNA. New sequencing technologies and ongoing efforts for whole genome sequencing will provide exciting opportunities for avian numt research in the near future., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Mitogenomic trees for Bivalvia have proved problematic in the past, but several highly divergent lineages were missing from these analyses and increased representation of these groups may yet improve resolution. Here, we add seven new sequences from the Anomalodesmata and one unidentified semelid species (Bryopa lata, Euciroa cf. queenslandica, Laternula elliptica, Laternula truncata, Lyonsia norwegica, Myadora brevis, Tropidomya abbreviata, "Abra" sp.). We show that relationships in a mitogenomic tree for the Class are improved by the addition of seven anomalodesmatans from this highly divergent clade, but are still not completely consistent with relationships recovered in studies of nuclear genes. We suggest that some anomalous relationships (for instance the non-monophyly of Bivalvia) may be partially explained by compositional heterogeneity in the mitogenome and suggest that the addition of more taxa may help resolve both this effect and possible instances of long branch attraction. We also identify several curious features about anomalodesmatan mitogenomes. For example, many protein-coding gene boundaries are poorly defined in marine bivalves, but particularly so in anomalodesmatans, primarily due to non-conserved boundary sequences. The use of transcriptomic and genomic data together enabled better definition of gene boundaries, the identification of possible pseudogenes and suggests that most genes are translated monocistronically, which contrasts with many other studies. We also identified a possible case of gene duplication of ND5 in Myadora brevis (Myochamidae). Mitogenome size in the Anomalodesmata ranges from very small compact molecules, with the smallest for Laternula elliptica (Laternulidae) only 14,622bp, to Bryopa lata (Clavagellidae) which is at least 31,969bp long and may be >40,000bp. Finally, sampled species show a high degree of sequence divergence and variable gene order, although intraspecific variation in Laternula elliptica is very low., (Crown Copyright © 2017. Published by Elsevier Inc. All rights reserved.)

Here we report a comprehensive analysis of the vomeronasal receptor repertoire in pigs. We identified a total of 25 V1R sequences consisting of 10 functional genes, 3 pseudogenes, and 12 partial genes, while functional V2R and FPR genes were not present in the pig genome. Pig V1Rs were classified into three subfamilies, D, F, and J. Using direct high resolution sequencing-based typing of all functional V1Rs from 10 individuals of 5 different breeds, a total of 24 SNPs were identified, indicating that the allelic diversity of V1Rs is much lower than that of the olfactory receptors. A high expression level of V1Rs was detected in the vomeronasal organ (VNO) and testes, while a low expression level of V1Rs was observed in all other tissues examined. Our results showed that pigs could serve as an interesting large animal model system to study pheromone-related neurobiology because of their genetic simplicity., (Copyright © 2015. Published by Elsevier Inc.)

The deregulation of cellular long non-coding RNA (lncRNA) expression during a human adenovirus infection was studied by deep sequencing. Expression of lncRNAs increased substantially following the progression of the infection. Among 645 significantly expressed lncRNAs, the expression of 398 was changed more than 2-fold. More than 80% of them were up-regulated and 80% of them were detected during the late phase. Based on the genomic locations of the deregulated lncRNAs in relation to known mRNAs and miRNAs, they were predicted to be involved in growth, structure, apoptosis and wound healing in the early phase, cell proliferation in the intermediate phase and protein synthesis, modification and transport in the late phase. The most significant functions of cellular RNA-binding proteins, previously shown to interact with the deregulated lncRNAs identified here, are involved in RNA splicing, nuclear export and translation events. We hypothesize that adenoviruses exploit the lncRNA network to optimize their reproduction., (Copyright © 2016 Elsevier Inc. All rights reserved.)

The NADPH-dependent cytochrome P450 oxidoreductase (POR) is the obligate electron donor to eukaryotic microsomal cytochromes P450 enzymes. The number of PORs within plant species is limited to one to four isoforms, with the most common being two PORs per plant. These enzymes provide electrons to a huge number of different cytochromes P450s (from 50 to several hundred within one plant). Within the eudicotyledons, PORs can be divided into two major clades, POR 1 and POR 2. Based on our own sequencing analysis and publicly available data, we have identified 45 PORs from the angiosperm order Apiales. These were subjected to a phylogenetic analysis along with 237 other publicly available (NCBI and oneKP) POR sequences found within the clade Asterids. Here, we show that the order Apiales only harbor members of the POR 2 clade, which are further divided into two distinct subclades. This is in contrast to most other eudicotyledon orders that have both POR 1 and POR 2. This suggests that through gene duplications and one gene deletion, Apiales only contain members of the POR 2 clade. Three POR 2 isoforms from Thapsia garganica L., Apiaceae, were all full-length in an Illumina root transcriptome dataset (available from the SRA at NCBI). All three genes were shown to be functional upon reconstitution into nanodiscs, confirming that none of the isoforms are pseudogenes., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Species of the fungus Fusarium collectively cause disease on almost all crop plants and produce numerous natural products (NPs), including some of the mycotoxins of greatest concern to agriculture. Many Fusarium NPs are derived from polyketide synthases (PKSs), large multi-domain enzymes that catalyze sequential condensation of simple carboxylic acids to form polyketides. To gain insight into the biosynthesis of polyketide-derived NPs in Fusarium, we retrieved 488 PKS gene sequences from genome sequences of 31 species of the fungus. In addition to these apparently functional PKS genes, the genomes collectively included 81 pseudogenized PKS genes. Phylogenetic analysis resolved the PKS genes into 67 clades, and based on multiple lines of evidence, we propose that homologs in each clade are responsible for synthesis of a polyketide that is distinct from those synthesized by PKSs in other clades. The presence and absence of PKS genes among the species examined indicated marked differences in distribution of PKS homologs. Comparisons of Fusarium PKS genes and genes flanking them to those from other Ascomycetes provided evidence that Fusarium has the genetic potential to synthesize multiple NPs that are the same or similar to those reported in other fungi, but that have not yet been reported in Fusarium. The results also highlight ways in which such analyses can help guide identification of novel Fusarium NPs and differences in NP biosynthetic capabilities that exist among fungi., (Published by Elsevier Inc.)

Previous studies have reported inactivated copies of six enamel-related genes (AMBN, AMEL, AMTN, ENAM, KLK4, MMP20) and one dentin-related gene (DSPP) in one or more toothless vertebrates and/or vertebrates with enamelless teeth, thereby providing evidence that these genes are enamel or tooth-specific with respect to their critical functions that are maintained by natural selection. Here, we employ available genome sequences for edentulous and enamelless mammals to evaluate the enamel specificity of four genes (WDR72, SLC24A4, FAM83H, C4orf26) that have been implicated in amelogenesis imperfecta, a condition in which proper enamel formation is abrogated during tooth development. Coding sequences for WDR72, SCL24A4, and FAM83H are intact in four edentulous taxa (Chinese pangolin, three baleen whales) and three taxa (aardvark, nine-banded armadillo, Hoffmann's two-toed sloth) with enamelless teeth, suggesting that these genes have critical functions beyond their involvement in tooth development. By contrast, genomic data for C4orf26 reveal inactivating mutations in pangolin and bowhead whale as well as evidence for deletion of this gene in two minke whale species. Hybridization capture of exonic regions and PCR screens provide evidence for inactivation of C4orf26 in eight additional baleen whale species. However, C4orf26 is intact in all three species with enamelless teeth that were surveyed, as well as in 95 additional mammalian species with enamel-capped teeth. Estimates of selection intensity suggest that dN/dS ratios on branches leading to taxa with enamelless teeth are similar to the dN/dS ratio on branches leading to taxa with enamel-capped teeth. Based on these results, we conclude that C4orf26 is tooth-specific, but not enamel-specific, with respect to its essential functions that are maintained by natural selection. A caveat is that an alternative splice site variant, which translates exon 3 in a different reading frame, is putatively functional in Catarrhini and may have evolved an additional role in this primate clade., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Low-copy nuclear genes (LCNGs) have complex genetic architectures and evolutionary dynamics. However, unlike multicopy nuclear genes, LCNGs are rarely subject to gene conversion or concerted evolution, and they have higher mutation rates than organellar or nuclear ribosomal DNA markers, so they have great potential for improving the robustness of phylogenetic reconstructions at all taxonomic levels. In this study, our first objective is to evaluate the evolutionary dynamics of the LCNG β-amylase by testing for potential pseudogenization, paralogy, homeology, recombination, and phylogenetic incongruence within a broad representation of the main Pooideae lineages. Our second objective is to determine whether β-amylase shows sufficient phylogenetic signal to reconstruct the evolutionary history of the Pooid grasses. A multigenic (ITS, matK, ndhF, trnTL, and trnLF) tree of the study group provided a framework for assessing the β-amylase phylogeny. Eight accessions showed complete absence of selection, suggesting putative pseudogenic copies or other relaxed selection pressures; resolution of Vulpia alopecuros 2x clones indicated its potential (semi) paralogy; and homeologous copies of allopolyploid species Festuca simensis, F. fenas, and F. arundinacea tracked their Mediterranean origin. Two recombination events were found within early-diverged Pooideae lineages, and five within the PACCMAD clade. The unexpected phylogenetic relationships of 37 grass species (26% of the sampled species) highlight the frequent occurrence of non-treelike evolutionary events, so this LCNG should be used with caution as a phylogenetic marker. However, once the pitfalls are identified and removed, the phylogenetic reconstruction of the grasses based on the β-amylase exon+intron positions is optimal at all taxonomic levels., (Copyright © 2015 Elsevier Inc. All rights reserved.)

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., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Wood decay mechanisms in Agaricomycotina have been traditionally separated in two categories termed white and brown rot. Recently the accuracy of such a dichotomy has been questioned. Here, we present the genome sequences of the white-rot fungus Cylindrobasidium torrendii and the brown-rot fungus Fistulina hepatica both members of Agaricales, combining comparative genomics and wood decay experiments. C. torrendii is closely related to the white-rot root pathogen Armillaria mellea, while F. hepatica is related to Schizophyllum commune, which has been reported to cause white rot. Our results suggest that C. torrendii and S. commune are intermediate between white-rot and brown-rot fungi, but at the same time they show characteristics of decay that resembles soft rot. Both species cause weak wood decay and degrade all wood components but leave the middle lamella intact. Their gene content related to lignin degradation is reduced, similar to brown-rot fungi, but both have maintained a rich array of genes related to carbohydrate degradation, similar to white-rot fungi. These characteristics appear to have evolved from white-rot ancestors with stronger ligninolytic ability. F. hepatica shows characteristics of brown rot both in terms of wood decay genes found in its genome and the decay that it causes. However, genes related to cellulose degradation are still present, which is a plesiomorphic characteristic shared with its white-rot ancestors. Four wood degradation-related genes, homologs of which are frequently lost in brown-rot fungi, show signs of pseudogenization in the genome of F. hepatica. These results suggest that transition toward a brown-rot lifestyle could be an ongoing process in F. hepatica. Our results reinforce the idea that wood decay mechanisms are more diverse than initially thought and that the dichotomous separation of wood decay mechanisms in Agaricomycotina into white rot and brown rot should be revisited., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Regressive evolution involves the degeneration of formerly useful structures in a lineage over time, and may be accompanied by the molecular decay of phenotype-specific genes. The mammalian eye has repeatedly undergone degeneration in taxa that occupy dim-light environments including subterranean habitats. Here we assess whether a decrease in the amount of light that reaches the retina is associated with increased regression of retinal genes, whether the phototransduction and visual cycle pathways degrade in a predictable pattern, and if the timing of retinal gene loss is associated with the entrance of mammalian lineages into subterranean environments. Sequence data were obtained from the publically available genomes of the Cape golden mole (Chrysochloris asiatica), naked mole-rat (Heterocephalus glaber) and star-nosed mole (Condylura cristata) for 65 genes associated with phototransduction, the visual cycle, and other retinal functions. Gene sequences were inspected for inactivating mutations and, when present, pseudogene sequences were compared to sequences from subaerial outgroup species. To test whether retinal degeneration is correlated with historical entrances into subterranean environments, estimated dates of retinal gene inactivation were compared to the fossil record and phylogenetic inferences of ancestral fossoriality. Our results show that (1) lower levels of light available to the retina correspond with an increase in the number of retinal pseudogenes, (2) retinal protein networks generally degrade in a predictable manner, although the extensive loss of cone phototransduction genes in Heterocephalus raises further questions regarding SWS1-cone monochromacy versus functional rod monochromacy in this species, and (3) inactivation dates of retinal genes usually post-date inferred entrances into subterranean habitats., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Flanders virus (FLAV) and Hart Park virus (HPV) are rhabdoviruses that circulate in mosquito-bird cycles in the eastern and western United States, respectively, and constitute the only two North American representatives of the Hart Park serogroup. Previously, it was suggested that FLAV is unique among the rhabdoviruses in that it contains two pseudogenes located between the P and M genes, while the cognate sequence for HPV has been lacking. Herein, we demonstrate that FLAV and HPV do not contain pseudogenes in this region, but encode three small functional proteins designated as U1-U3 that apparently arose by gene duplication. To further investigate the U1-U3 region, we conducted the first large-scale evolutionary analysis of a member of the Hart Park serogroup by analyzing over 100 spatially and temporally distinct FLAV isolates. Our phylogeographic analysis demonstrates that although FLAV appears to be slowly evolving, phylogenetically divergent lineages co-circulate sympatrically., (© 2013 Published by Elsevier Inc.)

Deficiency caused by mutations in the creatine transporter gene (SLC6A8/CT1) is an X-linked form of intellectual disability. The presence of highly homologous pseudogenes and high GC content of SLC6A8 genomic sequence complicates the molecular diagnosis of this disorder. To minimize the pseudogene interference, exons 2 to 13 of SLC6A8 were amplified as a single PCR product using gene-specific long-range PCR (LR-PCR) primers. The GC-rich exon 1 and its flanking intronic sequences were amplified separately in a short fragment under GC-rich conditions and a touchdown PCR program. Traditional Sanger sequence analysis of all coding exons of SLC6A8 from a 3-year-old boy with creatine transporter deficiency did not detect deleterious mutations. The long-range PCR product was used as template followed by massively parallel sequencing (MPS) on HiSeq2000. We were able to detect a tandem duplication involving part of exons 11 and 12 in the SLC6A8 gene. The deduced c.1592_1639dup133 mutation was confirmed to be a hemizygous insertion by targeted genomic DNA and cDNA Sanger sequencing. Combination of deep sequencing technology with long-range PCR revealed a novel intragenic duplication in the SLC6A8 gene, providing a definitive molecular diagnosis of creatine transporter deficiency in a male patient., (© 2013 Elsevier Inc. All rights reserved.)

The ecological model, Daphnia pulex (Cladocera: Daphniidae), is broadly distributed in Holarctic freshwater habitats and has been the subject of multidisciplinary study for over half a century, but never has a natural RNA virus infection been reported in daphnids. Here we report on a group of paleoviruses related to RNA dependent RNA polymerase in the genome of D. pulex. Phylogenetic analysis suggests that these paleoviruses are derived from a viral lineage within the genus Phlebovirus. Comparison of the genomic sequences flanking individual paleoviruses reveal that some are orthologous viral insertions having been present in the common ancestor of the D. pulex species complex, which is millions of years old. Still, we detected some sites that have the signature of purifying selection. In contrast, other paleoviruses in this group seem to be unique to specific host lineages and even contain undisrupted open reading frames, suggesting either more recent acquisition, or selective maintenance., (© 2013 Elsevier Inc. All rights reserved.)

There is emerging evidence that human solid tumor cells originate from cancer stem cells (CSCs). In cancer cell lines, tumor-initiating CSCs are mainly found in the side population (SP) that has the capacity to extrude dyes such as Hoechst 33342. We found that Nanog is expressed specifically in SP cells of human gastrointestinal (GI) cancer cells. Nucleotide sequencing revealed that NanogP8 but not Nanog was expressed in GI cancer cells. Transfection of NanogP8 into GI cancer cell lines promoted cell proliferation, while its inhibition by anti-Nanog siRNA suppressed the proliferation. Immunohistochemical staining of primary GI cancer tissues revealed NanogP8 protein to be strongly expressed in 3 out of 60 cases. In these cases, NanogP8 was found especially in an infiltrative part of the tumor, in proliferating cells with Ki67 expression. These data suggest that NanogP8 is involved in GI cancer development in a fraction of patients, in whom it presumably acts by supporting CSC proliferation., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Opsin gene sequences were first reported in the 1980s. The goal of that research was to test the hypothesis that human opsins were members of a single gene family and that variation in human color vision was mediated by mutations in these genes. While the new data supported both hypotheses, the greatest contribution of this work was, arguably, that it provided the data necessary for PCR-based surveys in a diversity of other species. Such studies, and recent whole genome sequencing projects, have uncovered exceptionally large opsin gene repertoires in ray-finned fishes (taxon, Actinopterygii). Guppies and zebrafish, for example, have 10 visual opsin genes each. Here we review the duplication and divergence events that have generated these gene collections. Phylogenetic analyses revealed that large opsin gene repertories in fish have been generated by gene duplication and divergence events that span the age of the ray-finned fishes. Data from whole genome sequencing projects and from large-insert clones show that tandem duplication is the primary mode of opsin gene family expansion in fishes. In some instances gene conversion between tandem duplicates has obscured evolutionary relationships among genes and generated unique key-site haplotypes. We mapped amino acid substitutions at so-called key-sites onto phylogenies and this exposed many examples of convergence. We found that dN/dS values were higher on the branches of our trees that followed gene duplication than on branches that followed speciation events, suggesting that duplication relaxes constraints on opsin sequence evolution. Though the focus of the review is opsin sequence evolution, we also note that there are few clear connections between opsin gene repertoires and variation in spectral environment, morphological traits, or life history traits., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Objectives: The PTEN pseudogene, PTENP1, was recently shown to play a role in cell proliferation in a prostate cancer model. In the present study, we sought to determine whether PTENP1 is expressed in endometrial cancer (EMCA) cell lines and primary tumors along with the microRNAs (miRNAs) that are predicted to regulate PTEN and PTENP1 transcript levels., Methods: RNA was prepared from six EMCA cell lines, three normal endometrial samples, and 61 primary tumors. TaqMan® RT-PCR was used to quantitate PTEN expression in all specimens and PTENP1 expression in cell lines, and normal endometrial (NE) samples. PTENP1 expression was evaluated using conventional RT-PCR in primary tumors. MicroRNA profiling was undertaken using NanoString(TM) technology in AN3CA and KLE cell lines. The relationship between PTEN transcript levels, PTENP1 expression, and PTEN mutation status was investigated., Results: All NE samples, cell lines, and primary tumors expressed PTEN. PTENP1 transcript was expressed in NE, cell lines, and 34/61 (56%) primary tumors. The median relative PTEN level was 2.9 arbitrary expression units in PTENP1-positive tumors and 2.3 in PTENP1-negative tumors (p=0.09). PTEN levels in wild-type and haploinsufficient tumors were variable compared to PTEN-null tumors (p=0.015). Four microRNAs predicted to bind PTEN/PTENP1 ranked in the top 20 most abundant microRNA subtypes in the AN3CA and KLE cell lines., Conclusions: PTENP1 is expressed in NE and EMCA cell lines, as are PTEN/PTENP1 targeting inhibitory miRNAs (cell lines). Further studies are needed to evaluate the impact of PTEN/PTENP1/miRNA interactions on tumorigenesis regulation in EMCA., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Gonadotropin releasing hormone (GnRH) is a key regulator of sexual development and reproduction in vertebrates. Fish have either two or three pre-pro-GnRH genes, encoding structurally distinct peptides. We identified three pre-pro-GnRH genes in Atlantic cod (Gadus morhua, gmGnRH) using RT-PCR, RACE-PCR and BAC DNA library clone sequencing based on synteny searching. Gene identity was confirmed by sequence alignment and subsequent phylogenetic analysis. The expression of these genes was measured by quantitative PCR in the brain and pituitary of female cod throughout their reproductive cycle and in peripheral tissues. All three gmGnRH genes have highly conserved deduced decapeptide sequences, but sequence and phylogenetic data for gmGnRH1 suggest that this is a pseudogene. gmGnRH1 shares low identity with all fish GnRH variants and grouped with the GnRH3 clade. Although gmGnRH1 is a putative pseudogene, it is transcribed in multiple tissues but at low levels in the brain, indicating the loss of conserved hypophysiotrophic function. Phylogenetic analysis reveals that gmGnRH2 and gmGnRH3 variants are located in variant-specific clades. Both gmGnRH2 and gmGnRH3 transcripts are most abundant in the brain, with lower expression in pituitaries and ovaries. Brain gmGnRH3 gene expression increases in spawning fish and is expressed in the pituitary during puberty. Brain gmGnRH2 transcripts are highly expressed relative to gmGnRH3 before and during spawning. Sequence and expression data suggest that gmGnRH1 is a pseudogene and that gmGnRH3 is likely the hypophysiotrophic form of GnRH in Atlantic cod., (Copyright © 2011 Elsevier Inc. All rights reserved.)

The phylogeny of Onychophora (velvet worms) is unresolved and even the monophyly of the two major onychophoran subgroups, Peripatidae and Peripatopsidae, is uncertain. Previous studies of complete mitochondrial genomes from two onychophoran species revealed two strikingly different gene arrangement patterns from highly conserved in a representative of Peripatopsidae to highly derived in a species of Peripatidae, suggesting that these data might be informative for clarifying the onychophoran phylogeny. In order to assess the diversity of mitochondrial genomes among onychophorans, we analyzed the complete mitochondrial genome of Metaperipatus inae, a second representative of Peripatopsidae from Chile. Compared to the proposed ancestral gene order in Onychophora, the mitochondrial genome of M. inae shows dramatic rearrangements, although all protein-coding and ribosomal RNA genes are encoded on the same strands as in the ancestral peripatopsid genome. The retained strand affiliation of all protein-coding and ribosomal RNA genes and the occurrence of three control regions and several pseudo-genes suggest that the derived mitochondrial gene arrangement pattern in M. inae evolved by partial genome duplications, followed by a subsequent loss of redundant genes. Our findings, thus, confirm the diversity of the mitochondrial gene arrangement patterns among onychophorans and support their utility for clarifying the phylogeography of Onychophora, in particular of the Peripatopsidae species from South Africa and Chile., (Copyright 2010 Elsevier Inc. All rights reserved.)

Chronic granulomatous Disease (CGD) is an immunodeficiency disorder affecting about 1 in 250,000 individuals. The disease is caused by mutations in the genes encoding the components of the leukocyte NADPH oxidase. This enzyme produces superoxide, which is essential in the process of intracellular pathogen killing by phagocytic leukocytes. Four of the five genes involved in CGD are autosomal; these are CYBA, encoding p22-phox, NCF2, encoding p67-phox, NCF1, encoding p47-phox, and NCF4, encoding p40-phox. This article lists all mutations identified in these genes in the autosomal forms of CGD. Moreover, polymorphisms in these genes are also given, which should facilitate the recognition of future disease-causing mutations.

Molecular studies of six species from the ancient extant seed plant Cycas, covering a wide range of its morphological diversity and all major areas of distribution, revealed a high level of intra-individual polymorphism of the internal transcribed spacer (ITS1, 5.8S, and ITS2) region, indicative of incomplete nrDNA concerted evolution. Through a range of comparisons of sequence characteristics to functional cDNA ITS copies, including sequence length and substitution variation, GC content, secondary structure stability, the presence of a conserved motif in the 5.8S gene, and evolutionary rates, the PCR amplified divergent genomic DNA ITS paralogs were identified as either putative pseudogenes, recombinants or functional paralogs. This incomplete ITS concerted evolution may be linked to the high number of nucleolar organizer regions in the Cycas genome, and the incomplete lineage sorting due to recent species divergence in the genus. Based on the distribution of a 14 bp deletion, an early evolutionary origin of the pseudogenes is indicated, possibly predating the diversification of Cycas. Due to their early origin combined with the unconstraint evolution of the ITS region in pseudogenes, they accumulate high levels of homoplastic mutations. This leads to random relationships among the pseudogenes due to long-branch attractions, whereas the phylogenetic relationships inferred from the functional ITS paralogs grouped the sequences in species specific clades (except for C. circinalis and C. rumphii). The findings of our extensive study will have a wide significance, for the evolution of these molecular sequences, and their utilization as a major marker for reconstructing phylogenies., ((c) 2009 Elsevier Inc. All rights reserved.)

We have previously reported that some pseudogenes are expressed in Mycobacterium leprae (M. leprae), the causative agent of leprosy, and that their expression levels alter upon infection of macrophages. We attempted to further examine the expression of pseudogene and non-coding genomic region in M. leprae, in this study. 19 Pseudogenes, 17 non-coding genomic regions, and 21 coding genes expression in M. leprae maintained in the footpads of the hypertensive nude rat (SHR/NCrj-rnu) were examined by reverse transcriptase polymerase chain reaction (RT-PCR). The expression of some of these pseudogenes, non-coding genomic regions and coding genes were also examined in M. leprae from skin smear specimens obtained from patients with lepromatous leprosy by RT-PCR. Transcripts from pseudogenes, non-coding genomic regions and coding genes examined in this study were clearly observed in M. leprae. The expression patterns of some of these transcripts vary greatly among different leprosy patients. These results indicate that some of pseudogenes and non-coding genomic regions are transcribed in M. leprae and analysis of RNA expression patterns including pseudogene and non-coding genomic region in M. leprae may be useful in understanding the pathological states of infected patients.

C/D box small nucleolar RNAs (snoRNAs) guide site-specific 2'-O-methylation of RNAs. Nearly all C/D box snoRNAs with known targets are involved in rRNA modification. In vertebrates, snoRNAs are encoded in introns of various genes and their processing is coupled with splicing of host gene pre-mRNA. Here, the genes encoding C/D box snoRNAs that guide 2'-O-methylation of rRNA were identified and analyzed in vertebrate genomes. The number of copies of most C/D box snoRNA genes proved to be lower in placental mammals compared to other vertebrates. This can be due to smaller oocytes and accordingly lower number of ribosomes in them in eutherians. The targets of snoRNAs encoded by single-copy and multiple-copy genes proved to have different distribution in rRNAs. The causes of this difference are discussed. In some cases, the transcripts of homologous C/D box RNA genes were shown to guide the modification of neighboring nucleotides in rRNA. C/D box snoRNA pseudogenes were found in all vertebrate classes. Three novel C/D box snoRNAs were found in Xenopus tropicalis that may guide 2'-O-methylation of Xenopus-specific rRNA sites. A list of 922 annotated C/D box snoRNA genes is presented.

Molecular study of Stone flounder Kareius bicoloratus revealed a high level of intra-individual polymorphism of ribosomal DNA (18S rDNA, ITS1 and 5.8S rDNA). This polymorphism includes not only multiple functional genes but also putative pseudogenes and recombinants. Three major types (Type I, II and III) of 18S rDNA differing in lengths coexisted in the genome of each studied individual. Type II and Type III have a 15-bp-deletion and Type III has an additional 7-bp-deletion relative to Type I. Expression analyses gave clear evidence that Type I 18S rDNA was abundantly expressed. In contrast, transcription of Type II and III 18S rDNA were not detected. Further sequence analysis revealed that Type I was the expected and functional 18S rDNA in K. bicoloratus, while Type II and Type III were pseudogenes. Phylogeny of the ITS1 sequences formed three groups (Group I, II and III), of which Group I linking with Type I 18S rDNA were functional, whereas the remaining two Groups (II and III) linking with Type II and Type III 18S rDNA were pseudogenes. Groups II and III ITS1 sequences were shorter (685-717bp in length) than Group I ITS1 (762-776bp), but had much higher sequence divergence (9.8%) than the later ones did (4.9%). Group II sequences were recombinants between Group I and Group III sequences, providing the possibility that the rDNA of K. bicoloratus is in a transition stage of concerted evolution. In addition, by comparing standard PCR and high-stringency PCR, we concluded that DMSO might prevent functional ITS1 sequences from being amplified instead of enhancing its amplification efficiency. The analysis suggests a non-concerted evolution of rDNA in K. bicoloratus. Features of rDNA functional genes and pseudogenes are described. Possible causes for the lack of homogenization of rDNA paralogues within individuals are discussed.

The phylogeny and phylogeography of the Old World wood mice (subgenus Sylvaemus, genus Apodemus, Muridae) are well-documented. Nevertheless, the distributions of species, such as A. fulvipectus and A. ponticus remain dubious, as well as their phylogenetic relationships with A. sylvaticus. We analysed samples of Apodemus spp. across Europe using the mitochondrial cytochrome-b gene (cyt-b) and compared the DNA and amino-acid compositions of previously published sequences. The main result stemming from this study is the presence of a well-differentiated lineage of Sylvaemus including samples of various species (A. sylvaticus, A. fulvipectus, A. ponticus) from distant locations, which were revealed to be nuclear copies of the mitochondrial cyt-b. The presence of this cryptic pseudogene in published sequences is supported by different pathways. This has led to important errors in previous molecular trees and hence to partial misinterpretations in the phylogeny of Apodemus.

Demonstrating an association between a polymorphism and a disease phenotype through case-control studies requires reliable large-scale genotyping, but accurate measurement of copy number variation has proven to be technically challenging. Here we build on our previous experience with Paralogue Ratio Tests (PRT) to develop PRT copy number determination at the CCL3L1/CCL4L1 copy number variant. A multiplex PRT assay based on four independent comparative PCRs results in a convenient, accurate and robust method of multiallelic copy number measurement suitable for use in large-scale case-control studies, which can unambiguously assign virtually all samples tested to discrete copy number classes.

A pseudogene with 94% similarity to mitochondrial cytochrome c oxidase subunit I (COI) was identified and localized to chromosome 4 of Drosophila ananassae. Because this chromosome is believed to have reduced recombination, its history can be traced using the pseudo-COI sequence. Pseudo-COI sequences were obtained from 27 iso-female lines of six taxa belonging to the D. ananassae species cluster in which reproductive isolation is incomplete. The phylogenetic network constructed from seven recognized haplotypes (#0-#6) indicated that different taxa inhabiting the same geographic area share the haplotypes: #1 from Papua New Guinean populations of D. ananassae and pallidosa-like-Wau; #2 from Papua New Guinean populations of D. ananassae, pallidosa-like, and papuensis-like; and #4 from South Pacific populations of D. ananassae and D. pallidosa. Taxon-K has a unique haplotype (#6), and 18 mutation steps separate it from the closest haplotype, #2. We discuss the possibility of chromosome 4 introgression beyond taxon boundaries.

Molecular studies of 19 species of the genus Pyrus L. revealed different degrees of intra-individual polymorphism of the internal transcribed spacer (ITS 1, 5.8S rDNA and ITS 2) region due to the existence of putative non-functional copies (pseudogenes), putative recombinants and non-concerted evolution among functional copies. Different types of ITS pseudogenes displaying lower GC content and unstable secondary structure were preferentially amplified under normal PCR conditions. Functional ITS copies were successfully obtained in all investigated accessions under the modified PCR conditions. All pseudogenes were highly divergent from their corresponding functional copies and formed a monophyletic group in the phylogenetic tree based on all paralogs, indicating they were of relatively early origin. Functional ITS copies led to confused and poorly resolved phylogeny as a result of low sequence divergence, existence of unidentified ancient recombinants and a high degree of intra-individual functional ITS polymorphism, while certain types of pseudogenes and some relict pseudogenes offered more credible clues for the evolutionary history of Pyrus species.

Currently, the spatial patterns of mitochondrial genes and how the genomic localization of (pseudo)genes originated from mitochondrial DNA remain largely unexplained. The aim of this study was to elucidate the organization of mitochondrial (pseudo)genes given their evolutionary origin. We used a keyword finding method and a bootstrapping method to estimate parameter values that represent the distribution pattern of mitochondrial genes in the nuclear genome. Almost half of mitochondrial genes showing physical clusters were located in the pericentromeric and subtelomeric regions of the chromosome. Most interestingly, the size of these clusters ranged from 0.085 to 3.2 Mb (average+/-SD 1.3+/-0.73 Mb), which coincides with the size of the evolutionary pocket, or the average size of evolutionary breakpoint regions. Our findings imply that the localization of mitochondrial genes in the human genome is determined independent of adaptation.

A substantial proportion of human genes contain tissue-specifically DNA-methylated regions (TDMRs). However, little is known about the evolutionary conservation of differentially methylated loci, how they evolve, and the signals that regulate them. We have studied TDMR conservation in the PLG and TBX gene families and in 32 pseudogene-parental gene pairs. Among the members of the recently evolved PLG gene family, 5'-UTR methylation is conserved and inversely correlated with the cognate gene expression, indicating as well a conserved regulatory role of DNA methylation. Conversely, many genes of the much older TBX family display complementary tissue-specific methylation, suggesting an epigenetic complementation in the evolution of this gene family. Similar to gene families, unprocessed pseudogenes arose from gene duplications and we found TDMR conservation in some pseudogene-parental gene pairs displaying short evolutionary distances. However, for the majority of unprocessed pseudogenes and for all processed pseudogenes examined, we found that tissue-specific methylation arose de novo after gene duplication.

The bee tribe Ctenoplectrini, with two genera, comprises nine species in tropical Africa and ten in Asia and Australia. Most of them collect floral oil, pollen, and nectar from Cucurbitaceae, but three species are thought to be cleptoparasites. The unusual morphology of Ctenoplectrini has made it difficult to infer their closest relatives, in turn preventing an understanding of these bees' geographic and temporal origin. We used two mitochondrial and two nuclear markers (4741 nucleotides) generated for most of the species to test the monophyly of the tribe, its relationships to other Apidae, and its biogeographic history. Ctenoplectrini are strongly supported as monophyletic and closest to the Long-horned bees, Eucerini. The presumably cleptoparasitic species form a clade (Ctenoplectrina) that is sister to the remaining species (Ctenoplectra), confirming the independent evolution of cleptoparasitism in this tribe. Tree topology and molecular dating together suggest that Ctenoplectrini originated in Africa in the Early Eocene and that Ctenoplectra dispersed twice from Africa to Asia, sometime in the Late Eocene, 30-40 my ago, from where one species reached the Australian continent via Indonesia and New Guinea in the mid-Miocene, c. 13 my ago. Dry and cool mid-Miocene climates also coincide with the divergence between Ctenoplectra bequaerti from West Africa and Ctenoplectra terminalis from East and South Africa, perhaps related to fragmentation of the equatorial African rainforest belt.

In the rodent species Microtus cabrerae, males as well as females present several copies of the SRY gene, a single-copy gene located on the Y chromosome in most mammals. Using different PCR approaches, we have characterized the sequence, structure, and organization of the SRY copies and their flanking regions distributed on the X and Y chromosomes of this species. All copies of SRY analyzed, including those from the Y chromosome, proved to be nonfunctional pseudogenes, as they have internal stop codons. In addition, we demonstrated the association of SRY pseudogenes with different fragments of L1 and LTR retroelements in both sex chromosomes of M. cabrerae. Examining the possible origin of SRY pseudogene and retroposons association, we propose that retroposons could have been involved in the mechanism of SRY gene amplification on the Y chromosome and in the transference of the Y-linked SRY copies to the X-chromosome heterochromatin.

A phylogeny of the fungal phylum Basidiomycota is presented based on a survey of 160 taxa and five nuclear genes. Two genes, rpb2, and tef1, are presented in detail. The rpb2 gene is more variable than tef1 and recovers well-supported clades at shallow and deep taxonomic levels. The tef1 gene recovers some deep and ordinal-level relationships but with greater branch support from nucleotides compared to amino acids. Intron placement is dynamic in tef1, often lineage-specific, and diagnostic for many clades. Introns are fewer in rpb2 and tend to be highly conserved by position. When both protein-coding loci are combined with sequences of nuclear ribosomal RNA genes, 18 inclusive clades of Basidiomycota are strongly supported by Bayesian posterior probabilities and 16 by parsimony bootstrapping. These numbers are greater than produced by single genes and combined ribosomal RNA gene regions. Combination of nrDNA with amino acid sequences, or exons with third codon positions removed, produces strong measures of support, particularly for deep internodes of Basidiomycota, which have been difficult to resolve with confidence using nrDNA data alone. This study produces strong boostrap support and significant posterior probabilities for the first time for the following monophyletic groups: (1) Ustilaginomycetes plus Hymenomycetes, (2) an inclusive cluster of hymenochaetoid, corticioid, polyporoid, Thelephorales, russuloid, athelioid, Boletales, and euagarics clades, (3) Thelephorales plus the polyporoid clade, (4) the polyporoid clade, and (5) the cantharelloid clade. Strong support is also recovered for the basal position of the Dacrymycetales in the Hymenomycetidae and paraphyly of the Exobasidiomycetidae.

We applied a comprehensive data-mining strategy to examine the repertoires of rat and mouse odorant receptors (ORs) and type 1 pheromone receptors (V1Rs) using the mm5 (mouse) and rn3 (rat) genomes. We identified 1576 rat OR genes, including 292 pseudogenes. The rat V1R repertoire is composed of 115 intact genes and 72 pseudogenes. The mouse OR and V1R databases were updated using the new assembly mm5, from which 1375 mouse ORs and 308 V1Rs were identified, with more than 100 putative pseudogenes from mm2 now identified as intact because of the higher sequence quality. With these new data we have conducted a series of genomic analyses of the OR and V1R genes from mouse and rat. Orthologous OR clusters were identified in mouse and rat and comparison analysis was performed at three incremental levels: families, coding sequences, and motifs. At the family level, we found that V1R genes have more species-specific families than OR genes. About 20% of intact V1R genes have no orthologous counterpart in the same family, whereas less than 1% of intact ORs are similarly isolated. At the coding sequence level, OR genes are more conserved between mouse and rat than V1R genes. OR genes share greater similarity with their orthologous counterparts than with their closest neighbor, whereas V1R genes show the opposite tendency. Motifs were identified to obtain biological insights. Motifs specific for species or families were found in OR and V1R genes, which may result in the differential pheromone-dependent behaviors and perception of odors between mouse and rat.

Chromosome deletions do abound in cancer and are detected in certain regions in a non-random manner. Although their relevance remains elusive, it is a general agreement that segmental losses provide the cell with selective growth advantage. Consequently these may contain genes and/or regulatory sequences that control normal growth and inhibit malignancy. We have developed a monochromosomal hybrid based experimental model for the generation and functional analysis of deletions, that is called "elimination test" (Et). Focused on human chromosome 3 - that was known to carry multiple 3p deletions - the Et was expected to restrict a 3p tumor suppressor region to a sufficiently small segment that permits the selection of a critically important candidate gene. Surprisingly, we detected three regions that were lost in all or majority of tumors: CER1 (3p21.3, Mb: 43.32-45.74), CER2 (3p22, Mb: 37.83-39.06) and FER (3p14.3-p21.2, Mb: 50.12-58.03). In contrast a 3q26-qter region (CRR) was regularly retained. CER1 - our main focus - contains multiple genes that may inhibit tumor growth, but 3 genes, RIS1, LF (LTF) and LIMD1 have already the necessary experimental support to be considered bona fide tumor suppressors. Tumor suppressor region borders display instability features including: (1) they break in evolution and in tumors, (2) they evolve horizontally, and (3) they are enriched with pseudogene insertions. The most remarkable features at the breakpoint cluster regions were segmental duplications that drive horizontal evolution and contribute to cancer associated instability.

Molecular studies of 21 species of the large Cactaceae genus Mammillaria representing a variety of intrageneric taxonomic levels revealed a high degree of intra-individual polymorphism of the internal transcribed spacer region (ITS1, 5.8S rDNA, ITS2). Only a few of these ITS copies belong to apparently functional genes, whereas most are probably non-functional (pseudogenes). As a multiple gene family, the ITS region is subjected to concerted evolution. However, the high degree of intra-individual polymorphism of up to 36% in ITS1 and up to 35% in ITS2 suggests a non-concerted evolution of these loci in Mammillaria. Conserved angiosperm motifs of ITS1 and ITS2 were compared between genomic and cDNA ITS clones of Mammillaria. Some of these motifs (e.g., ITS1 motif 1, 'TGGT' within ITS2) in combination with the determination of GC-content, length comparisons of the spacers and ITS2 secondary structure (helices II and III) are helpful in the identification of pseudogene rDNA regions.

We have discovered a family of small secreted proteins in Homo sapiens and Mus musculus. The IGF-like (IGFL) genes encode proteins of approximately 100 amino acids that contain 11 conserved cysteine residues at fixed positions, including two CC motifs. In H. sapiens, the family is composed of four genes and two pseudogenes that are referred as IGFL1 to IGFL4 and IGFL1P1 and IGFL1P2, respectively. Human IGFL genes are clustered together on chromosome 19 within a 35-kb interval. M. musculus has a single IGFL family member that is located on chromosome 7. Further, evolutionary analysis shows a lack of direct orthology between any of the four human members and the mouse gene. This relationship between the mouse and the human family members suggests that the multiple members in the human complement have arisen from recent duplication events that appear limited to the primate lineage. Structural considerations and sequence comparisons would suggest that IGFL proteins are distantly related to the IGF superfamily of growth factors. IGFL mRNAs display specific expression patterns; they are expressed in fetal tissues, breast, and prostate, and in many cancers as well, and this pattern is consistent with that of the IGF family members.