Package: RAIDS 1.5.0

Pascal Belleau

RAIDS: Accurate Inference of Genetic Ancestry from Cancer Sequences

This package implements specialized algorithms that enable genetic ancestry inference from various cancer sequences sources (RNA, Exome and Whole-Genome sequences). This package also implements a simulation algorithm that generates synthetic cancer-derived data. This code and analysis pipeline was designed and developed for the following publication: Belleau, P et al. Genetic Ancestry Inference from Cancer-Derived Molecular Data across Genomic and Transcriptomic Platforms. Cancer Res 1 January 2023; 83 (1): 49–58.

Authors:Pascal Belleau [cre, aut], Astrid Deschênes [aut], David A. Tuveson [aut], Alexander Krasnitz [aut]

RAIDS_1.5.0.tar.gz
RAIDS_1.5.0.zip(r-4.5)RAIDS_1.5.0.zip(r-4.4)RAIDS_1.5.0.zip(r-4.3)
RAIDS_1.5.0.tgz(r-4.4-any)RAIDS_1.5.0.tgz(r-4.3-any)
RAIDS_1.5.0.tar.gz(r-4.5-noble)RAIDS_1.5.0.tar.gz(r-4.4-noble)
RAIDS_1.5.0.tgz(r-4.4-emscripten)RAIDS_1.5.0.tgz(r-4.3-emscripten)
RAIDS.pdf |RAIDS.html
RAIDS/json (API)
NEWS

# Install 'RAIDS' in R:
install.packages('RAIDS', repos = c('https://bioc.r-universe.dev', 'https://cloud.r-project.org'))

Peer review:

Bug tracker:https://github.com/krasnitzlab/raids/issues

Datasets:
  • demoKnownSuperPop1KG - The known super population ancestry of the demo 1KG reference profiles.
  • demoPCA1KG - The PCA results of the demo 1KG reference dataset for demonstration purpose. Beware that the PCA has been run on a very small subset of the 1KG reference dataset and should not be used to call ancestry inference on a real profile.
  • demoPCASyntheticProfiles - The PCA result of demo synthetic profiles projected on the demo subset 1KG reference PCA.
  • demoPedigreeEx1 - The pedigree information about a demo profile called 'ex1'.
  • matKNNSynthetic - A small 'data.frame' containing the inferred ancestry on the synthetic profiles.
  • pedSynthetic - A small 'data.frame' containing the information related to synthetic profiles. The ancestry of the profiles used to generate the synthetic profiles must be present.
  • snpPositionDemo - A small 'data.frame' containing the SNV information.

On BioConductor:RAIDS-1.5.0(bioc 3.21)RAIDS-1.4.0(bioc 3.20)

geneticssoftwaresequencingwholegenomeprincipalcomponentgeneticvariabilitydimensionreductionbiocviewsancestrycancer-genomicsexome-sequencinggenomicsinferencer-languagerna-seqrna-sequencingwhole-genome-sequencing

6.30 score 5 stars 22 scripts 104 downloads 39 exports 155 dependencies

Last updated 23 days agofrom:bb243c99fc. Checks:OK: 1 NOTE: 6. Indexed: yes.

TargetResultDate
Doc / VignettesOKNov 19 2024
R-4.5-winNOTENov 19 2024
R-4.5-linuxNOTENov 19 2024
R-4.4-winNOTENov 19 2024
R-4.4-macNOTENov 19 2024
R-4.3-winNOTENov 19 2024
R-4.3-macNOTENov 19 2024

