Compensated Hypogonadism Identified in Males with Cluster Headache: A Prospective Case-Controlled Study.

OBJECTIVE: Androgens have been hypothesized to be involved in the pathophysiology of cluster headache due to the male predominance, but whether androgens are altered in patients with cluster headache remains unclear. METHODS: We performed a prospective, case-controlled study in adult males with cluster headache. Sera were measured for hormones including testosterone, luteinizing hormone (LH), and sex hormone-binding globulin in 60 participants with episodic cluster headache (during a bout and in remission), 60 participants with chronic cluster headache, and 60 age- and sex-matched healthy controls. Free testosterone (fT) was calculated according to the Vermeulen equation. Shared genetic risk variants were assessed between cluster headache and testosterone concentrations. RESULTS: The mean fT/LH ratio was reduced by 35% (95% confidence interval [CI]: 21%-47%, p < 0.0001) in patients with chronic cluster headache and by 24% (95% CI: 9%-37%, p = 0.004) in patients with episodic cluster headache compared to controls after adjusting for age, sleep duration, and use of acute medication. Androgen concentrations did not differ between bouts and remissions. Furthermore, a shared genetic risk allele, rs112572874 (located in the intron of the microtubule associated protein tau (MAPT) gene on chromosome 17), between fT and cluster headache was identified. INTERPRETATION: Our results demonstrate that the male endocrine system is altered in patients with cluster headache to a state of compensated hypogonadism, and this is not an epiphenomenon associated with sleep or the use of acute medication. Together with the identified shared genetic risk allele, this may suggest a pathophysiological link between cluster headache and fT. ANN NEUROL 2024;95:1149-1161.

The PACAP pathway is independent of CGRP in mouse models of migraine: possible new drug target?

Calcitonin gene-related peptide (CGRP)-antagonizing drugs represent a major advance in migraine treatment. However, up to 50% of patients do not benefit from monoclonal antibodies against CGRP or its receptor. Here, we test the hypothesis that a closely related peptide, pituitary adenylate cyclase-activating peptide (PACAP-38), works independently of CGRP and thus might represent a new, alternative drug target. To understand differences in CGRP- and PACAP-mediated migraine pain, we used mouse models of provoked migraine-like pain based on multiple stimulations and subsequent measurement of tactile sensitivity response with von Frey filaments. Genetically modified mice lacking either functional CGRP receptors (Ramp1 knockout) or TRPA1 channels (Trpa1 knockout) were used together with CGRP-targeting antibodies and chemical inhibitors in wild-type mice (ntotal = 299). Ex vivo myograph studies were used to measure dilatory responses to CGRP and PACAP-38 in mouse carotid arteries. PACAP-38 provoked significant hypersensitivity and dilated the carotid arteries independently of CGRP. In contrast, glyceryl trinitrate-induced hypersensitivity is dependent on CGRP. Contrary to previous results with the migraine-inducing substances glyceryl trinitrate, cilostazol and levcromakalim, PACAP-38-induced hypersensitivity worked only partially through inhibition of ATP-sensitive potassium channels. Using multiple migraine-relevant models, these findings establish the PACAP-38 pathway as distinct from other migraine provoking pathways such as CGRP and glyceryl trinitrate. PACAP antagonism may therefore be a novel therapeutic target of particular interest in patients unresponsive to CGRP-antagonizing drugs.

Smooth muscle ATP-sensitive potassium channels mediate migraine-relevant hypersensitivity in mouse models.

