Search

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17β-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

Several species in the oomycete genus Peronosclerospora cause downy mildew on maize and can result in significant yield losses in Asia. Bio-surveillance of these pathogens is a high priority to prevent epidemics on maize in the United States and consequent damage to the US economy. The unresolved taxonomy and dearth of molecular resources for Peronosclerospora spp. hinder these efforts. P. sorghi is a pathogen of sorghum and maize with a global distribution, for which limited diversity has been detected in the southern USA. We characterized the genome, transcriptome, and mitogenome of an isolate, representing the US pathotype 6. The highly homozygous genome was assembled using 10× Genomics linked reads and scaffolded using Hi-C into 13 chromosomes. The total assembled length was 303.2 Mb, larger than any other oomycete previously assembled. The mitogenome was 38 kb, similar in size to other oomycetes, although it had a unique gene order. Nearly 20,000 genes were annotated in the nuclear genome, more than described for other downy mildew causing oomycetes. The 13 chromosomes of P. sorghi were highly syntenic with the 17 chromosomes of Peronospora effusa with conserved centromeric regions and distinct chromosomal fusions. The increased assembly size and gene count of P. sorghi is due to extensive retrotransposition, resulting in putative pseudogenization. Ancestral genes had higher transcript abundance and were enriched for differential expression. This study provides foundational resources for analysis of Peronosclerospora and comparisons to other oomycete genera. Further genomic studies of global Peronosclerospora spp. will determine the suitability of the mitogenome, ancestral genes, and putative pseudogenes for marker development and taxonomic relationships.

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

The cricetine rodents Peromyscus leucopus and P. maniculatus are key reservoirs for several zoonotic diseases in North America. We determined the complete circular mitochondrial genome sequences of representatives of 3 different stock colonies of P. leucopus, one stock colony of P. maniculatus and two wild populations of P. leucopus. The genomes were syntenic with that of the murids Mus musculus and Rattus norvegicus. Phylogenetic analysis confirmed that these two Peromyscus species are sister taxa in a clade with P. polionotus and also uncovered a distinction between P. leucopus populations in the eastern and the central United States. In one P. leucopus lineage four extended regions of mitochondrial pseudogenes were identified in the nuclear genome. RNA-seq analysis revealed transcription of the entire genome and differences from controls in the expression profiles of mitochondrial genes in the blood, but not in liver or brain, of animals infected with the zoonotic pathogen Borrelia hermsii. PCR and sequencing of the D-loop of the mitochondrion identified 32 different haplotypes among 118 wild P. leucopus at a Connecticut field site. These findings help to further establish P. leucopus as a model organism for studies of emerging infectious diseases, ecology, and in other disciplines.

Statins are the most commonly prescribed cardiovascular disease drug, but their inter-individual efficacy varies considerably. Genetic factors uncovered to date have only explained a small proportion of variation in low-density lipoprotein cholesterol (LDLC) lowering. To identify novel markers and determinants of statin response, we used whole transcriptome sequence data collected from simvastatin and control incubated lymphoblastoid cell lines (LCLs) established from participants of the Cholesterol and Pharmacogenetics (CAP) simvastatin clinical trial. We looked for genes whose statin-induced expression changes were most different between LCLs derived from individuals with high versus low plasma LDLC statin response during the CAP trial. We created a classification model of 82 "signature" gene expression changes that distinguished high versus low LDLC statin response. One of the most differentially changing genes was zinc finger protein 542 pseudogene (ZNF542P), the signature gene with changes most correlated with statin-induced change in cellular cholesterol ester, an in vitro marker of statin response. ZNF542P knock-down in a human hepatoma cell line increased intracellular cholesterol ester levels upon simvastatin treatment. Together, these findings imply a role for ZNF542P in LDLC response to simvastatin and, importantly, highlight the potential significance of noncoding RNAs as a contributing factor to variation in drug response.

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

BackgroundPseudogenes are non-functional sequences in the genome with homologous sequences that are functional (i.e. genes). They are abundant in eukaryotes where they have been extensively investigated, while in prokaryotes they are significantly scarcer and less well studied. Here we conduct a comprehensive analysis of the evolution of orthologs of Mycobacterium leprae pseudogenes in prokaryotes. The leprosy pathogen M. leprae is of particular interest since it contains an unusually large number of pseudogenes, comprising approximately 40% of its entire genome. The analysis is conducted in both broad and narrow phylogenetic ranges.ResultsWe have developed an informatics-based approach to characterize the evolution of pseudogenes. This approach combines tools from phylogenomics, genomics, and transcriptomics. The results we obtain are used to assess the contributions of two mechanisms for pseudogene formation: failed horizontal gene transfer events and disruption of native genes.ConclusionsWe conclude that, although it was reported that in most bacteria the former is most likely responsible for the majority of pseudogenization events, in mycobacteria, and in particular in M. leprae with its exceptionally high pseudogene numbers, the latter predominates. We believe that our study sheds new light on the evolution of pseudogenes in bacteria, by utilizing new methodologies that are applied to the unusually abundant M. leprae pseudogenes and their orthologs.

Statins are the most commonly prescribed cardiovascular disease drug, but their inter-individual efficacy varies considerably. Genetic factors uncovered to date have only explained a small proportion of variation in low-density lipoprotein cholesterol (LDLC) lowering. To identify novel markers and determinants of statin response, we used whole transcriptome sequence data collected from simvastatin and control incubated lymphoblastoid cell lines (LCLs) established from participants of the Cholesterol and Pharmacogenetics (CAP) simvastatin clinical trial. We looked for genes whose statin-induced expression changes were most different between LCLs derived from individuals with high versus low plasma LDLC statin response during the CAP trial. We created a classification model of 82 "signature" gene expression changes that distinguished high versus low LDLC statin response. One of the most differentially changing genes was zinc finger protein 542 pseudogene (ZNF542P), the signature gene with changes most correlated with statin-induced change in cellular cholesterol ester, an in vitro marker of statin response. ZNF542P knock-down in a human hepatoma cell line increased intracellular cholesterol ester levels upon simvastatin treatment. Together, these findings imply a role for ZNF542P in LDLC response to simvastatin and, importantly, highlight the potential significance of noncoding RNAs as a contributing factor to variation in drug response.

