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

Mapping adult plant stem rust resistance in barley accessions Hietpas-5 and GAW-79

Key message Major stem rust resistance QTLs proposed to be Rpg2 from Hietpas-5 and Rpg3 from GAW-79 were identified in chromosomes 2H and 5H, respectively, and will enhance the diversity of stem rust resistance in barley improvement programs. Stem rust is a devastating disease of cereal crops worldwide. In barley (Hordeum vulgare ssp. vulgare), the disease is caused by two pathogens: Puccinia graminis f. sp. secalis (Pgs) and Puccinia graminis f. sp. tritici (Pgt). In North America, the stem rust resistance gene Rpg1 has protected barley from serious losses for more than 60 years; however, widely virulent Pgt races from Africa in the Ug99 group threaten the crop. The accessions Hietpas-5 (CIho 7124) and GAW-79 (PI 382313) both possess moderate-to-high levels of adult plant resistance to stem rust and are the sources of the resistance genes Rpg2 and Rpg3, respectively. To identify quantitative trait loci (QTL) for stem rust resistance in Hietpas-5 and GAW-79, two biparental populations were developed with Hiproly (PI 60693), a stem rust-susceptible accession. Both populations were phenotyped to the North American Pgt races of MCCFC, QCCJB, and HKHJC in St. Paul, Minnesota, and to African Pgt races (predominately TTKSK in the Ug99 group) in Njoro, Kenya. In the Hietpas-5/Hiproly population, a major effect QTL was identified in chromosome 2H, which is proposed as the location for Rpg2. In the GAW-79/Hiproly population, a major effect QTL was identified in chromosome 5H and is the proposed location for Rpg3. These QTLs will enhance the diversity of stem rust resistance in barley improvement programs.

PGD: Pineapple Genomics Database

