Comparative genomic analysis of eutherian Mas-related G protein-coupled receptor genes

The present study made attempts to update comprehensive eutherian Mas-related G protein-coupled receptor gene data sets, using public eutherian genomic sequence data sets and new genomics and molecular evolution tests. Among 254 potential coding sequences, the most comprehensive gene data set of eutherian Mas-related G protein-coupled receptor genes included 119 complete coding sequences that described eight major gene clusters. The present analysis integrated gene annotations, phylogenetic analysis and protein molecular evolution analysis and first explained differential gene expansion patterns of eutherian Mas-related G protein-coupled receptor genes. The updated classification and nomenclature of eutherian Mas-related G protein-coupled receptor genes were proposed as new framework of future experiments.     

Comparative genomic analysis of eutherian interferon-γ-inducible GTPases

The interferon-γ-inducible GTPases, IFGGs, are intracellular proteins involved in immune response against pathogens. A comprehensive comparative genomic review and analysis of eutherian IFGGs was carried out using public genomic sequences. The 64 eutherian IFGG genes were examined in detail and annotated. The eutherian IFGG promoter types were first catalogued followed by a phylogenetic analysis of eutherian IFGGs, which described five major IFGG clusters. The patterns of differential gene expansions and protein regions that may regulate IFGG catalytic features suggested a new classification of eutherian IFGGs. This mini-review has also provided new tests of reliability of public genomic sequences as well as tests of protein molecular evolution.     

Comparative genomic analysis of soybean flowering genes

Flowering is an important agronomic trait that determines crop yield. Soybean is a major oilseed legume crop used for human and animal feed. Legumes have unique vegetative and floral complexities. Our understanding of the molecular basis of flower initiation and development in legumes is limited. Here, we address this by using a computational approach to examine flowering regulatory genes in the soybean genome in comparison to the most studied model plant, Arabidopsis. For this comparison, a genome-wide analysis of orthologue groups was performed, followed by an in silico gene expression analysis of the identified soybean flowering genes. Phylogenetic analyses of the gene families highlighted the evolutionary relationships among these candidates. Our study identified key flowering genes in soybean and indicates that the vernalisation and the ambient-temperature pathways seem to be the most variant in soybean. A comparison of the orthologue groups containing flowering genes indicated that, on average, each Arabidopsis flowering gene has 2-3 orthologous copies in soybean. Our analysis highlighted that the CDF3, VRN1, SVP, AP3 and PIF3 genes are paralogue-rich genes in soybean. Furthermore, the genome mapping of the soybean flowering genes showed that these genes are scattered randomly across the genome. A paralogue comparison indicated that the soybean genes comprising the largest orthologue group are clustered in a 1.4 Mb region on chromosome 16 of soybean. Furthermore, a comparison with the undomesticated soybean (Glycine soja) revealed that there are hundreds of SNPs that are associated with putative soybean flowering genes and that there are structural variants that may affect the genes of the light-signalling and ambient-temperature pathways in soybean. Our study provides a framework for the soybean flowering pathway and insights into the relationship and evolution of flowering genes between a short-day soybean and the long-day plant, Arabidopsis.     

Housekeeping genes for phylogenetic analysis of eutherian relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.     

Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity

Background

Ralstonia solanacearum is a vascular soil-borne plant pathogen with an unusually broad host range. This economically destructive and globally distributed bacterium has thousands of distinct lineages within a heterogeneous and taxonomically disputed species complex. Some lineages include highly host-adapted strains (ecotypes), such as the banana Moko disease-causing strains, the cold-tolerant potato brown rot strains (also known as R3bv2) and the recently emerged Not Pathogenic to Banana (NPB) strains.

Results

These distinct ecotypes offer a robust model to study host adaptation and the emergence of ecotypes because the polyphyletic Moko strains include lineages that are phylogenetically close to the monophyletic brown rot and NPB strains. Draft genomes of eight new strains belonging to these three model ecotypes were produced to complement the eleven publicly available R. solanacearum genomes. Using a suite of bioinformatics methods, we searched for genetic and evolutionary features that distinguish ecotypes and propose specific hypotheses concerning mechanisms of host adaptation in the R. solanacearum species complex. Genome-wide, few differences were identified, but gene loss events, non-synonymous polymorphisms, and horizontal gene transfer were identified among type III effectors and were associated with host range differences.