BackgroundPseudogenes are non-functional sequences in the genome with homologous sequences that are functional (i.e. genes). They are abundant in eukaryotes where they have been extensively investigated, while in prokaryotes they are significantly scarcer and less well studied. Here we conduct a comprehensive analysis of the evolution of orthologs of Mycobacterium leprae pseudogenes in prokaryotes. The leprosy pathogen M. leprae is of particular interest since it contains an unusually large number of pseudogenes, comprising approximately 40% of its entire genome. The analysis is conducted in both broad and narrow phylogenetic ranges.ResultsWe have developed an informatics-based approach to characterize the evolution of pseudogenes. This approach combines tools from phylogenomics, genomics, and transcriptomics. The results we obtain are used to assess the contributions of two mechanisms for pseudogene formation: failed horizontal gene transfer events and disruption of native genes.ConclusionsWe conclude that, although it was reported that in most bacteria the former is most likely responsible for the majority of pseudogenization events, in mycobacteria, and in particular in M. leprae with its exceptionally high pseudogene numbers, the latter predominates. We believe that our study sheds new light on the evolution of pseudogenes in bacteria, by utilizing new methodologies that are applied to the unusually abundant M. leprae pseudogenes and their orthologs.

BackgroundAnalyzing close species with diverse developmental modes is instrumental for investigating the evolutionary significance of physiological, anatomical and behavioral features at a molecular level. Many examples of trait loss are known in metazoan populations living in dark environments. Tunicates are the closest living relatives of vertebrates and typically present a lifecycle with distinct motile larval and sessile adult stages. The nervous system of the motile larva contains melanized cells associated with geotactic and light-sensing organs. It has been suggested that these are homologous to vertebrate neural crest-derived melanocytes. Probably due to ecological adaptation to distinct habitats, several species of tunicates in the Molgulidae family have tailless (anural) larvae that fail to develop sensory organ-associated melanocytes. Here we studied the evolution of Tyrosinase family genes, indispensible for melanogenesis, in the anural, unpigmented Molgula occulta and in the tailed, pigmented Molgula oculata by using phylogenetic, developmental and molecular approaches.ResultsWe performed an evolutionary reconstruction of the tunicate Tyrosinase gene family: in particular, we found that M. oculata possesses genes predicted to encode one Tyrosinase (Tyr) and three Tyrosinase-related proteins (Tyrps) while M. occulta has only Tyr and Tyrp.a pseudogenes that are not likely to encode functional proteins. Analysis of Tyr sequences from various M. occulta individuals indicates that different alleles independently acquired frameshifting short indels and/or larger mobile genetic element insertions, resulting in pseudogenization of the Tyr locus. In M. oculata, Tyr is expressed in presumptive pigment cell precursors as in the model tunicate Ciona robusta. Furthermore, a M. oculata Tyr reporter gene construct was active in the pigment cell precursors of C. robusta embryos, hinting at conservation of the regulatory network underlying Tyr expression in tunicates. In contrast, we did not observe any expression of the Tyr pseudogene in M. occulta embryos. Similarly, M. occulta Tyr allele expression was not rescued in pigmented interspecific M. occulta × M. oculata hybrid embryos, suggesting deleterious mutations also to its cis-regulatory sequences. However, in situ hybridization for transcripts from the M. occulta Tyrp.a pseudogene revealed its expression in vestigial pigment cell precursors in this species.ConclusionsWe reveal a complex evolutionary history of the melanogenesis pathway in tunicates, characterized by distinct gene duplication and loss events. Our expression and molecular data support a tight correlation between pseudogenization of Tyrosinase family members and the absence of pigmentation in the immotile larvae of M. occulta. These results suggest that relaxation of purifying selection has resulted in the loss of sensory organ-associated melanocytes and core genes in the melanogenesis biosynthetic pathway in M. occulta.

BackgroundGenome duplication has played a pivotal role in the evolution of many eukaryotic lineages, including the vertebrates. A relatively recent vertebrate genome duplication is that in Xenopus laevis, which resulted from the hybridization of two closely related species about 17 million years ago. However, little is known about the consequences of this duplication at the level of the genome, the epigenome, and gene expression.ResultsThe X. laevis genome consists of two subgenomes, referred to as L (long chromosomes) and S (short chromosomes), that originated from distinct diploid progenitors. Of the parental subgenomes, S chromosomes have degraded faster than L chromosomes from the point of genome duplication until the present day. Deletions appear to have the largest effect on pseudogene formation and loss of regulatory regions. Deleted regions are enriched for long DNA repeats and the flanking regions have high alignment scores, suggesting that non-allelic homologous recombination has played a significant role in the loss of DNA. To assess innovations in the X. laevis subgenomes we examined p300-bound enhancer peaks that are unique to one subgenome and absent from X. tropicalis. A large majority of new enhancers comprise transposable elements. Finally, to dissect early and late events following interspecific hybridization, we examined the epigenome and the enhancer landscape in X. tropicalis × X. laevis hybrid embryos. Strikingly, young X. tropicalis DNA transposons are derepressed and recruit p300 in hybrid embryos.ConclusionsThe results show that erosion of X. laevis genes and functional regulatory elements is associated with repeats and non-allelic homologous recombination and furthermore that young repeats have also contributed to the p300-bound regulatory landscape following hybridization and whole-genome duplication.