BACKGROUND: Opening of KATP channels by systemic levcromakalim treatment triggers attacks in migraine patients and hypersensitivity to von Frey stimulation in a mouse model. Blocking of these channels is effective in several preclinical migraine models. It is unknown in what tissue and cell type KATP-induced migraine attacks are initiated and which KATP channel subtype is targeted. METHODS: In mouse models, we administered levcromakalim intracerebroventricularly, intraperitoneally and intraplantarily and compared the nociceptive responses by von Frey and hotplate tests. Mice with a conditional loss-of-function mutation in the smooth muscle KATP channel subunit Kir6.1 were given levcromakalim and GTN and examined with von Frey filaments. Arteries were tested for their ability to dilate ex vivo. mRNA expression, western blotting and immunohistochemical stainings were made to identify relevant target tissue for migraine induced by KATP channel opening. RESULTS: Systemic administration of levcromakalim induced hypersensitivity but central and local administration provided antinociception respectively no effect. The Kir6.1 smooth muscle knockout mouse was protected from both GTN and levcromakalim induced hypersensitivity, and their arteries had impaired dilatory response to the latter. mRNA and protein expression studies showed that trigeminal ganglia did not have significant KATP channel expression of any subtype, whereas brain arteries and dura mater primarily expressed the Kir6.1 + SUR2B subtype. CONCLUSION: Hypersensitivity provoked by GTN and levcromakalim in mice is dependent on functional smooth muscle KATP channels of extracerebral origin. These results suggest a vascular contribution to hypersensitivity induced by migraine triggers.

Microbial genome analysis: the COG approach.

For the past 20 years, the Clusters of Orthologous Genes (COG) database had been a popular tool for microbial genome annotation and comparative genomics. Initially created for the purpose of evolutionary classification of protein families, the COG have been used, apart from straightforward functional annotation of sequenced genomes, for such tasks as (i) unification of genome annotation in groups of related organisms; (ii) identification of missing and/or undetected genes in complete microbial genomes; (iii) analysis of genomic neighborhoods, in many cases allowing prediction of novel functional systems; (iv) analysis of metabolic pathways and prediction of alternative forms of enzymes; (v) comparison of organisms by COG functional categories; and (vi) prioritization of targets for structural and functional characterization. Here we review the principles of the COG approach and discuss its key advantages and drawbacks in microbial genome analysis.

No additive effect of combining sumatriptan and olcegepant in the GTN mouse model of migraine.

INTRODUCTION: Despite recent advances in migraine treatment there is a need for therapies with higher clinical efficacy and/or fewer side effects. Triptans (5-HT1B/1D/1F agonists) are essential in the present treatment regime and gepants (CGRP-receptor antagonists) are recognized as effective in acute migraine treatment. Triptans and gepants have different mechanisms of action and here we tested the hypothesis that a combination of these drugs (sumatriptan and olcegepant) would result in an additive effect. METHODS: Using the validated glyceryl trinitrate mouse model of migraine, we initially tested dose-response relationships of sumatriptan (0.1, 0.3, and 0.6 mg/kg IP) and olcegepant (0.25, 0.50, and 1.0 mg/kg IP) to find suitable high and low doses. Subsequently, we performed a combination study of the two drugs with a low and a high dose. All experiments were vehicle (placebo) controlled and blinded. RESULTS: Sumatriptan significantly reduced glyceryl trinitrate-induced allodynia (F(4,54) = 13.51, p < 0.0001) at all doses. Olcegepant also reduced glyceryl trinitrate-induced allodynia (F(4,53) = 16.11, p < 0.0001) with the two higher doses being significantly effective. Combining 0.50 mg/kg olcegepant with 0.1 or 0.6 mg/kg sumatriptan did not have any improved effect compared to either drug alone (p > 0.50 on all days) in our mouse model. CONCLUSION: Combining olcegepant and sumatriptan did not have an additive effect compared to single-drug treatment in this study. Triptan-gepant combinations will therefore most likely not improve migraine treatment. Nevertheless, further studies are necessary, and combinations should also be examined in patients with migraine.

CGRP-dependent signalling pathways involved in mouse models of GTN- cilostazol- and levcromakalim-induced migraine.

