Search

Background Kuruma shrimp, a major commercial shrimp species in the world, has two cryptic or sibling species, Marsupenaeus japonicus and Marsupenaeus pulchricaudatus. Codon usage analysis would contribute to our understanding of the genetic and evolutionary characteristics of the two Marsupenaeus species. In this study, we analyzed codon usage and related indices using coding sequences (CDSs) from RNA-seq data. Results Using CodonW 1.4.2 software, we performed the codon bias analysis of transcriptomes obtained from hepatopancreas tissues, which indicated weak codon bias. Almost all parameters had similar correlations for both species. The gene expression level (FPKM) was negatively correlated with A/T3s. We determined 12 and 14 optimal codons for M. japonicus and M. pulchricaudatus, respectively, and all optimal codons have a C/G-ending. The two Marsupenaeus species had different usage frequencies of codon pairs, which contributed to further analysis of transcriptional differences between them. Orthologous genes that underwent positive selection (ω > 1) had a higher correlation coefficient than that of experienced purifying selection (ω Conclusions This study provides a relatively comprehensive understanding of the correlations among codon usage bias, gene expression, and selection pressures of CDSs for M. japonicus and M. pulchricaudatus. The genetic evolution was driven by mutations and selection pressure. Moreover, the results point out new insights into the specificities and evolutionary characteristics of the two Marsupenaeus species.

Proanthocyanidins (PAs) are a group of oligomers or polymers composed of monomeric flavanols. They offer many benefits for human fitness, such as antioxidant, anticancer, and anti-inflammatory activities. To date, three types of PA have been observed in nature: procyanidins, propelargonidins, and prodelphinidins. These are synthesized as some of the end-products of the flavonoid pathway by different consecutive enzymatic activities, from the same precursor—naringenin. Although the general biosynthetic pathways of PAs have been reported in a few model plant species, little is known about the species-specific pathways in major crops containing different types of PA. In the present study, we identified the species-specific pathways in 10 major crops, based on the presence/absence of flavanol-based intermediates in the metabolic pathway, and found 202 orthologous genes in the reference genomic database of each species, which may encode for key enzymes involved in the biosynthetic pathways of PAs. Parallel enzymatic reactions in the pathway are responsible for the ratio between PAs and anthocyanins, as well as among the three types of PAs. Our study suggests a promising strategy for molecular breeding, to regulate the content of PAs and anthocyanins and improve the nutritional quality of food sources globally.

Background Auxins play an important role in plant growth and development; the auxins responsive gene; auxin/indole-3-acetic acid (Aux/IAA), small auxin-up RNAs (SAUR) and Gretchen Hagen3 (GH3) control their mechanisms. The GH3 genes function in homeostasis by the catalytic activities in auxin conjugation and bounding free indole-3-acetic acid (IAA) to amino acids. Results In our study, we identified the GH3 genes in three cotton species; Gossypium hirsutum, Gossypium arboreum and Gossypium raimondii, analyzed their chromosomal distribution, phylogenetic relationships, cis-regulatory element function and performed virus induced gene silencing of the novel Gh_A08G1120 (GH3.5) gene. The phylogenetic tree showed four clusters of genes with clade 1, 3 and 4 having mainly members of the GH3 of the cotton species while clade 2 was mainly members belonging to Arabidopsis. There were no paralogous genes, and few orthologous genes were observed between Gossypium and other species. All the GO terms were detected, but only 14 genes were found to have described GO terms in upland cotton, more biological functions were detected, as compared to the other functions. The GH3.17 subfamily harbored the highest number of the cis-regulatory elements, most having promoters towards dehydration-responsiveness. The RNA expression analysis revealed that 10 and 8 genes in drought and salinity stress conditions respectively were upregulated in G. hirsutum. All the genes that were upregulated in plants under salt stress conditions were also upregulated in drought stress; moreover, Gh_A08G1120 (GH3.5) exhibited a significant upregulation across the two stress factors. Functional characterization of Gh_A08G1120 (GH3.5) through virus-induced gene silencing (VIGS) revealed that the VIGS plants ability to tolerate drought and salt stresses was significantly reduced compared to the wild types. The chlorophyll content, relative leaf water content (RLWC), and superoxide dismutase (SOD) concentration level were reduced significantly while malondialdehyde concentration and ion leakage as a measure of cell membrane stability (CMS) increased in VIGS plants under drought and salt stress conditions. Conclusion This study revealed the significance of the GH3 genes in enabling the plant’s adaptation to drought and salt stress conditions as evidenced by the VIGS results and RT-qPCR analysis. Electronic supplementary material The online version of this article (10.1186/s12863-019-0756-6) contains supplementary material, which is available to authorized users.

Made available in DSpace on 2019-10-06T16:35:19Z (GMT). No. of bitstreams: 0 Previous issue date: 2019-01-01. Added 1 bitstream(s) on 2019-10-09T18:35:28Z : No. of bitstreams: 1 S1415-47572019000100145.pdf: 1807906 bytes, checksum: 429de26cfda2d5142dade2f2716528e9 (MD5) Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) The filter cake from sugar cane processing is rich in organic matter and nutrients, which favors the proliferation of microorganisms with potential to deconstruct plant biomass. From the metagenomic data of this material, we assembled a draft genome that was phylogenetically related to Thermomonospora curvata DSM 43183, which shows the functional and ecological importance of this bacterium in the filter cake. Thermomonospora is a gram-positive bacterium that produces cellulases in compost, and it can survive temperatures of 60 ºC. We identified a complete set of biomass depolymerizing enzymes in the draft genome of Thermomonospora sp. CIT 1, such as α-amylase, catalase-peroxidases, β-mannanase, and arabinanase, demonstrating the potential of this bacterium to deconstruct the components of starch, lignin, and hemicellulose. In addition, the draft genome of Thermomonospora sp. CIT 1 contains 18 genes that do not share identity with five other species of Thermomonospora, suggesting that this bacterium has different genetic characteristics than those present in genomes reported so far for this genus. These findings add a new dimension to the current understanding of the functional profile of this microorganism that inhabits agro-industrial waste, which may boost new gene discoveries and be of importance for application in the production of bioethanol. Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP) Departamento de Tecnologia Laboratório de Bioinformática Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP) Laboratório Multiusuário Centralizado para Sequenciamento de DNA em Larga Escala e Análise de Expressão Gênica (LMSeq) Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP) Departamento de Biologia Aplicada à Agropecuária Laboratório de Genética Aplicada Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP) Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP) Departamento de Tecnologia Laboratório de Bioinformática Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP) Laboratório Multiusuário Centralizado para Sequenciamento de DNA em Larga Escala e Análise de Expressão Gênica (LMSeq) Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP) Departamento de Biologia Aplicada à Agropecuária Laboratório de Genética Aplicada Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista (UNESP)

