Package 'Rgff'

Test consistency and order of a GFF file

Description

This function tests the consistency and order of a GFF file.

Usage

check_gff(inFile, fileType = c("AUTO", "GFF3", "GTF"))

Arguments

inFile	Path to the input GFF file
fileType	Version of the input file (GTF/GFF3). Default AUTO: determined from the file name.

Details

The following list indicates the code and description of the issues detected in GFF3 files

NCOLUMNS_EXCEEDED

Input file contains lines with more than 9 fields

NCOLUMNS_INFERIOR

Input file contains lines with less than 9 fields

TOO_MANY_FEATURE_TYPES

Input file contains too many (more than 100) different feature types

NO_IDs

ID attribute not found in any feature

DUPLICATED_IDs

There are duplicated IDs

ID_IN_MULTIPLE_CHR

The same ID has been found in more than one chromosome

NO_PARENTs

Parent attribute not found in any feature

MISSING_PARENT_IDs

There are missing Parent IDs

PARENT_IN DIFFERENT CHR

There are features whose Parent is located in a different chromosome

PARENT_DEFINED_BEFORE_ID

Feature ids referenced in Parent attribute before being defined as ID

NOT_GROUPED_BY_CHR

Features are not grouped by chromosome

NOT_SORTED_BY_COORDINATE

Features are not sorted by start coordinate

NOT_VALID_WARNING

File cannot be recognized as valid GFF3. Parsing warnings.

NOT_VALID_ERROR

File cannot be recognized as valid GFF3. Parsing errors.

The following list indicates the code and description of the issues detected in GTF files

NCOLUMNS_EXCEEDED

Input file contains lines with more than 9 fields

NCOLUMNS_INFERIOR

Input file contains lines with less than 9 fields

TOO_MANY_FEATURE_TYPES

Input file contains too many (more than 100) different feature types

NO_GENE_ID_ATTRIBUTE

gene_id attribute not found in any feature

MISSING_GENE_IDs

There are features without gene_id attribute

NO_GENE_FEATURES

Gene features are not included in this GTF file

DUPLICATED_GENE_IDs

There are duplicated gene_ids

GENE_ID_IN_MULTIPLE_CHR

The same gene_id has been found in more than one chromosome

NO_TRANSCRIPT_ID_ATTRIBUTE

transcript_id attribute not found in any feature There are no elements with transcript_id attribute

MISSING_TRANSCRIPT_IDs

There are features without transcript_id attribute

NO_TRANSCRIPT_FEATURES

Transcript features are not included in this GTF file

DUPLICATED_TRANSCRIPT_IDs

There are duplicated transcript_ids

TRANSCRIPT_ID_IN_MULTIPLE_CHR

The same transcript_id has been found in more than one chromosome

DUPLICATED_GENE_AND_TRANSCRIPT_IDs

Same id has been defined as gene_id and transcript_id

NOT_GROUPED_BY_CHR

Features are not grouped by chromosome

NOT_SORTED_BY_COORDINATE

Features are not sorted by start coordinate

NOT_VALID_WARNING

File cannot be recognized as valid GTF. Parsing warnings.

NOT_VALID_ERROR

File cannot be recognized as valid GTF. Parsing errors.

Value

A data frame of detected issues, including a short code name, a description and estimated severity each. In no issues are detected the function will return an empty data frame.

Examples

test_gff3<-system.file("extdata", "eden.gff3", package="Rgff")
check_gff(test_gff3)

---

Analyses the feature type hierarchy of a GFF file

Description

Based on the feature type hierarchy a GFF file, this function creates and returns a feature tree or a feature dependency table.

Usage

get_features(
  inFile,
  includeCounts = FALSE,
  outFormat = c("tree", "data.frame", "JSON"),
  fileType = c("AUTO", "GFF3", "GTF")
)

Arguments

inFile	Path to the input GTF/GFF3 features file
includeCounts	Include number of occurrences of each feature and subfeature
outFormat	Output format of the function. Available formats are: tree (DEFAULT), data.frame and JSON.
fileType	Version of the input file (GTF/GFF3). Default AUTO: determined from the file name.

Value

Depending on the outFormat selected returns a feature tree (tree), a feature dependency table as data.frame (data.frame) or a feature dependency table as JSON object (JSON)

Examples

test_gff3<-system.file("extdata", "AthSmall.gff3", package="Rgff")
get_features(test_gff3)

---

Summarizes the number of features of each type in a GFF file

Description

This function summarizes the number of features of each type in a GFF file and returns the statistics

Usage

gff_stats(inFile)

Arguments

inFile

Path to the input GFF file

Value

A tibble with the summary data

Examples

