msreport 0.0.29__tar.gz → 0.0.30__tar.gz

{msreport-0.0.29 → msreport-0.0.30}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: msreport
-Version: 0.0.29
+Version: 0.0.30
 Summary: Post processing and analysis of quantitative proteomics data
 Author-email: "David M. Hollenstein" <hollenstein.david@gmail.com>
 License-Expression: Apache-2.0
@@ -29,7 +29,7 @@ Requires-Dist: seaborn>=0.12.0
 Requires-Dist: statsmodels>=0.13.2
 Requires-Dist: typing_extensions>=4
 Provides-Extra: r
-Requires-Dist: rpy2!=3.5.13,>=3.5.3; extra == "r"
+Requires-Dist: rpy2<3.5.13,>=3.5.3; extra == "r"
 Provides-Extra: dev
 Requires-Dist: mypy>=1.15.0; extra == "dev"
 Requires-Dist: pytest>=8.3.5; extra == "dev"
@@ -40,6 +40,7 @@ Dynamic: license-file
 # MsReport
 
 [![Project Status: WIP – Initial development is in progress, but there has not yet been a stable, usable release suitable for the public.](https://www.repostatus.org/badges/latest/wip.svg)](https://www.repostatus.org/#wip)
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.15309090.svg)](https://doi.org/10.5281/zenodo.15309090)
 ![Python Version from PEP 621 TOML](https://img.shields.io/python/required-version-toml?tomlFilePath=https%3A%2F%2Fraw.githubusercontent.com%2Fhollenstein%2Fmsreport%2Fmain%2Fpyproject.toml)
 [![Run tests](https://github.com/hollenstein/msreport/actions/workflows/run-tests.yml/badge.svg)](https://github.com/hollenstein/msreport/actions/workflows/run-tests.yml)
 
@@ -55,6 +56,7 @@ bottom-up mass spectrometry experiments.
     - [Additional requirements](#additional-requirements)
     - [Optional Dependencies](#optional-dependencies)
 - [Development status](#development-status)
+- [How to cite](#how-to-cite)
 
 ## What is MsReport?
 
@@ -134,3 +136,9 @@ For example, the R home directory might look like this on Windows: `C:\Program F
 ## Development status
 
 MsReport is a stable and reliable library that has been used on a daily basis for over two years in the Mass Spectrometry Facility at the Max Perutz Labs and the Mass Spectrometry Facility of IMP/IMBA/GMI. While the current interface of MsReport is stable, the library is still under active development, with new features being added regularly. Please note that a major rewrite is planned, which may introduce changes to the API in the future.
+
+## How to cite
+
+If you use MsReport for your research or publications, please include the following citation and consider giving the project a star on GitHub.
+
+> Hollenstein, D. M., & Hartl, M. (2025). hollenstein/msreport: v0.0.29 (0.0.29). Zenodo. https://doi.org/10.5281/zenodo.15309090

{msreport-0.0.29 → msreport-0.0.30}/README.md

@@ -1,6 +1,7 @@
 # MsReport
 
 [![Project Status: WIP – Initial development is in progress, but there has not yet been a stable, usable release suitable for the public.](https://www.repostatus.org/badges/latest/wip.svg)](https://www.repostatus.org/#wip)
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.15309090.svg)](https://doi.org/10.5281/zenodo.15309090)
 ![Python Version from PEP 621 TOML](https://img.shields.io/python/required-version-toml?tomlFilePath=https%3A%2F%2Fraw.githubusercontent.com%2Fhollenstein%2Fmsreport%2Fmain%2Fpyproject.toml)
 [![Run tests](https://github.com/hollenstein/msreport/actions/workflows/run-tests.yml/badge.svg)](https://github.com/hollenstein/msreport/actions/workflows/run-tests.yml)
 
@@ -16,6 +17,7 @@ bottom-up mass spectrometry experiments.
     - [Additional requirements](#additional-requirements)
     - [Optional Dependencies](#optional-dependencies)
 - [Development status](#development-status)
+- [How to cite](#how-to-cite)
 
 ## What is MsReport?
 
@@ -95,3 +97,9 @@ For example, the R home directory might look like this on Windows: `C:\Program F
 ## Development status
 
 MsReport is a stable and reliable library that has been used on a daily basis for over two years in the Mass Spectrometry Facility at the Max Perutz Labs and the Mass Spectrometry Facility of IMP/IMBA/GMI. While the current interface of MsReport is stable, the library is still under active development, with new features being added regularly. Please note that a major rewrite is planned, which may introduce changes to the API in the future.
+
+## How to cite
+
+If you use MsReport for your research or publications, please include the following citation and consider giving the project a star on GitHub.
+
+> Hollenstein, D. M., & Hartl, M. (2025). hollenstein/msreport: v0.0.29 (0.0.29). Zenodo. https://doi.org/10.5281/zenodo.15309090

{msreport-0.0.29 → msreport-0.0.30}/msreport/__init__.py

@@ -8,4 +8,4 @@ from msreport.fasta import import_protein_database
 from msreport.qtable import Qtable
 from msreport.reader import FragPipeReader, MaxQuantReader, SpectronautReader
 
-__version__ = "0.0.29"
+__version__ = "0.0.30"

{msreport-0.0.29 → msreport-0.0.30}/msreport/plot/distribution.py

@@ -204,7 +204,8 @@ def experiment_ratios(
     mask = np.all([(qtable.data[f"Events {exp}"] > 0) for exp in experiments], axis=0)
     if exclude_invalid:
         mask = mask & qtable["Valid"]
-    experiment_data = experiment_data[mask]
+    # Use `mask.to_numpy` to solve issue with different indices of mask and dataframe
+    experiment_data = experiment_data[mask.to_numpy()]
     pseudo_reference = np.nanmean(experiment_data, axis=1)
     ratio_data = experiment_data.subtract(pseudo_reference, axis=0)
 

{msreport-0.0.29 → msreport-0.0.30}/msreport/qtable.py

@@ -27,13 +27,11 @@ class Qtable:
         design: A pandas.DataFrame describing the experimental design.
     """
 
-    _default_id_column = "Representative protein"
-
     def __init__(
         self,
         data: pd.DataFrame,
-        design: Optional[pd.DataFrame] = None,
-        id_column: str = "Representative protein",
+        design: pd.DataFrame,
+        id_column: str,
     ):
         """Initializes the Qtable.
 
@@ -42,12 +40,13 @@ class Qtable:
 
         Args:
             data: A dataframe containing quantitative proteomics data in a wide format.
+                The index of the dataframe must contain unique values.
             design: A dataframe describing the experimental design that must at least
                 contain the columns "Sample" and "Experiment". The "Sample" entries
                 should correspond to the Sample names present in the quantitative
                 columns of the data.
             id_column: The name of the column that contains the unique identifiers for
-                the entries in the data table. Default is "Representative protein".
+                the entries in the data table.
 
         Raises:
             KeyError: If the specified id_column is not found in data.
@@ -76,8 +75,7 @@ class Qtable:
         self._id_column = id_column
         if "Valid" not in self.data.columns:
             self.data["Valid"] = True
-        if design is not None:
-            self.add_design(design)
+        self.add_design(design)
 
         self._expression_columns: list[str] = []
         self._expression_features: list[str] = []
@@ -438,6 +436,11 @@ class Qtable:
 
         Returns:
             An instance of Qtable loaded from the specified files.
+
+        Raises:
+            ValueError: If the loaded config file does not contain the
+                "Unique ID column" key. This is due to the qtable being saved with a
+                version of msreport <= 0.0.27.
         """
         filepaths = _get_qtable_export_filepaths(directory, basename)
         with open(filepaths["config"]) as openfile:
@@ -458,13 +461,20 @@ class Qtable:
                 filepaths["design"], sep="\t", index_col=0, keep_default_na=True
             )
 
-        qtable = Qtable(data, design)
+        if "Unique ID column" not in config_data:
+            # Mention that the qtable was likely saved with a version of msreport <= 0.0.27
+            raise ValueError(
+                "The qtable config file does not contain the 'Unique ID column' key. "
+                "This is likely due to the qtable being saved with a version of "
+                "msreport <= 0.0.27."
+            )
+        id_column = config_data["Unique ID column"]
+
+        qtable = Qtable(data, design, id_column)
         qtable._expression_columns = config_data["Expression columns"]
         qtable._expression_features = config_data["Expression features"]
         qtable._expression_sample_mapping = config_data["Expression sample mapping"]
         # This check is required for backwards compatibility with msreport <= 0.0.27
-        if "Unique ID column" in config_data:
-            qtable._id_column = config_data["Unique ID column"]
         return qtable
 
     def to_tsv(self, path: str, index: bool = False):
@@ -570,7 +580,7 @@ class Qtable:
         self._expression_sample_mapping = {}
 
     def __copy__(self) -> Qtable:
-        new_instance = Qtable(self.data, self.design)
+        new_instance = Qtable(self.data, self.design, self.id_column)
         # Copy all private attributes
         for attr in dir(self):
             if (

{msreport-0.0.29 → msreport-0.0.30}/msreport/reader.py

@@ -541,6 +541,8 @@ class FragPipeReader(ResultReader):
     """FragPipe result reader.
 
     Methods:
+        import_design: Reads a "fragpipe-files.fp-manifest" file and returns a
+            processed design dataframe.
         import_proteins: Reads a "combined_protein.tsv" or "protein.tsv" file and
             returns a processed dataframe, conforming to the MsReport naming
             convention.
@@ -583,12 +585,19 @@ class FragPipeReader(ResultReader):
         "ions": "combined_ion.tsv",
         "ion_evidence": "ion.tsv",
         "psm_evidence": "psm.tsv",
+        "design": "fragpipe-files.fp-manifest",
     }
     isobar_filenames: dict[str, str] = {
         "proteins": "protein.tsv",
         "peptides": "peptide.tsv",
         "ions": "ion.tsv",
     }
+    sil_filenames: dict[str, str] = {
+        "proteins": "combined_protein_label_quant.tsv",
+        "peptides": "combined_modified_peptide_label_quant.tsv",
+        "ions": "combined_ion_label_quant.tsv",
+    }
+
     protected_columns: list[str] = []
     sample_column_tags: list[str] = [
         "Spectral Count",
@@ -609,6 +618,7 @@ class FragPipeReader(ResultReader):
         "Modified Sequence": "Modified sequence",  # Modified peptide and ion
         "Start": "Start position",  # Peptide and ion
         "End": "End position",  # Peptide and ion
+        "Mapped Proteins": "Mapped proteins",  # All PSM, ion, and peptide tables
         "Combined Total Peptides": "Total peptides",  # From LFQ
         "Total Peptides": "Total peptides",  # From TMT
         "Description": "Protein name",
@@ -638,7 +648,11 @@ class FragPipeReader(ResultReader):
     protein_info_tags: list[str] = []
 
     def __init__(
-        self, directory: str, isobar: bool = False, contaminant_tag: str = "contam_"
+        self,
+        directory: str,
+        isobar: bool = False,
+        sil: bool = False,
+        contaminant_tag: str = "contam_",
     ) -> None:
         """Initializes the FragPipeReader.
 
@@ -646,16 +660,69 @@ class FragPipeReader(ResultReader):
             directory: Location of the FragPipe result folder
             isobar: Set to True if quantification strategy was TMT, iTRAQ or similar;
                 default False.
+            sil: Set to True if the FragPipe result files are from a stable isotope
+                labeling experiment, such as SILAC; default False.
             contaminant_tag: Prefix of Protein ID entries to identify contaminants;
                 default "contam_".
         """
+        if sil and isobar:
+            raise ValueError("Cannot set both 'isobar' and 'sil' to True.")
         self._add_data_directory(directory)
         self._isobar: bool = isobar
+        self._sil: bool = sil
         self._contaminant_tag: str = contaminant_tag
-        if not isobar:
+        if isobar:
+            self.filenames = self.isobar_filenames
+        elif sil:
+            self.filenames = self.sil_filenames
+        else:
             self.filenames = self.default_filenames
+
+    def import_design(
+        self, filename: Optional[str] = None, sort: bool = False
+    ) -> pd.DataFrame:
+        """Reads a 'fp-manifest' file and returns a processed design dataframe.
+
+        Args:
+            filename: Allows specifying an alternative filename, otherwise the default
+                filename is used.
+            sort: If True, the design dataframe is sorted by "Experiment" and
+                "Replicate"; default False.
+
+        Returns:
+            A dataframe containing the processed design table with columns:
+            "Sample", "Experiment", "Replicate", "Rawfile".
+
+        Raises:
+            FileNotFoundError: If the specified manifest file does not exist.
+        """
+        if filename is None:
+            filepath = os.path.join(self.data_directory, self.filenames["design"])
         else:
-            self.filenames = self.isobar_filenames
+            filepath = os.path.join(self.data_directory, filename)
+        if not os.path.exists(filepath):
+            raise FileNotFoundError(
+                f"File '{filepath}' does not exist. Please check the file path."
+            )
+        fp_manifest = pd.read_csv(filepath, sep="\t", header=None, dtype=str)
+        fp_manifest.columns = ["Path", "Experiment", "Bioreplicate", "Data type"]
+
+        design = pd.DataFrame(
+            {
+                "Sample": fp_manifest["Experiment"] + "_" + fp_manifest["Bioreplicate"],
+                "Experiment": fp_manifest["Experiment"],
+                "Replicate": fp_manifest["Bioreplicate"],
+                "Rawfile": fp_manifest["Path"].apply(
+                    # Required to handle Windows and Unix style paths on either system
+                    lambda x: x.replace("\\", "/").split("/")[-1]
+                ),
+            }
+        )
+
+        if sort:
+            design.sort_values(by=["Experiment", "Replicate"], inplace=True)
+            design.reset_index(drop=True, inplace=True)
+        return design
 
     def import_proteins(
         self,
@@ -737,6 +804,7 @@ class FragPipeReader(ResultReader):
         df = self._read_file("peptides" if filename is None else filename)
         df["Protein reported by software"] = _extract_protein_ids(df["Protein"])
         df["Representative protein"] = df["Protein reported by software"]
+        df["Mapped Proteins"] = self._collect_mapped_proteins(df)
         # Note that _add_protein_entries would need to be adapted for the peptide table.
         # df = self._add_protein_entries(df)
         if rename_columns:
@@ -793,6 +861,8 @@ class FragPipeReader(ResultReader):
         #         'Indistinguishable Proteins' to the ion table.
         df["Protein reported by software"] = _extract_protein_ids(df["Protein"])
         df["Representative protein"] = df["Protein reported by software"]
+        df["Mapped Proteins"] = self._collect_mapped_proteins(df)
+
         if rename_columns:
             df = self._rename_columns(df, prefix_column_tags)
         if rewrite_modifications and rename_columns:
@@ -879,6 +949,8 @@ class FragPipeReader(ResultReader):
         #         'Indistinguishable Proteins' to the ion table.
         df["Protein reported by software"] = _extract_protein_ids(df["Protein"])
         df["Representative protein"] = df["Protein reported by software"]
+        df["Mapped Proteins"] = self._collect_mapped_proteins(df)
+
         if rename_columns:
             df = self._rename_columns(df, prefix_column_tags)
         if rewrite_modifications and rename_columns:
@@ -928,23 +1000,7 @@ class FragPipeReader(ResultReader):
 
         df["Protein reported by software"] = _extract_protein_ids(df["Protein"])
         df["Representative protein"] = df["Protein reported by software"]
-        df["Mapped Proteins"] = df["Mapped Proteins"].astype(str).replace("nan", "")
-
-        # FP only lists additional mapped proteins in the "Mapped Proteins" column
-        # MsReport reports all matching proteins in the "Mapped proteins" column
-        mapped_proteins_entries = []
-        for protein, mapped_protein_fp in zip(
-            df["Representative protein"], df["Mapped Proteins"], strict=True
-        ):
-            if mapped_protein_fp == "":
-                mapped_proteins = [protein]
-            else:
-                additional_mapped_proteins = msreport.reader._extract_protein_ids(
-                    mapped_protein_fp.split(", ")
-                )
-                mapped_proteins = [protein] + additional_mapped_proteins
-            mapped_proteins_entries.append(";".join(mapped_proteins))
-        df["Mapped proteins"] = mapped_proteins_entries
+        df["Mapped Proteins"] = self._collect_mapped_proteins(df)
 
         if rename_columns:
             df = self._rename_columns(df, prefix_tag=True)
@@ -980,6 +1036,35 @@ class FragPipeReader(ResultReader):
             df[key] = protein_entry_table[key]
         return df
 
+    def _collect_mapped_proteins(self, df: pd.DataFrame) -> list[str]:
+        """Generates a list of mapped proteins entries.
+
+        This method extracts protein IDs from the 'Representative protein' and the
+        'Mapped Proteins' column and combines them into a single string for each row,
+        where multiple protein IDs are separated by semicolons.
+
+        Args:
+            df: DataFrame containing the 'Mapped Proteins' column.
+
+        Returns:
+            A list of mapped proteins entries.
+        """
+        mapped_proteins_entries = []
+        for protein, mapped_protein_fp in zip(
+            df["Representative protein"],
+            df["Mapped Proteins"].astype(str).replace("nan", ""),
+            strict=True,
+        ):
+            if mapped_protein_fp == "":
+                mapped_proteins = [protein]
+            else:
+                additional_mapped_proteins = msreport.reader._extract_protein_ids(
+                    mapped_protein_fp.split(", ")
+                )
+                mapped_proteins = [protein] + additional_mapped_proteins
+            mapped_proteins_entries.append(";".join(mapped_proteins))
+        return mapped_proteins_entries
+
     def _collect_leading_protein_entries(self, df: pd.DataFrame) -> list[list[str]]:
         """Generates a list of leading protein entries.
 
@@ -995,6 +1080,9 @@ class FragPipeReader(ResultReader):
             A list of the same length as the input dataframe. Each position contains a
             list of leading protein entries, which a minimum of one entry.
         """
+        if self._sil:  # No "Indistinguishable Proteins" columns in 'SIL' data
+            return [[p] for p in df["Protein"]]
+
         leading_protein_entries = []
         for protein_entry, indist_protein_entry in zip(
             df["Protein"], df["Indistinguishable Proteins"].fillna("").astype(str)

{msreport-0.0.29 → msreport-0.0.30}/msreport.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: msreport
-Version: 0.0.29
+Version: 0.0.30
 Summary: Post processing and analysis of quantitative proteomics data
 Author-email: "David M. Hollenstein" <hollenstein.david@gmail.com>
 License-Expression: Apache-2.0
@@ -29,7 +29,7 @@ Requires-Dist: seaborn>=0.12.0
 Requires-Dist: statsmodels>=0.13.2
 Requires-Dist: typing_extensions>=4
 Provides-Extra: r
-Requires-Dist: rpy2!=3.5.13,>=3.5.3; extra == "r"
+Requires-Dist: rpy2<3.5.13,>=3.5.3; extra == "r"
 Provides-Extra: dev
 Requires-Dist: mypy>=1.15.0; extra == "dev"
 Requires-Dist: pytest>=8.3.5; extra == "dev"
@@ -40,6 +40,7 @@ Dynamic: license-file
 # MsReport
 
 [![Project Status: WIP – Initial development is in progress, but there has not yet been a stable, usable release suitable for the public.](https://www.repostatus.org/badges/latest/wip.svg)](https://www.repostatus.org/#wip)
+[![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.15309090.svg)](https://doi.org/10.5281/zenodo.15309090)
 ![Python Version from PEP 621 TOML](https://img.shields.io/python/required-version-toml?tomlFilePath=https%3A%2F%2Fraw.githubusercontent.com%2Fhollenstein%2Fmsreport%2Fmain%2Fpyproject.toml)
 [![Run tests](https://github.com/hollenstein/msreport/actions/workflows/run-tests.yml/badge.svg)](https://github.com/hollenstein/msreport/actions/workflows/run-tests.yml)
 
@@ -55,6 +56,7 @@ bottom-up mass spectrometry experiments.
     - [Additional requirements](#additional-requirements)
     - [Optional Dependencies](#optional-dependencies)
 - [Development status](#development-status)
+- [How to cite](#how-to-cite)
 
 ## What is MsReport?
 
@@ -134,3 +136,9 @@ For example, the R home directory might look like this on Windows: `C:\Program F
 ## Development status
 
 MsReport is a stable and reliable library that has been used on a daily basis for over two years in the Mass Spectrometry Facility at the Max Perutz Labs and the Mass Spectrometry Facility of IMP/IMBA/GMI. While the current interface of MsReport is stable, the library is still under active development, with new features being added regularly. Please note that a major rewrite is planned, which may introduce changes to the API in the future.
+
+## How to cite
+
+If you use MsReport for your research or publications, please include the following citation and consider giving the project a star on GitHub.
+
+> Hollenstein, D. M., & Hartl, M. (2025). hollenstein/msreport: v0.0.29 (0.0.29). Zenodo. https://doi.org/10.5281/zenodo.15309090

{msreport-0.0.29 → msreport-0.0.30}/msreport.egg-info/requires.txt

@@ -12,7 +12,7 @@ statsmodels>=0.13.2
 typing_extensions>=4
 
 [R]
-rpy2!=3.5.13,>=3.5.3
+rpy2<3.5.13,>=3.5.3
 
 [dev]
 mypy>=1.15.0

{msreport-0.0.29 → msreport-0.0.30}/pyproject.toml

@@ -47,7 +47,7 @@ changelog = "https://github.com/hollenstein/msreport/blob/main/CHANGELOG.md"
 
 [project.optional-dependencies]
 R = [
-    "rpy2>=3.5.3,!=3.5.13",
+    "rpy2>=3.5.3,<3.5.13",
 ]
 dev = [
     "mypy>=1.15.0",

{msreport-0.0.29 → msreport-0.0.30}/tests/test_analyze.py

@@ -46,7 +46,7 @@ def example_data():
 
 @pytest.fixture
 def example_qtable(example_data):
-    qtable = msreport.qtable.Qtable(example_data["data"], design=example_data["design"])
+    qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
     qtable.set_expression_by_tag("Intensity")
     return qtable
 

{msreport-0.0.29 → msreport-0.0.30}/tests/test_plot.py

@@ -49,11 +49,21 @@ def example_data():
 
 @pytest.fixture
 def example_qtable(example_data):
-    qtable = msreport.qtable.Qtable(example_data["data"], design=example_data["design"])
+    qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
     qtable.set_expression_by_tag("Intensity")
     return qtable
 
 
+class TestExperimentRatios:
+    @pytest.fixture(autouse=True)
+    def _init_qtable(self, example_qtable):
+        self.qtable = example_qtable
+
+    def test_non_default_df_index_does_not_raise_error(self):
+        self.qtable.data.index = range(2, len(self.qtable.data) + 2)
+        fig, axes = msreport.plot.experiment_ratios(self.qtable)
+
+
 class TestVolcanoMa:
     @pytest.fixture(autouse=True)
     def _init_qtable(self, example_qtable):

{msreport-0.0.29 → msreport-0.0.30}/tests/test_qtable.py

@@ -68,43 +68,41 @@ def example_data():
 
 @pytest.fixture
 def example_qtable(example_data):
-    qtable = msreport.qtable.Qtable(
-        example_data["data"], design=example_data["design"], id_column="id"
-    )
+    qtable = msreport.qtable.Qtable(example_data["data"], design=example_data["design"], id_column="id")  # fmt: skip
     qtable.set_expression_by_tag("Intensity")
     return qtable
 
 
 class TestQtableInitialization:
     def test_data_is_added_to_qtable(self, example_data):
-        qtable = msreport.qtable.Qtable(example_data["data"])
+        qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
         assert qtable.data.equals(example_data["data"])
 
     def test_design_is_added_to_qtable(self, example_data):
-        qtable = msreport.qtable.Qtable(example_data["data"], design=example_data["design"])  # fmt: skip
+        qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
         assert qtable.design.equals(example_data["design"])
 
     def test_id_column_is_added_to_qtable(self, example_data):
-        qtable = msreport.qtable.Qtable(example_data["data"], id_column="id")
-        assert qtable.id_column == "id"
+        qtable = msreport.qtable.Qtable(example_data["data"], design=example_data["design"], id_column="id")  # fmt: skip
+        assert qtable._id_column == "id"
 
     def test_non_unique_data_index_raises_error(self, example_data):
         example_data["data"].index = [0 for _ in range(len(example_data["data"]))]
         with pytest.raises(ValueError):
-            msreport.qtable.Qtable(example_data["data"])
+            msreport.qtable.Qtable(example_data["data"], design=example_data["design"], id_column="Representative protein")  # fmt: skip
 
     def test_non_existing_id_column_raises_error(self, example_data):
         with pytest.raises(KeyError):
-            msreport.qtable.Qtable(example_data["data"], id_column="non_existing_column")  # fmt: skip
+            msreport.qtable.Qtable(example_data["data"], design=example_data["design"], id_column="non_existing_column")  # fmt: skip
 
     def test_id_column_containing_non_unique_values_raises_error(self, example_data):
         example_data["data"]["id"] = "1"
         with pytest.raises(ValueError):
-            msreport.qtable.Qtable(example_data["data"], id_column="id")
+            msreport.qtable.Qtable(example_data["data"], design=example_data["design"], id_column="id")  # fmt: skip
 
 
 def test_qtable_add_design(example_data):
-    qtable = msreport.qtable.Qtable(example_data["data"])
+    qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
     qtable.add_design(example_data["design"])
     assert qtable.design.equals(example_data["design"])
 
@@ -192,9 +190,7 @@ class TestMatchSamplesToTagColumns:
 class TestQtableGetData:
     @pytest.fixture(autouse=True)
     def _init_qtable(self, example_data):
-        self.qtable = msreport.qtable.Qtable(
-            example_data["data"], design=example_data["design"]
-        )
+        self.qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
 
     def test_get_data(self, example_data):
         assert self.qtable.get_data().equals(example_data["data"])
@@ -229,7 +225,7 @@ def test_qtable_contains(example_qtable, key, is_present):
 
 
 def test_qtable_get_design(example_data):
-    qtable = msreport.qtable.Qtable(example_data["data"], design=example_data["design"])
+    qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
     assert qtable.get_design().equals(example_data["design"])
 
 
@@ -265,9 +261,7 @@ class TestQtableGetExperiments:
 class TestQtableResetExpression:
     @pytest.fixture(autouse=True)
     def _init_qtable(self, example_data):
-        self.qtable = msreport.qtable.Qtable(
-            example_data["data"], design=example_data["design"]
-        )
+        self.qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
 
     def test_reset_of_parameters(self):
         self.qtable._expression_columns = ["test"]
@@ -285,7 +279,7 @@ class TestQtableResetExpression:
         self.qtable._reset_expression()
         data_columns = self.qtable.data.columns
         all_expression_columns_absent_in_data = not any(
-            [c in data_columns for c in example_data["expression_columns"]]
+            c in data_columns for c in example_data["expression_columns"]
         )
         assert all_expression_columns_absent_in_data
 
@@ -303,9 +297,7 @@ class TestQtableResetExpression:
 class TestQtableSetExpression:
     @pytest.fixture(autouse=True)
     def _init_qtable(self, example_data):
-        self.qtable = msreport.qtable.Qtable(
-            example_data["data"], design=example_data["design"]
-        )
+        self.qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
 
     def test_correct_setting_of_private_variables(self, example_data):
         self.qtable._set_expression(example_data["intensity_cols_to_samples"])
@@ -377,9 +369,7 @@ class TestQtableSetExpression:
 class TestQtableSetExpressionByTag:
     @pytest.fixture(autouse=True)
     def _init_qtable(self, example_data):
-        self.qtable = msreport.qtable.Qtable(
-            example_data["data"], design=example_data["design"]
-        )
+        self.qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
 
     def test_set_expression_by_tag(self, example_data):
         self.qtable.set_expression_by_tag(example_data["expression_tag"])
@@ -413,9 +403,7 @@ class TestQtableSetExpressionByTag:
 class TestQtableSetExpressionByColumn:
     @pytest.fixture(autouse=True)
     def _init_qtable(self, example_data):
-        self.qtable = msreport.qtable.Qtable(
-            example_data["data"], design=example_data["design"]
-        )
+        self.qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
 
     def test_set_expression_by_column(self, example_data):
         self.qtable.set_expression_by_column(example_data["intensity_cols_to_samples"])
@@ -450,9 +438,7 @@ class TestQtableSetExpressionByColumn:
 class TestQtableAddExpressionFeature:
     @pytest.fixture(autouse=True)
     def _init_qtable(self, example_data):
-        self.qtable = msreport.qtable.Qtable(
-            example_data["data"], design=example_data["design"]
-        )
+        self.qtable = msreport.qtable.Qtable(example_data["data"], example_data["design"], id_column="Representative protein")  # fmt: skip
 
     def test_with_series(self):
         new_data = self.qtable.data["id"].copy()