ViroBLAST: a stand-alone BLAST web server for flexible queries of multiple databases and user's datasets

Development of a wheat single gene FISH map for analyzing homoeologous relationship and chromosomal rearrangements within the Triticeae

A cytogenetic map of wheat was constructed using FISH with cDNA probes. FISH markers detected homoeology and chromosomal rearrangements of wild relatives, an important source of genes for wheat improvement. To transfer agronomically important genes from wild relatives to bread wheat (Triticum aestivum L., 2n = 6 x = 42, AABBDD) by induced homoeologous recombination, it is important to know the chromosomal relationships of the species involved. Fluorescence in situ hybridization (FISH) can be used to study chromosome structure. The genomes of allohexaploid bread wheat and other species from the Triticeae tribe are colinear to some extent, i.e., composed of homoeoloci at similar positions along the chromosomes, and with genic regions being highly conserved. To develop cytogenetic markers specific for genic regions of wheat homoeologs, we selected more than 60 full-length wheat cDNAs using BLAST against mapped expressed sequence tags and used them as FISH probes. Most probes produced signals on all three homoeologous chromosomes at the expected positions. We developed a wheat physical map with several cDNA markers located on each of the 14 homoeologous chromosome arms. The FISH markers confirmed chromosome rearrangements within wheat genomes and were successfully used to study chromosome structure and homoeology in wild Triticeae species. FISH analysis detected 1 U-6 U chromosome translocation in the genome of Aegilops umbellulata, showed colinearity between chromosome A of Ae. caudata and group-1 wheat chromosomes, and between chromosome arm 7S#3 L of Thinopyrum intermedium and the long arm of the group-7 wheat chromosomes.

VH1-44 gene usage defines a subset of canine B-cell lymphomas associated with better patient survival

The use of specific immunoglobulin heavy chain variable region (VH) genes has been associated with increased patient survival in human B-cell lymphomas (hBCL). Given the similarity of human and canine BCL (cBCL) in morphology and clinical treatment, we examined the choice of VH in cBCL and determined whether VH gene selection was a distinct feature associated with survival time in dogs. VH gene selection and mutational status in 52 cBCL, including 29 diffuse large B-cell lymphomas (cDLBCL, the most common subtype of cBCL), were analyzed by comparison with the 80 published canine germline VH gene sequences. We further examined the prognostic impact of the subgroups defined by these features on canine survival. We found that VH1-44 was preferentially expressed in the majority of the 52 cBCLs (60%) as well as in the majority of the cDLBCL subset (59%). VH1-44 gene expression was associated with a statistically better overall survival (p=0.039) in cBCL patients, as well as in the cDLBCL subset of patients (p=0.038). These findings suggest that VH gene selection in cBCL is not random and may therefore have functional implications for cBCL lymphomagenesis, in addition to being a useful prognostic biomarker.

Changes in the ghrelin hormone pathway maybe part of an unusual gastric system in monotremes

Ghrelin is a growth hormone (GH)-releasing and appetite-regulating peptide predominately released from the stomach. Ghrelin is evolutionarily highly conserved and known to have a wide range of functions including the regulation of metabolism by maintaining an insulin-glucose balance. The peptide is produced as a single proprotein, which is later proteolytically cleaved. Ghrelin exerts its biological function after O-n-octanoylation at residue serine 3, which is catalyzed by ghrelin O-acyl transferase (GOAT) and allows binding to the growth hormone secretagogue receptor (GHS-R 1a). Genes involved in the ghrelin pathway have been identified in a broad range of vertebrate species, however, little is known about this pathway in the basal mammalian lineage of monotremes (platypus and echidna). Monotremes are particularly interesting in this context, as they have undergone massive changes in stomach anatomy and physiology, accompanied by a striking loss of genes involved in gastric function. In this study, we investigated genes in the ghrelin pathway in monotremes. Using degenerate PCR, database searches and synteny analysis we found that genes encoding ghrelin and GOAT are missing in the platypus genome, whilst, as has been reported in other species, the GHSR is present and expressed in brain, pancreas, kidney, intestine, heart and stomach. This is the first report suggesting the loss of ghrelin in a mammal. The loss of this gene may be related to changes to the platypus digestive system and raises questions about the control of blood glucose levels and insulin response in monotreme mammals. In addition, the conservation of the ghrelin receptor gene in platypus indicates that another ligand(s) maybe acting via this receptor in monotremes.

Two independent pathways for self-recognition in Proteus mirabilis are linked by type VI-dependent export

Swarming colonies of the bacterium Proteus mirabilis are capable of self-recognition and territorial behavior. Swarms of independent P. mirabilis isolates can recognize each other as foreign and establish a visible boundary where they meet; in contrast, genetically identical swarms merge. The ids genes, which encode self-identity proteins, are necessary but not sufficient for this territorial behavior. Here we have identified two new gene clusters: one (idr) encodes rhs-related products, and another (tss) encodes a putative type VI secretion (T6S) apparatus. The Ids and Idr proteins function independently of each other in extracellular transport and in territorial behaviors; however, these self-recognition systems are linked via this type VI secretion system. The T6S system is required for export of select Ids and Idr proteins. Our results provide a mechanistic and physiological basis for the fundamental behaviors of self-recognition and territoriality in a bacterial model system.

IMPORTANCE

Our results support a model in which self-recognition in P. mirabilis is achieved by the combined action of two independent pathways linked by a shared machinery for export of encoded self-recognition elements. These proteins together form a mechanistic network for self-recognition that can serve as a foundation for examining the prevalent biological phenomena of territorial behaviors and self-recognition in a simple, bacterial model system.

Application of a new informatics tool for contamination screening in the HIV sequencing laboratory

BACKGROUND

Current HIV-1 sequencing-based methods for detecting drug resistance-associated mutations are open and susceptible to contamination. Informatic identification of clinical sequences that are nearly identical to one another may indicate specimen-to-specimen contamination or another laboratory-associated issue.

OBJECTIVES