BACKGROUND: Knowledge of exact signalling events during migraine attacks is lacking. Various substances are known to trigger migraine attacks in patients and calcitonin gene-related peptide antagonising drugs are effective against migraine pain. Here, we investigated the signalling pathways involved in three different mouse models of provoked migraine and relate them to calcitonin gene-related peptide and other migraine-relevant targets. METHODS: In vivo mouse models of glyceryl trinitrate-, cilostazol- and levcromakalim-induced migraine were applied utilising tactile sensitivity to von Frey filaments as measuring readout. Signalling pathways involved in the three models were dissected by use of specific knockout mice and chemical inhibitors. In vivo results were supported by ex vivo wire myograph experiments measuring arterial dilatory responses and ex vivo calcitonin gene-related peptide release from trigeminal ganglion and trigeminal nucleus caudalis from mice. RESULTS: Glyceryl trinitrate-induced hypersensitivity was dependent on both prostaglandins and transient receptor potential cation channel, subfamily A, member 1, whereas cilostazol- and levcromakalim-induced hypersensitivity were independent of both. All three migraine triggers activated calcitonin gene-related peptide signalling, as both receptor antagonism and antibody neutralisation of calcitonin gene-related peptide were effective inhibitors of hypersensitivity in all three models. Stimulation of trigeminal ganglia and brain stem tissue samples with cilostazol and levcromakalim did not result in release of calcitonin gene-related peptide, and vasodilation following levcromakalim stimulation was independent of CGRP receptor antagonism. CONCLUSION: The mouse models of glyceryl trinitrate-, cilostazol- and levcromakalim- induced migraine all involve calcitonin gene-related peptide signalling in a complex interplay between different cell/tissue types. These models are useful in the study of migraine mechanisms.

Functional gene networks reveal distinct mechanisms segregating in migraine families.

Migraine is the most common neurological disorder worldwide and it has been shown to have complex polygenic origins with a heritability of estimated 40-70%. Both common and rare genetic variants are believed to underlie the pathophysiology of the prevalent types of migraine, migraine with typical aura and migraine without aura. However, only common variants have been identified so far. Here we identify for the first time a gene module with rare mutations through a systems genetics approach integrating RNA sequencing data from brain and vascular tissues likely to be involved in migraine pathology in combination with whole genome sequencing of 117 migraine families. We found a gene module in the visual cortex, based on single nuclei RNA sequencing data, that had increased rare mutations in the migraine families and replicated this in a second independent cohort of 1930 patients. This module was mainly expressed by interneurons, pyramidal CA1, and pyramidal SS cells, and pathway analysis showed association with hormonal signalling (thyrotropin-releasing hormone receptor and oxytocin receptor signalling pathways), Alzheimer's disease pathway, serotonin receptor pathway and general heterotrimeric G-protein signalling pathways. Our results demonstrate that rare functional gene variants are strongly implicated in the pathophysiology of migraine. Furthermore, we anticipate that the results can be used to explain the critical mechanisms behind migraine and potentially improving the treatment regime for migraine patients.

No central action of CGRP antagonising drugs in the GTN mouse model of migraine.

INTRODUCTION: Clinically, calcitonin gene-related peptide antagonising drugs are recognized as effective in migraine treatment, but their site of action is debated. Only a small fraction of these compounds pass the blood-brain barrier and accesses the central nervous system. Regardless, it has been argued that the central nervous system is the site of action. Here, we test this hypothesis by bypassing the blood-brain barrier through intracerebroventricular injection of calcitonin gene-related peptide antagonising drugs. METHODS: We used the glyceryl trinitrate (GTN) mouse model, which is well validated by its response to specific migraine drugs. The calcitonin gene-related peptide receptor antagonist olcegepant and the calcitonin gene-related peptide monoclonal antibody ALD405 were administered either intraperitoneally or intracerebroventricularly. The outcome measure was cutaneous mechanical allodynia. RESULTS: Mice given olcegepant intraperitoneally + GTN on day 1 had a mean 50% withdrawal threshold of 1.2 g in contrast to mice receiving placebo + GTN, which had a threshold of 0.3 g (p < 0.001). Similarly, in the ALD405 + GTN group, mice had thresholds of 1.2 g versus 0.2 g in the placebo + GTN group (p < 0.001). However, both drugs were ineffective when delivered intracerebroventricularly, as control and active groups had identical mechanical sensitivity thresholds, 0.2 g versus 0.1 g and 0.1 g versus 0.1 g for olcegepant and ALD405, respectively (p > 0.99 in both cases). DISCUSSION: The site of action of olcegepant and of the monoclonal antibody ALD405 is outside the blood-brain barrier in this mouse model of migraine. It is likely that these results can be generalised to all gepants and all antibodies and that the results are relevant for human migraine.