BackgroundGenome duplication has played a pivotal role in the evolution of many eukaryotic lineages, including the vertebrates. A relatively recent vertebrate genome duplication is that in Xenopus laevis, which resulted from the hybridization of two closely related species about 17 million years ago. However, little is known about the consequences of this duplication at the level of the genome, the epigenome, and gene expression.ResultsThe X. laevis genome consists of two subgenomes, referred to as L (long chromosomes) and S (short chromosomes), that originated from distinct diploid progenitors. Of the parental subgenomes, S chromosomes have degraded faster than L chromosomes from the point of genome duplication until the present day. Deletions appear to have the largest effect on pseudogene formation and loss of regulatory regions. Deleted regions are enriched for long DNA repeats and the flanking regions have high alignment scores, suggesting that non-allelic homologous recombination has played a significant role in the loss of DNA. To assess innovations in the X. laevis subgenomes we examined p300-bound enhancer peaks that are unique to one subgenome and absent from X. tropicalis. A large majority of new enhancers comprise transposable elements. Finally, to dissect early and late events following interspecific hybridization, we examined the epigenome and the enhancer landscape in X. tropicalis × X. laevis hybrid embryos. Strikingly, young X. tropicalis DNA transposons are derepressed and recruit p300 in hybrid embryos.ConclusionsThe results show that erosion of X. laevis genes and functional regulatory elements is associated with repeats and non-allelic homologous recombination and furthermore that young repeats have also contributed to the p300-bound regulatory landscape following hybridization and whole-genome duplication.

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

Extracellular factors belonging to the TGF-β family play pivotal roles in the formation and patterning of germ layers during early Xenopus embryogenesis. Here, we show that the vg1 and nodal3 genes of Xenopus laevis are present in gene clusters on chromosomes XLA1L and XLA3L, respectively, and that both gene clusters have been completely lost from the syntenic S chromosome regions. The presence of gene clusters and chromosome-specific gene loss were confirmed by cDNA FISH analyses. Sequence and expression analyses revealed that paralogous genes in the vg1 and nodal3 clusters on the L chromosomes were also altered compared to their Xenopus tropicalis orthologs. X. laevis vg1 and nodal3 paralogs have potentially become pseudogenes or sub-functionalized genes and are expressed at different levels. As X. tropicalis has a single vg1 gene on chromosome XTR1, the ancestral vg1 gene in X. laevis appears to have been expanded on XLA1L. Of note, two reported vg1 genes, vg1(S20) and vg1(P20), reside in the cluster on XLA1L. The nodal3 gene cluster is also present on X. tropicalis chromosome XTR3, but phylogenetic analysis indicates that nodal3 genes in X. laevis and X. tropicalis were independently expanded and/or evolved in concert within each cluster by gene conversion. These findings provide insights into the function and molecular evolution of TGF-β family genes in response to allotetraploidization.

Extracellular factors belonging to the TGF-β family play pivotal roles in the formation and patterning of germ layers during early Xenopus embryogenesis. Here, we show that the vg1 and nodal3 genes of Xenopus laevis are present in gene clusters on chromosomes XLA1L and XLA3L, respectively, and that both gene clusters have been completely lost from the syntenic S chromosome regions. The presence of gene clusters and chromosome-specific gene loss were confirmed by cDNA FISH analyses. Sequence and expression analyses revealed that paralogous genes in the vg1 and nodal3 clusters on the L chromosomes were also altered compared to their Xenopus tropicalis orthologs. X. laevis vg1 and nodal3 paralogs have potentially become pseudogenes or sub-functionalized genes and are expressed at different levels. As X. tropicalis has a single vg1 gene on chromosome XTR1, the ancestral vg1 gene in X. laevis appears to have been expanded on XLA1L. Of note, two reported vg1 genes, vg1(S20) and vg1(P20), reside in the cluster on XLA1L. The nodal3 gene cluster is also present on X. tropicalis chromosome XTR3, but phylogenetic analysis indicates that nodal3 genes in X. laevis and X. tropicalis were independently expanded and/or evolved in concert within each cluster by gene conversion. These findings provide insights into the function and molecular evolution of TGF-β family genes in response to allotetraploidization.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

BackgroundDNA cytosine methylation is an important epigenetic modification that has significant effects on a variety of biological processes in animals. Avian species hold a crucial position in evolutionary history. In this study, we used whole-genome bisulfite sequencing (MethylC-seq) to generate single base methylation profiles of lungs in two genetically distinct and highly inbred chicken lines (Fayoumi and Leghorn) that differ in genetic resistance to multiple pathogens, and we explored the potential regulatory role of DNA methylation associated with immune response differences between the two chicken lines.MethodsThe MethylC-seq was used to generate single base DNA methylation profiles of Fayoumi and Leghorn birds. In addition, transcriptome profiling using RNA-seq from the same chickens and tissues were obtained to interrogate how DNA methylation regulates gene transcription on a genome-wide scale.ResultsThe general DNA methylation pattern across different regions of genes was conserved compared to other species except for hyper-methylation of repeat elements, which was not observed in chicken. The methylation level of miRNA and pseudogene promoters was high, which indicates that silencing of these genes may be partially due to promoter hyper-methylation. Interestingly, the promoter regions of more recently evolved genes tended to be more highly methylated, whereas the gene body regions of evolutionarily conserved genes were more highly methylated than those of more recently evolved genes. Immune-related GO (Gene Ontology) terms were significantly enriched from genes within the differentially methylated regions (DMR) between Fayoumi and Leghorn, which implicates DNA methylation as one of the regulatory mechanisms modulating immune response differences between these lines.ConclusionsThis study establishes a single-base resolution DNA methylation profile of chicken lung and suggests a regulatory role of DNA methylation in controlling gene expression and maintaining genome transcription stability. Furthermore, profiling the DNA methylomes of two genetic lines that differ in disease resistance provides a unique opportunity to investigate the potential role of DNA methylation in host disease resistance. Our study provides a foundation for future studies on epigenetic modulation of host immune response to pathogens in chickens.

