bridgeHomologues: Use bridge markers to cluster homologues into linkage groups

In mdavy86/polymapR: Linkage Analysis in Outcrossing Polyploids

Description

Clustering at high LOD scores results in marker clusters representing homologues. bridgeHomologues clusters these (pseudo)homologues to linkage groups using linkage information between 1.0 and bridge markers within a parent (e.g. 2.0 for a tetraploid). If parent-specific bridge markers (e.g. 2.0) cannot be used, biparental markers can also be used (e.g. 1.1, 1.2, 2.1, 2.2 and 1.3 markers). The linkage information between 1.0 and biparental markers can be combined.

Usage

bridgeHomologues(cluster_stack, cluster_stack2 = NULL, linkage_df,
  linkage_df2 = NULL, LOD_threshold = 10, automatic_clustering = TRUE,
  LG_number = 5, parentname = "", min_bridges = 2, log = NULL)

Arguments

cluster_stack	A data.frame with a column "marker" specifying markernames, and a column "cluster" specifying marker cluster
cluster_stack2	Optional. A cluster_stack for the other parent. Use this argument if cross-parent markers are used (e.g. when using 1.1 markers).
linkage_df	A linkage data.frame as output of linkage between bridge (e.g. 1.0 and 2.0) markers.
linkage_df2	Optional. A linkage_df specifying linkages between 1.0 and cross-parent markers in the other parent. Use this argument if cross-parent markers are used (e.g. when using 1.1, 2.1, 1.2 and/or 2.2 markers). The use of multiple types of cross-parent markers is allowed.
LOD_threshold	Integer. The LOD threshold specifying at which LOD score a link between 1.0 and bridge (e.g. 2.0) markers is used for clustering homologues.
automatic_clustering	Logical. Should clustering be executed without user input?
LG_number	Integer. Expected number of chromosomes (linkage groups)
parentname	Name of the parent. Used in the main title of the plot.
min_bridges	The minimum number of cross-parent bridges for a link to be considered. Make this number higher if there are a lot of spurious links. This argument is only used for biparental bridges (N.B. if cluster_stack2 and linkage_df2 are specified).
log	Character string specifying the log filename to which standard output should be written. If NULL log is send to stdout.

Value

A data.frame with markers classified by homologue and linkage group.

Examples

data("P1_homologues", "P2_homologues", "SN_DN_P1", "SN_SS_P1", "SN_SS_P2")
ChHomDf<-bridgeHomologues(cluster_stack = P1_homologues[["5"]],
                 linkage_df=SN_DN_P1,
                 LOD_threshold=4,
                 automatic_clustering=TRUE,
                 LG_number=5,
                 parentname="P1")

ChHomDf<-bridgeHomologues(cluster_stack = P1_homologues[["5"]],
                           cluster_stack2 = P2_homologues[["5"]],
                 linkage_df=SN_SS_P1,
                 linkage_df2=SN_SS_P2,
                 LOD_threshold=4,
                 automatic_clustering=TRUE,
                 LG_number=5,
                 parentname="P1")

mdavy86/polymapR documentation built on May 25, 2019, 9:35 p.m.