The micronutrients iron (Fe) and zinc (Zn) are not only essential for plant survival and proliferation but are crucial for human health. Increasing Fe and Zn levels in edible parts of plants, known as biofortification, is seen a sustainable approach to alleviate micronutrient deficiency in humans. Wheat, as one of the leading staple foods worldwide, is recognized as a prioritized choice for Fe and Zn biofortification. However, to date, limited molecular and physiological mechanisms have been elucidated for Fe and Zn homeostasis in wheat. The expanding molecular understanding of Fe and Zn homeostasis in model plants is providing invaluable resources to biofortify wheat. Recent advancements in NGS (next generation sequencing) technologies coupled with improved wheat genome assembly and high-throughput genotyping platforms have initiated a revolution in resources and approaches for wheat genetic investigations and breeding. Here, we summarize molecular processes and genes involved in Fe and Zn homeostasis in the model plants Arabidopsis and rice, identify their orthologs in the wheat genome, and relate them to known wheat Fe/Zn QTL (quantitative trait locus/loci) based on physical positions. The current study provides the first inventory of the genes regulating grain Fe and Zn homeostasis in wheat, which will benefit gene discovery and breeding, and thereby accelerate the release of Fe- and Zn-enriched wheats.

Background The circadian clock drives endogenous 24-h rhythms that allow organisms to adapt and prepare for predictable and repeated changes in their environment throughout the day-night (diurnal) cycle. Many components of the circadian clock in Arabidopsis thaliana have been functionally characterized, but comparatively little is known about circadian clocks in grass species including major crops like maize and sorghum. Results Comparative research based on protein homology and diurnal gene expression patterns suggests the function of some predicted clock components in grasses is conserved with their Arabidopsis counterparts, while others have diverged in function. Our analysis of diurnal gene expression in three panicoid grasses sorghum, maize, and foxtail millet revealed conserved and divergent evolution of expression for core circadian clock genes and for the overall transcriptome. We find that several classes of core circadian clock genes in these grasses differ in copy number compared to Arabidopsis, but mostly exhibit conservation of both protein sequence and diurnal expression pattern with the notable exception of maize paralogous genes. We predict conserved cis-regulatory motifs shared between maize, sorghum, and foxtail millet through identification of diurnal co-expression clusters for a subset of 27,196 orthologous syntenic genes. In this analysis, a Cochran–Mantel–Haenszel based method to control for background variation identified significant enrichment for both expected and novel 6–8 nucleotide motifs in the promoter regions of genes with shared diurnal regulation predicted to function in common physiological activities. Conclusions This study illustrates the divergence and conservation of circadian clocks and diurnal regulatory networks across syntenic orthologous genes in panacoid grass species. Further, conserved local regulatory sequences contribute to the architecture of these diurnal regulatory networks that produce conserved patterns of diurnal gene expression.