BackgroundDNA cytosine methylation is an important epigenetic modification that has significant effects on a variety of biological processes in animals. Avian species hold a crucial position in evolutionary history. In this study, we used whole-genome bisulfite sequencing (MethylC-seq) to generate single base methylation profiles of lungs in two genetically distinct and highly inbred chicken lines (Fayoumi and Leghorn) that differ in genetic resistance to multiple pathogens, and we explored the potential regulatory role of DNA methylation associated with immune response differences between the two chicken lines.MethodsThe MethylC-seq was used to generate single base DNA methylation profiles of Fayoumi and Leghorn birds. In addition, transcriptome profiling using RNA-seq from the same chickens and tissues were obtained to interrogate how DNA methylation regulates gene transcription on a genome-wide scale.ResultsThe general DNA methylation pattern across different regions of genes was conserved compared to other species except for hyper-methylation of repeat elements, which was not observed in chicken. The methylation level of miRNA and pseudogene promoters was high, which indicates that silencing of these genes may be partially due to promoter hyper-methylation. Interestingly, the promoter regions of more recently evolved genes tended to be more highly methylated, whereas the gene body regions of evolutionarily conserved genes were more highly methylated than those of more recently evolved genes. Immune-related GO (Gene Ontology) terms were significantly enriched from genes within the differentially methylated regions (DMR) between Fayoumi and Leghorn, which implicates DNA methylation as one of the regulatory mechanisms modulating immune response differences between these lines.ConclusionsThis study establishes a single-base resolution DNA methylation profile of chicken lung and suggests a regulatory role of DNA methylation in controlling gene expression and maintaining genome transcription stability. Furthermore, profiling the DNA methylomes of two genetic lines that differ in disease resistance provides a unique opportunity to investigate the potential role of DNA methylation in host disease resistance. Our study provides a foundation for future studies on epigenetic modulation of host immune response to pathogens in chickens.

Not all pseudogenes are transcriptionally silent as previously thought. Pseudogene transcripts, although not translated, contribute to the non-coding RNA pool of the cell that regulates the expression of other genes. Pseudogene transcripts can also directly compete with the parent gene transcripts for mRNA stability and other cell factors, modulating their expression levels. Tissue-specific and cancer-specific differential expression of these "functional" pseudogenes has been reported. To ascertain potential pseudogene:gene interactions in HIV-1 infection, we analyzed transcriptomes from infected and uninfected T-cells and found that 21 pseudogenes are differentially expressed in HIV-1 infection. This is interesting because parent genes of one-third of these differentially-expressed pseudogenes are implicated in HIV-1 life cycle, and parent genes of half of these pseudogenes are involved in different viral infections. Our bioinformatics analysis identifies candidate pseudogene:gene interactions that may be of significance in HIV-1 infection. Experimental validation of these interactions would establish that retroviruses exploit this newly-discovered layer of host gene expression regulation for their own benefit.

Not all pseudogenes are transcriptionally silent as previously thought. Pseudogene transcripts, although not translated, contribute to the non-coding RNA pool of the cell that regulates the expression of other genes. Pseudogene transcripts can also directly compete with the parent gene transcripts for mRNA stability and other cell factors, modulating their expression levels. Tissue-specific and cancer-specific differential expression of these "functional" pseudogenes has been reported. To ascertain potential pseudogene:gene interactions in HIV-1 infection, we analyzed transcriptomes from infected and uninfected T-cells and found that 21 pseudogenes are differentially expressed in HIV-1 infection. This is interesting because parent genes of one-third of these differentially-expressed pseudogenes are implicated in HIV-1 life cycle, and parent genes of half of these pseudogenes are involved in different viral infections. Our bioinformatics analysis identifies candidate pseudogene:gene interactions that may be of significance in HIV-1 infection. Experimental validation of these interactions would establish that retroviruses exploit this newly-discovered layer of host gene expression regulation for their own benefit.

Rod monochromacy is a rare condition in vertebrates characterized by the absence of cone photoreceptor cells. The resulting phenotype is colourblindness and low acuity vision in dim-light and blindness in bright-light conditions. Early reports of xenarthrans (armadillos, sloths and anteaters) suggest that they are rod monochromats, but this has not been tested with genomic data. We searched the genomes of Dasypus novemcinctus (nine-banded armadillo), Choloepus hoffmanni (Hoffmann's two-toed sloth) and Mylodon darwinii (extinct ground sloth) for retinal photoreceptor genes and examined them for inactivating mutations. We performed PCR and Sanger sequencing on cone phototransduction genes of 10 additional xenarthrans to test for shared inactivating mutations and estimated the timing of inactivation for photoreceptor pseudogenes. We concluded that a stem xenarthran became an long-wavelength sensitive-cone monochromat following a missense mutation at a critical residue in SWS1, and a stem cingulate (armadillos, glyptodonts and pampatheres) and stem pilosan (sloths and anteaters) independently acquired rod monochromacy early in their evolutionary history following the inactivation of LWS and PDE6C, respectively. We hypothesize that rod monochromacy in armadillos and pilosans evolved as an adaptation to a subterranean habitat in the early history of Xenarthra. The presence of rod monochromacy has major implications for understanding xenarthran behavioural ecology and evolution.

Here we present the results of a large-scale bioinformatics annotation of non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models of hand-curated families from the Rfam database to infer conserved RNA families within each avian genome. We supplement these annotations with predictions from the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We identify 34 lncRNA-associated loci that are conserved between birds and mammals and validate 12 of these in chicken. We report several intriguing cases where a reported mammalian lncRNA, but not its function, is conserved. We also demonstrate extensive conservation of classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g., snoRNAs and miRNAs) in birds. Furthermore, we describe numerous "losses" of several RNA families, and attribute these to either genuine loss, divergence or missing data. In particular, we show that many of these losses are due to the challenges associated with assembling avian microchromosomes. These combined results illustrate the utility of applying homology-based methods for annotating novel vertebrate genomes.