To design an informatic tool to rapidly identify potential contamination in the clinical laboratory using sequence analysis and to establish reference ranges for sequence variation in the HIV-1 protease and reverse transcriptase regions among a U.S. patient population.

STUDY DESIGN

We developed an open-source tool named HIV Contamination Detection (HIVCD). HIVCD was utilized to make pairwise comparisons of nearly 8000 partial HIV-1 pol gene sequences from patients across the United States and to calculate percent identities (PIDs) for each pair. ROC analysis and standard deviations of PID data were used to determine reference ranges for between-patient and within-patient comparisons and to guide selection of a threshold for identifying abnormally high PID between two unrelated sequences.

RESULTS

The PID reference range for between-patient comparisons ranged from 83.8 to 95.7% while within-patient comparisons ranged from 96 to 100%. Interestingly, 48% of between-patient sequence pairs with a PID>96.5 were geographically related. The selected threshold for abnormally high PIDs was 96 (AUC=0.993, sensitivity=0.980, specificity=0.999). During routine use, HIVCD identified a specimen mix-up and the source of contamination of a negative control.

CONCLUSIONS

In our experience, HIVCD is easily incorporated into laboratory workflow, useful for identifying potential laboratory errors, and contributes to quality testing. This type of analysis should be incorporated into routine laboratory practice.

HVP10 encoding V-PPase is a prime candidate for the barley HvNax3 sodium exclusion gene: evidence from fine mapping and expression analysis

In cereals, a common salinity tolerance mechanism is to limit accumulation of Na(+) in the shoot. In a cross between the barley variety Barque-73 (Hordeum vulgare ssp. vulgare) and the accession CPI-71284 of wild barley (H. vulgare ssp. spontaneum), the HvNax3 locus on chromosome 7H was found to determine a ~10-25 % difference in leaf Na(+) accumulation in seedlings grown in saline hydroponics, with the beneficial exclusion trait originating from the wild parent. The Na(+) exclusion allele was also associated with a 13-21 % increase in shoot fresh weight. The HvNax3 locus was delimited to a 0.4 cM genetic interval, where it cosegregated with the HVP10 gene for vacuolar H(+)-pyrophosphatase (V-PPase). Sequencing revealed that the mapping parents encoded identical HVP10 proteins, but salinity-induced mRNA expression of HVP10 was higher in CPI-71284 than in Barque-73, in both roots and shoots. By contrast, the expression of several other genes predicted by comparative mapping to be located in the HvNax3 interval was similar in the two parent lines. Previous work demonstrated roles for V-PPase in ion transport and salinity tolerance. We therefore considered transcription levels of HVP10 to be a possible basis for variation in shoot Na(+) accumulation and biomass production controlled by the HvNax3 locus under saline conditions. Potential mechanisms linking HVP10 expression patterns to the observed phenotypes are discussed.

Genome-wide identification and characterization of microRNA genes and their targets in flax (Linum usitatissimum): Characterization of flax miRNA genes

MicroRNAs (miRNAs) are small (20-24 nucleotide long) endogenous regulatory RNAs that play important roles in plant growth and development. They regulate gene expression at the post-transcriptional level by translational repression or target degradation and gene silencing. In this study, we identified 116 conserved miRNAs belonging to 23 families from the flax (Linum usitatissimum L.) genome using a computational approach. The precursor miRNAs varied in length; while most of the mature miRNAs were 21 nucleotide long, intergenic and showed conserved signatures of RNA polymerase II transcripts in their upstream regions. Promoter region analysis of the flax miRNA genes indicated prevalence of MYB transcription factor binding sites. Four miRNA gene clusters containing members of three phylogenetic groups were identified. Further, 142 target genes were predicted for these miRNAs and most of these represent transcriptional regulators. The miRNA encoding genes were expressed in diverse tissues as determined by digital expression analysis as well as real-time PCR. The expression of fourteen miRNAs and nine target genes was independently validated using the quantitative reverse transcription PCR (qRT-PCR). This study suggests that a large number of conserved plant miRNAs are also found in flax and these may play important roles in growth and development of flax.

PcarnBase: development of a transcriptomic database for the brain coral Platygyra carnosus

The aims of this study were to sequence the transcriptome and organize the sequence data into a searchable database for the brain coral Platygyra carnosus, a structure-forming dominant species along the coast of southern China. We collected healthy and tumorous coral tissues from two locations, extracted RNA from each tissue sample, pooled the RNA from all tissue samples, generated a cDNA library from the pooled samples, and conducted paired-end sequencing of the cDNA library using the Illumina platform to produce 59.6 M clean sequences with a read length of 90 bp. De novo assembly of the sequence data resulted in 162,468 unigenes with an average length of 606 bp (range, 201 to 23,923 bp). This is the largest transcriptome dataset for a species of coral whose genome has not been sequenced. A BLASTx search against the NCBI protein database showed that 55,355 of the unigenes matched at least a sequence with an E-value of < 0.00001; 59 % of the matched sequences are from Metazoa, 13 % are from Alveolata to which the symbiont Symbiodinium belongs, and 7 % are from bacteria. A database (PcarnBase) was constructed to provide easy access to the unigenes with attributes such as NCBI protein annotation, GO annotation, and KEGG pathway. It will facilitate functional genomic studies of P. carnosus, such as biomarker discovery for bleaching, tumor formation, and disease development at the gene or protein level, involvement of coral symbiotic algae in the host coral's stress responses, and genetic basis of stress resistance.

prfectBLAST: a platform-independent portable front end for the command terminal BLAST+ stand-alone suite

prfectBLAST is a multiplatform graphical user interface (GUI) for the stand-alone BLAST+ suite of applications. It allows researchers to do nucleotide or amino acid sequence similarity searches against public (or user-customized) databases that are locally stored. It does not require any dependencies or installation and can be used from a portable flash drive. prfectBLAST is implemented in Java version 6 (SUN) and runs on all platforms that support Java and for which National Center for Biotechnology Information has made available stand-alone BLAST executables, including MS Windows, Mac OS X, and Linux. It is free and open source software, made available under the GNU General Public License version 3 (GPLv3) and can be downloaded at www.cicy.mx/sitios/jramirez or http://code.google.com/p/prfectblast/.