Exports:add1KG2SampleGDSaddBlockFromDetFileaddGeneBlockGDSRefAnnotaddGeneBlockRefAnnotaddRef2GDS1KGaddStudy1KgcomputeAncestryFromSyntheticcomputeAncestryFromSyntheticFilecomputeKNNRefSamplecomputeKNNRefSyntheticcomputePCAMultiSyntheticcomputePCARefSamplecomputePoolSyntheticAncestryGrcomputeSyntheticROCcreateAccuracyGraphcreateAUROCGraphcreateStudy2GDS1KGestimateAllelicFractiongenerateGDS1KGgenerateMapSnvSelgeneratePhase1KG2GDSgeneratePhaseRefgetRef1KGPopgetRefSuperPopgroupChr1KGSNVidentifyRelativeidentifyRelativeRefinferAncestryinferAncestryGeneAwareprepPed1KGprepSyntheticpruningSamplerunExomeAncestryrunRNAAncestryselect1KGPopselect1KGPopForSyntheticsnvListVCFsplitSelectByPopsyntheticGeno

Dependencies:abindAnnotationDbiAnnotationFilteraskpassbackportsBHBiobaseBiocGenericsBiocIOBiocParallelBiostringsbitbit64bitopsblobbootbroomBSgenomecachemclassclicliprcodetoolscolorspacecpp11crayoncurldata.tableDBIDelayedArrayDNAcopydplyrensembldbfansifarverfastmapforcatsforeachformatRformula.toolsfutile.loggerfutile.optionsgdsfmtgenericsGENESISGenomeInfoDbGenomeInfoDbDataGenomicAlignmentsGenomicFeaturesGenomicRangesggplot2glmnetgluegtableGWASExactHWGWASToolshavenhmshttrigraphIRangesisobanditeratorsjomojsonliteKEGGRESTlabelinglambda.rlatticelazyevallifecyclelme4lmtestlogistfmagrittrMASSMatrixMatrixGenericsMatrixModelsmatrixStatsmemoisemgcvmicemimeminqamitmlmunsellnlmenloptrnnetnumDerivopenssloperator.toolsordinalpanpillarpkgconfigplogrplyrpngprettyunitspROCprogressProtGenericspurrrquantregquantsmoothR6RColorBrewerRcppRcppEigenRCurlreadrreshape2restfulrRhtslibrjsonrlangrpartRsamtoolsRSQLitertracklayerS4ArraysS4VectorssandwichscalesSeqArraySeqVarToolsshapesnowSNPRelateSparseArraySparseMstringistringrSummarizedExperimentsurvivalsystibbletidyrtidyselecttzdbucminfUCSC.utilsutf8VariantAnnotationvctrsviridisLitevroomwithrXMLXVectoryamlzlibbioczoo

Accurate Inference of Genetic Ancestry from Cancer-derived Sequences

Rendered fromRAIDS.Rmdusingknitr::rmarkdownon Nov 19 2024.

Last update: 2024-10-22
Started: 2023-03-30

Population reference dataset GDS files

Rendered fromCreate_Reference_GDS_File.Rmdusingknitr::rmarkdownon Nov 19 2024.

Last update: 2023-09-29
Started: 2023-08-14

Using wrappper functions

Rendered fromWrappers.Rmdusingknitr::rmarkdownon Nov 19 2024.

Last update: 2024-10-16
Started: 2024-10-16

Readme and manuals

Help Manual