Interferon lambda 4 gene (IFNL4) encodes IFN-λ4, a new member of the IFN-λ family with antiviral activity. In humans IFNL4 open reading frame is truncated by a polymorphic frame-shift insertion that eliminates IFN-λ4 and turns IFNL4 into a polymorphic pseudogene. Functional IFN-λ4 has antiviral activity but the elimination of IFN-λ4 through pseudogenization is strongly associated with improved clearance of hepatitis C virus (HCV) infection. We show that functional IFN-λ4 is conserved and evolutionarily constrained in mammals and thus functionally relevant. However, the pseudogene has reached moderately high frequency in Africa, America, and Europe, and near fixation in East Asia. In fact, the pseudogenizing variant is among the 0.8% most differentiated SNPs between Africa and East Asia genome-wide. Its raise in frequency is associated with additional evidence of positive selection, which is strongest in East Asia, where this variant falls in the 0.5% tail of SNPs with strongest signatures of recent positive selection genome-wide. Using a new Approximate Bayesian Computation (ABC) approach we infer that the pseudogenizing allele appeared just before the out-of-Africa migration and was immediately targeted by moderate positive selection; selection subsequently strengthened in European and Asian populations resulting in the high frequency observed today. This provides evidence for a changing adaptive process that, by favoring IFN-λ4 inactivation, has shaped present-day phenotypic diversity and susceptibility to disease.

Mapping reads to a reference sequence is a common step when analyzing allele effects in high-throughput sequencing data. The choice of reference is critical because its effect on quantitative sequence analysis is non-negligible. Recent studies suggest aligning to a single standard reference sequence, as is common practice, can lead to an underlying bias depending on the genetic distances of the target sequences from the reference. To avoid this bias, researchers have resorted to using modified reference sequences. Even with this improvement, various limitations and problems remain unsolved, which include reduced mapping ratios, shifts in read mappings and the selection of which variants to include to remove biases. To address these issues, we propose a novel and generic multi-alignment pipeline. Our pipeline integrates the genomic variations from known or suspected founders into separate reference sequences and performs alignments to each one. By mapping reads to multiple reference sequences and merging them afterward, we are able to rescue more reads and diminish the bias caused by using a single common reference. Moreover, the genomic origin of each read is determined and annotated during the merging process, providing a better source of information to assess differential expression than simple allele queries at known variant positions. Using RNA-seq of a diallel cross, we compare our pipeline with the single-reference pipeline and demonstrate our advantages of more aligned reads and a higher percentage of reads with assigned origins. Database URL: http://csbio.unc.edu/CCstatus/index.py?run=Pseudo.

Mapping reads to a reference sequence is a common step when analyzing allele effects in high-throughput sequencing data. The choice of reference is critical because its effect on quantitative sequence analysis is non-negligible. Recent studies suggest aligning to a single standard reference sequence, as is common practice, can lead to an underlying bias depending on the genetic distances of the target sequences from the reference. To avoid this bias, researchers have resorted to using modified reference sequences. Even with this improvement, various limitations and problems remain unsolved, which include reduced mapping ratios, shifts in read mappings and the selection of which variants to include to remove biases. To address these issues, we propose a novel and generic multi-alignment pipeline. Our pipeline integrates the genomic variations from known or suspected founders into separate reference sequences and performs alignments to each one. By mapping reads to multiple reference sequences and merging them afterward, we are able to rescue more reads and diminish the bias caused by using a single common reference. Moreover, the genomic origin of each read is determined and annotated during the merging process, providing a better source of information to assess differential expression than simple allele queries at known variant positions. Using RNA-seq of a diallel cross, we compare our pipeline with the single-reference pipeline and demonstrate our advantages of more aligned reads and a higher percentage of reads with assigned origins. Database URL: http://csbio.unc.edu/CCstatus/index.py?run=Pseudo.

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the β-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the β-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

Interferon lambda 4 gene (IFNL4) encodes IFN-λ4, a new member of the IFN-λ family with antiviral activity. In humans IFNL4 open reading frame is truncated by a polymorphic frame-shift insertion that eliminates IFN-λ4 and turns IFNL4 into a polymorphic pseudogene. Functional IFN-λ4 has antiviral activity but the elimination of IFN-λ4 through pseudogenization is strongly associated with improved clearance of hepatitis C virus (HCV) infection. We show that functional IFN-λ4 is conserved and evolutionarily constrained in mammals and thus functionally relevant. However, the pseudogene has reached moderately high frequency in Africa, America, and Europe, and near fixation in East Asia. In fact, the pseudogenizing variant is among the 0.8% most differentiated SNPs between Africa and East Asia genome-wide. Its raise in frequency is associated with additional evidence of positive selection, which is strongest in East Asia, where this variant falls in the 0.5% tail of SNPs with strongest signatures of recent positive selection genome-wide. Using a new Approximate Bayesian Computation (ABC) approach we infer that the pseudogenizing allele appeared just before the out-of-Africa migration and was immediately targeted by moderate positive selection; selection subsequently strengthened in European and Asian populations resulting in the high frequency observed today. This provides evidence for a changing adaptive process that, by favoring IFN-λ4 inactivation, has shaped present-day phenotypic diversity and susceptibility to disease.

BackgroundWe describe the genome of the western painted turtle, Chrysemys picta bellii, one of the most widespread, abundant, and well-studied turtles. We place the genome into a comparative evolutionary context, and focus on genomic features associated with tooth loss, immune function, longevity, sex differentiation and determination, and the species' physiological capacities to withstand extreme anoxia and tissue freezing.ResultsOur phylogenetic analyses confirm that turtles are the sister group to living archosaurs, and demonstrate an extraordinarily slow rate of sequence evolution in the painted turtle. The ability of the painted turtle to withstand complete anoxia and partial freezing appears to be associated with common vertebrate gene networks, and we identify candidate genes for future functional analyses. Tooth loss shares a common pattern of pseudogenization and degradation of tooth-specific genes with birds, although the rate of accumulation of mutations is much slower in the painted turtle. Genes associated with sex differentiation generally reflect phylogeny rather than convergence in sex determination functionality. Among gene families that demonstrate exceptional expansions or show signatures of strong natural selection, immune function and musculoskeletal patterning genes are consistently over-represented.ConclusionsOur comparative genomic analyses indicate that common vertebrate regulatory networks, some of which have analogs in human diseases, are often involved in the western painted turtle's extraordinary physiological capacities. As these regulatory pathways are analyzed at the functional level, the painted turtle may offer important insights into the management of a number of human health disorders.