ATP sensitive potassium (KATP) channel inhibition: A promising new drug target for migraine.

BACKGROUND: Recently, the adenosine triphosphate (ATP) sensitive potassium channel opener levcromakalim was shown to induce migraine attacks with a far higher incidence than any previous provoking agent such as calcitonin gene-related peptide. Here, we show efficacy of ATP sensitive potassium channel inhibitors in two validated rodent models of migraine. METHODS: In female spontaneous trigeminal allodynic rats, the sensitivity of the frontal region of the head was tested by an electronic von Frey filament device. In mice, cutaneous hypersensitivity was induced by repeated glyceryl trinitrate or levcromakalim injections over nine days, as measured with von Frey filaments in the hindpaw. Release of calcitonin gene-related peptide from dura mater and trigeminal ganglion was studied ex vivo. RESULTS: The ATP sensitive potassium channel inhibitor glibenclamide attenuated the spontaneous cephalic hypersensitivity in spontaneous trigeminal allodynic rats and glyceryl trinitrate-induced hypersensitivity of the hindpaw in mice. It also inhibited CGRP release from dura mater and the trigeminal ganglion isolated from spontaneous trigeminal allodynic rats. The hypersensitivity was also diminished by the structurally different ATP sensitive potassium channel inhibitor gliquidone. Mice injected with the ATP sensitive potassium channel opener levcromakalim developed a progressive hypersensitivity that was completely blocked by glibenclamide, confirming target engagement. CONCLUSION: The results suggest that ATP sensitive potassium channel inhibitors could be novel and highly effective drugs in the treatment of migraine.

Prokaryotic Virus Orthologous Groups (pVOGs): a resource for comparative genomics and protein family annotation.

Viruses are the most abundant and diverse biological entities on earth, and while most of this diversity remains completely unexplored, advances in genome sequencing have provided unprecedented glimpses into the virosphere. The Prokaryotic Virus Orthologous Groups (pVOGs, formerly called Phage Orthologous Groups, POGs) resource has aided in this task over the past decade by using automated methods to keep pace with the rapid increase in genomic data. The uses of pVOGs include functional annotation of viral proteins, identification of genes and viruses in uncharacterized DNA samples, phylogenetic analysis, large-scale comparative genomics projects, and more. The pVOGs database represents a comprehensive set of orthologous gene families shared across multiple complete genomes of viruses that infect bacterial or archaeal hosts (viruses of eukaryotes will be added at a future date). The pVOGs are constructed within the Clusters of Orthologous Groups (COGs) framework that is widely used for orthology identification in prokaryotes. Since the previous release of the POGs, the size has tripled to nearly 3000 genomes and 300 000 proteins, and the number of conserved orthologous groups doubled to 9518. User-friendly webpages are available, including multiple sequence alignments and HMM profiles for each VOG. These changes provide major improvements to the pVOGs database, at a time of rapid advances in virus genomics. The pVOGs database is hosted jointly at the University of Iowa at http://dmk-brain.ecn.uiowa.edu/pVOGs and the NCBI at ftp://ftp.ncbi.nlm.nih.gov/pub/kristensen/pVOGs/home.html.

ATGC database and ATGC-COGs: an updated resource for micro- and macro-evolutionary studies of prokaryotic genomes and protein family annotation.

The Alignable Tight Genomic Clusters (ATGCs) database is a collection of closely related bacterial and archaeal genomes that provides several tools to aid research into evolutionary processes in the microbial world. Each ATGC is a taxonomy-independent cluster of 2 or more completely sequenced genomes that meet the objective criteria of a high degree of local gene order (synteny) and a small number of synonymous substitutions in the protein-coding genes. As such, each ATGC is suited for analysis of microevolutionary variations within a cohesive group of organisms (e.g. species), whereas the entire collection of ATGCs is useful for macroevolutionary studies. The ATGC database includes many forms of pre-computed data, in particular ATGC-COGs (Clusters of Orthologous Genes), multiple sequence alignments, a set of 'index' orthologs representing the most well-conserved members of each ATGC-COG, the phylogenetic tree of the organisms within each ATGC, etc. Although the ATGC database contains several million proteins from thousands of genomes organized into hundreds of clusters (roughly a 4-fold increase since the last version of the ATGC database), it is now built with completely automated methods and will be regularly updated following new releases of the NCBI RefSeq database. The ATGC database is hosted jointly at the University of Iowa at dmk-brain.ecn.uiowa.edu/ATGC/ and the NCBI at ftp.ncbi.nlm.nih.gov/pub/kristensen/ATGC/atgc_home.html.