Help pageTopics
Add the genotype information for the list of pruned SNVs into the Profile GDS fileadd1KG2SampleGDS
Append information associated to ld blocks, as indexes, into the Population Reference SNV Annotation GDS fileaddBlockFromDetFile
Append information associated to blocks, as indexes, into the Population Reference SNV Annotation GDS fileaddGeneBlockGDSRefAnnot
Append information associated to blocks, as indexes, into the Population Reference SNV Annotation GDS fileaddGeneBlockRefAnnot
Add the information about the unrelated patients to the Reference GDS fileaddRef2GDS1KG
Append information about the 1KG samples into the Profile GDS fileaddStudy1Kg
Select the optimal K and D parameters using the synthetic data and infer the ancestry of a specific profilecomputeAncestryFromSynthetic
Select the optimal K and D parameters using the synthetic data and infer the ancestry of a specific profilecomputeAncestryFromSyntheticFile
Run a k-nearest neighbors analysis on one specific profilecomputeKNNRefSample
Run a k-nearest neighbors analysis on a subset of the synthetic datasetcomputeKNNRefSynthetic
Project synthetic profiles onto existing principal component axes generated using the reference 1KG profilescomputePCAMultiSynthetic
Project specified profile onto PCA axes generated using known reference profilescomputePCARefSample
Run a PCA analysis and a K-nearest neighbors analysis on a small set of synthetic data using all 1KG profiles except the ones used to generate the synthetic profilescomputePoolSyntheticAncestryGr
Calculate the AUROC of the inferences for specific values of D and K using the inferred ancestry results from the synthetic profiles.computeSyntheticROC
Generate accuracy graphcreateAccuracyGraph
Generate accuracy graphcreateAUROCGraph
Create the Profile GDS file(s) for one or multiple specific profiles using the information from a RDS Sample description file and the 1KG GDS filecreateStudy2GDS1KG
The known super population ancestry of the demo 1KG reference profiles.demoKnownSuperPop1KG
The PCA results of the demo 1KG reference dataset for demonstration purpose. Beware that the PCA has been run on a very small subset of the 1KG reference dataset and should not be used to call ancestry inference on a real profile.demoPCA1KG
The PCA result of demo synthetic profiles projected on the demo subset 1KG reference PCA.demoPCASyntheticProfiles
The pedigree information about a demo profile called 'ex1'.demoPedigreeEx1
Estimate the allelic fraction of the pruned SNVs for a specific profileestimateAllelicFraction
Generate the GDS file that will contain the information from Reference data set (reference data set)generateGDS1KG
Generate the filter SNP information file in RDS formatgenerateMapSnvSel
Adding the phase information into the Reference GDS filegeneratePhase1KG2GDS
Adding the phase information into the Reference GDS filegeneratePhaseRef
Extract the specified column from the 1KG GDS 'sample.ref' node for the reference profiles (real ancestry assignation)getRef1KGPop
Extract the from the 1KG GDS 'sample.ref' node for the reference profiles (real ancestry assignation)getRefSuperPop
Merge the genotyping files per chromosome into one filegroupChr1KGSNV
Identify genetically unrelated patients in GDS Reference fileidentifyRelative
Identify genetically unrelated patients in GDS Reference fileidentifyRelativeRef
Run most steps leading to the ancestry inference call on a specific DNA profileinferAncestry
Run most steps leading to the ancestry inference call on a specific RNA profileinferAncestryGeneAware
A small 'data.frame' containing the inferred ancestry on the synthetic profiles.matKNNSynthetic
A small 'data.frame' containing the information related to synthetic profiles. The ancestry of the profiles used to generate the synthetic profiles must be present.pedSynthetic
Prepare the pedigree file using pedigree information from ReferenceprepPed1KG
Add information related to the synthetic profiles (study and synthetic reference profiles information) into a Profile GDS fileprepSynthetic
Compute the list of pruned SNVs for a specific profile using the information from the Reference GDS file and a linkage disequilibrium analysispruningSample
Run most steps leading to the ancestry inference call on a specific exome profilerunExomeAncestry
Run most steps leading to the ancestry inference call on a specific RNA profilerunRNAAncestry
Random selection of a specific number of reference profiles in each subcontinental population present in the 1KG GDS fileselect1KGPop
Random selection of a specific number of reference profiles in each subcontinental population present in the 1KG GDS file ( same as select1KGPop but the function doesn't need gds object as parameters but the file name of the referenceGDS )select1KGPopForSynthetic
A small 'data.frame' containing the SNV information.snpPositionDemo
Generate a VCF with the information from the SNPs that pass a cut-off thresholdsnvListVCF
Group samples per subcontinental populationsplitSelectByPop
Generate synthetic profiles for each cancer profile and 1KG reference profile combination and add them to the Profile GDS filesyntheticGeno