Pineapple occupies an important phylogenetic position as its reference genome is a model for studying the evolution the Bromeliaceae family and the crassulacean acid metabolism (CAM) photosynthesis. Here, we developed a pineapple genomics database (PGD, http://pineapple.angiosperms.org/pineapple/html/index.html) as a central online platform for storing and integrating genomic, transcriptomic, function annotation and genetic marker data for pineapple (Ananas comosus (L.) Merr.). The PGD currently hosts significant search tools and available datasets for researchers to study comparative genomics, gene expression, gene co-expression molecular marker, and gene annotation of A. comosus (L). PGD also performed a series of additional pages for a genomic browser that visualizes genomic data interactively, bulk data download, a detailed user manual, and data integration information. PGD was developed with the capacity to integrate future data resources, and will be used as a long-term and open access database to facilitate the study of the biology, distribution, and the evolution of pineapple and the relative plant species. An email-based helpdesk is also available to offer support with the website and requests of specific datasets from the research community.

What lies beneath? Fungal diversity at the bottom of Lake Michigan and Lake Superior

Fungi are phylogenetically diverse organisms found in nearly every environment as key contributors to the processes of nutrient cycling and decomposition. To date, most fungal diversity has been documented from terrestrial habitats leaving aquatic habitats underexplored. In particular, comparatively little is known about fungi inhabiting freshwater lakes, particularly the benthic zone, which may serve as an untapped resource for fungal biodiversity. Advances in technology allowing for direct sequencing of DNA from environmental samples provide a new opportunity to investigate freshwater benthic fungi. In this study, we employed both culture-dependent and culture-independent methods to evaluate the diversity of fungi in one of the largest freshwater systems on Earth, the North American Laurentian Great Lakes. This study presents the first comprehensive survey of fungi from sediment from Lake Michigan and Lake Superior, resulting in 465 fungal taxa with only 7% of sequence overlap between these two methods. Additionally, culture-independent analyses of the ITS1 and ITS2 regions revealed 49% and 72%, respectively, of the OTUs did not match a described fungal taxonomic group below kingdom Fungi. The low level of sequence overlap between methods and high percentage of fungal taxa that can only be classified at the kingdom level suggests an immense amount of fungal diversity remains to be studied in these aquatic fungal communities.

AmpuBase: a transcriptome database for eight species of apple snails (Gastropoda: Ampullariidae)

BACKGROUND

Gastropoda, with approximately 80,000 living species, is the largest class of Mollusca. Among gastropods, apple snails (family Ampullariidae) are globally distributed in tropical and subtropical freshwater ecosystems and many species are ecologically and economically important. Ampullariids exhibit various morphological and physiological adaptations to their respective habitats, which make them ideal candidates for studying adaptation, population divergence, speciation, and larger-scale patterns of diversity, including the biogeography of native and invasive populations. The limited availability of genomic data, however, hinders in-depth ecological and evolutionary studies of these non-model organisms.

RESULTS

Using Illumina Hiseq platforms, we sequenced 1220 million reads for seven species of apple snails. Together with the previously published RNA-Seq data of two apple snails, we conducted de novo transcriptome assembly of eight species that belong to five genera of Ampullariidae, two of which represent Old World lineages and the other three New World lineages. There were 20,730 to 35,828 unigenes with predicted open reading frames for the eight species, with N50 (shortest sequence length at 50% of the unigenes) ranging from 1320 to 1803 bp. 69.7% to 80.2% of these unigenes were functionally annotated by searching against NCBI's non-redundant, Gene Ontology database and the Kyoto Encyclopaedia of Genes and Genomes. With these data we developed AmpuBase, a relational database that features online BLAST functionality for DNA/protein sequences, keyword searching for unigenes/functional terms, and download functions for sequences and whole transcriptomes.

CONCLUSIONS

In summary, we have generated comprehensive transcriptome data for multiple ampullariid genera and species, and created a publicly accessible database with a user-friendly interface to facilitate future basic and applied studies on ampullariids, and comparative molecular studies with other invertebrates.

Identification and Validation of a New Source of Low Grain Cadmium Accumulation in Durum Wheat

Cadmium (Cd) is a heavy metal that has no known biological function and is toxic for many living organisms. The maximum level of Cd concentration allowed in the international market for wheat grain is 0.2 mg kg⁻¹. Because phenotyping for Cd uptake is expensive and time consuming, molecular markers associated with genes conferring low Cd uptake would expedite selection and lead to the development of durum cultivars with reduced Cd concentrations. Here, we identified single nucleotide polymorphisms (SNPs) associated with a novel low Cd uptake locus in the durum experimental line D041735, which has hexaploid common wheat in its pedigree. Genetic analysis revealed a single major QTL for Cd uptake on chromosome arm 5BL within a 0.3 cM interval flanked by SNP markers. Analysis of the intervening sequence revealed a gene with homology to an aluminum-induced protein as a candidate gene. Validation and allelism tests revealed that the low Cd uptake gene identified in this study is different from the closely linked Cdu1-B gene, which also resides on 5BL. This study therefore showed that the durum experimental line D041735 contains a novel low Cd uptake gene that was likely acquired from hexaploid wheat.

Functional regulation of Q by microRNA172 and transcriptional co-repressor TOPLESS in controlling bread wheat spikelet density

Bread wheat (Triticum aestivum) spike architecture is an important agronomic trait. The Q gene plays a key role in the domestication of bread wheat spike architecture. However, the regulatory mechanisms of Q expression and transcriptional activity remain largely unknown. In this study, we show that overexpression of bread wheat tae-miR172 caused a speltoid-like spike phenotype, reminiscent of that in wheat plants with the q gene. The reduction in Q transcript levels in the tae-miR172 overexpression transgenic bread wheat lines suggests that the Q expression can be suppressed by tae-miR172 in bread wheat. Indeed, our RACE analyses confirmed that the Q mRNA is targeted by tae-miR172 for cleavage. According to our analyses, the Q protein is localized in nucleus and confers transcriptional repression activity. Meanwhile, the Q protein could physically interact with the bread wheat transcriptional co-repressor TOPLESS (TaTPL). Specifically, the N-terminal ethylene-responsive element binding factor-associated amphiphilic repression (EAR) (LDLNVE) motif but not the C-terminal EAR (LDLDLR) motif of Q protein mediates its interaction with the CTLH motif of TaTPL. Moreover, we show that the N-terminal EAR motif of Q protein is also essentially required for the transcriptional repression activity of Q protein. Taken together, we reveal the functional regulation of Q protein by tae-miR172 and transcriptional co-repressor TaTPL in controlling the bread wheat spike architecture.

ASGDB: a specialised genomic resource for interpreting Anopheles sinensis insecticide resistance

BACKGROUND

Anopheles sinensis is an important malaria vector in Southeast Asia. The widespread emergence of insecticide resistance in this mosquito species poses a serious threat to the efficacy of malaria control measures, particularly in China. Recently, the whole-genome sequencing and de novo assembly of An. sinensis (China strain) has been finished. A series of insecticide-resistant studies in An. sinensis have also been reported. There is a growing need to integrate these valuable data to provide a comprehensive database for further studies on insecticide-resistant management of An. sinensis.

RESULTS

A bioinformatics database named An. sinensis genome database (ASGDB) was built. In addition to being a searchable database of published An. sinensis genome sequences and annotation, ASGDB provides in-depth analytical platforms for further understanding of the genomic and genetic data, including visualization of genomic data, orthologous relationship analysis, GO analysis, pathway analysis, expression analysis and resistance-related gene analysis. Moreover, ASGDB provides a panoramic view of insecticide resistance studies in An. sinensis in China. In total, 551 insecticide-resistant phenotypic and genotypic reports on An. sinensis distributed in Chinese malaria-endemic areas since the mid-1980s have been collected, manually edited in the same format and integrated into OpenLayers map-based interface, which allows the international community to assess and exploit the high volume of scattered data much easier. The database has been given the URL: http://www.asgdb.org /.

CONCLUSIONS

ASGDB was built to help users mine data from the genome sequence of An. sinensis easily and effectively, especially with its advantages in insecticide resistance surveillance and control.

Association mapping of drought tolerance-related traits in barley to complement a traditional biparental QTL mapping study

Association mapping of drought-related traits in barley was used to increase the density of existing QTL maps without recreating mapping populations. We used 109 spring barley genotypes exhibiting high or low drought tolerance to elucidate the associations between diversity array technology sequencing (DArTseq) and single nucleotide polymorphism (SNP) markers and various physiological parameters related to plant responses to drought conditions. The study was performed in controlled conditions (growth chambers), drought tolerance was phenotyped in the four-leaf seedlings. We identified 58 associations including 34 new markers (i.e., 16 DArTseq and 18 SNP markers). The results for three markers were consistent with the data obtained in an earlier traditional biparental QTL mapping study. The regions neighboring markers on linkage group 2H contained the highest number of significant marker-trait associations. Five markers related to the photosynthetic activity of photosystem II were detected on chromosome 4H. The lowest number of associations were observed for the sequences neighboring DArT markers on linkage group 6H. A chromosome 3H region related to water use efficiency and net photosynthesis rate in both biparental QTL, and association study, may be particularly valuable, as these parameters correspond to the ability of plants to remain highly productive under water deficit stress. Our findings confirm that association mapping can increase the density of existing QTL maps without recreating mapping populations.

Genome sequence and comparative analysis of Jiangella alba YIM 61503T isolated from a medicinal plant Maytenus austroyunnanensis

A draft genome sequence of Jiangella alba YIM 61503^T revealed a genome size of 7,664,864 bp arranged in 33 scaffolds. The genome was predicted to contain 7196 predicted genes, including 51 coding for RNA. Phylogenetic and comparative analyses of the draft genome of J. alba YIM 61503^T with the available genomes of other Jiangella species suggested a proximal similarity between strains J. alba YIM 61503^T and J. muralis DSM 45357^T, while indicating a high divergence between J. gansuensis YIM 002^T and other Jiangella species. The genome of J. alba YIM 61503^T also revealed genes involved in indole-3-acetic acid biosynthesis and an alkylresorcinols gene cluster. Further, detection of phosphotransferase genes in the genome of all Jiangella species indicated that they can uptake and phosphorylate sugars. The presences of TreX-Z, TreS and OtsA-OtsB genes in some of the Jiangella strains also indicated a possible mechanism for their tolerance of high salinity. Besides providing new insights into its genetic features, our results suggested that J. alba YIM 61503^T could be a potential strain for further genome mining studies. The release of this genome may, therefore, provide a better prospect for understanding "evolutionary taxonomy" about this genus in future.

Identification and characterization of regulatory network components for anthocyanin synthesis in barley aleurone

BACKGROUND

Among natural populations, there are different colours of barley (Hordeum vulgare L.). The colour of barley grains is directly related to the accumulation of different pigments in the aleurone layer, pericarp and lemma. Blue grain colour is due to the accumulation of anthocyanins in the aleurone layer, which is dependent on the presence of five Blx genes that are not sequenced yet (Blx1, Blx3 and Blx4 genes clustering on chromosome 4HL and Blx2 and Blx5 on 7HL). Due to the health benefits of anthocyanins, blue-grained barley can be considered as a source of dietary food. The goal of the current study was to identify and characterize components of the anthocyanin synthesis regulatory network for the aleurone layer in barley.

RESULTS

The candidate genes for components of the regulatory complex MBW (consisting of transcription factors MYB, bHLH/MYC and WD40) for anthocyanin synthesis in barley aleurone were identified. These genes were designated HvMyc2 (4HL), HvMpc2 (4HL), and HvWD40 (6HL). HvMyc2 was expressed in aleurone cells only. A loss-of-function (frame shift) mutation in HvMyc2 of non-coloured compared to blue-grained barley was revealed. Unlike aleurone-specific HvMyc2, the HvMpc2 gene was expressed in different tissues; however, its activity was not detected in non-coloured aleurone in contrast to a coloured aleurone, and allele-specific mutations in its promoter region were found. The single-copy gene HvWD40, which encodes the required component of the regulatory MBW complex, was expressed constantly in coloured and non-coloured tissues and had no allelic differences. HvMyc2 and HvMpc2 were genetically mapped using allele-specific developed CAPS markers developed. HvMyc2 was mapped in position between SSR loci XGBS0875-4H (3.4 cM distal) and XGBM1048-4H (3.4 cM proximal) matching the region chromosome 4HL where the Blx-cluster was found. In this position, one of the anthocyanin biosynthesis structural genes (HvF3'5'H) was also mapped using an allele-specific CAPS-marker developed in the current study.

CONCLUSIONS

The genes involved in anthocyanin synthesis in the barley aleurone layer were identified and characterized, including components of the regulatory complex MBW, from which the MYC-encoding gene (HvMyc2) appeared to be the main factor underlying variation of barley by aleurone colour.

Patterns and rates of viral evolution in HIV-1 subtype B infected females and males

Biological sex differences affect the course of HIV infection, with untreated women having lower viral loads compared to their male counterparts but, for a given viral load, women have a higher rate of progression to AIDS. However, the vast majority of data on viral evolution, a process that is clearly impacted by host immunity and could be impacted by sex differences, has been derived from men. We conducted an intensive analysis of HIV-1 gag and env-gp120 evolution taken over the first 6–11 years of infection from 8 Women’s Interagency HIV Study (WIHS) participants who had not received combination antiretroviral therapy (ART). This was compared to similar data previously collected from men, with both groups infected with HIV-1 subtype B. Early virus populations in men and women were generally homogenous with no differences in diversity between sexes. No differences in ensuing nucleotide substitution rates were found between the female and male cohorts studied herein. As previously reported for men, time to peak diversity in env-gp120 in women was positively associated with time to CD4+ cell count below 200 (P = 0.017), and the number of predicted N-linked glycosylation sites generally increased over time, followed by a plateau or decline, with the majority of changes localized to the V1-V2 region. These findings strongly suggest that the sex differences in HIV-1 disease progression attributed to immune system composition and sensitivities are not revealed by, nor do they impact, global patterns of viral evolution, the latter of which proceeds similarly in women and men.

HIV-1 transmission networks in high risk fishing communities on the shores of Lake Victoria in Uganda: A phylogenetic and epidemiological approach

BACKGROUND

Fishing communities around Lake Victoria in sub-Saharan Africa have been characterised as a population at high risk of HIV-infection.

METHODS

Using data from a cohort of HIV-positive individuals aged 13-49 years, enrolled from 5 fishing communities on Lake Victoria between 2009-2011, we sought to identify factors contributing to the epidemic and to understand the underlying structure of HIV transmission networks. Clinical and socio-demographic data were combined with HIV-1 phylogenetic analyses. HIV-1 gag-p24 and env-gp-41 sub-genomic fragments were amplified and sequenced from 283 HIV-1-infected participants. Phylogenetic clusters with ≥2 highly related sequences were defined as transmission clusters. Logistic regression models were used to determine factors associated with clustering.

RESULTS

Altogether, 24% (n = 67/283) of HIV positive individuals with sequences fell within 34 phylogenetically distinct clusters in at least one gene region (either gag or env). Of these, 83% occurred either within households or within community; 8/34 (24%) occurred within household partnerships, and 20/34 (59%) within community. 7/12 couples (58%) within households clustered together. Individuals in clusters with potential recent transmission (11/34) were more likely to be younger 71% (15/21) versus 46% (21/46) in un-clustered individuals and had recently become resident in the community 67% (14/21) vs 48% (22/46). Four of 11 (36%) potential transmission clusters included incident-incident transmissions. Independently, clustering was less likely in HIV subtype D (adjusted Odds Ratio, aOR = 0.51 [95% CI 0.26-1.00]) than A and more likely in those living with an HIV-infected individual in the household (aOR = 6.30 [95% CI 3.40-11.68]).

CONCLUSIONS

A large proportion of HIV sexual transmissions occur within house-holds and within communities even in this key mobile population. The findings suggest localized HIV transmissions and hence a potential benefit for the test and treat approach even at a community level, coupled with intensified HIV counselling to identify early infections.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Exploring new alleles for frost tolerance in winter rye

Rye genetic resources provide a valuable source of new alleles for the improvement of frost tolerance in rye breeding programs. Frost tolerance is a must-have trait for winter cereal production in northern and continental cropping areas. Genetic resources should harbor promising alleles for the improvement of frost tolerance of winter rye elite lines. For frost tolerance breeding, the identification of quantitative trait loci (QTL) and the choice of optimum genome-based selection methods are essential. We identified genomic regions involved in frost tolerance of winter rye by QTL mapping in a biparental population derived from a highly frost tolerant selection from the Canadian cultivar Puma and the European elite line Lo157. Lines per se and their testcrosses were phenotyped in a controlled freeze test and in multi-location field trials in Russia and Canada. Three QTL on chromosomes 4R, 5R, and 7R were consistently detected across environments. The QTL on 5R is congruent with the genomic region harboring the Frost resistance locus 2 (Fr-2) in Triticeae. The Puma allele at the Fr-R2 locus was found to significantly increase frost tolerance. A comparison of predictive ability obtained from the QTL-based model with different whole-genome prediction models revealed that besides a few large, also small QTL effects contribute to the genomic variance of frost tolerance in rye. Genomic prediction models assigning a high weight to the Fr-R2 locus allow increasing the selection intensity for frost tolerance by genome-based pre-selection of promising candidates.

Complete Genome Analysis of Thermus parvatiensis and Comparative Genomics of Thermus spp. Provide Insights into Genetic Variability and Evolution of Natural Competence as Strategic Survival Attributes

Thermophilic environments represent an interesting niche. Among thermophiles, the genus Thermus is among the most studied genera. In this study, we have sequenced the genome of Thermus parvatiensis strain RL, a thermophile isolated from Himalayan hot water springs (temperature >96°C) using PacBio RSII SMRT technique. The small genome (2.01 Mbp) comprises a chromosome (1.87 Mbp) and a plasmid (143 Kbp), designated in this study as pTP143. Annotation revealed a high number of repair genes, a squeezed genome but containing highly plastic plasmid with transposases, integrases, mobile elements and hypothetical proteins (44%). We performed a comparative genomic study of the group Thermus with an aim of analysing the phylogenetic relatedness as well as niche specific attributes prevalent among the group. We compared the reference genome RL with 16 Thermus genomes to assess their phylogenetic relationships based on 16S rRNA gene sequences, average nucleotide identity (ANI), conserved marker genes (31 and 400), pan genome and tetranucleotide frequency. The core genome of the analyzed genomes contained 1,177 core genes and many singleton genes were detected in individual genomes, reflecting a conserved core but adaptive pan repertoire. We demonstrated the presence of metagenomic islands (chromosome:5, plasmid:5) by recruiting raw metagenomic data (from the same niche) against the genomic replicons of T. parvatiensis. We also dissected the CRISPR loci wide all genomes and found widespread presence of this system across Thermus genomes. Additionally, we performed a comparative analysis of competence loci wide Thermus genomes and found evidence for recent horizontal acquisition of the locus and continued dispersal among members reflecting that natural competence is a beneficial survival trait among Thermus members and its acquisition depicts unending evolution in order to accomplish optimal fitness.

Bioinformatics in the plant genomic and phenomic domain: The German contribution to resources, services and perspectives

Plant genetic resources are a substantial opportunity for plant breeding, preservation and maintenance of biological diversity. As part of the German Network for Bioinformatics Infrastructure (de.NBI) the German Crop BioGreenformatics Network (GCBN) focuses mainly on crop plants and provides both data and software infrastructure which are tailored to the needs of the plant research community. Our mission and key objectives include: (1) provision of transparent access to germplasm seeds, (2) the delivery of improved workflows for plant gene annotation, and (3) implementation of bioinformatics services that link genotypes and phenotypes. This review introduces the GCBN's spectrum of web-services and integrated data resources that address common research problems in the plant genomics community.

Proteomic comparison of Chelidonium majus L. latex in different phases of plant development

Chelidonium majus L. (Papaveraceae) latex is used in traditinonal folk medicine to treat papillae, warts, condylomas, which are visible effects of human papilloma virus (HPV) infections. The aim of this work was to provide new insights into the biology and medicinal use of C. majus milky sap in the flowering and fruit ripening period of the plant by comparing the protein content between samples collected on respective developmental stages using LC-MS-based label-free proteome approach. For quantification, the multiplexed LC-MS data were processed using comparative chemometric approach. Progenesis LC-MS results showed that in green fruit phase (stage IV), comparing to flowering phase (stage III) of plant development, a range of proteins with higher abundance were identified as stress- and defense-related. On the other hand at stage III very intense protein synthesis, processes of transcription, protein folding and active transport of molecules (ABC transporters) are well represented. 2-DE protein maps showed an abundant set of spots with similar MWs (about 30-35 kDa) and pIs (ca. 5.5-6.5), which were identified as major latex proteins (MLPs). Therefore we suggest that biological activity of C. majus latex could be related to its protein content, which shifts during plant development from intense biosynthetic processes (biosynthesis and transport of small molecules, like alkaloids) to plant defense mechanisms against pathogens. Further studies will help to elucidate if these defense-related and pathogenesis-related proteins, like MLP, together with small-molecule compounds, could inhibit viral infection, what could be a step to fully understand the medicinal activity of C. majus latex.

Towards a whole-genome sequence for rye (Secale cereale L.)

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Breast milk and in utero transmission of HIV-1 select for envelope variants with unique molecular signatures

Background

Mother-to-child transmission of human immunodeficiency virus-type 1 (HIV-1) poses a serious health threat in developing countries, and adequate interventions are as yet unrealized. HIV-1 infection is frequently initiated by a single founder viral variant, but the factors that influence particular variant selection are poorly understood.

Results

Our analysis of 647 full-length HIV-1 subtype C and G viral envelope sequences from 22 mother–infant pairs reveals unique genotypic and phenotypic signatures that depend upon transmission route. Relative to maternal strains, intrauterine HIV transmission selects infant variants that have shorter, less-glycosylated V1 loops that are more resistant to soluble CD4 (sCD4) neutralization. Transmission through breastfeeding selects for variants with fewer potential glycosylation sites in gp41, are more sensitive to the broadly neutralizing antibodies PG9 and PG16, and that bind sCD4 with reduced cooperativity. Furthermore, experiments with Affinofile cells indicate that infant viruses, regardless of transmission route, require increased levels of surface CD4 receptor for productive infection.

Conclusions

These data provide the first evidence for transmission route-specific selection of HIV-1 variants, potentially informing therapeutic strategies and vaccine designs that can be tailored to specific modes of vertical HIV transmission.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-017-0331-z) contains supplementary material, which is available to authorized users.

