Search

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.

The Ca2+-sensor protein calmodulin (CaM) is a major regulator of multiple cell functions. A unique and puzzling feature of human, and all so far investigated mammals, is the presence of three distinct CaM genes on different chromosomes, which code for identical proteins. How this case of apparent genetic redundancy evolved and why it could be to the advantage of the mammalian organisms is not well established. With a main focus on humans, this article aims to review existing literature addressing how the genes nonetheless differ in function. Clearly, the three CaM genes are differentially expressed in different tissues, during development, in response to different stimuli, and other factors including environmental conditions. As shown in hippocampal neurons, different mRNAs from the CAM genes may even localize differently within the same cell. Regulation of CaM gene expression is achieved by a variety of regulatory elements present in the three genes, including different promotor/insulator elements and 3′- and 5′-noncoding regions differing in length and sequence, as well as regulation by epigenetic factors and miRNAs. Here, we hypothesize that predicted differences in mRNA stability and translational efficiency due to divergent codon usage could play an additional regulatory role as the three genes differ markedly in their use of synonymous codons. CALM3, predicted to produce a relatively stable mRNA may be important where the transcription level is low or transiently absent, e.g. during spermatogenesis. In contrast, CALM2 with a predicted much shorter mRNA half-life, may provide better temporal control of CaM levels. Deciphering the underlying mechanisms responsible for all this complexity may help to understand why this unique multigenic arrangement may be an advantage for the optimal spatio-temporal expression of CaM in higher eukaryotes. Finally, we discuss the expression of the CaM genes in selected human pathologies, and how mutations in these genes are respo

The Ca2+-sensor protein calmodulin (CaM) is a major regulator of multiple cell functions. A unique and puzzling feature of human, and all so far investigated mammals, is the presence of three distinct CaM genes on different chromosomes, which code for identical proteins. How this case of apparent genetic redundancy evolved and why it could be to the advantage of the mammalian organisms is not well established. With a main focus on humans, this article aims to review existing literature addressing how the genes nonetheless differ in function. Clearly, the three CaM genes are differentially expressed in different tissues, during development, in response to different stimuli, and other factors including environmental conditions. As shown in hippocampal neurons, different mRNAs from the CAM genes may even localize differently within the same cell. Regulation of CaM gene expression is achieved by a variety of regulatory elements present in the three genes, including different promotor/insulator elements and 3′- and 5′-noncoding regions differing in length and sequence, as well as regulation by epigenetic factors and miRNAs. Here, we hypothesize that predicted differences in mRNA stability and translational efficiency due to divergent codon usage could play an additional regulatory role as the three genes differ markedly in their use of synonymous codons. CALM3, predicted to produce a relatively stable mRNA may be important where the transcription level is low or transiently absent, e.g. during spermatogenesis. In contrast, CALM2 with a predicted much shorter mRNA half-life, may provide better temporal control of CaM levels. Deciphering the underlying mechanisms responsible for all this complexity may help to understand why this unique multigenic arrangement may be an advantage for the optimal spatio-temporal expression of CaM in higher eukaryotes. Finally, we discuss the expression of the CaM genes in selected human pathologies, and how mutations in these g

Hepatocellular carcinoma (HCC), the main form of primary liver cancer, is the second leading cause of cancer-related deaths worldwide after lung cancer. Multiple aetiologies have been associated with the development of HCC, which arises in most cases in the context of a chronically inflamed liver. HCC is in fact an inflammation-driven cancer, with the TNF and IL6 families of cytokines playing key roles in maintaining a chronic inflammatory state, promoting hepatocarcinogenesis. IL6 signals mainly through the JAK1/STAT3 signal transduction pathway and is known to play key roles in liver physiology and disease. In the interest of identifying novel players and downstream effectors of the IL6/JAK1/STAT3 signalling pathway that may contribute to the signal transduction of IL6 in liver-derived cells, we have been investigating the expression of long non-coding RNAs (lncRNAs) in response to treatment with the designer cytokine Hyper-IL6. Indeed, lncRNAs have recently emerged as a key layer of biological regulation and have been shown to be differentially expressed in cancer, including HCC. Upon analysis of time series transcriptomics data, we have identified hundreds of lncRNAs to be differentially expressed in HepG2, HuH7, and Hep3B hepatoma cells upon cytokine stimulation; 26 of which are common to the three cell lines tested. qPCR validation experiments have been performed for several lncRNAs, such as for the liver-specific lncRNA linc-ELL2. By functionally characterising identified clusters of IL6-regulated coding and non-coding genes in hepatoma cells, we propose, based on a guilt-by-association hypothesis, novel functions for previously poorly characterized lncRNAs and pseudogenes such as AL391422.4 or TUBA5P. Several lncRNA genes seem to be co-regulated with a protein-coding gene localized in their vicinity. For example, Hyper-IL-6 increases the mRNA and protein levels of XBP1, a well-known regulator of the unfolded protein response. At the same time, the expression