OPTIMAS-DW: a comprehensive transcriptomics, metabolomics, ionomics, proteomics and phenomics data resource for maize

BACKGROUND

Maize is a major crop plant, grown for human and animal nutrition, as well as a renewable resource for bioenergy. When looking at the problems of limited fossil fuels, the growth of the world's population or the world's climate change, it is important to find ways to increase the yield and biomass of maize and to study how it reacts to specific abiotic and biotic stress situations. Within the OPTIMAS systems biology project maize plants were grown under a large set of controlled stress conditions, phenotypically characterised and plant material was harvested to analyse the effect of specific environmental conditions or developmental stages. Transcriptomic, metabolomic, ionomic and proteomic parameters were measured from the same plant material allowing the comparison of results across different omics domains. A data warehouse was developed to store experimental data as well as analysis results of the performed experiments.

DESCRIPTION

The OPTIMAS Data Warehouse (OPTIMAS-DW) is a comprehensive data collection for maize and integrates data from different data domains such as transcriptomics, metabolomics, ionomics, proteomics and phenomics. Within the OPTIMAS project, a 44K oligo chip was designed and annotated to describe the functions of the selected unigenes. Several treatment- and plant growth stage experiments were performed and measured data were filled into data templates and imported into the data warehouse by a Java based import tool. A web interface allows users to browse through all stored experiment data in OPTIMAS-DW including all data domains. Furthermore, the user can filter the data to extract information of particular interest. All data can be exported into different file formats for further data analysis and visualisation. The data analysis integrates data from different data domains and enables the user to find answers to different systems biology questions. Finally, maize specific pathway information is provided.

CONCLUSIONS

With OPTIMAS-DW a data warehouse for maize was established, which is able to handle different data domains, comprises several analysis results that will support researchers within their work and supports systems biological research in particular. The system is available at http://www.optimas-bioenergy.org/optimas_dw.

Genome-wide identification and characterization of nucleotide binding site leucine-rich repeat genes in linseed reveal distinct patterns of gene structure

Plants employ different disease-resistance genes to detect pathogens and to induce defense responses. The largest class of these genes encodes proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains. To identify the putative NBS-LRR encoding genes from linseed, we analyzed the recently published linseed genome sequence and identified 147 NBS-LRR genes. The NBS domain was used for phylogeny construction and these genes were classified into two well-known families, non-TIR (CNL) and TIR related (TNL), and formed eight clades in the neighbor-joining bootstrap tree. Eight different gene structures were observed among these genes. An unusual domain arrangement was observed in the TNL family members, predominantly in the TNL-5 clade members belonging to class D. About 12% of the genes observed were linseed specific. The study indicated that the linseed genes probably have an ancient origin with few progenitor genes. Quantitative expression analysis of five genes showed inducible expression. The in silico expression evidence was obtained for a few of these genes, and the expression was not correlated with the presence of any particular regulatory element or with unusual domain arrangement in those genes. This study will help in understanding the evolution of these genes, the development of disease resistant varieties, and the mechanism of disease resistance in linseed.

Proteome profiling of flax (Linum usitatissimum) seed: characterization of functional metabolic pathways operating during seed development

Flax (Linum usitatissimum L.) seeds are an important source of food and feed due to the presence of various health promoting compounds, making it a nutritionally and economically important plant. An in-depth analysis of the proteome of developing flax seed is expected to provide significant information with respect to the regulation and accumulation of such storage compounds. Therefore, a proteomic analysis of seven seed developmental stages (4, 8, 12, 16, 22, 30, and 48 days after anthesis) in a flax variety, NL-97 was carried out using a combination of 1D-SDS-PAGE and LC-MSE methods. A total 1716 proteins were identified and their functional annotation revealed that a majority of them were involved in primary metabolism, protein destination, storage and energy. Three carbon assimilatory pathways appeared to operate in flax seeds. Reverse transcription quantitative PCR of selected 19 genes was carried out to understand their roles during seed development. Besides storage proteins, methionine synthase, RuBisCO and S-adenosylmethionine synthetase were highly expressed transcripts, highlighting their importance in flax seed development. Further, the identified proteins were mapped onto developmental seed specific expressed sequence tag (EST) libraries of flax to obtain transcriptional evidence and 81% of them had detectable expression at the mRNA level. This study provides new insights into the complex seed developmental processes operating in flax.

Phylogenetic insights into regional HIV transmission

OBJECTIVES

Despite prevention efforts, new HIV diagnoses continue in the southern United States, where the epidemic is characterized by significant racial/ethnic disparities. We integrated phylogenetic analyses with clinical data to reveal trends in local HIV transmission.

DESIGN

Cross-sectional analysis of 1671 HIV-infected individuals each with one B-subtype pol sequence obtained during chronic (82%; UNC Center for AIDS Research Clinical Cohort) or acute/recent (18%; Duke/UNC Acute HIV Consortium) infection.

METHODS

Phylogenies were inferred using neighbor joining to select related sequences then confirmed with Bayesian methods. We characterized transmission clusters (clades n ≥ 3 sequences supported by posterior probabilities = 1) by factors including race/ethnicity and transmission risk. Factors associated with cluster membership were evaluated for newly diagnosed patients.

RESULTS

Overall, 72% were men, 59% black and 39% men who have sex with men (MSM). A total of 557 (33%) sequences grouped in either 108 pairs (n = 216) or 67 clusters (n = 341). Clusters ranged from three to 36 (median 4) members. Composition was delineated primarily by race, with 28% exclusively black, and to a lesser extent by risk group. Both MSM and heterosexuals formed discrete clusters, although substantial mixing was observed. In multivariable analysis, patients with age 30 years or less (P = 0.009), acute infection (P = 0.02), local residence (P = 0.002) and transmitted drug resistance (P = 0.02) were more likely to be cluster members, whereas Latinos were less likely (P < 0.001).

CONCLUSION

