Skip to main content

MetaboLights open metabolomics data repository utility tools and API.

Project description

Metabolights MetaboLights Utils Library

MetaboLights MetaboLights Github ISA-Tab version Licence

Python Coverage


metaboligts-utils is a lightweight API and command line interface (CLI) to use or search MetaboLights public studies, update MetaboLights submitted studies and ISA-Tab metadata files.


CLI Features

  • Download and list MetaboLights public study metadata files.
  • Search MetaboLights public studies [beta release].
  • Download, upload and validate MetaboLights submitted studies
  • Create assay file templates.

API Features

  • Read and update ISA files with minimum pyhton package dependency.
  • Read MetaboLights study metadata files
  • Calculate hash value of each ISA metadata file or study ISA metadata folder
  • Read ISA table files (s_.txt, a_.txt, m_*.txt) with Pagination support.
  • Multi-column filters and sort options on ISA table files.
  • Apply actions on ISA table files to manipulate them.
  • Json serializable models with pydantic library.

Installation


The following command installs metabolights-utils from the Python Package Index. You will need Python 3.8+ and pip3 on your operating system.

  • (If python3 is not installed) Download and install Python
  • (If pip3 is not installed) Install pip3
  • Install metabolights-utils library
cd <directory to create a virtual environment>

# install metabolights-utils on new a virtual environment named mtbls-venv (you may change it)
python3 -m venv mtbls-venv
source mtbls-venv/bin/activate
pip install --upgrade pip
pip3 install -U metabolights-utils

# test mtbls command
mtbls --version

CLI (mtbls) Usage

After installation of metabolights-utils, mtbls command will be enabled to use or search MetaboLights studies.

Commands to use MetaboLights Public Studies

cd <directory to installed virtual environment.>
# activate your virtual environment
source mtbls-venv/bin/activate

# prints mtbls public command usage
mtbls --help

    # Usage: mtbls [OPTIONS] COMMAND [ARGS]...

    # Options:
    #   --help  Show this message and exit.

    # Commands:
    #   model       Commands to explain MetaboLights study data model.
    #   public      Commands to use MetaboLights public study data and ISA...
    #   submission  Commands to use MetaboLights study submission REST API.

# prints mtbls public command usage
mtbls public --help

        # Usage: mtbls public [OPTIONS] COMMAND [ARGS]...

        # Options:
        # --help  Show this message and exit.

        # Commands:
        # describe  View summary of any public study content.
        # download  Download study data and metadata files from MetaboLights FTP...
        # list      List studies and study folder content.
        # remove    Delete local study data and metadata files.

# prints mtbls public list command
mtbls public list --help

# lists all public studies listed on FTP server
mtbls public list

# lists all public studies on local storage
mtbls public list -l

# list root directory content of study
mtbls public list MTBLS3

# lists the content of FILES folder
mtbls public list MTBLS3 FILES

# downloads study metadata files from MetaboLights FTP server.
mtbls public download MTBLS3

    # DOWNLOADED      a_MTBLS3_live_metabolite_profiling_mass_spectrometry.txt
    # DOWNLOADED      s_MTBLS3.txt
    # DOWNLOADED      m_MTBLS3_live_metabolite_profiling_mass_spectrometry_v2_maf.tsv
    # DOWNLOADED      i_Investigation.txt

# uses local copies of study metadata files
mtbls public download MTBLS3
    # SKIPPED m_MTBLS3_live_metabolite_profiling_mass_spectrometry_v2_maf.tsv
    # SKIPPED a_MTBLS3_live_metabolite_profiling_mass_spectrometry.txt
    # SKIPPED i_Investigation.txt
    # SKIPPED s_MTBLS3.txt

# force to download study metadata files
mtbls public download MTBLS3 -o


# downloads study FILES folder (data files) from MetaboLights FTP server.
mtbls public download MTBLS3 FILES

# downloads a study data file from MetaboLights FTP server.
mtbls public download MTBLS3 FILES/Cecilia_AA_rerun45.raw



# help for mtbls study describe command
mtbls public describe --help

# prints summary of MTBLS3 study.
mtbls public describe MTBLS3

# prints MTBLS3 study title.
mtbls public describe MTBLS3 "$.investigation.studies[0].title"