Conclusions

This extensive comparative genomics analysis uncovered relatively few divergent features among closely related strains with contrasting biological characteristics; however, several virulence factors were associated with the emergence of Moko, NPB and brown rot and could explain host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1474-8) contains supplementary material, which is available to authorized users.     

Evolution of the interferon alpha gene family in eutherian mammals

Interferon alpha (IFNA) genes code for proteins with important signaling roles during the innate immune response. Phylogenetically, IFNA family members in eutherians (placental mammals) cluster together in a species-specific manner except for closely related species (i.e. Homo sapiens and Pan troglodytes) where gene-specific clustering is evident. Previous research has been unable to clarify whether gene conversion or recent gene duplication accounts for gene-specific clustering, partly because the similarity of members of the IFNA family within species has made it historically difficult to identify the exact composition of IFNA gene families. IFNA gene families were fully characterized in recently available genomes from Canis familiaris, Macaca mulatta, P. troglodytes and Rattus norvegicus, and combined with previously characterized IFNA gene families from H. sapiens and Mus musculus, for the analysis of both whole and partial gene conversion events using a variety of statistical methods. Gene conversion was inferred in every eutherian species analyzed and comparison of the IFNA gene family locus between primate species revealed independent gene duplication in M. mulatta. Thus, both gene conversion and gene duplication have shaped the evolution of the IFNA gene family in eutherian species. Scenarios may be envisaged whereby the increased production of a specific IFN-alpha protein would be beneficial against a particular pathogenic infection. Gene conversion, similar to duplication, provides a mechanism by which the protein product of a specific IFNA gene can be increased.     

Comparative genomic hybridization analysis of Enterococcus faecalis: identification of genes absent from food strains

Enterococcus faecalis, a member of the natural microbiota of animal and human intestinal tracts, is also present as a natural contaminant in a variety of fermented foods. Over the last decade, E. faecalis has emerged as a major cause of nosocomial infections. We investigated the genetic diversity in 30 clinical and food isolates, including strains V583 and MMH594, in order to determine whether clinical and food isolates could be distinguished. Data were obtained using comparative genomic hybridization and specific PCR with a total of 202 probes of E. faecalis, selected using the available V583 genome sequence and part of the MMH594 pathogenicity island. The cognate genes encoded mainly exported proteins. Hybridization data were analyzed by a two-component mixture model that estimates the probability of any given gene to be either present or absent in the strains. A total of 78 genes were found to be variable, as they were absent in at least one isolate. Most of the variable genes were clustered in regions that, in the published V583 sequence, related to prophages or mobile genetic elements. The variable genes were distributed in three main groups: (i) genes equally distributed between clinical and dairy food isolates, (ii) genes absent from dairy food-related isolates, and (iii) genes present in MMH594 and V583 strains only. Further analysis of the distribution of the last gene group in 70 other isolates confirmed that six of the probed genes were always absent in dairy food-related isolates, whereas they were detected in clinical and/or commensal isolates. Two of them corresponded to prophages that were not detected in the cognate isolates, thus possibly extending the number of genes absent from dairy food isolates. Genes specifically detected in clinical isolates may prove valuable for the development of new risk assessment markers for food safety studies and for identification of new factors that may contribute to host colonization or infection.     

Comparative genomic analysis of dha regulon and related genes for anaerobic glycerol metabolism in bacteria