Hepatocellular carcinoma (HCC), the main form of primary liver cancer, is the second leading cause of cancer-related deaths worldwide after lung cancer. Multiple aetiologies have been associated with the development of HCC, which arises in most cases in the context of a chronically inflamed liver. HCC is in fact an inflammation-driven cancer, with the TNF and IL6 families of cytokines playing key roles in maintaining a chronic inflammatory state, promoting hepatocarcinogenesis. IL6 signals mainly through the JAK1/STAT3 signal transduction pathway and is known to play key roles in liver physiology and disease. In the interest of identifying novel players and downstream effectors of the IL6/JAK1/STAT3 signalling pathway that may contribute to the signal transduction of IL6 in liver-derived cells, we have been investigating the expression of long non-coding RNAs (lncRNAs) in response to treatment with the designer cytokine Hyper-IL6. Indeed, lncRNAs have recently emerged as a key layer of biological regulation and have been shown to be differentially expressed in cancer, including HCC. Upon analysis of time series transcriptomics data, we have identified hundreds of lncRNAs to be differentially expressed in HepG2, HuH7, and Hep3B hepatoma cells upon cytokine stimulation; 26 of which are common to the three cell lines tested. qPCR validation experiments have been performed for several lncRNAs, such as for the liver-specific lncRNA linc-ELL2. By functionally characterising identified clusters of IL6-regulated coding and non-coding genes in hepatoma cells, we propose, based on a guilt-by-association hypothesis, novel functions for previously poorly characterized lncRNAs and pseudogenes such as AL391422.4 or TUBA5P. Several lncRNA genes seem to be co-regulated with a protein-coding gene localized in their vicinity. For example, Hyper-IL-6 increases the mRNA and protein levels of XBP1, a well-known regulator of the unfolded protein response. At the same time, the expression

The Ca2+-sensor protein calmodulin (CaM) is a major regulator of multiple cell functions. A unique and puzzling feature of human, and all so far investigated mammals, is the presence of three distinct CaM genes on different chromosomes, which code for identical proteins. How this case of apparent genetic redundancy evolved and why it could be to the advantage of the mammalian organisms is not well established. With a main focus on humans, this article aims to review existing literature addressing how the genes nonetheless differ in function. Clearly, the three CaM genes are differentially expressed in different tissues, during development, in response to different stimuli, and other factors including environmental conditions. As shown in hippocampal neurons, different mRNAs from the CAM genes may even localize differently within the same cell. Regulation of CaM gene expression is achieved by a variety of regulatory elements present in the three genes, including different promotor/insulator elements and 3′- and 5′-noncoding regions differing in length and sequence, as well as regulation by epigenetic factors and miRNAs. Here, we hypothesize that predicted differences in mRNA stability and translational efficiency due to divergent codon usage could play an additional regulatory role as the three genes differ markedly in their use of synonymous codons. CALM3, predicted to produce a relatively stable mRNA may be important where the transcription level is low or transiently absent, e.g. during spermatogenesis. In contrast, CALM2 with a predicted much shorter mRNA half-life, may provide better temporal control of CaM levels. Deciphering the underlying mechanisms responsible for all this complexity may help to understand why this unique multigenic arrangement may be an advantage for the optimal spatio-temporal expression of CaM in higher eukaryotes. Finally, we discuss the expression of the CaM genes in selected human pathologies, and how mutations in these g