# prints MTBLS3 study title.
mtbls public describe MTBLS3 "$.investigation.studies[0].study_assays.assays[*].assay_technique.name"

# prints MTBLS3 study protocol names.
mtbls public describe MTBLS3 "$.investigation.studies[0].study_protocols[*].protocols[*].name"

# prints study model descriptions.
mtbls model explain

# prints descriptions of property: model -> investigation.
mtbls model explain investigation

# prints descriptions of property: model -> assays.
mtbls model explain assays

# prints descriptions of property: model -> investigation -> studies.
mtbls model explain investigation.studies

# deletes local data and metadata files of MTBLS1 study.
mtbls public remove MTBLS1

Commands to search MetaboLights Public Studies

# prints help for public study search command.
mtbls public search

    # Usage: mtbls public search [OPTIONS] [QUERY]

    #   Search public studies with query keywords. If there are multiple search
    #   keywords and no join operator (+, |) defined, results are merged with the
    #   selected query join operator (and, or)

    #   query: query terms that will be in search. e.g. cancer, (mus musculus)

    # Options:
    #   -u, --search_rest_api_url TEXT  MetaboLights search API URL.
    #   -s, --skip INTEGER              Skip n items from the matched items.
    #   -l, --limit INTEGER             Maximum number items in response. Maximum
    #                                   return size is 100 items.
    #   -j, --query_join_operator TEXT  If multiple keywords are defined and there
    #                                   is no join operator (+, |) in query, One of
    #                                   the 'and' (default) or 'or' operator will be
    #                                   used.
    #   -b, --body TEXT                 Advanced filter options in json format.
    #                                   Please read the API documentation.
    #   --id, --study_ids               Shows only MetaboLights accession numbers.
    #   --raw                           Shows raw result in json format.
    #   --help                          Show this message and exit.

# prints summary of indexed public study on MetaboLights
mtbls public search MTBLS1

# prints indexed data of a public study on MetaboLights in json format
mtbls public search MTBLS1 --raw


# prints only number of search hits for query (targeted and lipid)
mtbls public search "targeted + lipid" --limit=0


# prints summary of first 10 studies in search results 
mtbls public search "targeted + lipid" --skip=0 --limit=10


# prints summary of 10 studies in search results 
mtbls public search "targeted + lipid" --skip=10 --limit=10

# prints summary of the first public study in search results  
mtbls public search "(mus musculus) + cancer + (lipidomics | lipidomic)" --skip=0 --limit=1


# prints number of the public studies in search results (search cancer term if there are organism samples from mus musculus **and** homo sapiens)
mtbls public search "cancer" --body '{"ontologyFilters": {"organism":{ "joinOperator": "and", "values": ["mus musculus", "homo sapiens"]}}}' --limit=0

# prints public study ids of the search result
mtbls public search "cancer" --body '{"ontologyFilters": {"organism":{ "joinOperator": "and", "values": ["mus musculus", "homo sapiens"]}}}' --limit=100 --id

# prints number of the public studies in search results (search cancer term if there is an organism sample from mus musculus **or** homo sapiens)
mtbls public search "cancer" --body '{"ontologyFilters": {"organism":{ "joinOperator": "and", "values": ["mus musculus", "homo sapiens"]}}}' --limit=0

# prints only study ids of search result (list studies that contain goose or chicken term and have LC-MS or GC-MS assays) 
mtbls public search "(goose | chicken)" --body '{"assayTechniqueNameFilter": { "joinOperator": "or", "values": ["LC-MS", "GC-MS"] }}' --limit=100 --id

# prints studies both NMR and LC-MS assays (no query)
mtbls public search --body '{"assayTechniqueNameFilter": { "joinOperator": "and", "values": ["NMR", "LC-MS"] }}' --limit=100 --id

# select term
mtbls public search "(cow + -sheep + -(blood serum | milk))" --body '{"assayTechniqueNameFilter": { "joinOperator": "or", "values": ["LC-MS", "GC-MS"] }}'

# prints only study ids 
mtbls public search "(cow | -sheep)" --body '{"assayTechniqueNameFilter": { "joinOperator": "or", "values": ["LC-MS", "GC-MS"] }}' --limit=100 --id