FARME DB: a functional antibiotic resistance element database

Antibiotic resistance (AR) is a major global public health threat but few resources exist that catalog AR genes outside of a clinical context. Current AR sequence databases are assembled almost exclusively from genomic sequences derived from clinical bacterial isolates and thus do not include many microbial sequences derived from environmental samples that confer resistance in functional metagenomic studies. These environmental metagenomic sequences often show little or no similarity to AR sequences from clinical isolates using standard classification criteria. In addition, existing AR databases provide no information about flanking sequences containing regulatory or mobile genetic elements. To help address this issue, we created an annotated database of DNA and protein sequences derived exclusively from environmental metagenomic sequences showing AR in laboratory experiments. Our Functional Antibiotic Resistant Metagenomic Element (FARME) database is a compilation of publically available DNA sequences and predicted protein sequences conferring AR as well as regulatory elements, mobile genetic elements and predicted proteins flanking antibiotic resistant genes. FARME is the first database to focus on functional metagenomic AR gene elements and provides a resource to better understand AR in the 99% of bacteria which cannot be cultured and the relationship between environmental AR sequences and antibiotic resistant genes derived from cultured isolates.Database URL: http://staff.washington.edu/jwallace/farme.

viruSITE-integrated database for viral genomics

Viruses are the most abundant biological entities and the reservoir of most of the genetic diversity in the Earth's biosphere. Viral genomes are very diverse, generally short in length and compared to other organisms carry only few genes. viruSITE is a novel database which brings together high-value information compiled from various resources. viruSITE covers the whole universe of viruses and focuses on viral genomes, genes and proteins. The database contains information on virus taxonomy, host range, genome features, sequential relatedness as well as the properties and functions of viral genes and proteins. All entries in the database are linked to numerous information resources. The above-mentioned features make viruSITE a comprehensive knowledge hub in the field of viral genomics.