The Ca2+-sensor protein calmodulin (CaM) is a major regulator of multiple cell functions. A unique and puzzling feature of human, and all so far investigated mammals, is the presence of three distinct CaM genes on different chromosomes, which code for identical proteins. How this case of apparent genetic redundancy evolved and why it could be to the advantage of the mammalian organisms is not well established. With a main focus on humans, this article aims to review existing literature addressing how the genes nonetheless differ in function. Clearly, the three CaM genes are differentially expressed in different tissues, during development, in response to different stimuli, and other factors including environmental conditions. As shown in hippocampal neurons, different mRNAs from the CAM genes may even localize differently within the same cell. Regulation of CaM gene expression is achieved by a variety of regulatory elements present in the three genes, including different promotor/insulator elements and 3′- and 5′-noncoding regions differing in length and sequence, as well as regulation by epigenetic factors and miRNAs. Here, we hypothesize that predicted differences in mRNA stability and translational efficiency due to divergent codon usage could play an additional regulatory role as the three genes differ markedly in their use of synonymous codons. CALM3, predicted to produce a relatively stable mRNA may be important where the transcription level is low or transiently absent, e.g. during spermatogenesis. In contrast, CALM2 with a predicted much shorter mRNA half-life, may provide better temporal control of CaM levels. Deciphering the underlying mechanisms responsible for all this complexity may help to understand why this unique multigenic arrangement may be an advantage for the optimal spatio-temporal expression of CaM in higher eukaryotes. Finally, we discuss the expression of the CaM genes in selected human pathologies, and how mutations in these genes are respo

G-protein coupled receptors (GPCRs) play a key role in a broad range of biological processes by binding to a wide variety of signaling molecules, which have resulted in 34% of all FDA-approved drugs which target GPCRs. The human genome encodes for approximately 800 GPCR members of which about 140 non-olfactory receptors remain orphans with an unknown function and endogenous ligand. Despite prolonged efforts to deorphanize the unresolved receptors, they remain orphans until this day. By studying scientific publications, this thesis has clarified the challenges with the deorphanization of GPCRs to explain why there are still so many orphan GPCRs when they have confirmed involvement in so many human disorders.

G-protein coupled receptors (GPCRs) play a key role in a broad range of biological processes by binding to a wide variety of signaling molecules, which have resulted in 34% of all FDA-approved drugs which target GPCRs. The human genome encodes for approximately 800 GPCR members of which about 140 non-olfactory receptors remain orphans with an unknown function and endogenous ligand. Despite prolonged efforts to deorphanize the unresolved receptors, they remain orphans until this day. By studying scientific publications, this thesis has clarified the challenges with the deorphanization of GPCRs to explain why there are still so many orphan GPCRs when they have confirmed involvement in so many human disorders.

Missoula, MT

Mobile genetic elements often encode integrases which catalyze the site-specific insertion of their genetic information into the host genome and the reverse reaction of excision. Hyperthermophilic archaea harbor integrases belonging to the SSV-family which carry the MGE recombination site within their open reading frame. Upon integration into the host genome, SSV integrases disrupt their own gene into two inactive pseudogenes and are termed suicidal for this reason. The evolutionary maintenance of suicidal integrases, concurring with the high prevalence and multiples recruitments of these recombinases by archaeal MGEs, is highly paradoxical. To elucidate this phenomenon, we analyzed the wide phylogenomic distribution of a prominent class of suicidal integrases which revealed a highly variable integration site specificity. Our results highlighted the remarkable hybrid nature of these enzymes encoded from the assembly of inactive pseudogenes of different origins. The characterization of the biological properties of one of these integrases, IntpT26-2 showed that this enzyme was active over a wide range of temperatures up to 99 °C and displayed a less stringent site specificity requirement than comparable integrases. These observations concurred in explaining the pervasiveness of these suicidal integrases in the most hyperthermophilic organisms. The biochemical and phylogenomic data presented here revealed a target site switching system operating on highly thermostable integrases and suggested a new model for split gene reconstitution. By generating fast-evolving pseudogenes at high frequency, suicidal integrases constitute a powerful model to approach the molecular mechanisms involved in the generation of active genes variants by the recombination of proto-genes.