# prints only study ids if study has 'SCIEX QTRAP' term in mass spectrometry assay files
mtbls public search "(SCIEX QTRAP)"  --body '{"assayMainTechniqueFilter": { "joinOperator": "or", "values": ["MS"] }}' --limit 100 --id

# prints study ids which have both targeted and untargeted design descriptor terms
mtbls public search --body '{"ontologyFilters": {"design_descriptor":{ "joinOperator": "and", "values": ["targeted", "untargeted"]}}}' --limit=100 --id

# prints study ids which have cancer or covid design descriptor terms
mtbls public search --body '{"ontologyFilters": {"design_descriptor":{ "joinOperator": "or", "values": ["cancer", "disease"]}}}' --limit=100 --id

# prints study ids which have multiomics design descriptor term
mtbls public search --body '{"ontologyFilters": {"design_descriptor":{ "joinOperator": "or", "values": ["multiomics", "multi-omics"]}}}' --limit=100 --id

# prints number of studies aggregated by organism  (filters aggregations less than 5 item and shows only top 50 organisms)
mtbls public search --body '{"aggregations": [{ "aggregationName": "organisms", "fieldName": "organisms.term", "maxItemCount": 50, "minItemCount": 5}]}' --limit=0 --raw

# prints number of studies by assay techniques
mtbls public search --body '{"aggregations": [{ "aggregationName": "assay_techniques", "fieldName": "assayTechniques.name", "maxItemCount": 50, "minItemCount": 1}]}' --limit=0 --raw

Commands to update MetaboLights Submitted (private) Studies

# prints help for submitted studies.
mtbls submission

    # Usage: mtbls submission [OPTIONS] COMMAND [ARGS]...

    #   Commands to use MetaboLights study submission REST API.

    # Options:
    #   --help  Show this message and exit.

    # Commands:
    #   create-assay  Creates a study assay and maf file.
    #   describe      View summary of any user submitted study content.
    #   download      Download submission study metadata files.
    #   list          List submitted studies and study folder content.
    #   login         Creates a file path to use connect private FTP server and...
    #   upload        Uploads local metadata files to private FTP and start...
    #   validate      Validate submitted study and save validation report on...


# saves credentials to connect private FTP and use MetaboLights API
mtbls submission login

# lists studies submitted by user
mtbls submission list

# lists root folder content of submitted/public study
mtbls submission list <MTBLSXXX>

# lists subfoler content of submitted/public study
mtbls submission list <MTBLSXXX> FILES

# summarizes the submitted/public study
mtbls submission describe <MTBLSXXX>


# downloads the submitted/public study
mtbls submission download <MTBLSXXX>


# uploads all metadata files, starts private FTP folder-Study Folder sync task for the submitted study
mtbls submission upload <MTBLSXXX> -o



# uploads all metadata files, starts private FTP folder-Study Folder sync task for the submitted study
mtbls submission upload <MTBLSXXX> -o


mtbls submission create-assay
    # Usage: mtbls submission create-assay [OPTIONS] STUDY_ID ASSAY_TECHNIQUE

    #   Creates a study assay and maf file.

    #   study_id: MetaboLights study accession number (MTBLSxxxx). assay_technique:
    #   Valid assay techniques:  NMR, LC-MS, LC-DAD, GC-MS, GCxGC-MS, GC-FID, DI-MS,
    #   FIA-MS, CE-MS, MALDI-MS, MSImaging

    #   Acronyms: Diode array detection (LC-DAD), Tandem MS (GCxGC-MS), Flame
    #   ionisation detector (GC-FID), Direct infusion (DI-MS), Flow injection
    #   analysis (FIA-MS), Capillary electrophoresis (CE-MS), Matrix-assisted laser
    #   desorption-ionisation imaging mass spectrometry (MALDI-MS), Nuclear magnetic
    #   resonance (NMR), Mass spec spectrometry (MSImaging)

    #   scan_polarity (optional): Valid only for LC-MS, LC-DAD, GC-MS, GCxGC-MS, GC-
    #   FID. Valid values: positive, negative, alternating

    #   column_type (optional): Valid only for LC-MS, LC-DAD. Valid values: hilic,
    #   reverse phase, direct infusion

    # Options:
    #   -p, --local_path TEXT           Local storage root path. Folder will be
    #                                   created if it does not exist.
    #   -u, --rest_api_base_url TEXT    MetaboLights study submission API base URL.
    #   -x, --local_cache_path TEXT     Path to store cache files of study
    #                                   submission file indices, study models, etc.
    #   -c, --credentials_file_path TEXT
    #                                   Path to store cache files of study
    #                                   submission file indices, study models, etc.
    #   --scan_polarity, --polarity TEXT
    #                                   Scan polarity of the assay.
    #   --column_type TEXT              Column type of the assay.
    #   --help                          Show this message and exit.