Significance This work represents a large, order-wide evolutionary analysis of the order Chiroptera (bats). Our pipeline for assembling sequence data and curating orthologous multiple sequence alignments includes methods for improving results when combining genomic and transcriptomic data sources. The resulting phylogenetic tree divides the order Chiroptera into Yinpterochiroptera and Yangochiroptera, in disagreement with the previous division into Megachiroptera and Microchiroptera and in agreement with some other recent molecular studies, and also provides evidence for other contested branch placements. We also performed a genome-wide analysis of positive selection and found 181 genes with signatures of positive selection. Enrichment analysis shows these positively selected genes to be primarily related to immune responses but also, surprisingly, collagen formation., Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus, the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus, the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.

International audience; Plant cell walls play a fundamental role in several plant traits and also influence crop use as livestock nutrition or biofuel production. The Glycosyltransferase family 61 (GT61) is involved in the synthesis of cell wall xylans. In grasses (Poaceae), a copy number expansion was reported for the GT61 family, and raised the question of the evolutionary history of this gene family in a broader taxonomic context. A phylogenetic study was performed on GT61 members from 13 species representing the major angiosperm clades, in order to classify the genes, reconstruct the evolutionary history of this gene family and study its expansion in monocots. Four orthogroups (OG) were identified in angiosperms with two of them displaying a copy number expansion in monocots. These copy number expansions resulted from both tandem and segmental duplications during the genome evolution of monocot lineages. Positive selection footprints were detected on the ancestral branch leading to one of the orthogroups suggesting that the gene number expansion was accompanied by functional diversification, at least partially. We propose an OG-based classification framework for the GT61 genes at different taxonomic levels of the angiosperm useful for any further functional or translational biology study.