Evolutionary switches between two serine codon sets are driven by selection.

Serine is the only amino acid that is encoded by two disjoint codon sets so that a tandem substitution of two nucleotides is required to switch between the two sets. Previously published evidence suggests that, for the most evolutionarily conserved serines, the codon set switch occurs by simultaneous substitution of two nucleotides. Here we report a genome-wide reconstruction of the evolution of serine codons in triplets of closely related species from diverse prokaryotes and eukaryotes. The results indicate that the great majority of codon set switches proceed by two consecutive nucleotide substitutions, via a threonine or cysteine intermediate, and are driven by selection. These findings imply a strong pressure of purifying selection in protein evolution, which in the case of serine codon set switches occurs via an initial deleterious substitution quickly followed by a second, compensatory substitution. The result is frequent reversal of amino acid replacements and, at short evolutionary distances, pervasive homoplasy.

Reconstruction of the evolution of microbial defense systems.

BACKGROUND: Evolution of bacterial and archaeal genomes is a highly dynamic process that involves intensive loss of genes as well as gene gain via horizontal transfer, with a lesser contribution from gene duplication. The rates of these processes can be estimated by comparing genomes that are linked by an evolutionary tree. These estimated rates of genome dynamics events substantially differ for different functional classes of genes. The genes involved in defense against viruses and other invading DNA are among those that are gained and lost at the highest rates. RESULTS: We employed a stochastic birth-and-death model to obtain maximum likelihood estimates of the rates of gain and loss of defense genes in 35 groups of closely related bacterial genomes and one group of archaeal genomes. We find that on average, the defense genes experience 1.4 fold higher flux than the rest of microbial genes. This excessive flux of defense genes over the genomic mean is consistent across diverse microbial groups. The few exceptions include intracellular parasites with small, degraded genomes that possess few defense systems which are more stable than in other microbes. Generally, defense genes follow the previously established pattern of genome dynamics, with gene family loss being about 3 times more common than gain and an order of magnitude more common than expansion or contraction of gene families. Case by case analysis of the evolutionary dynamics of defense genes indicates frequent multiple events in the same locus and widespread involvement of mobile elements in the gain and loss of defense genes. CONCLUSIONS: Evolution of microbial defense systems is highly dynamic but, notwithstanding the host-parasite arms race, generally follows the same trends that have been established for the rest of the genes. Apart from the paucity and the low flux of defense genes in parasitic bacteria with deteriorating genomes, there is no clear connection between the evolutionary regime of defense systems and microbial life style.

Mammary alveolar epithelial cells convert to brown adipocytes in post-lactating mice.

During pregnancy and lactation, subcutaneous white adipocytes in the mouse mammary gland transdifferentiate reversibly to milk-secreting epithelial cells. In this study, we demonstrate by transmission electron microscopy that in the post-lactating mammary gland interscapular multilocular adipocytes found close to the mammary alveoli contain milk protein granules. Use of the Cre-loxP recombination system allowed showing that the involuting mammary gland of whey acidic protein-Cre/R26R mice, whose secretory alveolar cells express the lacZ gene during pregnancy, contains some X-Gal-stained and uncoupling protein 1-positive interscapular multilocular adipocytes. These data suggest that during mammary gland involution some milk-secreting epithelial cells in the anterior subcutaneous depot may transdifferentiate to brown adipocytes, highlighting a hitherto unappreciated feature of mouse adipose organ plasticity.

A census of α-helical membrane proteins in double-stranded DNA viruses infecting bacteria and archaea.