# creates and downloads LC-MS assay and linked maf files for a study
mtbls submission create-assay <MTBLSXXX> LC-MS --polarity positive --column_type hilic


# creates and downloads NMR assay and linked maf files for a study
mtbls submission create-assay <MTBLSXXX> NMR


mtbls submission validate
    # Usage: mtbls submission validate [OPTIONS] STUDY_ID

    #   Validate submitted study and save validation report on local storage.

    #   study_id: MetaboLights study accession number (MTBLSxxxx).

    # Options:
    #   -u, --rest_api_base_url TEXT    MetaboLights study submission API base URL.
    #   -x, --local_cache_path TEXT     Path to store cache files of study
    #                                   submission file indices, study models, etc.
    #   -c, --credentials_file_path TEXT
    #                                   Path to store cache files of study
    #                                   submission file indices, study models, etc.
    #   -v, --validation_file_path TEXT
    #                                   Path to store validation file.
    #   --help                          Show this message and exit.


# runs study validation and creates tsv validation report for errors. 
mtbls submission validate <MTBLSXXX>

# runs study validation and creates tsv validation report with given filename on current working directory. 
mtbls submission validate <MTBLSXXX> --validation_file_path errors.tsv

API Usage

Load MetaboLights study model from a directory


Read i_Investigation.txt, s_,txt, a_.txt and m_*.tsv files in a study folder.

    provider = MetabolightsStudyProvider()
    model, messages = provider.load_study(
        study_id,
        study_path,
        load_assay_files=True,
        load_sample_file=True,
        load_maf_files=True
    )

Read and update Investigation files


Read and update an investigation file. Results are json serializable, so you can use them with REST APIs.

import os
import pathlib
import uuid
from typing import List

from metabolights_utils.isatab import Reader, Writer
from metabolights_utils.isatab.reader import InvestigationFileReader, InvestigationFileReaderResult
from metabolights_utils.isatab.writer import InvestigationFileWriter
from metabolights_utils.models.isa.common import OntologyItem
from metabolights_utils.models.isa.investigation_file import Assay, Study


def test_investigation_file_write_01():
    file_path = pathlib.Path("tests/test-data/MTBLS398/i_Investigation.txt")
    reader: InvestigationFileReader = Reader.get_investigation_file_reader()
    result: InvestigationFileReaderResult = reader.read(file_buffer_or_path=file_path)

    assay: Assay = result.investigation.studies[0].study_assays.assays[0]
    assay.measurement_type = OntologyItem(
        term="test",
        term_source_ref="test source",
        term_accession_number="test accesion",
    )

    tmp_file_name = uuid.uuid4().hex
    tmp_path = pathlib.Path(f"test-temp/test_{tmp_file_name}.txt")
    writer: InvestigationFileWriter = Writer.get_investigation_file_writer()
    writer.write(result.investigation, file_buffer_or_path=tmp_path)

    result: InvestigationFileReaderResult = reader.read(file_buffer_or_path=tmp_path)
    assay_read: Assay = result.investigation.studies[0].study_assays.assays[0]
    assert assay_read.measurement_type.term == "test"
    assert assay_read.measurement_type.term_source_ref == "test source"
    assert assay_read.measurement_type.term_accession_number == "test accesion"
    os.remove(tmp_path)

Read and Update ISA Table Files with Pagination Support


  • Select page size (number of rows in a page) and read results with the selected page size.
  • Define custom row offset and read only selected rows. Actual row indices in a result may be unordered after filter and sort operations.
  • Read the selected columns you defined. If a selected column has additional columns (Term Source REF, etc) and these columns are not defined by user, they will also be in the result. Column names may be different than header if there are multiple columns with same header.
  • If columns are not selected, all table columns will be returned in result. If columns are selected, result will contain only these columns in selected order.

