seqCNA: an R package for DNA copy number analysis in cancer using high-throughput sequencing

Background

Deviations in the amount of genomic content that arise during tumorigenesis, called copy number alterations, are structural rearrangements that can critically affect gene expression patterns. Additionally, copy number alteration profiles allow insight into cancer discrimination, progression and complexity. On data obtained from high-throughput sequencing, improving quality through GC bias correction and keeping false positives to a minimum help build reliable copy number alteration profiles.

Results

We introduce seqCNA, a parallelized R package for an integral copy number analysis of high-throughput sequencing cancer data. The package includes novel methodology on (i) filtering, reducing false positives, and (ii) GC content correction, improving copy number profile quality, especially under great read coverage and high correlation between GC content and copy number. Adequate analysis steps are automatically chosen based on availability of paired-end mapping, matched normal samples and genome annotation.

Conclusions

seqCNA, available through Bioconductor, provides accurate copy number predictions in tumoural data, thanks to the extensive filtering and better GC bias correction, while providing an integrated and parallelized workflow.

Electronic supplementary material

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

Whole-genome sequencing of Mesorhizobium huakuii 7653R provides molecular insights into host specificity and symbiosis island dynamics

Background

Evidence based on genomic sequences is urgently needed to confirm the phylogenetic relationship between Mesorhizobium strain MAFF303099 and M. huakuii. To define underlying causes for the rather striking difference in host specificity between M. huakuii strain 7653R and MAFF303099, several probable determinants also require comparison at the genomic level. An improved understanding of mobile genetic elements that can be integrated into the main chromosomes of Mesorhizobium to form genomic islands would enrich our knowledge of how genome dynamics may contribute to Mesorhizobium evolution in general.

Results

In this study, we sequenced the complete genome of 7653R and compared it with five other Mesorhizobium genomes. Genomes of 7653R and MAFF303099 were found to share a large set of orthologs and, most importantly, a conserved chromosomal backbone and even larger perfectly conserved synteny blocks. We also identified candidate molecular differences responsible for the different host specificities of these two strains. Finally, we reconstructed an ancestral Mesorhizobium genomic island that has evolved into diverse forms in different Mesorhizobium species.

Conclusions

Our ortholog and synteny analyses firmly establish MAFF303099 as a strain of M. huakuii. Differences in nodulation factors and secretion systems T3SS, T4SS, and T6SS may be responsible for the unique host specificities of 7653R and MAFF303099 strains. The plasmids of 7653R may have arisen by excision of the original genomic island from the 7653R chromosome.

Electronic supplementary material

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

Identification,isolation and expression analysis of eight stress-related <Emphasis Type="Italic">R2R3</Emphasis>-<Emphasis Type="Italic">MYB</Emphasis> genes in tartary buckwheat (<Emphasis Type="Italic">Fagopyrum tataricum</Emphasis>)

Key message

Eight R2R3 - MYB genes in tartary buckwheat were identified, and their expression patterns were comprehensively analyzed, which reveals role in plant response to abiotic stresses.

Abstract

The proteins of the R2R3-MYB superfamily play key roles in the growth and development processes as well as defense responses in plants. However, their characteristics and functions have not been fully investigated in tartary buckwheat (Fagopyrum tataricum), a strongly abiotic resistant coarse cereal. In this article, eight tartary buckwheat R2R3-MYB genes were isolated with full-length cDNA and DNA sequences. Phylogenetic analysis of the members of the R2R3-MYB superfamily between Arabidopsis and tartary buckwheat revealed that the assumed functions of the eight tartary buckwheat R2R3-MYB proteins are divided into five Arabidopsis functional subgroups that are involved in abiotic stress. Expression analysis during abiotic stress and exogenous phytohormone treatments identified that the eight R2R3-MYB genes responded to one or more treatments. This study is the first comprehensive analysis of the R2R3-MYB gene family in tartary buckwheat under abiotic stress.

     

MPAgenomics: an R package for multi-patient analysis of genomic markers

Background

Last generations of Single Nucleotide Polymorphism (SNP) arrays allow to study copy-number variations in addition to genotyping measures.

Results

MPAgenomics, standing for multi-patient analysis (MPA) of genomic markers, is an R-package devoted to: (i) efficient segmentation and (ii) selection of genomic markers from multi-patient copy number and SNP data profiles. It provides wrappers from commonly used packages to streamline their repeated (sometimes difficult) manipulation, offering an easy-to-use pipeline for beginners in R.The segmentation of successive multiple profiles (finding losses and gains) is performed with an automatic choice of parameters involved in the wrapped packages. Considering multiple profiles in the same time, MPAgenomics wraps efficient penalized regression methods to select relevant markers associated with a given outcome.

Conclusions

MPAgenomics provides an easy tool to analyze data from SNP arrays in R. The R-package MPAgenomics is available on CRAN.     