This contribution contains data associated to the research article entitled “Exploring the genes of yerba mate (Ilex paraguariensis A. St.-Hil.) by NGS and de novo transcriptome assembly” (Debat et al., 2014) [1]. By means of a bioinformatic approach involving extensive NGS data analyses, we provide a resource encompassing the full transcriptome assembly of yerba mate, the first available reference for the Ilex L. genus. This dataset (Supplementary files 1 and 2) consolidates the transcriptome-wide assembled sequences of I. paraguariensis with further comprehensive annotation of the protein coding genes of yerba mate via the integration of Arabidopsis thaliana databases. The generated data is pivotal for the characterization of agronomical relevant genes in the tree crop yerba mate -a non-model species- and related taxa in Ilex. The raw sequencing data dissected here is available at DDBJ/ENA/GenBank (NCBI Resource Coordinators, 2016) [2] Sequence Read Archive (SRA) under the accession SRP043293 and the assembled sequences have been deposited at the Transcriptome Shotgun Assembly Sequence Database (TSA) under the accession GFHV00000000. Fil: Aguilera, Patricia Mabel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Biología Subtropical. Universidad Nacional de Misiones. Instituto de Biología Subtropical; Argentina Fil: Debat, Humberto Julio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Biología Subtropical. Universidad Nacional de Misiones. Instituto de Biología Subtropical; Argentina Fil: Grabiele, Mauro. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Biología Subtropical. Universidad Nacional de Misiones. Instituto de Biología Subtropical; Argentina. Instituto Nacional de Tecnología Agropecuaria. Centro de Investigaciones Agropecuarias. Instituto de Patología Vegetal; Argentina