Integration of molecular, clinical and demographic data offers a unique view into the structure of local transmission networks. Clustering by black race, youth and transmitted drug resistance and inability to identify Latino clusters will inform prevention, testing and linkage to care strategies.

High-resolution genome-wide mapping of AHR and ARNT binding sites by ChIP-Seq

The aryl hydrocarbon receptor (AHR) and AHR nuclear translocator (ARNT) activated complex regulates genes in response to the environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). AHR has also emerged as a potential therapeutic target for the treatment of human diseases and different cancers, including breast cancer. To better understand AHR and ARNT signaling in breast cancer cells, we used chromatin immunoprecipitation linked to high-throughput sequencing to identify AHR- and ARNT-binding sites across the genome in TCDD-treated MCF-7 cells. We identified 2594 AHR-bound, 1352 ARNT-bound, and 882 AHR/ARNT cobound regions. No significant differences in the genomic distribution of AHR and ARNT were observed. Approximately 60% of the cobound regions contained at least one core an aryl hydrocarbon response element (AHRE), 5'-GCGTG-3'. AHR/ARNT peak density was the highest within 1 kb of transcription start sites (TSS); however, a number of AHR/ARNT cobound regions were located as far as 100 kb from TSS. De novo motif discovery identified a symmetrical variation of the AHRE (5'-GTGCGTG-3'), as well as FOXA1 and SP1 binding motifs. Microarray analysis identified 104 TCDD-responsive genes where 98 genes were upregulated by TCDD. Of the 104 regulated genes, 69 (66.3%) were associated with an AHR- or ARNT-bound region within 100 kb of their TSS. Overall our study identified AHR/ARNT cobound regions across the genome, revealed the importance but not absolute requirement for an AHRE in AHR/ARNT interactions with DNA, and identified a modified AHRE motif, thereby increasing our understanding of AHR/ARNT signaling pathway.

Induced mutations in the starch branching enzyme II (SBEII) genes increase amylose and resistant starch content in durum wheat

Starch is the largest component of the wheat (Triticum aestivum L.) grain and consists of approximately 70-80% amylopectin and 20-30% amylose. Amylopectin is a highly-branched, readily digested polysaccharide, whereas amylose has few branches and forms complexes that resist digestion and mimic dietary fiber (resistant starch). Down-regulation of the starch branching enzyme II (SBEII) gene by RNA interference (RNAi) was previously shown to increase amylose content in both hexaploid and tetraploid wheat. We generated ethyl methane sulphonate (EMS) mutants for the SBEIIa-A and SBEIIa-B homoeologs in the tetraploid durum wheat variety Kronos (T. turgidum ssp. durum L.). Single-gene mutants showed non-significant increases in amylose and resistant starch content, but a double mutant combining a SBEIIa-A knock-out mutation with a SBEIIa-B splice-site mutation showed a 22% increase in amylose content (P<0.0001) and a 115% increase in resistant starch content (P<0.0001). In addition, we obtained mutants for the A and B genome copies of the paralogous SBEIIb gene, mapped them 1-2 cM from SBEIIa, and generated double SBEIIa-SBEIIb mutants to study the effect of the SBEIIb gene in the absence of SBEIIa. These mutants are available to those interested in increasing amylose content and resistant starch in durum wheat.

The barley EST DNA Replication and Repair Database (bEST-DRRD) as a tool for the identification of the genes involved in DNA replication and repair

BACKGROUND

The high level of conservation of genes that regulate DNA replication and repair indicates that they may serve as a source of information on the origin and evolution of the species and makes them a reliable system for the identification of cross-species homologs. Studies that had been conducted to date shed light on the processes of DNA replication and repair in bacteria, yeast and mammals. However, there is still much to be learned about the process of DNA damage repair in plants.

DESCRIPTION

These studies, which were conducted mainly using bioinformatics tools, enabled the list of genes that participate in various pathways of DNA repair in Arabidopsis thaliana (L.) Heynh to be outlined; however, information regarding these mechanisms in crop plants is still very limited. A similar, functional approach is particularly difficult for a species whose complete genomic sequences are still unavailable. One of the solutions is to apply ESTs (Expressed Sequence Tags) as the basis for gene identification. For the construction of the barley EST DNA Replication and Repair Database (bEST-DRRD), presented here, the Arabidopsis nucleotide and protein sequences involved in DNA replication and repair were used to browse for and retrieve the deposited sequences, derived from four barley (Hordeum vulgare L.) sequence databases, including the "Barley Genome version 0.05" database (encompassing ca. 90% of barley coding sequences) and from two databases covering the complete genomes of two monocot models: Oryza sativa L. and Brachypodium distachyon L. in order to identify homologous genes. Sequences of the categorised Arabidopsis queries are used for browsing the repositories, which are located on the ViroBLAST platform. The bEST-DRRD is currently used in our project during the identification and validation of the barley genes involved in DNA repair.

CONCLUSIONS

The presented database provides information about the Arabidopsis genes involved in DNA replication and repair, their expression patterns and models of protein interactions. It was designed and established to provide an open-access tool for the identification of monocot homologs of known Arabidopsis genes that are responsible for DNA-related processes. The barley genes identified in the project are currently being analysed to validate their function.

RCDB: Renal Cancer Gene Database

BACKGROUND

Renal cell carcinoma or RCC is one of the common and most lethal urological cancers, with 40% of the patients succumbing to death because of metastatic progression of the disease. Treatment of metastatic RCC remains highly challenging because of its resistance to chemotherapy as well as radiotherapy, besides surgical resection. Whereas RCC comprises tumors with differing histological types, clear cell RCC remains the most common. A major problem in the clinical management of patients presenting with localized ccRCC is the inability to determine tumor aggressiveness and accurately predict the risk of metastasis following surgery. As a measure to improve the diagnosis and prognosis of RCC, researchers have identified several molecular markers through a number of techniques. However the wealth of information available is scattered in literature and not easily amenable to data-mining. To reduce this gap, this work describes a comprehensive repository called Renal Cancer Gene Database, as an integrated gateway to study renal cancer related data.

