The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

Vertical and Horizontal Inheritance Consistence Analysis

Description:

The package implements the VHICA method described in Wallau et al. (2016). The purpose of the method is to detect horizontal transfers of transposable elements, by contrasting the divergence of transposable element sequences with that of regular genes. Two files should be provided, for both a set of reference genes and transposable element sequences: (i) pairwise divergence across species (e.g., dS), (ii) codon usage bias for all genes and elements in all species.

Details:

This package contains three main functions. * ‘read.vhica’: reads sequence files and generates an object of class ‘vhica’ that will be used for further analysis. * ‘plot.vhica’: plots the VHICA regression between two species, and displays how far transposable elements (or any kind of other sequences) are from the reference genes.

In addition, it provides tools to calculate divergence (‘div’) and codon usage bias (‘CUB’), which are necessary to apply the VHICA method.

Author(s):

Implementation: Arnaud Le Rouzic

Scientists who designed the method: Gabriel Wallau, Aurélie Hua-Van, Arnaud Le~Rouzic.

Maintainer: Arnaud Le Rouzic arnaud.le-rouzic@universite-paris-saclay.fr

References:

Wallau, G. L., Capy, P., Loreto, E., Le Rouzic, A., & Hua-Van, A. (2016). VHICA, a new method to discriminate between vertical and horizontal transposon transfer: Application to the mariner family within Drosophila. Molecular biology and evolution, 33(4), 1094-1109.

Examples:

 file.cb <- system.file("extdata", "mini-cbias.txt", package="vhica")
 file.div <- system.file("extdata", "mini-div.txt", package="vhica")
 file.tree <- if(require("ape")) system.file("extdata", "phylo.nwk", package="vhica") else NULL
 vc <- read.vhica(cb.filename=file.cb, div.filename=file.div)
 plot(vc, "dere", "dana")
 im <- image(vc, "mellifera:6", treefile=file.tree, skip.void=TRUE)
 summary(im)
 

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.