Pseudogenes have been shown to acquire unique regulatory roles from more and more organisms. We report the observation of a cluster of siRNAs derived from pseudogenes of African Trypanosoma brucei using high through-put analysis. We show that these pseudogene-derived siRNAs suppress gene expression through RNA interference. The discovery that siRNAs may originate from pseudogenes and regulate gene expression in a unicellular eukaryote provides insights into the functional roles of pseudogenes and into the origin of noncoding small RNAs.

Pseudogenes have been shown to acquire unique regulatory roles from more and more organisms. We report the observation of a cluster of siRNAs derived from pseudogenes of African Trypanosoma brucei using high through-put analysis. We show that these pseudogene-derived siRNAs suppress gene expression through RNA interference. The discovery that siRNAs may originate from pseudogenes and regulate gene expression in a unicellular eukaryote provides insights into the functional roles of pseudogenes and into the origin of noncoding small RNAs.

Anaplasma phagocytophilum is a tick-transmitted bacterial pathogen of humans and other animals, and is an obligate intracellular parasite. Throughout the course of infection, hosts acquire temporary resistance to granulocytic anaplasmosis as they develop immunity specific for the major antigen, major surface protein 2 (Msp2). However, the bacterium then utilizes a novel recombination mechanism shuffling functional pseudogenes sequentially into an expression cassette with conserved 5' and 3' ends, bypassing host immunity. Approximately 100 pseudogenes are present in the only fully sequenced human-origin HZ genome, representing the possibility for almost unlimited antigenic diversity. In the present study, we identified a select group of 20% of the A. phagocytophilum HZ msp2 pseudogenes that have matched preferentially to human, canine, and equine expression cassettes. Pseudogenes cluster predominantly in one spatial run limited to a single genomic island in less than 50% of the genome but phylogenetically related pseudogenes are neither necessarily located in close proximity on the genome nor share similar percent identity with expression cassettes. Pseudogenes near the expression cassette (and the origin) are more likely to be expressed than those farther away. Taken together, these findings suggest that there may be natural selection pressure to retain pseudogenes in one cluster near the putative origin of replication, even though global recombination shuffles pseudogenes around the genome, separating pseudogenes that share genetic origins as well as those with similar identities.

Anaplasma phagocytophilum is a tick-transmitted bacterial pathogen of humans and other animals, and is an obligate intracellular parasite. Throughout the course of infection, hosts acquire temporary resistance to granulocytic anaplasmosis as they develop immunity specific for the major antigen, major surface protein 2 (Msp2). However, the bacterium then utilizes a novel recombination mechanism shuffling functional pseudogenes sequentially into an expression cassette with conserved 5' and 3' ends, bypassing host immunity. Approximately 100 pseudogenes are present in the only fully sequenced human-origin HZ genome, representing the possibility for almost unlimited antigenic diversity. In the present study, we identified a select group of 20% of the A. phagocytophilum HZ msp2 pseudogenes that have matched preferentially to human, canine, and equine expression cassettes. Pseudogenes cluster predominantly in one spatial run limited to a single genomic island in less than 50% of the genome but phylogenetically related pseudogenes are neither necessarily located in close proximity on the genome nor share similar percent identity with expression cassettes. Pseudogenes near the expression cassette (and the origin) are more likely to be expressed than those farther away. Taken together, these findings suggest that there may be natural selection pressure to retain pseudogenes in one cluster near the putative origin of replication, even though global recombination shuffles pseudogenes around the genome, separating pseudogenes that share genetic origins as well as those with similar identities.

BackgroundIntegrons are mechanisms that facilitate horizontal gene transfer, allowing bacteria to integrate and express foreign DNA. These are important in the exchange of antibiotic resistance determinants, but can also transfer a diverse suite of genes unrelated to pathogenicity. Here, we provide a systematic analysis of the distribution and diversity of integron intI genes and integron-containing bacteria.ResultsWe found integrons in 103 different pathogenic and non-pathogenic bacteria, in six major phyla. Integrons were widely scattered, and their presence was not confined to specific clades within bacterial orders. Nearly 1/3 of the intI genes that we identified were pseudogenes, containing either an internal stop codon or a frameshift mutation that would render the protein product non-functional. Additionally, 20% of bacteria contained more than one integrase gene. dN/dS ratios revealed mutational hotspots in clades of Vibrio and Shewanella intI genes. Finally, we characterized the gene cassettes associated with integrons in Methylobacillus flagellatus KT and Dechloromonas aromatica RCB, and found a heavy metal efflux gene as well as genes involved in protein folding and stability.ConclusionOur analysis suggests that the present distribution of integrons is due to multiple losses and gene transfer events. While, in some cases, the ability to integrate and excise foreign DNA may be selectively advantageous, the gain, loss, or rearrangment of gene cassettes could also be deleterious, selecting against functional integrases. Thus, such a high fraction of pseudogenes may suggest that the selective impact of integrons on genomes is variable, oscillating between beneficial and deleterious, possibly depending on environmental conditions.

BackgroundIntegrons are mechanisms that facilitate horizontal gene transfer, allowing bacteria to integrate and express foreign DNA. These are important in the exchange of antibiotic resistance determinants, but can also transfer a diverse suite of genes unrelated to pathogenicity. Here, we provide a systematic analysis of the distribution and diversity of integron intI genes and integron-containing bacteria.ResultsWe found integrons in 103 different pathogenic and non-pathogenic bacteria, in six major phyla. Integrons were widely scattered, and their presence was not confined to specific clades within bacterial orders. Nearly 1/3 of the intI genes that we identified were pseudogenes, containing either an internal stop codon or a frameshift mutation that would render the protein product non-functional. Additionally, 20% of bacteria contained more than one integrase gene. dN/dS ratios revealed mutational hotspots in clades of Vibrio and Shewanella intI genes. Finally, we characterized the gene cassettes associated with integrons in Methylobacillus flagellatus KT and Dechloromonas aromatica RCB, and found a heavy metal efflux gene as well as genes involved in protein folding and stability.ConclusionOur analysis suggests that the present distribution of integrons is due to multiple losses and gene transfer events. While, in some cases, the ability to integrate and excise foreign DNA may be selectively advantageous, the gain, loss, or rearrangment of gene cassettes could also be deleterious, selecting against functional integrases. Thus, such a high fraction of pseudogenes may suggest that the selective impact of integrons on genomes is variable, oscillating between beneficial and deleterious, possibly depending on environmental conditions.