FINDINGS

Renal Cancer Gene Database is a manually curated compendium of 240 protein-coding and 269 miRNA genes contributing to the etiology and pathogenesis of various forms of renal cell carcinomas. The protein coding genes have been classified according to the kind of gene alteration observed in RCC. RCDB also includes the miRNAsdysregulated in RCC, along with the corresponding information regarding the type of RCC and/or metastatic or prognostic significance. While some of the miRNA genes showed an association with other types of cancers few were unique to RCC. Users can query the database using keywords, category and chromosomal location of the genes. The knowledgebase can be freely accessed via a user-friendly web interface at http://www.juit.ac.in/attachments/jsr/rcdb/homenew.html.

CONCLUSIONS

It is hoped that this database would serve as a useful complement to the existing public resources and as a good starting point for researchers and physicians interested in RCC genetics.

Enchytraeus albidus microarray: enrichment, design, annotation and database (EnchyBASE)

Enchytraeus albidus (Oligochaeta) is an ecologically relevant species used as standard test organisms for risk assessment. Effects of stressors in this species are commonly determined at the population level using reproduction and survival as endpoints. The assessment of transcriptomic responses can be very useful e.g. to understand underlying mechanisms of toxicity with gene expression fingerprinting. In the present paper the following is being addressed: 1) development of suppressive subtractive hybridization (SSH) libraries enriched for differentially expressed genes after metal and pesticide exposures; 2) sequencing and characterization of all generated cDNA inserts; 3) development of a publicly available genomic database on E. albidus. A total of 2100 Expressed Sequence Tags (ESTs) were isolated, sequenced and assembled into 1124 clusters (947 singletons and 177 contigs). From these sequences, 41% matched known proteins in GenBank (BLASTX, e-value ≤ 10(-5)) and 37% had at least one Gene Ontology (GO) term assigned. In total, 5.5% of the sequences were assigned to a metabolic pathway, based on KEGG. With this new sequencing information, an Agilent custom oligonucleotide microarray was designed, representing a potential tool for transcriptomic studies. EnchyBASE (http://bioinformatics.ua.pt/enchybase/) was developed as a web freely available database containing genomic information on E. albidus and will be further extended in the near future for other enchytraeid species. The database so far includes all ESTs generated for E. albidus from three cDNA libraries. This information can be downloaded and applied in functional genomics and transcription studies.

Profile of T cell recognition of HIV type 1 consensus group M Gag and Nef peptides in a clade A1- and D-infected Ugandan population

Reagents for evaluating non-clade B HIV-specific T cell responses are uncommon. Peptides based on highly conserved HIV-1 consensus group M sequences that are phylogenetically closer to most circulating strains may provide potential alternative reagents in populations with diverse infections, and may be relevant for vaccine design. Recognition of such reagents in clade A1-and D-infected populations has not been previously evaluated. Interferon (IFN)-γ ELISpot assay was used to evaluate T cell recognition of Gag and Nef peptides based on consensus group M sequences in 50 treatment-naive adults predominantly infected with HIV-1 clades A1 and D. Gag-induced T cell responses were correlated with gag sequence diversity. Infecting clades were determined from gag sequences for 45 of the 50 subjects as 40% clade A1 (18/45), 45% clade D (20/45), 2% clade C (1/45), 2% A1/C recombinant (1/45), 2% A1/D (1/45), 7% CRF10_CD (3/45), and 2% U (unclassifiable) (1/45). The mean genetic divergence and diversity of clade A and D gag region compared to group M consensus sequences at synonymous and nonsynonymous nucleotide and amino acid levels were not always significant. Gag peptides were targeted at significantly higher frequency [88% (44/50)] than Nef [64% (32/50)]; p=0.014, although their mean IFN-γ magnitudes were comparable ([3703 (95% CI 2567-4839)] vs. [2120 (95% CI 478-3762)]), respectively. Measurable virus-induced IFN-γ responses were detected in 96% (48/50) individuals, primarily targeting the more conserved Gag p24 and Nef central core regions. Use of these reagents to screen for HIV-specific IFN-γ responses may mitigate the challenge of viral diversity; although this targeting is apparently biased toward a few highly conserved epitopes.

Community annotation and bioinformatics workforce development in concert--Little Skate Genome Annotation Workshops and Jamborees

Recent advances in high-throughput DNA sequencing technologies have equipped biologists with a powerful new set of tools for advancing research goals. The resulting flood of sequence data has made it critically important to train the next generation of scientists to handle the inherent bioinformatic challenges. The North East Bioinformatics Collaborative (NEBC) is undertaking the genome sequencing and annotation of the little skate (Leucoraja erinacea) to promote advancement of bioinformatics infrastructure in our region, with an emphasis on practical education to create a critical mass of informatically savvy life scientists. In support of the Little Skate Genome Project, the NEBC members have developed several annotation workshops and jamborees to provide training in genome sequencing, annotation and analysis. Acting as a nexus for both curation activities and dissemination of project data, a project web portal, SkateBase (http://skatebase.org) has been developed. As a case study to illustrate effective coupling of community annotation with workforce development, we report the results of the Mitochondrial Genome Annotation Jamborees organized to annotate the first completely assembled element of the Little Skate Genome Project, as a culminating experience for participants from our three prior annotation workshops. We are applying the physical/virtual infrastructure and lessons learned from these activities to enhance and streamline the genome annotation workflow, as we look toward our continuing efforts for larger-scale functional and structural community annotation of the L. erinacea genome.

Viral linkage in HIV-1 seroconverters and their partners in an HIV-1 prevention clinical trial

BACKGROUND

Characterization of viruses in HIV-1 transmission pairs will help identify biological determinants of infectiousness and evaluate candidate interventions to reduce transmission. Although HIV-1 sequencing is frequently used to substantiate linkage between newly HIV-1 infected individuals and their sexual partners in epidemiologic and forensic studies, viral sequencing is seldom applied in HIV-1 prevention trials. The Partners in Prevention HSV/HIV Transmission Study (ClinicalTrials.gov #NCT00194519) was a prospective randomized placebo-controlled trial that enrolled serodiscordant heterosexual couples to determine the efficacy of genital herpes suppression in reducing HIV-1 transmission; as part of the study analysis, HIV-1 sequences were examined for genetic linkage between seroconverters and their enrolled partners.

METHODOLOGY/PRINCIPAL FINDINGS

We obtained partial consensus HIV-1 env and gag sequences from blood plasma for 151 transmission pairs and performed deep sequencing of env in some cases. We analyzed sequences with phylogenetic techniques and developed a Bayesian algorithm to evaluate the probability of linkage. For linkage, we required monophyletic clustering between enrolled partners' sequences and a Bayesian posterior probability of ≥ 50%. Adjudicators classified each seroconversion, finding 108 (71.5%) linked, 40 (26.5%) unlinked, and 3 (2.0%) indeterminate transmissions, with linkage determined by consensus env sequencing in 91 (84%). Male seroconverters had a higher frequency of unlinked transmissions than female seroconverters. The likelihood of transmission from the enrolled partner was related to time on study, with increasing numbers of unlinked transmissions occurring after longer observation periods. Finally, baseline viral load was found to be significantly higher among linked transmitters.

CONCLUSIONS/SIGNIFICANCE

In this first use of HIV-1 sequencing to establish endpoints in a large clinical trial, more than one-fourth of transmissions were unlinked to the enrolled partner, illustrating the relevance of these methods in the design of future HIV-1 prevention trials in serodiscordant couples. A hierarchy of sequencing techniques, analysis methods, and expert adjudication contributed to the linkage determination process.

Comment on "The origins of sexually transmitted HIV among men who have sex with men"

Whether HIV from seminal cells or free HIV in semen is the origin of transmitted virus has important implications for the design of transmission prevention strategies. We found that a recent claim that HIV originates from seminal plasma and not from seminal cells was erroneous, because it was based on biological specimens that had been mislabeled, mixed-up, or contaminated. The origin of transmitted virus from semen therefore remains an open question.

Mutation of HIV-1 genomes in a clinical population treated with the mutagenic nucleoside KP1461

The deoxycytidine analog KP1212, and its prodrug KP1461, are prototypes of a new class of antiretroviral drugs designed to increase viral mutation rates, with the goal of eventually causing the collapse of the viral population. Here we present an extensive analysis of viral sequences from HIV-1 infected volunteers from the first “mechanism validation” phase II clinical trial of a mutagenic base analog in which individuals previously treated with antiviral drugs received 1600 mg of KP1461 twice per day for 124 days. Plasma viral loads were not reduced, and overall levels of viral mutation were not increased during this short-term study, however, the mutation spectrum of HIV was altered. A large number (N = 105 per sample) of sequences were analyzed, each derived from individual HIV-1 RNA templates, after 0, 56 and 124 days of therapy from 10 treated and 10 untreated control individuals (>7.1 million base pairs of unique viral templates were sequenced). We found that private mutations, those not found in more than one viral sequence and likely to have occurred in the most recent rounds of replication, increased in treated individuals relative to controls after 56 (p = 0.038) and 124 (p = 0.002) days of drug treatment. The spectrum of mutations observed in the treated group showed an excess of A to G and G to A mutations (p = 0.01), and to a lesser extent T to C and C to T mutations (p = 0.09), as predicted by the mechanism of action of the drug. These results validate the proposed mechanism of action in humans and should spur development of this novel antiretroviral approach.

MULTBLAST: A web application for multiple BLAST searches

Basic Local Alignment Search Tool, (BLAST) allows the comparison of a query sequence/s to a database of sequences and identifies those sequences that are similar to the query above a user-defined threshold. We have developed a user friendly web application, MULTBLAST that runs a series of BLAST searches on a user-supplied list of proteins against one or more target protein or nucleotide databases. The application pre-processes the data, launches each individual BLAST search on the University of Nevada, Reno's-TimeLogic DeCypher® system (available from Active Motif, Inc.) and retrieves and combines all the results into a simple, easy to read output file. The output file presents the list of the query proteins, followed by the BLAST results for the matching sequences from each target database in consecutive columns. This format is especially useful for either comparing the results from the different target databases, or analyzing the results while keeping the identification of each target database separate.

YeastWeb: a workset-centric web resource for gene family analysis in yeast

Background

Currently, a number of yeast genomes with different physiological features have been sequenced and annotated, which provides invaluable information to investigate yeast genetics, evolutionary mechanism, structure and function of gene families.

Description

YeastWeb is a novel database created to provide access to gene families derived from the available yeast genomes by assigning the genes into putative families. It has many useful features that complement existing databases, such as SGD, CYGD and Génolevures: 1) Detailed computational annotation was conducted with each entry with InterProScan, EMBOSS and functional/pathway databases, such as GO, COG and KEGG; 2) A well established user-friendly environment was created to allow users to retrieve the annotated genes and gene families using functional classification browser, keyword search or similarity-based search; 3) Workset offers users many powerful functions to manage the retrieved data efficiently, associate the individual items easily and save the intermediate results conveniently; 4) A series of comparative genomics and molecular evolution analysis tools are neatly implemented to allow users to view multiple sequence alignments and phylogenetic tree of gene families. At present, YeastWeb holds the gene families clustered from various MCL inflation values from a total of 13 available yeast genomes.