The dihydroxyacetone (dha) regulon of bacteria encodes genes for the anaerobic metabolism of glycerol. In this work, genomic data are used to analyze and compare the dha regulon and related genes in different organisms in silico with respect to gene organization, sequence similarity, and possible functions. Database searches showed that among the organisms, the genomes of which have been sequenced so far, only two, i.e., Klebsiella pneumoniae MGH 78578 and Clostridium perfringens contain a complete dha regulon bearing all known enzymes. The components and their organization in the dha regulon of these two organisms differ considerably from each other and also from the previously partially sequenced dha regulons in Citrobacter freundii, Clostridium pasteurianum, and Clostridium butyricum. Unlike all of the other organisms, genes for the oxidative and reductive pathways of anaerobic glycerol metabolism in C. perfringens are located in two separate organization units on the chromosome. Comparisons of deduced protein sequences of genes with similar functions showed that the dha regulon components in K. pneumoniae and C. freundii have high similarities (80-95%) but lower similarities to those of the Clostridium species (30-80%). Interestingly, the protein sequence similarities among the dha genes of the Clostridium species are in many cases even lower than those between the Clostridium species and K. pneumoniae or C. freundii, suggesting two different types of dha regulon in the Clostridium species studied. The in silico reconstruction and comparison of dha regulons revealed several new genes in the microorganisms studied. In particular, a novel dha kinase that is phosphoenolpyruvate-dependent is identified and experimentally confirmed for K. pneumoniae in addition to the known ATP-dependent dha kinase. This finding gives new insights into the regulation of glycerol metabolism in K. pneumoniae and explains some hitherto not well understood experimental observations.     

Comparative genomic analysis of human and chimpanzee proteases

Proteolytic enzymes are implicated in multiple physiological and pathological processes. The availability of the sequence of the chimpanzee genome has allowed us to determine that the chimpanzee degradome-the repertoire of protease genes from this organism-is composed of at least 559 protease and protease-like genes and is virtually identical to that of human, containing 561 genes. Despite the high degree of conservation between both genomes, we have identified important differences that vary from deletion of whole genes to small insertion/deletion events or single nucleotide changes that lead to the specific gene inactivation in one species, mostly affecting immune system genes. For example, the genes encoding PRSS33/EOS, a macrophage serine protease conserved in most mammals, and GGTLA1 are absent in chimpanzee, while the gene for metalloprotease MMP23A, located in chromosome 1p36, has been specifically duplicated in the human genome together with its neighbor gene CDC2L1. Other differences arise from single nucleotide changes in protease genes, such as NAPSB and CASP12, resulting in the presence of functional genes in chimpanzee and pseudogenes in human. Finally, we have confirmed that the Trypanosoma lytic factor HPR is inactive in chimpanzee, likely contributing to the susceptibility of chimpanzees to T. brucei infection. This study provides the first analysis of the chimpanzee degradome and might contribute to the understanding of the molecular bases underlying variations in host defense mechanisms between human and chimpanzee.     

Comparative genomic analysis of hyperthermophilic archaeal Fuselloviridae viruses

The complete genome sequences of two Sulfolobus spindle-shaped viruses (SSVs) from acidic hot springs in Kamchatka (Russia) and Yellowstone National Park (United States) have been determined. These nonlytic temperate viruses were isolated from hyperthermophilic Sulfolobus hosts, and both viruses share the spindle-shaped morphology characteristic of the Fuselloviridae family. These two genomes, in combination with the previously determined SSV1 genome from Japan and the SSV2 genome from Iceland, have allowed us to carry out a phylogenetic comparison of these geographically distributed hyperthermal viruses. Each virus contains a circular double-stranded DNA genome of approximately 15 kbp with approximately 34 open reading frames (ORFs). These Fusellovirus ORFs show little or no similarity to genes in the public databases. In contrast, 18 ORFs are common to all four isolates and may represent the minimal gene set defining this viral group. In general, ORFs on one half of the genome are colinear and highly conserved, while ORFs on the other half are not. One shared ORF among all four genomes is an integrase of the tyrosine recombinase family. All four viral genomes integrate into their host tRNA genes. The specific tRNA gene used for integration varies, and one genome integrates into multiple loci. Several unique ORFs are found in the genome of each isolate.     

Comparative genomic analysis of bacteriocin-producing Weissella cibaria 110

Weissella cibaria 110 was isolated from plaa-som, a Thai fermented fish product, and known to produce the weissellicin 110 bacteriocin. We carried out comprehensive comparative genomic analysis of W. cibaria 110 with four other non-bacteriocin-producing W. cibaria strains and identified potential antibiotic-resistant genes. We further identified a type III restriction-modification system, a TA system, and a bacteriocin gene cluster that are unique in W. cibaria 110. Genes related to bacteriocin biosynthesis are organized in clusters and are encoded with minimum genetic machinery consisting of structural cognate immunity genes, including ABC transporter and immunity protein. Finally, we predicted W. cibaria 110 to produce a class IId bacteriocin, weissellicin 110, which is 31 amino acids in length and contains a 21-amino-acid N-terminal leader peptide. This is the first bacteriocin-producing sequencing genome in W. cibaria, and we describe the difference between the bacteriocin-producing and non bacteriocin-producing strains from genome point of view.

     