Mobile genetic elements often encode integrases which catalyze the site-specific insertion of their genetic information into the host genome and the reverse reaction of excision. Hyperthermophilic archaea harbor integrases belonging to the SSV-family which carry the MGE recombination site within their open reading frame. Upon integration into the host genome, SSV integrases disrupt their own gene into two inactive pseudogenes and are termed suicidal for this reason. The evolutionary maintenance of suicidal integrases, concurring with the high prevalence and multiples recruitments of these recombinases by archaeal MGEs, is highly paradoxical. To elucidate this phenomenon, we analyzed the wide phylogenomic distribution of a prominent class of suicidal integrases which revealed a highly variable integration site specificity. Our results highlighted the remarkable hybrid nature of these enzymes encoded from the assembly of inactive pseudogenes of different origins. The characterization of the biological properties of one of these integrases, IntpT26-2 showed that this enzyme was active over a wide range of temperatures up to 99 °C and displayed a less stringent site specificity requirement than comparable integrases. These observations concurred in explaining the pervasiveness of these suicidal integrases in the most hyperthermophilic organisms. The biochemical and phylogenomic data presented here revealed a target site switching system operating on highly thermostable integrases and suggested a new model for split gene reconstitution. By generating fast-evolving pseudogenes at high frequency, suicidal integrases constitute a powerful model to approach the molecular mechanisms involved in the generation of active genes variants by the recombination of proto-genes.

Mobile genetic elements often encode integrases which catalyze the site-specific insertion of their genetic information into the host genome and the reverse reaction of excision. Hyperthermophilic archaea harbor integrases belonging to the SSV-family which carry the MGE recombination site within their open reading frame. Upon integration into the host genome, SSV integrases disrupt their own gene into two inactive pseudogenes and are termed suicidal for this reason. The evolutionary maintenance of suicidal integrases, concurring with the high prevalence and multiples recruitments of these recombinases by archaeal MGEs, is highly paradoxical. To elucidate this phenomenon, we analyzed the wide phylogenomic distribution of a prominent class of suicidal integrases which revealed a highly variable integration site specificity. Our results highlighted the remarkable hybrid nature of these enzymes encoded from the assembly of inactive pseudogenes of different origins. The characterization of the biological properties of one of these integrases, IntpT26-2 showed that this enzyme was active over a wide range of temperatures up to 99 °C and displayed a less stringent site specificity requirement than comparable integrases. These observations concurred in explaining the pervasiveness of these suicidal integrases in the most hyperthermophilic organisms. The biochemical and phylogenomic data presented here revealed a target site switching system operating on highly thermostable integrases and suggested a new model for split gene reconstitution. By generating fast-evolving pseudogenes at high frequency, suicidal integrases constitute a powerful model to approach the molecular mechanisms involved in the generation of active genes variants by the recombination of proto-genes.

Mobile genetic elements often encode integrases which catalyze the site-specific insertion of their genetic information into the host genome and the reverse reaction of excision. Hyperthermophilic archaea harbor integrases belonging to the SSV-family which carry the MGE recombination site within their open reading frame. Upon integration into the host genome, SSV integrases disrupt their own gene into two inactive pseudogenes and are termed suicidal for this reason. The evolutionary maintenance of suicidal integrases, concurring with the high prevalence and multiples recruitments of these recombinases by archaeal MGEs, is highly paradoxical. To elucidate this phenomenon, we analyzed the wide phylogenomic distribution of a prominent class of suicidal integrases which revealed a highly variable integration site specificity. Our results highlighted the remarkable hybrid nature of these enzymes encoded from the assembly of inactive pseudogenes of different origins. The characterization of the biological properties of one of these integrases, IntpT26-2 showed that this enzyme was active over a wide range of temperatures up to 99 °C and displayed a less stringent site specificity requirement than comparable integrases. These observations concurred in explaining the pervasiveness of these suicidal integrases in the most hyperthermophilic organisms. The biochemical and phylogenomic data presented here revealed a target site switching system operating on highly thermostable integrases and suggested a new model for split gene reconstitution. By generating fast-evolving pseudogenes at high frequency, suicidal integrases constitute a powerful model to approach the molecular mechanisms involved in the generation of active genes variants by the recombination of proto-genes.

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.