Conclusions

Given the great interest in yeast research, YeastWeb has the potential to become a useful resource for the scientific community of yeast biologists and related researchers investigating the evolutionary relationship of yeast gene families. YeastWeb is available at http://centre.bioinformatics.zj.cn/Yeast/.

Restriction of HIV-1 genotypes in breast milk does not account for the population transmission genetic bottleneck that occurs following transmission

Background

Breast milk transmission of HIV-1 remains a major route of pediatric infection. Defining the characteristics of viral variants to which breastfeeding infants are exposed is important for understanding the genetic bottleneck that occurs in the majority of mother-to-child transmissions. The blood-milk epithelial barrier markedly restricts the quantity of HIV-1 in breast milk, even in the absence of antiretroviral drugs. The basis of this restriction and the genetic relationship between breast milk and blood variants are not well established.

Methodology/Principal Findings

We compared 356 HIV-1 subtype C gp160 envelope (env) gene sequences from the plasma and breast milk of 13 breastfeeding women. A trend towards lower viral population diversity and divergence in breast milk was observed, potentially indicative of clonal expansion within the breast. No differences in potential N-linked glycosylation site numbers or in gp160 variable loop amino acid lengths were identified. Genetic compartmentalization was evident in only one out of six subjects in whom contemporaneously obtained samples were studied. However, in samples that were collected 10 or more days apart, six of seven subjects were classified as having compartmentalized viral populations, highlighting the necessity of contemporaneous sampling for genetic compartmentalization studies. We found evidence of CXCR4 co-receptor using viruses in breast milk and blood in nine out of the thirteen subjects, but no evidence of preferential localization of these variants in either tissue.