The web interface of the database was designed so as to offer an easy-to-navigate, intuitive and user-friendly environment. It provides sophisticated text searching and a taxonomy-based browsing system. viruSITE also allows for an alternative approach based on sequence search. A proprietary genome browser generates a graphical representation of viral genomes. In addition to retrieving and visualising data, users can perform comparative genomics analyses using a variety of tools.

Database URL: http://www.virusite.org/

Pharmacometabolomics Informs Viromics toward Precision Medicine

Nowadays, we are experiencing the big data era with the emerging challenge of single data interpretation. Although the advent of high-throughput technologies as well as chemo- and bio-informatics tools presents pan-omics data as the way forward to precision medicine, personalized health care and tailored-made therapeutics can be only envisaged when interindividual variability in response to/toxicity of xenobiotics can be interpreted and thus, predicted. We know that such variability is the net outcome of genetics (host and microbiota) and environmental factors (diet, lifestyle, polypharmacy, and microbiota) and for this, tremendous efforts have been made to clarify key-molecules from correlation to causality to clinical significance. Herein, we focus on the host–microbiome interplay and its direct and indirect impact on efficacy and toxicity of xenobiotics and we inevitably wonder about the role of viruses, as the least acknowledged ones. We present the emerging discipline of pharmacometabolomics-informed viromics, in which pre-dose metabotypes can assist modeling and prediction of interindividual response to/toxicity of xenobiotics. Such features, either alone or in combination with host genetics, can power biomarker discovery so long as the features are variable among patients, stable enough to be of predictive value, and better than pre-existing tools for predicting therapeutic efficacy/toxicity.

Down-regulation of the glucan synthase-like 6 gene (HvGsl6) in barley leads to decreased callose accumulation and increased cell wall penetration by Blumeria graminis f. sp. hordei

The recent characterization of the polysaccharide composition of papillae deposited at the barley cell wall during infection by the powdery mildew pathogen, Blumeria graminis f. sp. hordei (Bgh), has provided new targets for the generation of enhanced disease resistance. The role of callose in papilla-based penetration resistance of crop species is largely unknown because the genes involved in the observed callose accumulation have not been identified unequivocally. We have employed both comparative and functional genomics approaches to identify the functional orthologue of AtGsl5 in the barley genome. HvGsl6 (the barley glucan synthase-like 6 gene), which has the highest sequence identity to AtGsl5, is the only Bgh-induced gene among the HvGsls examined in this study. Through double-stranded RNA interference (dsRNAi)-mediated silencing of HvGsl6, we have shown that the down-regulation of HvGsl6 is associated with a lower accumulation of papillary and wound callose and a higher susceptibility to penetration of the papillae by Bgh, compared with control lines. The results indicate that the HvGsl6 gene is a functional orthologue of AtGsl5 and is involved in papillary callose accumulation in barley. The increased susceptibility of HvGsl6 dsRNAi transgenic lines to infection indicates that callose positively contributes to the barley fungal penetration resistance mechanism.

A Direct Comparison of Two Densely Sampled HIV Epidemics: The UK and Switzerland