[ntroduction] Gnetophytes, comprising the genera Ephedra, Gnetum and Welwitschia, are an understudied, enigmatic lineage of gymnosperms with a controversial phylogenetic relationship to other seed plants. Here we examined the organization of ribosomal DNA (rDNA) across representative species. [Methods] We applied high-throughput sequencing approaches to isolate and reconstruct rDNA units and to determine their intragenomic homogeneity. In addition, fluorescent in situ hybridization and Southern blot hybridization techniques were used to reveal the chromosome and genomic organization of rDNA. [Key results] The 5S and 35S rRNA genes were separate (S-type) in Gnetum montanum, Gnetum gnemon and Welwitschia mirabilis and linked (L-type) in Ephedra altissima. There was considerable variability in 5S rDNA abundance, ranging from as few as ~4000 (W. mirabilis) to >100 000 (G. montanum) copies. A similar large variation was also observed in 5S rDNA locus numbers (two to 16 sites per diploid cell). 5S rRNA pseudogenes were interspersed between functional genes forming a single unit in E. altissima and G. montanum. Their copy number was comparable or even higher than that of functional 5S rRNA genes. In E. altissima internal transcribed spacers of 35S rDNA were long and intrinsically repetitive while in G. montanum and W. mirabilis they were short without the subrepeats. [Conclusions] Gnetophytes are distinct from other gymnosperms and angiosperms as they display surprisingly large variability in rDNA organization and rDNA copy and locus numbers between genera, with no relationship between copy numbers and genome sizes apparent. Concerted evolution of 5S rDNA units seems to have led to the amplification of 5S pseudogenes in both G. montanum and E. altissima. Evolutionary patterns of rDNA show both gymnosperm and angiosperm features underlining the diversity of the group.

Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.

The wide implementation of next-generation sequencing (NGS) technologies has revolutionized the field of medical genetics. However, the short read lengths of currently used sequencing approaches pose a limitation for the identification of structural variants, sequencing repetitive regions, phasing of alleles and distinguishing highly homologous genomic regions. These limitations may significantly contribute to the diagnostic gap in patients with genetic disorders who have undergone standard NGS, like whole exome or even genome sequencing. Now, the emerging long-read sequencing (LRS) technologies may offer improvements in the characterization of genetic variation and regions that are difficult to assess with the prevailing NGS approaches. LRS has so far mainly been used to investigate genetic disorders with previously known or strongly suspected disease loci. While these targeted approaches already show the potential of LRS, it remains to be seen whether LRS technologies can soon enable true whole genome sequencing routinely. Ultimately, this could allow the de novo assembly of individual whole genomes used as a generic test for genetic disorders. In this article, we summarize the current LRS-based research on human genetic disorders and discuss the potential of these technologies to facilitate the next major advancements in medical genetics.

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.

The genomes of living organisms are composed of a multitude of functional units, which interact with each other and their environment in a highly regulated fashion, to facilitate the expression of an enduring (and evolving) phenotype. Several approaches have emerged in the effort to identify these functional units and explore their activities. In this thesis, I have taken a population genetics approach; evaluating how the distribution of genetic variation in the human genome has been shaped through the actions of natural selection on functional genomic elements. In the first paper, I interrogate a catalogue of elements derived from biochemical signatures for signals of selection; finding significant signals of purifying selection on regulatory elements, independent of linked-purifying selection. In the second paper, I explore the pseudogene class of genomic elements, and find that a large proportion of a particular subclass, transcribed duplicated pseudogenes, has experienced significant amounts of positive selection. In the third paper, I focus on protein coding genes and variants that disrupt their open reading frames. Specifically, I examine the distribution of loss-of-function variants in the Khoe-San population; gauging their functional significance and exploring the biological roles of affected genes. In the final paper, instead of using population genetics to uncover and explore genomic elements, I use a major genomic element – the Y chromosome – as an effective tool to study the evolutionary history of a human population.