Increased levels of soluble interleukin-6 receptor and CCL3 in COPD sputum

Background

COPD patients have increased numbers of macrophages and neutrophils in the lungs. Interleukin-6 (IL-6) trans-signaling via its soluble receptor sIL-6R, governs the influx of innate immune cells to inflammatory foci through regulation of the chemokine CCL3. We hypothesized that there would be enhanced levels of IL-6, sIL-6R and CCL3 in COPD sputum.

Methods

59 COPD patients, 15 HNS and 15 S underwent sputum induction and processing with phosphate buffered saline to obtain supernatants for IL-6, sIL-6R and CCL3 analysis. Cytoslides were produced for differential cell counting and immunocytochemistry (COPD; n = 3) to determine cell type surface expression of the CCL3 receptors CCR5 and CCR1.

Results

COPD patients expressed higher levels (p < 0.05) of sIL-6R and CCL3 compared to controls (sIL-6R medians pg/ml: COPD 166.4 vs S 101.1 vs HNS 96.4; CCL3 medians pg/ml: COPD 117.9 vs S 0 vs HNS 2.7). COPD sIL-6R levels were significantly correlated with sputum neutrophil (r = 0.5, p < 0.0001) and macrophage (r = 0.3, p = 0.01) counts. Immunocytochemical analysis revealed that CCR5 and CCR1 were exclusively expressed on airway macrophages.

Conclusion

Enhanced airway generation of sIL-6R may promote IL-6 trans-signaling in COPD. Associated upregulation of CCL3 may facilitate the recruitment of macrophages into the airways by ligation of CCR1 and CCR5.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0103-4) contains supplementary material, which is available to authorized users.     

BACA: bubble chArt to compare annotations

Background

DAVID is the most popular tool for interpreting large lists of gene/proteins classically produced in high-throughput experiments. However, the use of DAVID website becomes difficult when analyzing multiple gene lists, for it does not provide an adequate visualization tool to show/compare multiple enrichment results in a concise and informative manner.

Result

We implemented a new R-based graphical tool, BACA (Bubble chArt to Compare Annotations), which uses the DAVID web service for cross-comparing enrichment analysis results derived from multiple large gene lists. BACA is implemented in R and is freely available at the CRAN repository (http://cran.r-project.org/web/packages/BACA/).

Conclusion

The package BACA allows R users to combine multiple annotation charts into one output graph by passing DAVID website.

Electronic supplementary material

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

Beyond genomic variation - comparison and functional annotation of three Brassica rapa genomes: a turnip,a rapid cycling and a Chinese cabbage

Background

Brassica rapa is an economically important crop species. During its long breeding history, a large number of morphotypes have been generated, including leafy vegetables such as Chinese cabbage and pakchoi, turnip tuber crops and oil crops.

Results

To investigate the genetic variation underlying this morphological variation, we re-sequenced, assembled and annotated the genomes of two B. rapa subspecies, turnip crops (turnip) and a rapid cycling. We then analysed the two resulting genomes together with the Chinese cabbage Chiifu reference genome to obtain an impression of the B. rapa pan-genome. The number of genes with protein-coding changes between the three genotypes was lower than that among different accessions of Arabidopsis thaliana, which can be explained by the smaller effective population size of B. rapa due to its domestication. Based on orthology to a number of non-brassica species, we estimated the date of divergence among the three B. rapa morphotypes at approximately 250,000 YA, far predating Brassica domestication (5,000-10,000 YA).

Conclusions

By analysing genes unique to turnip we found evidence for copy number differences in peroxidases, pointing to a role for the phenylpropanoid biosynthesis pathway in the generation of morphological variation. The estimated date of divergence among three B. rapa morphotypes implies that prior to domestication there was already considerably divergence among B. rapa genotypes. Our study thus provides two new B. rapa reference genomes, delivers a set of computer tools to analyse the resulting pan-genome and uses these to shed light on genetic drivers behind the rich morphological variation found in B. rapa.

Electronic supplementary material

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

Identification and characterization of LIM gene family in Brassica rapa

Background

LIM (Lin-11, Isl-1 and Mec-3 domains) genes have been reported to trigger the formation of actin bundles, a major higher-order cytoskeletal assembly, in higher plants; however, the stress resistance related functions of these genes are still not well known. In this study, we collected 22 LIM genes designated as Brassica rapa LIM (BrLIM) from the Brassica database, analyzed the sequences, compared them with LIM genes of other plants and analyzed their expression after applying biotic and abiotic stresses in Chinese cabbage.

Results

Upon sequence analysis these genes were confirmed as LIM genes and found to have a high degree of homology with LIM genes of other species. These genes showed distinct clusters when compared to other recognized LIM proteins upon phylogenetic analysis. Additionally, organ specific expression of these genes was observed in Chinese cabbage plants, with BrPLIM2a, b, c, BrDAR1, BrPLIM2e, f and g only being expressed in flower buds. Furthermore, the expression of these genes (except for BrDAR1 and BrPLIM2e) was high in the early flowering stages. The remaining genes were expressed in almost all organs tested. All BrDAR genes showed higher expression in flower buds compared to other organs. These organ specific expressions were clearly correlated with the phylogenetic grouping. In addition, BrWLIM2c and BrDAR4 responded to Fusarium oxysporum f. sp. conglutinans infection, while commonly two BrDARs and eight BrLIMs responded to cold, ABA and pH (pH5, pH7 and pH9) stress treatments in Chinese cabbage plants.

Conclusion

Taken together, the results of this study indicate that BrLIM and BrDAR genes may be involved in resistance against biotic and abiotic stresses in Brassica.

Electronic supplementary material

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

Nutritional status modulates box C/D snoRNP biogenesis by regulated subcellular relocalization of the R2TP complex

Background

Box C/D snoRNPs, which are typically composed of box C/D snoRNA and the four core protein components Nop1, Nop56, Nop58, and Snu13, play an essential role in the modification and processing of pre-ribosomal RNA. The highly conserved R2TP complex, comprising the proteins Rvb1, Rvb2, Tah1, and Pih1, has been shown to be required for box C/D snoRNP biogenesis and assembly; however, the molecular basis of R2TP chaperone-like activity is not yet known.

Results

Here, we describe an unexpected finding in which the activity of the R2TP complex is required for Nop58 protein stability and is controlled by the dynamic subcellular redistribution of the complex in response to growth conditions and nutrient availability. In growing cells, the complex localizes to the nucleus and interacts with box C/D snoRNPs. This interaction is significantly reduced in poorly growing cells as R2TP predominantly relocalizes to the cytoplasm. The R2TP-snoRNP interaction is mainly mediated by Pih1.

Conclusions

The R2TP complex exerts a novel regulation on box C/D snoRNP biogenesis that affects their assembly and consequently pre-rRNA maturation in response to different growth conditions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0404-4) contains supplementary material, which is available to authorized users.     