Phylogenetic clustering approaches can elucidate HIV transmission dynamics. Comparisons across countries are essential for evaluating public health policies. Here, we used a standardised approach to compare the UK HIV Drug Resistance Database and the Swiss HIV Cohort Study while maintaining data-protection requirements. Clusters were identified in subtype A1, B and C pol phylogenies. We generated degree distributions for each risk group and compared distributions between countries using Kolmogorov-Smirnov (KS) tests, Degree Distribution Quantification and Comparison (DDQC) and bootstrapping. We used logistic regression to predict cluster membership based on country, sampling date, risk group, ethnicity and sex. We analysed >8,000 Swiss and >30,000 UK subtype B sequences. At 4.5% genetic distance, the UK was more clustered and MSM and heterosexual degree distributions differed significantly by the KS test. The KS test is sensitive to variation in network scale, and jackknifing the UK MSM dataset to the size of the Swiss dataset removed the difference. Only heterosexuals varied based on the DDQC, due to UK male heterosexuals who clustered exclusively with MSM. Their removal eliminated this difference. In conclusion, the UK and Swiss HIV epidemics have similar underlying dynamics and observed differences in clustering are mainly due to different population sizes.

The Triticeae Toolbox: Combining Phenotype and Genotype Data to Advance Small-Grains Breeding

The Triticeae Toolbox (; T3) is the database schema enabling plant breeders and researchers to combine, visualize, and interrogate the wealth of phenotype and genotype data generated by the Triticeae Coordinated Agricultural Project (TCAP). T3 enables users to define specific data sets for download in formats compatible with the external tools TASSEL, Flapjack, and R; or to use by software residing on the T3 server for operations such as Genome Wide Association and Genomic Prediction. New T3 tools to assist plant breeders include a Selection Index Generator, analytical tools to compare phenotype trials using common or user-defined indices, and a histogram generator for nursery reports, with applications using the Android OS, and a Field Plot Layout Designer in development. Researchers using T3 will soon enjoy the ability to design training sets, define core germplasm sets, and perform multivariate analysis. An increased collaboration with GrainGenes and integration with the small grains reference sequence resources will place T3 in a pivotal role for on-the-fly data analysis, with instant access to the knowledge databases for wheat and barley. T3 software is available under the GNU General Public License and is freely downloadable.

Reconstructing the HIV-1 CRF02_AG and CRF06_cpx epidemics in Burkina Faso and West Africa using early samples

BACKGROUND

HIV-1 circulating recombinant forms (CRFs) represent viral recombinant lineages that play a significant role in the global epidemic. Two of them dominate the epidemic in Burkina Faso: CRF06_cpx (first described in this country) and CRF02_AG. We reconstructed the phylodynamics of both recombinant viruses in Burkina Faso and throughout West Africa.

METHODS

We analysed CRF06_cpx and CRF02_AG sequences (protease/gp41) from early samples collected in Burkina Faso in 1986 together with other GenBank sequences (1984-2013) in 4 datasets: African CRF06_cpx (210/60); down-sampled CRF06_cpx (146/45); Burkina Faso CRF02_AG (130/39) and West/Central African CRF02_AG (691/298). For each dataset, we analysed both protease and gp41 jointly using the BEAST multilocus analysis and conducted phylogeographic analysis to reconstruct the early migration routes between countries.

RESULTS

The time to the most recent common ancestor (tMRCA) of CRF06_cpx was 1979 (1973-1983) for protease and 1981 (1978-1983) for gp41. The gp41 analysis inferred the origin of CRF06_cpx (or at least its parental subtype G lineage) in the Democratic Republic of Congo but migrated to Burkina Faso soon after (1982). Both genes showed that CRF06_cpx radiated to the rest of West Africa predominantly after around 1990. These results were robust to the oversampling of Burkina Faso sequences as they were confirmed in the down-sampled dataset. The tMRCA of the Burkina Faso CRF02_AG lineage was 1979 (1977-1983) for protease and 1980 (1978-1981) for gp41. However, we reconstructed its presence in West Africa much earlier (mid-1960s), with an initial origin in Cameroon and/or Nigeria, and its phylogeographic analysis revealed much interconnection within the region with a lack of country-specific phylogenetic patterns, which prevents tracking its exact migration routes.

CONCLUSIONS

Burkina Faso presents a relatively young HIV epidemic, with the diversification of the current in-country CRF02_AG and CRF06_cpx lineages taking place around 1980. This country represents the main source of CRF06_cpx in West Africa. The CRF02_AG epidemic started at least a decade earlier and showed much interchange between West African countries (especially involving coastal countries) suggesting great population mobility and an extensive viral spread in the region.

Putative nicotinic acetylcholine receptor subunits express differentially through the life cycle of codling moth, Cydia pomonella (Lepidoptera: Tortricidae)

Nicotinic acetylcholine receptors (nAChRs) are the targets of neonicotinoids and spinosads, two insecticides used in orchards to effectively control codling moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae). Orchardists in Washington State are concerned about the possibility of codling moth field populations developing resistance to these two insecticides. In an effort to help mitigate this issue, we initiated a project to identify and characterize codling moth nAChR subunits expressed in heads. This study had two main goals; (i) identify transcripts from a codling moth head transcriptome that encode for nAChR subunits, and (ii) determine nAChR subunit expression profiles in various life stages of codling moth. From a codling moth head transcriptome, 24 transcripts encoding for 12 putative nAChR subunit classes were identified and verified by PCR amplification, cloning, and sequence determination. Characterization of the deduced protein sequences encoded by putative nAChR transcripts revealed that they share the distinguishing features of the cys-loop ligand-gated ion channel superfamily with 9 α-type subunits and 3 β-type subunits identified. Phylogenetic analysis comparing these protein sequences to those of other insect nAChR subunits supports the identification of these proteins as nAChR subunits. Stage expression studies determined that there is clear differential expression of many of these subunits throughout the codling moth life cycle. The information from this study will be used in the future to monitor for potential target-site resistance mechanisms to neonicotinoids and spinosads in tolerant codling moth populations.

From FHB Resistance QTLs to Candidate Genes Identification in Triticum aestivum L

Fusarium head blight (FHB) caused by Fusarium graminearum is a worldwide destructive disease affecting cereals such as wheat. FHB resistance is a quantitative trait, and information for FHB resistance QTLs in wheat is available. However, little is known about genes underlying the FHB resistance QTL regions. Using a computational approach in this study, we have mined eight FHB resistance QTLs in wheat and predicted the candidate genes falling within these QTL intervals based on the available sequences and markers. A total of 18 genomic scaffolds located at chromosomes 2AL, 2DL, 3B and 4BS were prioritized to harbor FHB-resistant candidate genes. These genes are mainly involved in plant defense response, immune regulation and cellular detoxification. We believe that our results constitute a starting point for further validation to improve FHB-resistant bread wheat varieties.

Are PECTIN ESTERASE INHIBITOR Genes Involved in Mediating Resistance to Rhynchosporium commune in Barley?

A family of putative PECTIN ESTERASE INHIBITOR (PEI) genes, which were detected in the genomic region co-segregating with the resistance gene Rrs2 against scald caused by Rhynchosporium commune in barley, were characterized and tested for their possible involvement in mediating resistance to the pathogen by complementation and overexpression analysis. The sequences of the respective genes were derived from two BAC contigs originating from the susceptible cultivar ‘Morex’. For the genes HvPEI2, HvPEI3, HvPEI4 and HvPEI6, specific haplotypes for 18 resistant and 23 susceptible cultivars were detected after PCR-amplification and haplotype-specific CAPS-markers were developed. None of the tested candidate genes HvPEI2, HvPEI3 and HvPEI4 alone conferred a high resistance level in transgenic over-expression plants, though an improvement of the resistance level was observed especially with OE-lines for gene HvPEI4. These results do not confirm but also do not exclude an involvement of the PEI gene family in the response to the pathogen. A candidate for the resistance gene Rrs2 could not be identified yet. It is possible that Rrs2 is a PEI gene or another type of gene which has not been detected in the susceptible cultivar ‘Morex’ or the full resistance reaction requires the presence of several PEI genes.