test_gff3<-system.file("extdata", "AthSmall.gff3", package="Rgff")
gff_stats(test_gff3)

---

Summarizes the number of elements of each type in each chromosome of a GFF file

Description

This function summarizes the number of features of each type in each chromosome of a GFF file and returns the statistics

Usage

gff_stats_by_chr(inFile)

Arguments

inFile

Path to the input GFF file

Value

A tibble with the summary data

Examples

test_gff3<-system.file("extdata", "AthSmall.gff3", package="Rgff")
gff_stats_by_chr(test_gff3)

---

Converts a GTF file into a GFF3 file

Description

This function converts a GTF file into a GFF3 file mantaining the feature hierarchy defined by the gene_id and transcript_id attributes. The remaining attributes of each feature will be kept with the same name and value.

Usage

gtf_to_gff3(gtfFile, outFile, forceOverwrite = FALSE)

Arguments

gtfFile	Path to the input GTF file
outFile	Path to the output GFF3 file, inf not provided the output will be gtfFile.gff3
forceOverwrite	If output file exists, overwrite the existing file. (default FALSE)

Value

Path to the generated GFF3 file

Examples

## Not run: 
test_gtf<-system.file("extdata", "AthSmall.gtf", package="Rgff")
gtf_to_gff3(test_gtf)

## End(Not run)

---

Plots or exports an image of the feature tree from a GFF file

Description

This function plots the feature tree from a GFF file or, if an output file name is provided, exports an image of in the desired format ("png", "pdf" or "svg"). Packages "DiagrammeR", "DiagrammeRsvg" and "rsvg" must be installed to use this function.

Usage

plot_features(
  inFile,
  outFile,
  includeCounts = FALSE,
  fileType = c("AUTO", "GFF3", "GTF"),
  exportFormat = c("png", "pdf", "svg")
)

Arguments

inFile	Path to the input GFF file
outFile	Path to the output features image file, if not provided the tree will be plotted
includeCounts	Include number of occurrences of each subfeature
fileType	Version of the input file (GTF/GFF3). If not provided it will be determined from the file name.
exportFormat	Output image format when it is not possible to deduce it from the extension of outFile ("png", "pdf" or "svg"). Default, "png"

Value

Path of the output features image file

Examples

test_gff3<-system.file("extdata", "AthSmall.gff3", package="Rgff")
plot_features(test_gff3)

---

Creates a SAF file from a GTF/GFF3 features for the given pairs of blocks and features

Description

This function creates a SAF file from a GTF/GFF3 features for the given blocks and features

Usage

saf_from_gff(
  inFile,
  outFile,
  fileType = c("AUTO", "GFF3", "GTF"),
  forceOverwrite = FALSE,
  features = c("gene > exon"),
  sep = ">"
)

Arguments

inFile	Path to the input GFF file
outFile	Path to the output SAF file, if not provided the output path will be the input path with the suffix ".feature1-block1.feature2-block2(...).saf"
fileType	Version of the input file (GTF/GFF3). Default AUTO: determined from the file name.
forceOverwrite	If output file exists, overwrite the existing file. (default FALSE)
features	Vector of pairs of features/blocks, separated by '>' (see sep argument). In the case of features without defined blocks, only the feature is needed (see example)
sep	Separator of each "feature" and "block" provided in the feature argument (default '>')

Value

Path to the generated SAF file

Examples

test_gff3<-system.file("extdata", "AthSmall.gff3", package="Rgff")
## Default usage, extract gene features by exon blocks 
saf_from_gff(test_gff3)
## Define only feature without block to count reads within the whole genomic locus 
saf_from_gff(test_gff3, features=c("gene"))
## Define multiple features for counting readsoverlapping only in exonic regions 
saf_from_gff(test_gff3, features=c("gene > exon", "ncRNA_gene > exon"))

---

Sorts a GTF/GFF3 file

Description

This function produces a sorted GFF file from an unsorted GFF file. The default order is by Chromosome, Start, End (reverse) and feature (based on the precedency in feature tree)

Usage

sort_gff(
  inFile,
  outFile,
  fileType = c("AUTO", "GFF3", "GTF"),
  forceOverwrite = FALSE
)

Arguments

inFile	Path to the input GFF file
outFile	Path to the output sorted file, if not provided the output will be the input path (without extension) with the suffix sorted.gtf/gff3
fileType	Version of the input file (GTF/GFF3). Default AUTO: determined from the file name.
forceOverwrite	If output file exists, overwrite the existing file. (default FALSE)

Value

Path to the sorted feature file

Examples

test_gff3<-system.file("extdata", "eden.gff3", package="Rgff")
sort_gff(test_gff3)

---