Identification and distribution of the NBS-LRR gene family in the Cassava genome

Background

Plant resistance genes (R genes) exist in large families and usually contain both a nucleotide-binding site domain and a leucine-rich repeat domain, denoted NBS-LRR. The genome sequence of cassava (Manihot esculenta) is a valuable resource for analysing the genomic organization of resistance genes in this crop.

Results

With searches for Pfam domains and manual curation of the cassava gene annotations, we identified 228 NBS-LRR type genes and 99 partial NBS genes. These represent almost 1% of the total predicted genes and show high sequence similarity to proteins from other plant species. Furthermore, 34 contained an N-terminal toll/interleukin (TIR)-like domain, and 128 contained an N-terminal coiled-coil (CC) domain. 63% of the 327 R genes occurred in 39 clusters on the chromosomes. These clusters are mostly homogeneous, containing NBS-LRRs derived from a recent common ancestor.

Conclusions

This study provides insight into the evolution of NBS-LRR genes in the cassava genome; the phylogenetic and mapping information may aid efforts to further characterize the function of these predicted R genes.

Electronic supplementary material

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

mAPKL: R/ Bioconductor package for detecting gene exemplars and revealing their characteristics

Background

So far many algorithms have been proposed towards the detection of significant genes in microarray analysis problems. Several of those approaches are freely available as R-packages though their engagement in gene expression analysis by non-bioinformaticians is usually a frustrating task. Besides, only some of those packages offer a complete suite of tools starting from initial data import and ending to analysis report. Here we present an R/Bioconductor package that implements a hybrid gene selection method along with a bunch of functions to facilitate a thorough and convenient gene expression profiling analysis.

Results

mAPKL is an open-source R/Bioconductor package that implements the mAP-KL hybrid gene selection method. The advantage of this method is that selects a small number of gene exemplars while achieving comparable classification results to other well established algorithms on a variety of datasets and dataset sizes. The mAPKL package is accompanied with extra functionalities including (i) solid data import; (ii) data sampling following a user-defined proportion; (iii) preprocessing through several normalization and transformation alternatives; (iv) classification with the aid of SVM and performance evaluation; (v) network analysis of the significant genes (exemplars), including degree of centrality, closeness, betweeness, clustering coefficient as well as the construction of an edge list table; (vi) gene annotation analysis, (vii) pathway analysis and (viii) auto-generated analysis reporting.

Conclusions

Users are able to run a thorough gene expression analysis in a timely manner starting from raw data and concluding to network characteristics of the selected gene exemplars. Detailed instructions and example data are provided in the R package, which is freely available at Bioconductor under the GPL-2 or later license http://www.bioconductor.org/packages/3.1/bioc/html/mAPKL.html.

Electronic supplementary material

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