transPLANT Resources for Triticeae Genomic Data

The genome sequences of many important Triticeae species, including bread wheat ( L.) and barley ( L.), remained uncharacterized for a long time because their high repeat content, large sizes, and polyploidy. As a result of improvements in sequencing technologies and novel analyses strategies, several of these have recently been deciphered. These efforts have generated new insights into Triticeae biology and genome organization and have important implications for downstream usage by breeders, experimental biologists, and comparative genomicists. transPLANT () is an EU-funded project aimed at constructing hardware, software, and data infrastructure for genome-scale research in the life sciences. Since the Triticeae data are intrinsically complex, heterogenous, and distributed, the transPLANT consortium has undertaken efforts to develop common data formats and tools that enable the exchange and integration of data from distributed resources. Here we present an overview of the individual Triticeae genome resources hosted by transPLANT partners, introduce the objectives of transPLANT, and outline common developments and interfaces supporting integrated data access.

blastjs: a BLAST+ wrapper for Node.js

BACKGROUND

To cope with the ever-increasing amount of sequence data generated in the field of genomics, the demand for efficient and fast database searches that drive functional and structural annotation in both large- and small-scale genome projects is on the rise. The tools of the BLAST+ suite are the most widely employed bioinformatic method for these database searches. Recent trends in bioinformatics application development show an increasing number of JavaScript apps that are based on modern frameworks such as Node.js. Until now, there is no way of using database searches with the BLAST+ suite from a Node.js codebase.

RESULTS

We developed blastjs, a Node.js library that wraps the search tools of the BLAST+ suite and thus allows to easily add significant functionality to any Node.js-based application.

CONCLUSION

blastjs is a library that allows the incorporation of BLAST+ functionality into bioinformatics applications based on JavaScript and Node.js. The library was designed to be as user-friendly as possible and therefore requires only a minimal amount of code in the client application. The library is freely available under the MIT license at https://github.com/teammaclean/blastjs.

Transmission of Non-B HIV Subtypes in the United Kingdom Is Increasingly Driven by Large Non-Heterosexual Transmission Clusters

Background. The United Kingdom human immunodeficiency virus (HIV) epidemic was historically dominated by HIV subtype B transmission among men who have sex with men (MSM). Now 50% of diagnoses and prevalent infections are among heterosexual individuals and mainly involve non-B subtypes. Between 2002 and 2010, the prevalence of non-B diagnoses among MSM increased from 5.4% to 17%, and this study focused on the drivers of this change.

Methods. Growth between 2007 and 2009 in transmission clusters among 14 000 subtype A1, C, D, and G sequences from the United Kingdom HIV Drug Resistance Database was analysed by risk group.

Results. Of 1148 clusters containing at least 2 sequences in 2007, >75% were pairs and >90% were heterosexual. Most clusters (71.4%) did not grow during the study period. Growth was significantly lower for small clusters and higher for clusters of ≥7 sequences, with the highest growth observed for clusters comprising sequences from MSM and people who inject drugs (PWID). Risk group (P < .0001), cluster size (P < .0001), and subtype (P < .01) were predictive of growth in a generalized linear model.

Discussion. Despite the increase in non-B subtypes associated with heterosexual transmission, MSM and PWID are at risk for non-B infections. Crossover of subtype C from heterosexuals to MSM has led to the expansion of this subtype within the United Kingdom.

HIV Transmission Patterns Among Immigrant Latinos Illuminated by the Integration of Phylogenetic and Migration Data

Latinos represent a growing proportion of HIV cases in North Carolina (NC). Understanding how immigrants are involved in local HIV transmission is important to guide interventions. We used phylogenetics to characterize Latino involvement in local HIV transmission chains. Transmission clusters were identified from maximum-likelihood phylogenies constructed with HIV pol sequences from 177 Latinos and 1,496 non-Latinos receiving care in NC. Highly supported clusters involving one or more Latinos were characterized. Migration data were obtained from interviews and chart review. Factors associated with cluster membership were identified using log-binomial regression. Most Latinos were male (76%), immigrants (83%), and had HIV-1B (99%). Immigrants were more likely to report heterosexual risk (67% vs. 23%) than U.S.-born Latinos (p < 0.01). We identified 32 clusters that included one or more Latinos; these involved 53 Latinos (30%) and 41 non-Latinos. Immigrant and U.S.-born Latinos were equally likely to be in clusters, but immigrants were more likely to be in clusters with another Latino (78% vs. 29%; p = 0.02). Cluster composition by ethnicity and risk behavior varied by cluster size; larger clusters contained fewer immigrants and more men who have sex with men (MSM). Factors associated with immigrant membership in local transmission clusters included age <30 years [RR 2.34 (95% CI 1.47-3.73)], Mexican origin [RR 2.55 (95% CI 1.29-6.88)], and residing in the United States longer before diagnosis [RR 1.53 (95% CI 1.09-2.15), per 10 years]. While some Latinos immigrate with HIV infection, many immigrants are involved in transmission networks after arrival, particularly MSM. HIV testing and prevention interventions must consider this heterogeneity and may be better targeted by integrating phylogenetic analyses.

A high-throughput pipeline for detecting locus-specific polymorphism in hexaploid wheat (Triticum aestivum L.)

BACKGROUND

Bread wheat (Triticum aestivum L., 2n = 6x = 42) is an allohexaploid with a huge genome. Due to the presence of extensive homoeologs and paralogs, generating locus-specific sequences can be challenging, especially when a large number of sequences are required. Traditional methods of generating locus-specific sequences are rather strenuous and time-consuming if large numbers of sequences are to be handled.

RESULTS

To improve the efficiency of isolating sequences for targeted loci, a time-saving and high-throughput pipeline integrating orthologous sequence alignment, genomic sequence retrieving, and multiple sequence alignment was developed. This pipeline was successfully employed in retrieving and aligning homoeologous sequences and 83% of the primers designed based on the pipeline successfully amplified fragments from the targeted subgenomes.

CONCLUSIONS

The high-throughput pipeline developed in this study makes it feasible to efficiently identify locus-specific sequences for large numbers of sequences. It could find applications in all research projects where locus-specific sequences are required. In addition to generating locus-specific markers, the pipeline was also used in our laboratory to identify differentially expressed genes among the three subgenomes of bread wheat. Importantly, the pipeline is not only valuable for research in wheat but should also be applicable to other allopolyploid species.

Quantitative Trait Loci Associated with the Tocochromanol (Vitamin E) Pathway in Barley

The Genome-Wide Association Studies approach was used to detect Quantitative Trait Loci associated with tocochromanol concentrations using a panel of 1,466 barley accessions. All major tocochromanol types- α-, β-, δ-, γ-tocopherol and tocotrienol- were assayed. We found 13 single nucleotide polymorphisms associated with the concentration of one or more of these tocochromanol forms in barley, seven of which were within 2 cM of sequences homologous to cloned genes associated with tocochromanol production in barley and/or other plants. These associations confirmed a prior report based on bi-parental QTL mapping. This knowledge will aid future efforts to better understand the role of tocochromanols in barley, with specific reference to abiotic stress resistance. It will also be useful in developing barley varieties with higher tocochromanol concentrations, although at current recommended daily consumption amounts, barley would not be an effective sole source of vitamin E. However, it could be an important contributor in the context of whole grains in a balanced diet.