The genomes of living organisms are composed of a multitude of functional units, which interact with each other and their environment in a highly regulated fashion, to facilitate the expression of an enduring (and evolving) phenotype. Several approaches have emerged in the effort to identify these functional units and explore their activities. In this thesis, I have taken a population genetics approach; evaluating how the distribution of genetic variation in the human genome has been shaped through the actions of natural selection on functional genomic elements. In the first paper, I interrogate a catalogue of elements derived from biochemical signatures for signals of selection; finding significant signals of purifying selection on regulatory elements, independent of linked-purifying selection. In the second paper, I explore the pseudogene class of genomic elements, and find that a large proportion of a particular subclass, transcribed duplicated pseudogenes, has experienced significant amounts of positive selection. In the third paper, I focus on protein coding genes and variants that disrupt their open reading frames. Specifically, I examine the distribution of loss-of-function variants in the Khoe-San population; gauging their functional significance and exploring the biological roles of affected genes. In the final paper, instead of using population genetics to uncover and explore genomic elements, I use a major genomic element – the Y chromosome – as an effective tool to study the evolutionary history of a human population.

The genomes of living organisms are composed of a multitude of functional units, which interact with each other and their environment in a highly regulated fashion, to facilitate the expression of an enduring (and evolving) phenotype. Several approaches have emerged in the effort to identify these functional units and explore their activities. In this thesis, I have taken a population genetics approach; evaluating how the distribution of genetic variation in the human genome has been shaped through the actions of natural selection on functional genomic elements. In the first paper, I interrogate a catalogue of elements derived from biochemical signatures for signals of selection; finding significant signals of purifying selection on regulatory elements, independent of linked-purifying selection. In the second paper, I explore the pseudogene class of genomic elements, and find that a large proportion of a particular subclass, transcribed duplicated pseudogenes, has experienced significant amounts of positive selection. In the third paper, I focus on protein coding genes and variants that disrupt their open reading frames. Specifically, I examine the distribution of loss-of-function variants in the Khoe-San population; gauging their functional significance and exploring the biological roles of affected genes. In the final paper, instead of using population genetics to uncover and explore genomic elements, I use a major genomic element – the Y chromosome – as an effective tool to study the evolutionary history of a human population.

The genomes of living organisms are composed of a multitude of functional units, which interact with each other and their environment in a highly regulated fashion, to facilitate the expression of an enduring (and evolving) phenotype. Several approaches have emerged in the effort to identify these functional units and explore their activities. In this thesis, I have taken a population genetics approach; evaluating how the distribution of genetic variation in the human genome has been shaped through the actions of natural selection on functional genomic elements. In the first paper, I interrogate a catalogue of elements derived from biochemical signatures for signals of selection; finding significant signals of purifying selection on regulatory elements, independent of linked-purifying selection. In the second paper, I explore the pseudogene class of genomic elements, and find that a large proportion of a particular subclass, transcribed duplicated pseudogenes, has experienced significant amounts of positive selection. In the third paper, I focus on protein coding genes and variants that disrupt their open reading frames. Specifically, I examine the distribution of loss-of-function variants in the Khoe-San population; gauging their functional significance and exploring the biological roles of affected genes. In the final paper, instead of using population genetics to uncover and explore genomic elements, I use a major genomic element – the Y chromosome – as an effective tool to study the evolutionary history of a human population.

The genomes of living organisms are composed of a multitude of functional units, which interact with each other and their environment in a highly regulated fashion, to facilitate the expression of an enduring (and evolving) phenotype. Several approaches have emerged in the effort to identify these functional units and explore their activities. In this thesis, I have taken a population genetics approach; evaluating how the distribution of genetic variation in the human genome has been shaped through the actions of natural selection on functional genomic elements. In the first paper, I interrogate a catalogue of elements derived from biochemical signatures for signals of selection; finding significant signals of purifying selection on regulatory elements, independent of linked-purifying selection. In the second paper, I explore the pseudogene class of genomic elements, and find that a large proportion of a particular subclass, transcribed duplicated pseudogenes, has experienced significant amounts of positive selection. In the third paper, I focus on protein coding genes and variants that disrupt their open reading frames. Specifically, I examine the distribution of loss-of-function variants in the Khoe-San population; gauging their functional significance and exploring the biological roles of affected genes. In the final paper, instead of using population genetics to uncover and explore genomic elements, I use a major genomic element – the Y chromosome – as an effective tool to study the evolutionary history of a human population.

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. I

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.

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.

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.

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.

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.

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.