Example 1: Read ISA table file


Read selected assay file (a_.txt) rows and columns with pagination support. You can use same methods with sample (s_.txt) and metabolite assignment (m_*.tsv) files.

import pathlib

from metabolights_utils.isatab import Reader
from metabolights_utils.isatab.reader import (
    IsaTableFileReader,
    IsaTableFileReaderResult
)


def test_assay_file_success_01():
    file_path = pathlib.Path(
        "tests/test-data/MTBLS373/a_MTBLS373_sevinaskascreen_metabolite_profiling_mass_spectrometry.txt"
    )
    reader: IsaTableFileReader = Reader.get_assay_file_reader()

    # get page count. Default results per page is 100
    page_count: int = reader.get_total_pages(file_buffer_or_path=file_path)
    assert page_count == 147

    # get page count with custom page count.
    page_count = reader.get_total_pages(file_buffer_or_path=file_path, results_per_page=50)
    assert page_count == 294

    # get total row count
    total_rows_count = reader.get_total_row_count(file_buffer_or_path=file_path)
    assert total_rows_count == 14670

    # get isa table headers
    result: IsaTableFileReaderResult = reader.get_headers(file_buffer_or_path=file_path)
    assert len(result.parser_report.messages) == 0
    assert "Parameter Value[Column model]" in result.isa_table_file.table.columns

    # get isa table rows. Default offset is 0. Read 88 rows
    result: IsaTableFileReaderResult = reader.get_rows(file_buffer_or_path=file_path, limit=88)
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 88

    # get isa table rows. Read 70 rows from offset 14600
    result: IsaTableFileReaderResult = reader.get_rows(
        file_buffer_or_path=file_path, offset=14600, limit=70
    )
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 70
    assert result.isa_table_file.table.row_offset == 14600

    # get page 2 from isa table. Default page limit is 100. Read 100 items from offset 100
    result: IsaTableFileReaderResult = reader.get_page(page=2, file_buffer_or_path=file_path)
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 100

    # get page 2 from isa table. Read 50 items from offset 50
    result: IsaTableFileReaderResult = reader.get_page(
        page=2, results_per_page=50, file_buffer_or_path=file_path
    )
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 50

    # get last from isa table. Read 20 items from offset 14650
    result: IsaTableFileReaderResult = reader.get_page(
        page=294, results_per_page=50, file_buffer_or_path=file_path
    )
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 20

    # get page 2 with selected columns from isa table.
    # read 50 items from offset 50 (page 2)
    # Addition columns will be in result even if they are not selected.
    # Parameter Value[Autosampler model] is ontology column. So 2 new columns will be added to result.
    result: IsaTableFileReaderResult = reader.get_page(
        page=2,
        results_per_page=50,
        file_buffer_or_path=file_path,
        selected_columns=["Sample Name", "Parameter Value[Autosampler model]"],
    )
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 50
    assert len(result.isa_table_file.table.columns) == 4

Example 2: Update ISA table file


Load ISA table page and save after update table content

import os
import pathlib
import shutil

from metabolights_utils.isatab import Reader, Writer
from metabolights_utils.isatab.reader import (
    IsaTableFileReader,
    IsaTableFileReaderResult
)
from metabolights_utils.isatab.writer import IsaTableFileWriter


def test_assay_file_read_write():
    path_original = pathlib.Path(
        "tests/test-data/MTBLS1/a_MTBLS1_metabolite_profiling_NMR_spectroscopy.txt"
    )
    file_path = (
        "test-temp/test-data/MTBLS1/a_MTBLS1_metabolite_profiling_NMR_spectroscopy.txt"
    )
    os.makedirs(os.path.dirname(file_path), exist_ok=True)
    shutil.copy(path_original, file_path)
    helper: IsaTableFileReader = Reader.get_assay_file_reader()

    with open(file_path, "r", encoding="utf-8") as file_buffer:
        # read second page of assa file
        result: IsaTableFileReaderResult = helper.get_page(
            file_buffer_or_path=file_buffer,
            page=2,
            results_per_page=50,
            file_path=str(file_path),
            selected_columns=[
                "Sample Name",
                "Derived Spectral Data File",
                "Metabolite Assignment File",
            ],
        )
        assert len(result.parser_report.messages) == 0
        assert result.isa_table_file.table.row_count == 50
        assert len(result.isa_table_file.table.columns) == 3
    
    writer: IsaTableFileWriter = Writer.get_assay_file_writer()
    isa_table = result.isa_table_file.table

    sha256_hash = result.isa_table_file.sha256_hash
    # save same content without any update
    report = writer.save_isa_table(
        file_path=str(file_path), file_sha256_hash=sha256_hash, isa_table=isa_table
    )
    assert report.success

    first_column = result.isa_table_file.table.columns[0]
    result.isa_table_file.table.data[first_column][0] = "Updated Sample Name"
    # save updated content
    report = writer.save_isa_table(
        file_path=str(file_path), file_sha256_hash=sha256, isa_table=isa_table
    )
    assert report.success
    assert report.updated_file_sha256_hash != sha256
    assert not report.message