Conclusions/Significance

Despite marked restriction of HIV-1 quantities in milk, our data indicate intermixing of virus between blood and breast milk. Thus, we found no evidence that a restriction in viral genotype diversity in breast milk accounts for the genetic bottleneck observed following transmission. In addition, our results highlight the rapidity of HIV-1 env evolution and the importance of sample timing in analyses of gene flow.

4E10-resistant HIV-1 isolated from four subjects with rare membrane-proximal external region polymorphisms

Human antibody 4E10 targets the highly conserved membrane-proximal external region (MPER) of the HIV-1 transmembrane glycoprotein, gp41, and has extraordinarily broad neutralizing activity. It is considered by many to be a prototype for vaccine development. In this study, we describe four subjects infected with viruses carrying rare MPER polymorphisms associated with resistance to 4E10 neutralization. In one case resistant virus carrying a W680G substitution was transmitted from mother to infant. We used site-directed mutagenesis to demonstrate that the W680G substitution is necessary for conferring the 4E10-resistant phenotype, but that it is not sufficient to transfer the phenotype to a 4E10-sensitive Env. Our third subject carried Envs with a W680R substitution causing variable resistance to 4E10, indicating that residues outside the MPER are required to confer the phenotype. A fourth subject possessed a F673L substitution previously associated with 4E10 resistance. For all three subjects with W680 polymorphisms, we observed additional residues in the MPER that co-varied with position 680 and preserved charged distributions across this region. Our data provide important caveats for vaccine development targeting the MPER. Naturally occurring Env variants described in our study also represent unique tools for probing the structure-function of HIV-1 envelope.

LipocalinPred: a SVM-based method for prediction of lipocalins

BACKGROUND

Functional annotation of rapidly amassing nucleotide and protein sequences presents a challenging task for modern bioinformatics. This is particularly true for protein families sharing extremely low sequence identity, as for lipocalins, a family of proteins with varied functions and great diversity at the sequence level, yet conserved structures.

RESULTS

In the present study we propose a SVM based method for identification of lipocalin protein sequences. The SVM models were trained with the input features generated using amino acid, dipeptide and secondary structure compositions as well as PSSM profiles. The model derived using both PSSM and secondary structure emerged as the best model in the study. Apart from achieving a high prediction accuracy (>90% in leave-one-out), lipocalinpred correctly differentiates closely related fatty acid-binding proteins and triabins as non-lipocalins.

CONCLUSION

The method offers a promising approach as a lipocalin prediction tool, complementing PROSITE, Pfam and homology modelling methods.

Evidence for limited genetic compartmentalization of HIV-1 between lung and blood

Background

HIV-1 is frequently detected in the lungs of infected individuals and is likely important in the development of pulmonary opportunistic infections. The unique environment of the lung, rich in alveolar macrophages and with specialized local immune responses, may contribute to differential evolution or selection of HIV-1.

Methodology and Findings

We characterized HIV-1 in the lung in relation to contemporaneous viral populations in the blood. The C2-V5 region of HIV-1 env was sequenced from paired lung (induced sputum or bronchoalveolar lavage) and blood (plasma RNA and proviral DNA from sorted or unsorted PBMC) from 18 subjects. Compartmentalization between tissue pairs was assessed using 5 established tree or distance-based methods, including permutation tests to determine statistical significance. We found statistical evidence of compartmentalization between lung and blood in 10/18 subjects, although lung and blood sequences were intermingled on phylogenetic trees in all subjects. The subject showing the greatest compartmentalization contained many nearly identical sequences in BAL sample, suggesting clonal expansion may contribute to reduced viral diversity in the lung in some cases. However, HIV-1 sequences in lung were not more homogeneous overall, nor were we able to find a lung-specific genotype associated with macrophage tropism in V3. In all four subjects in whom predicted X4 genotypes were found in blood, predicted X4 genotypes were also found in lung.

Conclusions

Our results support a picture of continuous migration of HIV-1 between circulating blood and lung tissue, with perhaps a very limited degree of localized evolution or clonal replication.

Immunoinformatics: current trends and future directions

Immunoinformatics has recently emerged as a critical field for accelerating immunology research. Although still an evolving process, computational models now play instrumental roles, not only in directing the selection of key experiments, but also in the formulation of new testable hypotheses through detailed analysis of complex immunologic data that could not be achieved using traditional approaches alone. Immunomics, which combines traditional immunology with computer science, mathematics, chemistry, biochemistry, genomics and proteomics for the large-scale analysis of immune system function, offers new opportunities for future bench-to-bedside research. In this article, we review the latest trends and future directions of the field.

The Populus Genome Integrative Explorer (PopGenIE): a new resource for exploring the Populus genome