Lack of Effectiveness of Antiretroviral Therapy in Preventing HIV Infection in Serodiscordant Couples in Uganda: An Observational Study

Background

We examined the real-world effectiveness of ART as an HIV prevention tool among HIV serodiscordant couples in a programmatic setting in a low-income country.

Methods

We enrolled individuals from HIV serodiscordant couples aged ≥18 years of age in Jinja, Uganda from June 2009 – June 2011. In one group of couples the HIV positive partner was receiving ART as they met clinical eligibility criteria (a CD4 cell count ≤250 cells/ μL or WHO Stage III/IV disease). In the second group the infected partner was not yet ART-eligible. We measured HIV incidence by testing the uninfected partner every three months. We conducted genetic linkage studies to determine the source of new infections in seroconverting participants.

Results

A total of 586 couples were enrolled of which 249 (42%) of the HIV positive participants were receiving ART at enrollment, and an additional 99 (17%) initiated ART during the study. The median duration of follow-up was 1.5 years. We found 9 new infections among partners of participants who had been receiving ART for at least three months and 8 new infections in partners of participants who had not received ART or received it for less than three months, for incidence rates of 2.09 per 100 person-years (PYRs) and 2.30 per 100 PYRs, respectively. The incidence rate ratio for ART-use was 0.91 (95% confidence interval 0.31-2.70; p=0.999). The hazard ratio for HIV seroconversion associated with ART-use by the positive partner was 1.07 (95% CI 0.41-2.80). A total of 5/7 (71%) of the transmissions on ART and 6/7 (86%) of those not on ART were genetically linked.

Conclusion

Overall HIV incidence was low in comparison to previous studies of serodiscordant couples. However, ART-use was not associated with a reduced risk of HIV transmission in this study.

Evolutionary, Comparative and Functional Analyses of the Brassinosteroid Receptor Gene, BRI1, in Wheat and Its Relation to Other Plant Genomes

Brassinosteroids (BRs) are plant hormones, fundamental for the growth and development of plants. A trans-membrane protein receptor kinase, Brassinosteroid-Insensitive 1 (BRI1), is known to interact with BRs and be directly involved in plant development. This study investigates the structural organization of BRI1 orthologs in several taxa, with a specific interest in Triticum aestivum. True orthologs of Arabidopsis thaliana BRI1 (AtBRI1) from seven-plant species showed sequence identity ranging from 54% to 95% at the protein level. All gene sequences lacked introns, leading to speculation that post-transcriptional processing in TaBRI1 is similar to AtBRI1. Based on in silico analysis, a single copy of BRI1 was present in each of the three wheat genomes on the long arm of chromosome 3. Domain structure of BRI1 orthologs among different taxa showed multiple leucine rich repeats (LRRs), an island domain (ID), a juxtamembrane/transmembrane domain (JTMD), a catalytic kinase domain (KD), C and N-Terminal domains. The KD showed the highest level of conservation while the LRRs and JTMD were most variable. Phosphorylation of residues in the juxtamembrane domain, known to be involved in the activation of the KD, is conserved in TaBRI1. While TaBRI1 has well-defined differences in the ID and LRR domains, many residues involved in ligand binding are conserved. The activation loop present in the KD showed 100% conservation in all taxa. Despite residue differences, hydrophobicity was conserved in the BR binding pocket across taxa, suggesting that function may not differ as drastically as residue identity may suggest. Predicted 3D structure of AtBRI1 and TaBRI1 showed a conserved super helical assembly, a feature essential in protein-protein interactions. An unrooted phylogram showed TaBRI1 in the monocot clade to be distinct from that of dicots. New insight in the structure and functions of BRI1 may help in targeting BR pathway for crop improvement.

Analysis of the history and spread of HIV-1 in Uganda using phylodynamics

HIV prevalence has decreased in Uganda since the 1990s, but remains substantial within high-risk groups. Here, we reconstruct the history and spread of HIV subtypes A1 and D in Uganda and explore the transmission dynamics in high-risk populations. We analysed HIV pol sequences from female sex workers in Kampala (n = 42), Lake Victoria fisher-folk (n = 46) and a rural clinical cohort (n = 74), together with publicly available sequences from adjacent regions in Uganda (n = 412) and newly generated sequences from samples taken in Kampala in 1986 (n = 12). Of the sequences from the three Ugandan populations, 60 (37.1 %) were classified as subtype D, 54 (33.3 %) as subtype A1, 31 (19.1 %) as A1/D recombinants, six (3.7 %) as subtype C, one (0.6 %) as subtype G and 10 (6.2 %) as other recombinants. Among the A1/D recombinants we identified a new candidate circulating recombinant form. Phylodynamic and phylogeographic analyses using BEAST indicated that the Ugandan epidemics originated in 1960 (1950-1968) for subtype A1 and 1973 (1970-1977) for D, in rural south-western Uganda with subsequent spread to Kampala. They also showed extensive interconnection with adjacent countries. The sequence analysis shows both epidemics grew exponentially during the 1970s-1980s and decreased from 1992, which agrees with HIV prevalence reports in Uganda. Inclusion of sequences from the 1980s indicated the origin of both epidemics was more recent than expected and substantially narrowed the confidence intervals in comparison to previous estimates. We identified three transmission clusters and ten pairs, none of them including patients from different populations, suggesting active transmission within a structured transmission network.

Highly expressed captured genes and cross-kingdom domains present in Helitrons create novel diversity in Pleurotus ostreatus and other fungi

BACKGROUND

Helitrons are class-II eukaryotic transposons that transpose via a rolling circle mechanism. Due to their ability to capture and mobilize gene fragments, they play an important role in the evolution of their host genomes. We have used a bioinformatics approach for the identification of helitrons in two Pleurotus ostreatus genomes using de novo detection and homology-based searching. We have analyzed the presence of helitron-captured genes as well as the expansion of helitron-specific helicases in fungi and performed a phylogenetic analysis of their conserved domains with other representative eukaryotic species.

RESULTS

Our results show the presence of two helitron families in P. ostreatus that disrupt gene colinearity and cause a lack of synteny between their genomes. Both putative autonomous and non-autonomous helitrons were transcriptionally active, and some of them carried highly expressed captured genes of unknown origin and function. In addition, both families contained eukaryotic, bacterial and viral domains within the helitron's boundaries. A phylogenetic reconstruction of RepHel helicases using the Helitron-like and PIF1-like helicase conserved domains revealed a polyphyletic origin for eukaryotic helitrons.

CONCLUSION

P. ostreatus helitrons display features similar to other eukaryotic helitrons and do not tend to capture host genes or gene fragments. The occurrence of genes probably captured from other hosts inside the helitrons boundaries pose the hypothesis that an ancient horizontal transfer mechanism could have taken place. The viral domains found in some of these genes and the polyphyletic origin of RepHel helicases in the eukaryotic kingdom suggests that virus could have played a role in a putative lateral transfer of helitrons within the eukaryotic kingdom. The high similarity of some helitrons, along with the transcriptional activity of its RepHel helicases indicates that these elements are still active in the genome of P. ostreatus.

A mutant ('lab strain') of the hyperthermophilic archaeon Pyrococcus furiosus, lacking flagella, has unusual growth physiology