Previous research on the evolution of vision in mammals has typically focused on the duplication and loss of genes encoding opsins, the proteins that make up one moiety of the retinal photopigment complex. These gains and losses are typically attributed to invading new photic niches, with gains associated with bright light and losses with dim light niches. Additionally, the adaptive significance of opsin spectral tuning (i.e., preferential absorption wavelength) has largely remained a mystery for terrestrial mammals. For example, for almost 25 years researchers have been aware that some mammals possess ultraviolet-sensitive opsins whereas others possess violet-sensitive homologs. While scientists have proposed various explanations for this variation, there have been no formal tests of these hypotheses. In this dissertation, I expand upon this previous work in two ways: 1) exploring the loss of non-opsin visual genes in mammals occupying dim-light habitats, and 2) testing hypotheses of the evolution of an opsin’s spectral tuning in terrestrial environments. I used genome- and PCR-generated DNA sequence data in a comparative phylogenetic framework, and concluded the following: 1) subterranean mammals have variable degrees of visual gene decay, especially to the bright light visual system, as a result of their underground lifestyle, 2) armadillos, sloths and anteaters likely descended from a subterranean ancestor, as evidenced by the loss of genes encoding their bright light visual system, and 3) the ultraviolet/violet opsin in mammals has evolved in response to an increase in eye length and light exposure during mammalian history.

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 maintain

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.

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.

Background: Insertion sequences (ISs) are approximately 1 kbp long jumping genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (cut-and-paste) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. Results: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (local hopping) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. Conclusions: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wi, AuthorCount:4

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.

Background: Insertion sequences (ISs) are approximately 1 kbp long jumping genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (cut-and-paste) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. Results: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (local hopping) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. Conclusions: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wi, AuthorCount:4

Background: Insertion sequences (ISs) are approximately 1 kbp long "jumping" genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I ("cut-and-paste") Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. Results: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location ("local hopping") and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. Conclusions: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than gen

Background: Insertion sequences (ISs) are approximately 1 kbp long jumping genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (cut-and-paste) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. Results: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (local hopping) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. Conclusions: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wi, AuthorCount:4

Background: Insertion sequences (ISs) are approximately 1 kbp long jumping genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (cut-and-paste) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. Results: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (local hopping) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. Conclusions: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wi, AuthorCount:4

Background: Insertion sequences (ISs) are approximately 1 kbp long jumping genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (cut-and-paste) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. Results: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (local hopping) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. Conclusions: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wi, AuthorCount:4

We describe PGTools, an open source software suite for analysis and visualization of proteogenomic data. PGTools comprises applications, libraries, customized databases, and visualization tools for analysis of mass-spectrometry data using combined proteomic and genomic backgrounds. A single command is sufficient to search databases, calculate false discovery rates, group and annotate proteins, generate peptide databases from RNA-Seq transcripts, identify altered proteins associated with cancer, and visualize genome scale peptide data sets using sophisticated visualization tools. We experimentally confirm a subset of proteogenomic peptides in human PANC-1 cells and demonstrate the utility of PGTools using a colorectal cancer data set that led to the identification of 203 novel protein coding regions missed by conventional proteomic approaches. PGTools should be equally useful for individual proteogenomic investigations as well as international initiatives such as chromosome-centric Human Proteome Project (C-HPP). PGTools is available at http://qcmg.org/bioinformatics/PGTools.

Background: Insertion sequences (ISs) are approximately 1 kbp long jumping genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (cut-and-paste) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium 'Nostoc azollae' 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of 'Nostoc azollae' 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria. Results: Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (local hopping) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of 'Nostoc azollae' 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style. Conclusions: The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wi, AuthorCount:4