Populus has become an important model plant system. However, utilization of the increasingly extensive collection of genetics and genomics data created by the community is currently hindered by the lack of a central resource, such as a model organism database (MOD). Such MODs offer a single entry point to the collection of resources available within a model system, typically including tools for exploring and querying those resources. As a starting point to overcoming the lack of such an MOD for Populus, we present the Populus Genome Integrative Explorer (PopGenIE), an integrated set of tools for exploring the Populus genome and transcriptome. The resource includes genome, synteny and quantitative trait locus (QTL) browsers for exploring genetic data. Expression tools include an electronic fluorescent pictograph (eFP) browser, expression profile plots, co-regulation within collated transcriptomics data sets, and identification of over-represented functional categories and genomic hotspot locations. A number of collated transcriptomics data sets are made available in the eFP browser to facilitate functional exploration of gene function. Additional homology and data extraction tools are provided. PopGenIE significantly increases accessibility to Populus genomics resources and allows exploration of transcriptomics data without the need to learn or understand complex statistical analysis methods. PopGenIE is available at www.popgenie.org or via www.populusgenome.info.

FlyPhy: a phylogenomic analysis platform for Drosophila genes and gene families

Background

The availability of 12 fully sequenced Drosophila species genomes provides an excellent opportunity to explore the evolutionary mechanism, structure and function of gene families in Drosophila. Currently, several important resources, such as FlyBase, FlyMine and DroSpeGe, have been devoted to integrating genetic, genomic, and functional data of Drosophila into a well-organized form. However, all of these resources are gene-centric and lack the information of the gene families in Drosophila.

Description

FlyPhy is a comprehensive phylogenomic analysis platform devoted to analyzing the genes and gene families in Drosophila. Genes were classified into families using a graph-based Markov Clustering algorithm and extensively annotated by a number of bioinformatic tools, such as basic sequence features, functional category, gene ontology terms, domain organization and sequence homolog to other databases. FlyPhy provides a simple and user-friendly web interface to allow users to browse and retrieve the information at multiple levels. An outstanding feature of the FlyPhy is that all the retrieved results can be added to a workset for further data manipulation. For the data stored in the workset, multiple sequence alignment, phylogenetic tree construction and visualization can be easily performed to investigate the sequence variation of each given family and to explore its evolutionary mechanism.

Conclusion

With the above functionalities, FlyPhy will be a useful resource and convenient platform for the Drosophila research community. The FlyPhy is available at .

PlasmoGF: an integrated system for comparative genomics and phylogenetic analysis of Plasmodium gene families

Malaria, one of the world's most common diseases, is caused by the intracellular protozoan parasite known as Plasmodium. Recently, with the arrival of several malaria parasite genomes, we established an integrated system named PlasmoGF for comparative genomics and phylogenetic analysis of Plasmodium gene families. Gene families were clustered using the Markov Cluster algorithm implemented in TribeMCL program and could be searched using keywords, gene-family information, domain composition, Gene Ontology and BLAST. Moreover, a number of useful bioinformatics tools were implemented to facilitate the analysis of these putative Plasmodium gene families, including gene retrieval, annotation, sequence alignment, phylogeny construction and visualization. In the current version, PlasmoGF contained 8980 sets of gene families derived from six malaria parasite genomes: Plasmodium. falciparum, P. berghei, P. knowlesi, P. chabaudi, P. vivax and P. yoelii. The availability of such a highly integrated system would be of great interest for the community of researchers working on malaria parasite phylogenomics.

AVAILABILITY

PlasmoGF is freely available at http://bioinformatics.zj.cn/pgf/

The Personal Sequence Database: a suite of tools to create and maintain web-accessible sequence databases

Background

Large molecular sequence databases are fundamental resources for modern bioscientists. Whether for project-specific purposes or sharing data with colleagues, it is often advantageous to maintain smaller sequence databases. However, this is usually not an easy task for the average bench scientist.

Results

We present the Personal Sequence Database (PSD), a suite of tools to create and maintain small- to medium-sized web-accessible sequence databases. All interactions with PSD tools occur via the internet with a web browser. Users may define sequence groups within their database that can be maintained privately or published to the web for public use. A sequence group can be downloaded, browsed, searched by keyword or searched for sequence similarities using BLAST. Publishing a sequence group extends these capabilities to colleagues and collaborators. In addition to being able to manage their own sequence databases, users can enroll sequences in BLASTAgent, a BLAST hit tracking system, to monitor NCBI databases for new entries displaying a specified level of nucleotide or amino acid similarity.

Conclusion

The PSD offers a valuable set of resources unavailable elsewhere. In addition to managing sequence data and BLAST search results, it facilitates data sharing with colleagues, collaborators and public users. The PSD is hosted by the authors and is available at .

Recognition of HIV-1 peptides by host CTL is related to HIV-1 similarity to human proteins

Background

While human immunodeficiency virus type 1 (HIV-1)-specific cytotoxic T lymphocytes preferentially target specific regions of the viral proteome, HIV-1 features that contribute to immune recognition are not well understood. One hypothesis is that similarities between HIV and human proteins influence the host immune response, i.e., resemblance between viral and host peptides could preclude reactivity against certain HIV epitopes.

Methodology/Principal Findings

We analyzed the extent of similarity between HIV-1 and the human proteome. Proteins from the HIV-1 B consensus sequence from 2001 were dissected into overlapping k-mers, which were then probed against a non-redundant database of the human proteome in order to identify segments of high similarity. We tested the relationship between HIV-1 similarity to host encoded peptides and immune recognition in HIV-infected individuals, and found that HIV immunogenicity could be partially modulated by the sequence similarity to the host proteome. ELISpot responses to peptides spanning the entire viral proteome evaluated in 314 individuals showed a trend indicating an inverse relationship between the similarity to the host proteome and the frequency of recognition. In addition, analysis of responses by a group of 30 HIV-infected individuals against 944 overlapping peptides representing a broad range of individual HIV-1B Nef variants, affirmed that the degree of similarity to the host was significantly lower for peptides with reactive epitopes than for those that were not recognized.

Conclusions/Significance

Our results suggest that antigenic motifs that are scarcely represented in human proteins might represent more immunogenic CTL targets not selected against in the host. This observation could provide guidance in the design of more effective HIV immunogens, as sequences devoid of host-like features might afford superior immune reactivity.