BackgroundGene knockouts in a model organism such as mouse provide a valuable resource for the study of basic biology and human disease. Determining which gene has been inactivated by an untargeted gene trapping event poses a challenging annotation problem because gene trap sequence tags, which represent sequence near the vector insertion site of a trapped gene, are typically short and often contain unresolved residues. To understand better the localization of these sequences on the mouse genome, we compared stand-alone versions of the alignment programs BLAT, SSAHA, and MegaBLAST. A set of 3,369 sequence tags was aligned to build 34 of the mouse genome using default parameters for each algorithm. Known genome coordinates for the cognate set of full-length genes (1,659 sequences) were used to evaluate localization results.ResultsIn general, all three programs performed well in terms of localizing sequences to a general region of the genome, with only relatively subtle errors identified for a small proportion of the sequence tags. However, large differences in performance were noted with regard to correctly identifying exon boundaries. BLAT correctly identified the vast majority of exon boundaries, while SSAHA and MegaBLAST missed the majority of exon boundaries. SSAHA consistently reported the fewest false positives and is the fastest algorithm. MegaBLAST was comparable to BLAT in speed, but was the most susceptible to localizing sequence tags incorrectly to pseudogenes.ConclusionThe differences in performance for sequence tags and full-length reference sequences were surprisingly small. Characteristic variations in localization results for each program were noted that affect the localization of sequence at exon boundaries, in particular.

BackgroundGene knockouts in a model organism such as mouse provide a valuable resource for the study of basic biology and human disease. Determining which gene has been inactivated by an untargeted gene trapping event poses a challenging annotation problem because gene trap sequence tags, which represent sequence near the vector insertion site of a trapped gene, are typically short and often contain unresolved residues. To understand better the localization of these sequences on the mouse genome, we compared stand-alone versions of the alignment programs BLAT, SSAHA, and MegaBLAST. A set of 3,369 sequence tags was aligned to build 34 of the mouse genome using default parameters for each algorithm. Known genome coordinates for the cognate set of full-length genes (1,659 sequences) were used to evaluate localization results.ResultsIn general, all three programs performed well in terms of localizing sequences to a general region of the genome, with only relatively subtle errors identified for a small proportion of the sequence tags. However, large differences in performance were noted with regard to correctly identifying exon boundaries. BLAT correctly identified the vast majority of exon boundaries, while SSAHA and MegaBLAST missed the majority of exon boundaries. SSAHA consistently reported the fewest false positives and is the fastest algorithm. MegaBLAST was comparable to BLAT in speed, but was the most susceptible to localizing sequence tags incorrectly to pseudogenes.ConclusionThe differences in performance for sequence tags and full-length reference sequences were surprisingly small. Characteristic variations in localization results for each program were noted that affect the localization of sequence at exon boundaries, in particular.

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.

Pseudogenes have been defined as non-functional sequences of genomic DNA that are originally derived from functional genes, but exhibit degenerative features such as premature stop codons and frameshifts that prevent their expression. However, there is increasing evidence that pseudogenes are often evolutionarily conserved and may have retained some functional role or acquired new ones. Pseudogenes may exhibit non-functional features as well as functional ones. We investigate, as a model case, the beta-esterase gene cluster of Drosophila melanogaster that includes the Est-6 gene and the psiEst-6 putative pseudogene. We study four samples derived from natural populations of east Africa (Zimbabwe), Europe (Spain), North America (California), and South America (Venezuela). The level of nucleotide diversity is higher in Africa than in the non-African populations. There is twice more nucleotide diversity in psiEst-6 than in Est-6. Linkage disequilibrium within the beta-esterase gene cluster is strong in non-African samples, but much lower in Africa. The population recombination rate is the same for psiEst-6 and Est-6 in Africa, but significantly different in non-African samples. Intragenic gene conversion events are detected within Est-6 and, with much higher incidence, within psiEst-6; intergenic gene conversion events are rare. The extensive intragenic gene conversion within psiEst-6 can be explained by the invasion of retrotransposons that promote a form of homology-dependent gene conversion upon excision. Tests of neutrality with recombination are significant for the beta-esterase gene cluster in the non-African populations but not in Africa. The Est-6 gene sequences exhibit a well-known allozyme dimorphic structure. The sequences of psiEst-6 are also dimorphic in North and South America, but they do not correspond at all (South America) or only imperfectly (North America) to the Est-6 allozyme dimorphism. Sequence dimorphism is less pronounced in the European and African samples. We suggest that demographic history (bottleneck and admixture of genetically differentiated populations) is the major factor shaping the nucleotide pattern in the beta-esterase gene cluster. However, there are some clear indications of positive selection shaping the distribution of nucleotide polymorphism within the cluster. Intergenic epistatic selection may play an important role in the evolution of the beta-esterase gene cluster, preserving psiEst-6 from degenerative destruction and reflecting its functional interaction with Est-6. The Est-6 gene cluster of D. melanogaster represents an example of a functionally interacting complex ('intergene') in which two components (Est-6 and psiEst-6) or more are required to perform the final function.