Comparative genomic analysis of 18 Pseudomonas aeruginosa bacteriophages

A genomic analysis of 18 P. aeruginosa phages, including nine newly sequenced DNA genomes, indicates a tremendous reservoir of proteome diversity, with 55% of open reading frames (ORFs) being novel. Comparative sequence analysis and ORF map organization revealed that most of the phages analyzed displayed little relationship to each other.     

Comparative genomic analysis of Mycobacterium tuberculosis clinical isolates

Background

Due to excessive antibiotic use, drug-resistant Mycobacterium tuberculosis has become a serious public health threat and a major obstacle to disease control in many countries. To better understand the evolution of drug-resistant M. tuberculosis strains, we performed whole genome sequencing for 7 M. tuberculosis clinical isolates with different antibiotic resistance profiles and conducted comparative genomic analysis of gene variations among them.

Results

We observed that all 7 M. tuberculosis clinical isolates with different levels of drug resistance harbored similar numbers of SNPs, ranging from 1409–1464. The numbers of insertion/deletions (Indels) identified in the 7 isolates were also similar, ranging from 56 to 101. A total of 39 types of mutations were identified in drug resistance-associated loci, including 14 previously reported ones and 25 newly identified ones. Sixteen of the identified large Indels spanned PE-PPE-PGRS genes, which represents a major source of antigenic variability. Aside from SNPs and Indels, a CRISPR locus with varied spacers was observed in all 7 clinical isolates, suggesting that they might play an important role in plasticity of the M. tuberculosis genome. The nucleotide diversity (Л value) and selection intensity (dN/dS value) of the whole genome sequences of the 7 isolates were similar. The dN/dS values were less than 1 for all 7 isolates (range from 0.608885 to 0.637365), supporting the notion that M. tuberculosis genomes undergo purifying selection. The Л values and dN/dS values were comparable between drug-susceptible and drug-resistant strains.

Conclusions

In this study, we show that clinical M. tuberculosis isolates exhibit distinct variations in terms of the distribution of SNP, Indels, CRISPR-cas locus, as well as the nucleotide diversity and selection intensity, but there are no generalizable differences between drug-susceptible and drug-resistant isolates on the genomic scale. Our study provides evidence strengthening the notion that the evolution of drug resistance among clinical M. tuberculosis isolates is clearly a complex and diversified process.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-469) contains supplementary material, which is available to authorized users.     

Comparative mapping, genomic structure, and expression analysis of eight pseudo-response regulator genes in Brassica rapa

Circadian clocks regulate plant growth and development in response to environmental factors. In this function, clocks influence the adaptation of species to changes in location or climate. Circadian-clock genes have been subject of intense study in models such as Arabidopsis thaliana but the results may not necessarily reflect clock functions in species with polyploid genomes, such as Brassica species, that include multiple copies of clock-related genes. The triplicate genome of Brassica rapa retains high sequence-level co-linearity with Arabidopsis genomes. In B. rapa we had previously identified five orthologs of the five known Arabidopsis pseudo-response regulator (PRR) genes that are key regulators of the circadian clock in this species. Three of these B. rapa genes, BrPRR1, BrPPR5, and BrPPR7, are present in two copies each in the B. rapa genome, for a total of eight B. rapa PRR (BrPRR) orthologs. We have now determined sequences and expression characteristics of the eight BrPRR genes and mapped their positions in the B. rapa genome. Although both members of each paralogous pair exhibited the same expression pattern, some variation in their gene structures was apparent. The BrPRR genes are tightly linked to several flowering genes. The knowledge about genome location, copy number variation and structural diversity of these B. rapa clock genes will improve our understanding of clock-related functions in this important crop. This will facilitate the development of Brassica crops for optimal growth in new environments and under changing conditions.