Multi-column filters and sort options


  • Case sensitive or case insensitive multi-column sort is supported.

    • Multi-column sorts can be defined with combination of ascending and descending orders. For example; You can sort 'Parameter Value[Gender]' by ascending and Parameter Value[Age] by descending order.
    • Columns can be sorted as different data type. Supported sort data types are: str, int, float and datetime. datetime pattern can be defined.
    • Sort orders for invalid and empty values can also be defined. For example, If it is defined as VALID_EMPTY_INVALID, invalid values will be at the end. Empty values will follow valid values. This value order option is applicable for only int, datetime and float data types. All sort placement combinations are poossible for EMPTY, INVALID, VALID values.
    • You can define your custom sorters. A custom "enum-sorter" sorter has been already implemented. It sorts enums with given string values.
  • There are 10 different filters (plus (NOT) options of them). Any filter can be applied to any column. Multiple filters can be defined.

    • CONTAINS / NOT CONTAINS
    • EQUAL / NOT EQUAL
    • STARTSWITH / NOT STARTSWITH
    • ENDSWITH / NOT ENDSWITH
    • GREATER / NOT GREATER
    • GREATER_EQUAL / NOT GREATER_EQUAL
    • LESS / NOT LESS
    • LESS_EQUAL / NOT LESS_EQUAL
    • REGEX (regex match) / NOT REGEX (not regex match)
    • EMPTY / NOT EMPTY (None or empty)
  • You can define multiple filters. If one filter rejects row, row will not be selected (AND operation).

  • You can define one or more columns for a filter. If there are multiple columns for a filter. If any column matches, the filter selects the row (OR operation).

  • If you do not select any column for a filter, the filter will evaluate all columns. If filter matches with any column, it will select the row. Moreover, you can define some column names to skip them while filter is evaluating a row.

  • You can define your custom filters. Some custom filters have been already implemented.

    • "between-equal": Returns row if value between given min and max. Min and max inputs can be datetime, str, int or float.
    • "valid-datetime" Return row if value is valid datetime with given pattern. Default pattern is DD/MM/YYYY.
    • "valid-number": Return row if value is valid int or float.
    • "enum-contains": Gets a map to define a text for each enum value. Returns row if input parameter is in the enum text map. Enum data typese can be any allowed type (str, int, etc.).
      • Example: Enum values are 1, 2, 3, 4 (You store status values as in on database). Enum values are mapped to 1: "In Review", 2: "Published", 3: "In Curation", 4: "Public". If parameter is "Pub", all rows contain enum value 2 or 4 will be returned.

Example

Users can apply multiple filters and sort operations before retriving ISA table rows.

import pathlib

from metabolights_utils.isatab import Reader
from metabolights_utils.isatab.reader import (
    IsaTableFileReader,
    IsaTableFileReaderResult
)
from metabolights_utils.models.isa.common import (
    FilterOperation,
    SortType,
    TsvFileFilterOption,
    TsvFileSortOption,
    TsvFileSortValueOrder)