Pseudogenes have been defined as non-functional sequences of genomic DNA that are originally derived from functional genes, but exhibit degenerative features such as premature stop codons and frameshifts that prevent their expression. However, there is increasing evidence that pseudogenes are often evolutionarily conserved and may have retained some functional role or acquired new ones. Pseudogenes may exhibit non-functional features as well as functional ones. We investigate, as a model case, the beta-esterase gene cluster of Drosophila melanogaster that includes the Est-6 gene and the psiEst-6 putative pseudogene. We study four samples derived from natural populations of east Africa (Zimbabwe), Europe (Spain), North America (California), and South America (Venezuela). The level of nucleotide diversity is higher in Africa than in the non-African populations. There is twice more nucleotide diversity in psiEst-6 than in Est-6. Linkage disequilibrium within the beta-esterase gene cluster is strong in non-African samples, but much lower in Africa. The population recombination rate is the same for psiEst-6 and Est-6 in Africa, but significantly different in non-African samples. Intragenic gene conversion events are detected within Est-6 and, with much higher incidence, within psiEst-6; intergenic gene conversion events are rare. The extensive intragenic gene conversion within psiEst-6 can be explained by the invasion of retrotransposons that promote a form of homology-dependent gene conversion upon excision. Tests of neutrality with recombination are significant for the beta-esterase gene cluster in the non-African populations but not in Africa. The Est-6 gene sequences exhibit a well-known allozyme dimorphic structure. The sequences of psiEst-6 are also dimorphic in North and South America, but they do not correspond at all (South America) or only imperfectly (North America) to the Est-6 allozyme dimorphism. Sequence dimorphism is less pronounced in the European and African samples. We suggest that demographic history (bottleneck and admixture of genetically differentiated populations) is the major factor shaping the nucleotide pattern in the beta-esterase gene cluster. However, there are some clear indications of positive selection shaping the distribution of nucleotide polymorphism within the cluster. Intergenic epistatic selection may play an important role in the evolution of the beta-esterase gene cluster, preserving psiEst-6 from degenerative destruction and reflecting its functional interaction with Est-6. The Est-6 gene cluster of D. melanogaster represents an example of a functionally interacting complex ('intergene') in which two components (Est-6 and psiEst-6) or more are required to perform the final function.

To study genome evolution in wheat, we have sequenced and compared two large physical contigs of 285 and 142 kb covering orthologous low molecular weight (LMW) glutenin loci on chromosome 1AS of a diploid wheat species (Triticum monococcum subsp monococcum) and a tetraploid wheat species (Triticum turgidum subsp durum). Sequence conservation between the two species was restricted to small regions containing the orthologous LMW glutenin genes, whereas >90% of the compared sequences were not conserved. Dramatic sequence rearrangements occurred in the regions rich in repetitive elements. Dating of long terminal repeat retrotransposon insertions revealed different insertion events occurring during the last 5.5 million years in both species. These insertions are partially responsible for the lack of homology between the intergenic regions. In addition, the gene space was conserved only partially, because different predicted genes were identified on both contigs. Duplications and deletions of large fragments that might be attributable to illegitimate recombination also have contributed to the differentiation of this region in both species. The striking differences in the intergenic landscape between the A and A(m) genomes that diverged 1 to 3 million years ago provide evidence for a dynamic and rapid genome evolution in wheat species.

To study genome evolution in wheat, we have sequenced and compared two large physical contigs of 285 and 142 kb covering orthologous low molecular weight (LMW) glutenin loci on chromosome 1AS of a diploid wheat species (Triticum monococcum subsp monococcum) and a tetraploid wheat species (Triticum turgidum subsp durum). Sequence conservation between the two species was restricted to small regions containing the orthologous LMW glutenin genes, whereas >90% of the compared sequences were not conserved. Dramatic sequence rearrangements occurred in the regions rich in repetitive elements. Dating of long terminal repeat retrotransposon insertions revealed different insertion events occurring during the last 5.5 million years in both species. These insertions are partially responsible for the lack of homology between the intergenic regions. In addition, the gene space was conserved only partially, because different predicted genes were identified on both contigs. Duplications and deletions of large fragments that might be attributable to illegitimate recombination also have contributed to the differentiation of this region in both species. The striking differences in the intergenic landscape between the A and A(m) genomes that diverged 1 to 3 million years ago provide evidence for a dynamic and rapid genome evolution in wheat species.

Efforts to understand the genetic basis of evolutionary change have concentrated on proteins and their encoding DNA sequences. These studies have brought to light patterns and processes at the nucleotide level, yet the complex functional relationships between genetic variants and phenotypes remain poorly known. The realization that even a complete description of proteins and the effects of their activity will not suffice to understand the conditions under which they are time- and tissue-specifically expressed or repressed during development has refocused attention on cis-regulatory regions. In particular, promoter sequences are thought to hold the key for understanding the evolution of phenotypic differences between species. This is because of their complex organization into independent modules such that, unlike coding sequences in which mutations affect protein function every time the protein is expressed, mutations in cis-regulatory sequences may have minor or no pleiotropic effects. Complex information-encoding makes cis-regulatory regions poorly amenable to comparative methods designed for coding sequences. Some general conclusions are emerging as to how genetic variation is distributed across regulatory networks and the processes modulating the structure of this variation. We bring into this emerging scenario several recent findings pointing to different ways in which spliceosomal introns, pseudogenes and patterns of point mutation can be active players for the evolution of novel transcriptional profiles.

Efforts to understand the genetic basis of evolutionary change have concentrated on proteins and their encoding DNA sequences. These studies have brought to light patterns and processes at the nucleotide level, yet the complex functional relationships between genetic variants and phenotypes remain poorly known. The realization that even a complete description of proteins and the effects of their activity will not suffice to understand the conditions under which they are time- and tissue-specifically expressed or repressed during development has refocused attention on cis-regulatory regions. In particular, promoter sequences are thought to hold the key for understanding the evolution of phenotypic differences between species. This is because of their complex organization into independent modules such that, unlike coding sequences in which mutations affect protein function every time the protein is expressed, mutations in cis-regulatory sequences may have minor or no pleiotropic effects. Complex information-encoding makes cis-regulatory regions poorly amenable to comparative methods designed for coding sequences. Some general conclusions are emerging as to how genetic variation is distributed across regulatory networks and the processes modulating the structure of this variation. We bring into this emerging scenario several recent findings pointing to different ways in which spliceosomal introns, pseudogenes and patterns of point mutation can be active players for the evolution of novel transcriptional profiles.