BACKGROUND: Viruses are the most abundant and genetically diverse biological entities on earth, yet the repertoire of viral proteins remains poorly explored. As the number of sequenced virus genomes grows into the thousands, and the number of viral proteins into the hundreds of thousands, we report a systematic computational analysis of the point of first-contact between viruses and their hosts, namely viral transmembrane (TM) proteins. RESULTS: The complement of α-helical TM proteins in double-stranded DNA viruses infecting bacteria and archaea reveals large-scale trends that differ from those of their hosts. Viruses typically encode a substantially lower fraction of TM proteins than archaea or bacteria, with the notable exception of viruses with virions containing a lipid component such as a lipid envelope, internal lipid core, or inner membrane vesicle. Compared to bacteriophages, archaeal viruses are substantially enriched in membrane proteins. However, this feature is not always stable throughout the evolution of a viral lineage; for example, TM proteins are not part of the common heritage shared between Lipothrixviridae and Rudiviridae. In contrast to bacteria and archaea, viruses almost completely lack proteins with complicated membrane topologies composed of more than 4 TM segments, with the few detected exceptions being obvious cases of relatively recent horizontal transfer from the host. CONCLUSIONS: The dramatic differences between the membrane proteomes of cells and viruses stem from the fact that viruses do not depend on essential membranes for energy transformation, ion homeostasis, nutrient transport and signaling.

Indomethacin treatment prevents high fat diet-induced obesity and insulin resistance but not glucose intolerance in C57BL/6J mice.

Chronic low grade inflammation is closely linked to obesity-associated insulin resistance. To examine how administration of the anti-inflammatory compound indomethacin, a general cyclooxygenase inhibitor, affected obesity development and insulin sensitivity, we fed obesity-prone male C57BL/6J mice a high fat/high sucrose (HF/HS) diet or a regular diet supplemented or not with indomethacin (±INDO) for 7 weeks. Development of obesity, insulin resistance, and glucose intolerance was monitored, and the effect of indomethacin on glucose-stimulated insulin secretion (GSIS) was measured in vivo and in vitro using MIN6 ß-cells. We found that supplementation with indomethacin prevented HF/HS-induced obesity and diet-induced changes in systemic insulin sensitivity. Thus, HF/HS+INDO-fed mice remained insulin-sensitive. However, mice fed HF/HS+INDO exhibited pronounced glucose intolerance. Hepatic glucose output was significantly increased. Indomethacin had no effect on adipose tissue mass, glucose tolerance, or GSIS when included in a regular diet. Indomethacin administration to obese mice did not reduce adipose tissue mass, and the compensatory increase in GSIS observed in obese mice was not affected by treatment with indomethacin. We demonstrate that indomethacin did not inhibit GSIS per se, but activation of GPR40 in the presence of indomethacin inhibited glucose-dependent insulin secretion in MIN6 cells. We conclude that constitutive high hepatic glucose output combined with impaired GSIS in response to activation of GPR40-dependent signaling in the HF/HS+INDO-fed mice contributed to the impaired glucose clearance during a glucose challenge and that the resulting lower levels of plasma insulin prevented the obesogenic action of the HF/HS diet.

Genomes in turmoil: quantification of genome dynamics in prokaryote supergenomes.