def test_with_filter_and_sort_option_01():
    # Sample Name value does not start with 'control'  and
    # Parameter Value[Chromatography Instrument] value is 'Thermo Scientific TRACE GC Ultra'.
    # Both filters are applied in case insesitive mode.
    filter_options = [
        TsvFileFilterOption(
            search_columns=["Sample Name"],
            operation=FilterOperation.STARTSWITH,
            parameter="control",
            case_sensitive=False,
            negate_result=True,
        ),
        TsvFileFilterOption(
            search_columns=["Parameter Value[Chromatography Instrument]"],
            operation=FilterOperation.EQUAL,
            parameter="Thermo Scientific TRACE GC Ultra",
            case_sensitive=False,
        ),
    ]

    # sort by Sample Name and 'Parameter Value[Chromatography Instrument]'
    sort_options = [
        TsvFileSortOption(column_name="Sample Name", reverse=False),
        TsvFileSortOption(
            column_name="Parameter Value[Chromatography Instrument]",
            column_sort_type=SortType.STRING,
            reverse=True,
        ),
    ]

    file_path = pathlib.Path(
        "tests/test-data/MTBLS66/a_MTBLS66_GC_metabolite_profiling_mass_spectrometry.txt"
    )
    helper: IsaTableFileReader = Reader.get_assay_file_reader()

    # Result pages will contain second 111 rows (page 2) after filter and sort operations.
    # Selected columns are not set. Result will contain all columns.
    result: IsaTableFileReaderResult = helper.get_page(
        page=2,
        results_per_page=111,
        file_buffer_or_path=file_path,
        sort_options=sort_options,
        filter_options=filter_options,
        selected_columns=None,
    )
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 111

    selected_columns = [
        "Sample Name",
        "Parameter Value[Column model]",
        "Parameter Value[Column type]",
    ]

    # Result pages will contain third 50 rows (page 3) after filter and sort operations.
    # 3 columns are selected. Result will contain 3 selected columns.
    result: IsaTableFileReaderResult = helper.get_page(
        page=3,
        results_per_page=50,
        file_buffer_or_path=file_path,
        sort_options=sort_options,
        filter_options=filter_options,
        selected_columns=selected_columns,
    )
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 50

    # First filter applies regex epression math on Sample Name column in case insensitive mode
    # Second filter is exact match on Parameter Value[Chromatography Instrument]
    filter_options = [
        TsvFileFilterOption(
            search_columns=["Sample Name"],
            operation=FilterOperation.REGEX,
            parameter="^PG[\d]5.*_5$",
            case_sensitive=False,
        ),
        TsvFileFilterOption(
            search_columns=["Parameter Value[Chromatography Instrument]"],
            operation=FilterOperation.EQUAL,
            parameter="Thermo Scientific TRACE GC Ultra",
            case_sensitive=False,
        ),
    ]

    selected_columns = [
        "Sample Name",
        "Parameter Value[Column model]",
        "Parameter Value[Column type]",
        "Parameter Value[Autosampler model]",
    ]

    # Result pages will be 111 and read 2. page after filter operations (No sort options).
    # Selected columns are set. Result will contain 6 columns.
    # 4 selected columns + Term Source REF and Term Accession Number columns of Parameter Value[Autosampler model].
    result: IsaTableFileReaderResult = helper.get_page(
        page=1,
        results_per_page=111,
        file_buffer_or_path=file_path,
        selected_columns=selected_columns,
        filter_options=filter_options,
    )
    assert len(result.parser_report.messages) == 0
    assert result.isa_table_file.table.row_count == 60

Update ISA table files with actions

User can manuplate ISA table files in row, column or cell level with actions. Supported actions:

  • ADD_ROW: Insert rows to given indices
  • DELETE_ROW: delete selected rows
  • MOVE_ROW: move row to new index
  • ADD_COLUMN: Add new columns
  • DELETE_COLUMN: delete selected columns
  • MOVE_COLUMN: move column to new index
  • COPY_ROW: copy row data to other selected rows
  • COPY_COLUMN: copy column data to other selected columns
  • UPDATE_ROW_DATA: update selected row values
  • UPDATE_COLUMN_DATA: update selected column values
  • UPDATE_COLUMN_HEADER: update column headers
  • UPDATE_CELL_DATA: update cells given with row and column index

View model and action definitions on this file.

View examples on this file.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

metabolights_utils-1.1.3.tar.gz (249.6 kB view hashes)

Uploaded Source

Built Distribution

metabolights_utils-1.1.3-py3-none-any.whl (150.5 kB view hashes)

Uploaded Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page