Wide cross-species RNA-Seq comparison reveals convergent molecular mechanisms involved in nickel hyperaccumulation across dicotyledons.

The Anthropocene epoch is associated with the spreading of metals in the environment increasing oxidative and genotoxic stress on organisms. Interestingly, c. 520 plant species growing on metalliferous soils acquired the capacity to accumulate and tolerate a tremendous amount of nickel in their shoots. The wide phylogenetic distribution of these species suggests that nickel hyperaccumulation evolved multiple times independently. However, the exact nature of these mechanisms and whether they have been recruited convergently in distant species is not known. To address these questions, we have developed a cross-species RNA-Seq approach combining differential gene expression analysis and cluster of orthologous group annotation to identify genes linked to nickel hyperaccumulation in distant plant families. Our analysis reveals candidate orthologous genes encoding convergent function involved in nickel hyperaccumulation, including the biosynthesis of specialized metabolites and cell wall organization. Our data also point out that the high expression of IREG/Ferroportin transporters recurrently emerged as a mechanism involved in nickel hyperaccumulation in plants. We further provide genetic evidence in the hyperaccumulator Noccaea caerulescens for the role of the NcIREG2 transporter in nickel sequestration in vacuoles. Our results provide molecular tools to better understand the mechanisms of nickel hyperaccumulation and study their evolution in plants.

A regression model for estimating DNA copy number applied to capture sequencing data.

MOTIVATION: Target enrichment, also referred to as DNA capture, provides an effective way to focus sequencing efforts on a genomic region of interest. Capture data are typically used to detect single-nucleotide variants. It can also be used to detect copy number alterations, which is particularly useful in the context of cancer, where such changes occur frequently. In copy number analysis, it is a common practice to determine log-ratios between test and control samples, but this approach results in a loss of information as it disregards the total coverage or intensity at a locus. RESULTS: We modeled the coverage or intensity of the test sample as a linear function of the control sample. This regression approach is able to deal with regions that are completely deleted, which are problematic for methods that use log-ratios. To demonstrate the utility of our approach, we used capture data to determine copy number for a set of 600 genes in a panel of nine breast cancer cell lines. We found high concordance between our results and those generated using a single-nucleotide polymorphsim genotyping platform. When we compared our results with other log-ratio-based methods, including ExomeCNV, we found that our approach produced better overall correlation with SNP data. AVAILABILITY: The algorithm is implemented in C and R and the code can be downloaded from http://bioinformatics.nki.nl/ocs/ CONTACT: l.wessels@nki.nl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.