A mutant ('lab strain') of the hyperthermophilic archaeon Pyrococcus furiosus DSM3638 exhibited an extended exponential phase and atypical cell aggregation behavior. Genomic DNA from the mutant culture was sequenced and compared to wild-type (WT) DSM3638, revealing 145 genes with one or more insertions, deletions, or substitutions (12 silent, 33 amino acid substitutions, and 100 frame shifts). Approximately, half of the mutated genes were transposases or hypothetical proteins. The WT transcriptome revealed numerous changes in amino acid and pyrimidine biosynthesis pathways coincidental with growth phase transitions, unlike the mutant whose transcriptome reflected the observed prolonged exponential phase. Targeted gene deletions, based on frame-shifted ORFs in the mutant genome, in a genetically tractable strain of P. furiosus (COM1) could not generate the extended exponential phase behavior observed for the mutant. For example, a putative radical SAM family protein (PF2064) was the most highly up-regulated ORF (>25-fold) in the WT between exponential and stationary phase, although this ORF was unresponsive in the mutant; deletion of this gene in P. furiosus COM1 resulted in no apparent phenotype. On the other hand, frame-shifting mutations in the mutant genome negatively impacted transcription of a flagellar biosynthesis operon (PF0329-PF0338).Consequently, cells in the mutant culture lacked flagella and, unlike the WT, showed minimal evidence of exopolysaccharide-based cell aggregation in post-exponential phase. Electron microscopy of PF0331-PF0337 deletions in P. furiosus COM1 showed that absence of flagella impacted normal cell aggregation behavior and, furthermore, indicated that flagella play a key role, beyond motility, in the growth physiology of P. furiosus.

Unraveling the web of viroinformatics: computational tools and databases in virus research

The beginning of the second century of research in the field of virology (the first virus was discovered in 1898) was marked by its amalgamation with bioinformatics, resulting in the birth of a new domain--viroinformatics. The availability of more than 100 Web servers and databases embracing all or specific viruses (for example, dengue virus, influenza virus, hepatitis virus, human immunodeficiency virus [HIV], hemorrhagic fever virus [HFV], human papillomavirus [HPV], West Nile virus, etc.) as well as distinct applications (comparative/diversity analysis, viral recombination, small interfering RNA [siRNA]/short hairpin RNA [shRNA]/microRNA [miRNA] studies, RNA folding, protein-protein interaction, structural analysis, and phylotyping and genotyping) will definitely aid the development of effective drugs and vaccines. However, information about their access and utility is not available at any single source or on any single platform. Therefore, a compendium of various computational tools and resources dedicated specifically to virology is presented in this article.

Human cancer databases (review)

Cancer is one of the four major non‑communicable diseases (NCD), responsible for ~14.6% of all human deaths. Currently, there are >100 different known types of cancer and >500 genes involved in cancer. Ongoing research efforts have been focused on cancer etiology and therapy. As a result, there is an exponential growth of cancer‑associated data from diverse resources, such as scientific publications, genome‑wide association studies, gene expression experiments, gene‑gene or protein‑protein interaction data, enzymatic assays, epigenomics, immunomics and cytogenetics, stored in relevant repositories. These data are complex and heterogeneous, ranging from unprocessed, unstructured data in the form of raw sequences and polymorphisms to well‑annotated, structured data. Consequently, the storage, mining, retrieval and analysis of these data in an efficient and meaningful manner pose a major challenge to biomedical investigators. In the current review, we present the central, publicly accessible databases that contain data pertinent to cancer, the resources available for delivering and analyzing information from these databases, as well as databases dedicated to specific types of cancer. Examples for this wealth of cancer‑related information and bioinformatic tools have also been provided.

SkateBase, an elasmobranch genome project and collection of molecular resources for chondrichthyan fishes

Chondrichthyan fishes are a diverse class of gnathostomes that provide a valuable perspective on fundamental characteristics shared by all jawed and limbed vertebrates. Studies of phylogeny, species diversity, population structure, conservation, and physiology are accelerated by genomic, transcriptomic and protein sequence data. These data are widely available for many sarcopterygii (coelacanth, lungfish and tetrapods) and actinoptergii (ray-finned fish including teleosts) taxa, but limited for chondrichthyan fishes.  In this study, we summarize available data for chondrichthyes and describe resources for one of the largest projects to characterize one of these fish, Leucoraja erinacea, the little skate.  SkateBase ( http://skatebase.org) serves as the skate genome project portal linking data, research tools, and teaching resources.

De novo transcriptome sequencing and analysis of the cereal cyst nematode, Heterodera avenae

The cereal cyst nematode (CCN, Heterodera avenae) is a major pest of wheat (Triticum spp) that reduces crop yields in many countries. Cyst nematodes are obligate sedentary endoparasites that reproduce by amphimixis. Here, we report the first transcriptome analysis of two stages of H. avenae. After sequencing extracted RNA from pre parasitic infective juvenile and adult stages of the life cycle, 131 million Illumina high quality paired end reads were obtained which generated 27,765 contigs with N50 of 1,028 base pairs, of which 10,452 were annotated. Comparative analyses were undertaken to evaluate H. avenae sequences with those of other plant, animal and free living nematodes to identify differences in expressed genes. There were 4,431 transcripts common to H. avenae and the free living nematode Caenorhabditis elegans, and 9,462 in common with more closely related potato cyst nematode, Globodera pallida. Annotation of H. avenae carbohydrate active enzymes (CAZy) revealed fewer glycoside hydrolases (GHs) but more glycosyl transferases (GTs) and carbohydrate esterases (CEs) when compared to M. incognita. 1,280 transcripts were found to have secretory signature, presence of signal peptide and absence of transmembrane. In a comparison of genes expressed in the pre-parasitic juvenile and feeding female stages, expression levels of 30 genes with high RPKM (reads per base per kilo million) value, were analysed by qRT-PCR which confirmed the observed differences in their levels of expression levels. In addition, we have also developed a user-friendly resource, Heterodera transcriptome database (HATdb) for public access of the data generated in this study. The new data provided on the transcriptome of H. avenae adds to the genetic resources available to study plant parasitic nematodes and provides an opportunity to seek new effectors that are specifically involved in the H. avenae-cereal host interaction.

A customized Web portal for the genome of the ctenophore Mnemiopsis leidyi

BACKGROUND

Mnemiopsis leidyi is a ctenophore native to the coastal waters of the western Atlantic Ocean. A number of studies on Mnemiopsis have led to a better understanding of many key biological processes, and these studies have contributed to the emergence of Mnemiopsis as an important model for evolutionary and developmental studies. Recently, we sequenced, assembled, annotated, and performed a preliminary analysis on the 150-megabase genome of the ctenophore, Mnemiopsis. This sequencing effort has produced the first set of whole-genome sequencing data on any ctenophore species and is amongst the first wave of projects to sequence an animal genome de novo solely using next-generation sequencing technologies.

DESCRIPTION

The Mnemiopsis Genome Project Portal (http://research.nhgri.nih.gov/mnemiopsis/) is intended both as a resource for obtaining genomic information on Mnemiopsis through an intuitive and easy-to-use interface and as a model for developing customized Web portals that enable access to genomic data. The scope of data available through this Portal goes well beyond the sequence data available through GenBank, providing key biological information not available elsewhere, such as pathway and protein domain analyses; it also features a customized genome browser for data visualization.

CONCLUSIONS

We expect that the availability of these data will allow investigators to advance their own research projects aimed at understanding phylogenetic diversity and the evolution of proteins that play a fundamental role in metazoan development. The overall approach taken in the development of this Web site can serve as a viable model for disseminating data from whole-genome sequencing projects, framed in a way that best-serves the specific needs of the scientific community.

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.