BACKGROUND: Genomes of bacteria and archaea (collectively, prokaryotes) appear to exist in incessant flux, expanding via horizontal gene transfer and gene duplication, and contracting via gene loss. However, the actual rates of genome dynamics and relative contributions of different types of event across the diversity of prokaryotes are largely unknown, as are the sizes of microbial supergenomes, i.e. pools of genes that are accessible to the given microbial species. RESULTS: We performed a comprehensive analysis of the genome dynamics in 35 groups (34 bacterial and one archaeal) of closely related microbial genomes using a phylogenetic birth-and-death maximum likelihood model to quantify the rates of gene family gain and loss, as well as expansion and reduction. The results show that loss of gene families dominates the evolution of prokaryotes, occurring at approximately three times the rate of gain. The rates of gene family expansion and reduction are typically seven and twenty times less than the gain and loss rates, respectively. Thus, the prevailing mode of evolution in bacteria and archaea is genome contraction, which is partially compensated by the gain of new gene families via horizontal gene transfer. However, the rates of gene family gain, loss, expansion and reduction vary within wide ranges, with the most stable genomes showing rates about 25 times lower than the most dynamic genomes. For many groups, the supergenome estimated from the fraction of repetitive gene family gains includes about tenfold more gene families than the typical genome in the group although some groups appear to have vast, 'open' supergenomes. CONCLUSIONS: Reconstruction of evolution for groups of closely related bacteria and archaea reveals an extremely rapid and highly variable flux of genes in evolving microbial genomes, demonstrates that extensive gene loss and horizontal gene transfer leading to innovation are the two dominant evolutionary processes, and yields robust estimates of the supergenome size.

Orthologous gene clusters and taxon signature genes for viruses of prokaryotes.

Viruses are the most abundant biological entities on earth and encompass a vast amount of genetic diversity. The recent rapid increase in the number of sequenced viral genomes has created unprecedented opportunities for gaining new insight into the structure and evolution of the virosphere. Here, we present an update of the phage orthologous groups (POGs), a collection of 4,542 clusters of orthologous genes from bacteriophages that now also includes viruses infecting archaea and encompasses more than 1,000 distinct virus genomes. Analysis of this expanded data set shows that the number of POGs keeps growing without saturation and that a substantial majority of the POGs remain specific to viruses, lacking homologues in prokaryotic cells, outside known proviruses. Thus, the great majority of virus genes apparently remains to be discovered. A complementary observation is that numerous viral genomes remain poorly, if at all, covered by POGs. The genome coverage by POGs is expected to increase as more genomes are sequenced. Taxon-specific, single-copy signature genes that are not observed in prokaryotic genomes outside detected proviruses were identified for two-thirds of the 57 taxa (those with genomes available from at least 3 distinct viruses), with half of these present in all members of the respective taxon. These signatures can be used to specifically identify the presence and quantify the abundance of viruses from particular taxa in metagenomic samples and thus gain new insights into the ecology and evolution of viruses in relation to their hosts.

Contribution of phage-derived genomic islands to the virulence of facultative bacterial pathogens.

Facultative pathogens have extremely dynamic pan-genomes, to a large extent derived from bacteriophages and other mobile elements. We developed a simple approach to identify phage-derived genomic islands and apply it to show that pathogens from diverse bacterial genera are significantly enriched in clustered phage-derived genes compared with related benign strains. These findings show that genome expansion by integration of prophages containing virulence factors is a major route of evolution of facultative bacterial pathogens.

Computational methods for Gene Orthology inference.

Accurate inference of orthologous genes is a pre-requisite for most comparative genomics studies, and is also important for functional annotation of new genomes. Identification of orthologous gene sets typically involves phylogenetic tree analysis, heuristic algorithms based on sequence conservation, synteny analysis, or some combination of these approaches. The most direct tree-based methods typically rely on the comparison of an individual gene tree with a species tree. Once the two trees are accurately constructed, orthologs are straightforwardly identified by the definition of orthology as those homologs that are related by speciation, rather than gene duplication, at their most recent point of origin. Although ideal for the purpose of orthology identification in principle, phylogenetic trees are computationally expensive to construct for large numbers of genes and genomes, and they often contain errors, especially at large evolutionary distances. Moreover, in many organisms, in particular prokaryotes and viruses, evolution does not appear to have followed a simple 'tree-like' mode, which makes conventional tree reconciliation inapplicable. Other, heuristic methods identify probable orthologs as the closest homologous pairs or groups of genes in a set of organisms. These approaches are faster and easier to automate than tree-based methods, with efficient implementations provided by graph-theoretical algorithms enabling comparisons of thousands of genomes. Comparisons of these two approaches show that, despite conceptual differences, they produce similar sets of orthologs, especially at short evolutionary distances. Synteny also can aid in identification of orthologs. Often, tree-based, sequence similarity- and synteny-based approaches can be combined into flexible hybrid methods.