Background Polyploidy is considered a major driving force in genome expansion, yielding duplicated genes whose expression may be conserved or divergence as a consequence of polyploidization. Results We compared the genome sequences of tetraploid cotton (Gossypium hirsutum) and its two diploid progenitors, G. arboreum and G. raimondii, and found that the bHLH genes were conserved over the polyploidization. Oppositely, the expression of the homeolgous gene pairs was diversified. The biased homeologous proportion for bHLH family is significantly higher (64.6%) than the genome wide homeologous expression bias (40%). Compared with cacao (T. cacao), orthologous genes only accounted for a small proportion (41.7%) of whole cotton bHLHs family. The further Ks analysis indicated that bHLH genes underwent at least two distinct episodes of whole genome duplication: a recent duplication (1.0–60.0 million years ago, MYA, 0.005 60.0 MYA, 0.312

Gene duplication is one of the major evolutionary mechanisms providing raw material for the generation of genes with new or modified functions. The yeast Saccharomyces cerevisiae originated after an allopolyploidization event, which involved mating between two different ancestral yeast species. ScALT1 and ScALT2 codify proteins with 65% identity, which were proposed to be paralogous alanine transaminases. Further analysis of their physiological role showed that while ScALT1 encodes an alanine transaminase which constitutes the main pathway for alanine biosynthesis and the sole pathway for alanine catabolism, ScAlt2 does not display alanine transaminase activity and is not involved in alanine metabolism. Moreover, phylogenetic studies have suggested that ScALT1 and ScALT2 come from each one of the two parental strains which gave rise to the ancestral hybrid. The present work has been aimed to the understanding of the properties of the ancestral type Lacchancea kluyveri LkALT1 and Kluyveromyces lactis KlALT1, alanine transaminases in order to better understand the ScALT1 and ScALT2 evolutionary history. These ancestral -type species were chosen since they harbor ALT1 genes, which are related to ScALT2. Presented results show that, although LkALT1 and KlALT1 constitute ScALT1 orthologous genes, encoding alanine transaminases, both yeasts display LkAlt1 and KlAlt1 independent alanine transaminase activity and additional unidentified alanine biosynthetic and catabolic pathway(s). Furthermore, phenotypic analysis of null mutants uncovered the fact that KlAlt1 and LkAlt1 have an additional role, not related to alanine metabolism but is necessary to achieve wild type growth rate. Our study shows that the ancestral alanine transaminase function has been retained by the ScALT1 encoded enzyme, which has specialized its catabolic character, while losing the alanine independent role observed in the ancestral type enzymes. The fact that ScAlt2 conserves 64% identity with LkAlt1 and 66% with KlAlt1, suggests that ScAlt2 diversified after the ancestral hybrid was formed. ScALT2 functional diversification resulted in loss of both alanine transaminase activity and the additional alanine-independent LkAlt1 function, since ScALT2 did not complement the Lkalt1Δ phenotype. It can be concluded that LkALT1 and KlLALT1 functional role as alanine transaminases was delegated to ScALT1, while ScALT2 lost this role during diversification.

WRKY transcription factors have been extensively characterized in the past 20 years, but in wheat, studies on WRKY genes and their function are lagging behind many other species. To explore the function of wheat WRKY genes, we identified a TaWRKY68 gene from a common wheat cultivar. It encodes a protein comprising 313 amino acids which harbors 19 conserved motifs or active sites. Gene expression patterns were determined by analyzing microarray data of TaWRKY68 in wheat and of orthologous genes from maize, rice and barley using Genevestigator. TaWRKY68 orthologs were identified and clustered using DELTA-BLAST and COBALT programs available at NCBI. The results showed that these genes, which are expressed in all tissues tested, had relatively higher levels in the roots and were up-regulated in response to biotic stresses. Bioinformatics results were confirmed by RT-PCR experiments using wheat plants infected by Agrobacterium tumefaciens and Blumeria graminis, or treated with Deoxynivalenol, a Fusarium graminearum-induced mycotoxin in wheat or barley. In summary, TaWRKY68 functions differ during plant developmental stages and might be representing a hub gene function in wheat responses to various biotic stresses. It was also found that including data from major cereal genes in the bioinformatics analysis gave more accurate and comprehensive predictions of wheat gene functions.

Many important crops have received little attention by the scientific community, either because they are not considered economically important or due to their large and complex genomes. De novo transcriptome assembly, using next-generation sequencing data, is an attractive option for the study of these orphan crops. In spite of the large amount of sequencing data that can be generated, there is currently a lack of tools which can effectively help molecular breeders and biologists to mine this type of information. Our goal was to develop a tool that enables molecular breeders, without extensive bioinformatics knowledge, to efficiently study de novo transcriptome data from any orphan crop (http://www.bioinformatics.nl/denovobrowser/db/species/index). The Orphan Crops Browser has been designed to facilitate the following tasks (1) search and identification of candidate transcripts based on phylogenetic relationships between orthologous sequence data from a set of related species and (2) design specific and degenerate primers for expression studies in the orphan crop of interest. To demonstrate the usability and reliability of the browser, it was used to identify the putative orthologues of 17 known lignin biosynthetic genes from maize and sugarcane in the orphan crop Miscanthus sinensis. Expression studies in miscanthus stem internode tissue differing in maturation were subsequently carried out, to follow the expression of these genes during lignification. Our results showed a negative correlation between lignin content and gene expression. The present data are in agreement with recent findings in maize and other crops, and it is further discussed in this paper. Electronic supplementary material The online version of this article (doi:10.1007/s11032-015-0430-2) contains supplementary material, which is available to authorized users.

Background and aims Arachis batizocoi is a wild relative of cultivated peanut (A. hypogaea), an allotetraploid with an AABB genome. Arachis batizocoi was once considered the ancestral donor of the peanut B genome, but cytogenetics and DNA phylogenies have indicated a new genome classification, 'K'. These observations seem inconsistent with genetic studies and breeding that have shown that A. batizocoi can behave as a B genome. Methods The genetic behaviour, genome composition and phylogenetic position of A. batizocoi were studied using controlled hybridizations, induced tetraploidy, whole-genome in situ fluorescent hybridization (GISH) and molecular phylogenetics. Key results Sterile diploid hybrids containing AK genomes were obtained using A. batizocoi and the A genome species A. duranensis, A. stenosperma, A. correntina or A. villosa. From these, three types of AAKK allotetraploids were obtained, each in multiple independent polyploidy events. Induced allotetraploids were vigorous and fertile, and were hybridized to A. hypogaea to produce F1 hybrids. Even with the same parental combination, fertility of these F1 hybrids varied greatly, suggesting the influence of stochastic genetic or epigenetic events. Interestingly, hybrids with A. hypogaea ssp. hypogaea were significantly more fertile than those with the subspecies fastigiata. GISH in cultivated × induced allotetraploids hybrids (harbouring AABK genomes) and a molecular phylogeny using 16 intron sequences showed that the K genome is distinct, but more closely related to the B than to the A genome. Conclusions The K genome of A. batizocoi is more related to B than to the A genome, but is distinct. As such, when incorporated in an induced allotetraploid (AAKK) it can behave as a B genome in crosses with peanut. However, the fertility of hybrids and their progeny depends upon the compatibility of the A genome interactions. The genetic distinctness of A. batizocoi makes it an important source of allelic diversity in itself, especially in crosses involving A. hypogaea ssp. hypogaea.

A gene order database of 32 completely sequenced plastid genomes was developed. The data structure is formally identical to that of the feature tables in the major GenBank/EMBL/DDBJ databases. The quality of annotations was largely improved. A normalizing gene-labeling system across the complete plastid genomes was developed so that comparative studies are made available without having to go back to sequence analysis. Many incorrect coordinates of tRNA-encoding regions found in the major databases were corrected. We attempted to distinctively label tRNA genes with the anticodon sequence CAT, which encodes either the initiator tRNA, elongator tRNA, or Ile-tRNA. The database is available at http://www.rs.noda.tus.ac.jp/~kunisawa.

Diatoms are among the most diverse groups of phytoplankton in the ocean. Despite their widely recognized influence on ocean ecosystems and global biogeochemistry, little is known about the impact of this diversity on large-scale processes. Here, we examined the ramifications of between-species diversity by documenting the transcriptional response of three diatoms —Thalassiosira pseudonana, Fragilariopsis cylindrus, and Pseudo-nitzschia multiseries —to the onset of nitrate limitation of growth, a common limiting nutrient in the ocean. The three species shared 5583 clusters of orthologous genes based on OrthoMCL clustering of publically available diatom genomes. These clusters represent 30–54% of the predicted genes in each diatom genome. Less than 5% of genes within these core clusters displayed the same transcriptional responses across species when growth was limited by nitrate availability. Orthologs, such as those involved in nitrogen uptake and assimilation, as well as carbon metabolism, were differently expressed across the three species. The two pennate diatoms, F. cylindrus and P. multiseries, shared 3839 clusters without orthologs in the genome of the centric diatom T. pseudonana. A majority of these pennate-clustered genes, as well as the non-orthologous genes in each species, had minimal annotation information, but were often significantly differentially expressed under nitrate limitation, indicating their potential importance in the response to nitrogen availability. Despite these variations in the specific transcriptional response of each diatom, overall transcriptional patterns suggested that all three diatoms displayed a common physiological response to nitrate limitation that consisted of a general reduction in carbon fixation and carbohydrate and fatty acid metabolism and an increase in nitrogen recycling. Characterization of these finely tuned responses will help to better predict which types of diatoms will bloom under which sets of environmental factors.

Background With genomic sequences of many closely related bacterial strains made available by deep sequencing, it is now possible to investigate trends in prokaryotic microevolution. Positive selection is a sub-process of microevolution, in which a particular mutation is favored, causing the allele frequency to continuously shift in one direction. Wide scanning of prokaryotic genomes has shown that positive selection at the molecular level is much more frequent than expected. Genes with significant positive selection may play key roles in bacterial adaption to different environmental pressures. However, selection pressure analyses are computationally intensive and awkward to configure. Results Here we describe an open access web server, which is designated as PSP (Positive Selection analysis for Prokaryotic genomes) for performing evolutionary analysis on orthologous coding genes, specially designed for rapid comparison of dozens of closely related prokaryotic genomes. Remarkably, PSP facilitates functional exploration at the multiple levels by assignments and enrichments of KO, GO or COG terms. To illustrate this user-friendly tool, we analyzed Escherichia coli and Bacillus cereus genomes and found that several genes, which play key roles in human infection and antibiotic resistance, show significant evidence of positive selection. PSP is freely available to all users without any login requirement at: http://db-mml.sjtu.edu.cn/PSP/. Conclusions PSP ultimately allows researchers to do genome-scale analysis for evolutionary selection across multiple prokaryotic genomes rapidly and easily, and identify the genes undergoing positive selection, which may play key roles in the interactions of host-pathogen and/or environmental adaptation.

Background Aspartic proteases (APs) are a large family of proteolytic enzymes found in almost all organisms. In plants, they are involved in many biological processes, such as senescence, stress responses, programmed cell death, and reproduction. Prior to the present study, no grape AP gene(s) had been reported, and their research on woody species was very limited. Results In this study, a total of 50 AP genes (VvAP) were identified in the grape genome, among which 30 contained the complete ASP domain. Synteny analysis within grape indicated that segmental and tandem duplication events contributed to the expansion of the grape AP family. Additional analysis between grape and Arabidopsis demonstrated that several grape AP genes were found in the corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the divergence of grape and Arabidopsis. Phylogenetic relationships of the 30 VvAPs with the complete ASP domain and their Arabidopsis orthologs, as well as their gene and protein features were analyzed and their cellular localization was predicted. Moreover, expression profiles of VvAP genes in six different tissues were determined, and their transcript abundance under various stresses and hormone treatments were measured. Twenty-seven VvAP genes were expressed in at least one of the six tissues examined; nineteen VvAPs responded to at least one abiotic stress, 12 VvAPs responded to powdery mildew infection, and most of the VvAPs responded to SA and ABA treatments. Furthermore, integrated synteny and phylogenetic analysis identified orthologous AP genes between grape and Arabidopsis, providing a unique starting point for investigating the function of grape AP genes. Conclusions The genome-wide identification, evolutionary and expression analyses of grape AP genes provide a framework for future analysis of AP genes in defining their roles during stress response. Integrated synteny and phylogenetic analyses provide novel insight into the functions of less well-studied genes using information from their better understood orthologs.

The relationship between genome characteristics and several human diseases has been a central research goal in genomics. Many studies have shown that specific gene patterns, such as amino acid repetitions, are associated with human diseases. However, several open questions still remain, such as, how these tandem repeats appeared in the evolutionary path or how they have evolved in orthologous genes of related organisms. In this paper, we present a computational solution that facilitates comparative studies of orthologous genes from various organisms. The application uses various web services to gather gene sequence information, local algorithms for tandem repeats identification and similarity measures for gene clustering. published

Background A genome-wide comparative analysis of human and mouse gene expression patterns was performed in order to evaluate the evolutionary divergence of mammalian gene expression. Tissue-specific expression profiles were analyzed for 9,105 human-mouse orthologous gene pairs across 28 tissues. Expression profiles were resolved into species-specific coexpression networks, and the topological properties of the networks were compared between species. Results At the global level, the topological properties of the human and mouse gene coexpression networks are, essentially, identical. For instance, both networks have topologies with small-world and scale-free properties as well as closely similar average node degrees, clustering coefficients, and path lengths. However, the human and mouse coexpression networks are highly divergent at the local level: only a small fraction ( Conclusion The dissonance between global versus local network divergence suggests that the interspecies similarity of the global network properties is of limited biological significance, at best, and that the biologically relevant aspects of the architectures of gene coexpression are specific and particular, rather than universal. Nevertheless, there is substantial evolutionary conservation of the local network structure which is compatible with the notion that gene coexpression networks are subject to purifying selection.