A Quantitative Multivariate Model of Human Dendritic Cell-T Helper Cell Communication.

Cell-cell communication involves a large number of molecular signals that function as words of a complex language whose grammar remains mostly unknown. Here, we describe an integrative approach involving (1) protein-level measurement of multiple communication signals coupled to output responses in receiving cells and (2) mathematical modeling to uncover input-output relationships and interactions between signals. Using human dendritic cell (DC)-T helper (Th) cell communication as a model, we measured 36 DC-derived signals and 17 Th cytokines broadly covering Th diversity in 428 observations. We developed a data-driven, computationally validated model capturing 56 already described and 290 potentially novel mechanisms of Th cell specification. By predicting context-dependent behaviors, we demonstrate a new function for IL-12p70 as an inducer of Th17 in an IL-1 signaling context. This work provides a unique resource to decipher the complex combinatorial rules governing DC-Th cell communication and guide their manipulation for vaccine design and immunotherapies.

Live-cell chromosome dynamics and outcome of X chromosome pairing events during ES cell differentiation.

Random X inactivation represents a paradigm for monoallelic gene regulation during early ES cell differentiation. In mice, the choice of X chromosome to inactivate in XX cells is ensured by monoallelic regulation of Xist RNA via its antisense transcription unit Tsix/Xite. Homologous pairing events have been proposed to underlie asymmetric Tsix expression, but direct evidence has been lacking owing to their dynamic and transient nature. Here we investigate the live-cell dynamics and outcome of Tsix pairing in differentiating mouse ES cells. We find an overall increase in genome dynamics including the Xics during early differentiation. During pairing, however, Xic loci show markedly reduced movements. Upon separation, Tsix expression becomes transiently monoallelic, providing a window of opportunity for monoallelic Xist upregulation. Our findings reveal the spatiotemporal choreography of the X chromosomes during early differentiation and indicate a direct role for pairing in facilitating symmetry-breaking and monoallelic regulation of Xist during random X inactivation.

A critical function for transforming growth factor-beta, interleukin 23 and proinflammatory cytokines in driving and modulating human T(H)-17 responses.

Interleukin 17 (IL-17)-producing T helper 17 cells (T(H)-17 cells) have been described as a T helper cell subset distinct from T helper type 1 (T(H)1) and T(H)2 cells, with specific functions in antimicrobial defense and autoimmunity. The factors driving human T(H)-17 differentiation remain controversial. Using a systematic approach combining experimental and computational methods, we show here that transforming growth factor-beta, interleukin 23 (IL-23) and proinflammatory cytokines (IL-1beta and IL-6) were all essential for human T(H)-17 differentiation. However, individual T(H)-17 cell-derived cytokines, such as IL-17, IL-21, IL-22 and IL-6, as well as the global T(H)-17 cytokine profile, were differentially modulated by T(H)-17-promoting cytokines. Transforming growth factor-beta was critical, and its absence induced a shift from a T(H)-17 profile to a T(H)1-like profile. Our results shed new light on the regulation of human T(H)-17 differentiation and provide a framework for the global analysis of T helper responses.

Analytical Platforms 1: Use of Cultured Cells and Fluorescent Read-Out Coupled to NormaCurve Normalization in RPPA.

The analytic platform described in this chapter uses proteins extracted from cultured cells as an infinite source of material to set up, validate, and quality control an RPPA platform. Readout of the arrays uses near-infrared fluorescence labeling and data normalization is performed using the bioinformatics package NormaCurve.In the first part, we will describe the advantages, drawbacks, and different applications of cell line material for RPPA. In the second part, we will describe how the staining protocol, the method of readout, and the normalization method applied afterward are interconnected and should be considered together. Finally, we will describe the NormaCurve package, which is freely available, and its requirements for implementation.Four protocols are provided in this chapter: (1) Protein lysis of cell lines using a homemade Laemmli buffer, (2) RPPA staining for fluorescent readout including a signal amplification step, (3) total protein staining in the visible spectrum for normalization purposes, and (4) total protein staining in the near-infrared spectrum for normalization purposes.

Precision medicine in cancer: challenges and recommendations from an EU-funded cervical cancer biobanking study.

BACKGROUND: Cervical cancer (CC) remains a leading cause of gynaecological cancer-related mortality worldwide. CC pathogenesis is triggered when human papillomavirus (HPV) inserts into the genome, resulting in tumour suppressor gene inactivation and oncogene activation. Collecting tumour and blood samples is critical for identifying these genetic alterations. METHODS: BIO-RAIDs is the first prospective molecular profiling clinical study to include a substantial biobanking effort that used uniform high-quality standards and control of samples. In this European Union (EU)-funded study, we identified the challenges that were impeding the effective implementation of such a systematic and comprehensive biobanking effort. RESULTS: The challenges included a lack of uniform international legal and ethical standards, complexities in clinical and molecular data management, and difficulties in determining the best technical platforms and data analysis techniques. Some difficulties were encountered by all investigators, while others affected only certain institutions, regions, or countries. CONCLUSIONS: The results of the BIO-RAIDs programme highlight the need to facilitate and standardise regulatory procedures, and we feel that there is also a need for international working groups that make recommendations to regulatory bodies, governmental funding agencies, and academic institutions to achieve a proficient biobanking programme throughout EU countries. This represents the first step in precision medicine.

Histo-genomic stratification reveals the frequent amplification/overexpression of CCNE1 and BRD4 genes in non-BRCAness high grade ovarian carcinoma.

The treatment of epithelial ovarian cancer (EOC) is narrowly focused despite the heterogeneity of this disease in which outcomes remain poor. To stratify EOC patients for targeted therapy, we developed an approach integrating expression and genomic analyses including the BRCAness status. Gene expression and genomic profiling were used to identify genes recurrently (>5%) amplified and overexpressed in 105 EOC. The LST (Large-scale State Transition) genomic signature of BRCAness was applied to define molecular subgroups of EOC. Amplified/overexpressed genes clustered mainly in 3q, 8q, 19p and 19q. These changes were generally found mutually exclusive. In the 85 patients for which the genomic signature could be determined, genomic BRCAness was found in 52 cases (61.1%) and non-BRCAness in 33 (38.8%). A striking mutual exclusivity was observed between BRCAness and amplification/overexpression data. Whereas 3q and 8q alterations were preferentially observed in BRCAness EOC, most alterations on chromosome 19 were in non-BRCAness cases. CCNE1 (19q12) and BRD4 (19p13.1) amplification/overexpression was found in 19/33 (57.5%) of non-BRCAness cases. Such disequilibrium was also found in the TCGA EOC data set used for validation. Potential target genes are frequently amplified/overexpressed in non-BRCAness EOC. We report that BRD4, already identified as a target in several tumor models, is a new potential target in high grade non-BRCAness ovarian carcinoma.

Multi-factor data normalization enables the detection of copy number aberrations in amplicon sequencing data.

MOTIVATION: Because of its low cost, amplicon sequencing, also known as ultra-deep targeted sequencing, is now becoming widely used in oncology for detection of actionable mutations, i.e. mutations influencing cell sensitivity to targeted therapies. Amplicon sequencing is based on the polymerase chain reaction amplification of the regions of interest, a process that considerably distorts the information on copy numbers initially present in the tumor DNA. Therefore, additional experiments such as single nucleotide polymorphism (SNP) or comparative genomic hybridization (CGH) arrays often complement amplicon sequencing in clinics to identify copy number status of genes whose amplification or deletion has direct consequences on the efficacy of a particular cancer treatment. So far, there has been no proven method to extract the information on gene copy number aberrations based solely on amplicon sequencing. RESULTS: Here we present ONCOCNV, a method that includes a multifactor normalization and annotation technique enabling the detection of large copy number changes from amplicon sequencing data. We validated our approach on high and low amplicon density datasets and demonstrated that ONCOCNV can achieve a precision comparable with that of array CGH techniques in detecting copy number aberrations. Thus, ONCOCNV applied on amplicon sequencing data would make the use of additional array CGH or SNP array experiments unnecessary.

From prospective biobanking to precision medicine: BIO-RAIDs - an EU study protocol in cervical cancer.

BACKGROUND: Cervical cancer (CC) is -second to breast cancer- a dominant cause of gynecological cancer-related deaths worldwide. CC tumor biopsies and blood samples are of easy access and vital for the development of future precision medicine strategies. DESIGN: BIO-RAIDs is a prospective multicenter European study, presently recruiting patients in 6 EU countries. Tumor and liquid biopsies from patients with previously non-treated cervical cancer (stages IB2-IV) are collected at defined time points. Patients receive standard primary treatment according to the stage of their disease. 700 patients are planned to be enrolled. The main objectives are the discovery of -dominant molecular alterations, -signalling pathway activation, and -tumor micro-environment patterns that may predict response or resistance to treatment. An exhaustive molecular analysis is performed using 1° Next generation sequencing, 2° Reverse phase protein arrays and 3° Immuno-histochemistry. DISCUSSION: The clinical study BIO-RAIDs is activated in all planned countries, 170 patients have been recruited till now. This study will make an important contribution towards precision medicine treatments in cervical cancer. The results will support the development of clinical practice guidelines for cervical cancer patients to improve their prognosis and their quality of life. TRIAL REGISTRATION: Clinicaltrials.gov: NCT02428842 , registered 10 February 2015.

Spi-1/PU.1 activates transcription through clustered DNA occupancy in erythroleukemia.

Acute leukemias are characterized by deregulation of transcriptional networks that control the lineage specificity of gene expression. The aberrant overexpression of the Spi-1/PU.1 transcription factor leads to erythroleukemia. To determine how Spi-1 mechanistically influences the transcriptional program, we combined a ChIP-seq analysis with transcriptional profiling in cells from an erythroleukemic mouse model. We show that Spi-1 displays a selective DNA-binding that does not often cause transcriptional modulation. We report that Spi-1 controls transcriptional activation and repression partially through distinct Spi-1 recruitment to chromatin. We revealed several parameters impacting on Spi-1-mediated transcriptional activation. Gene activation is facilitated by Spi-1 occupancy close to transcriptional starting site of genes devoid of CGIs. Moreover, in those regions Spi-1 acts by binding to multiple motifs tightly clustered and with similar orientation. Finally, in contrast to the myeloid and lymphoid B cells in which Spi-1 exerts a physiological activity, in the erythroleukemic cells, lineage-specific cooperating factors do not play a prevalent role in Spi-1-mediated transcriptional activation. Thus, our work describes a new mechanism of gene activation through clustered site occupancy of Spi-1 particularly relevant in regard to the strong expression of Spi-1 in the erythroleukemic cells.

Serum and Salivary IgG and IgA Response After COVID-19 Messenger RNA Vaccination.

Importance: There is still considerable controversy in the literature regarding the capacity of intramuscular messenger RNA (mRNA) vaccination to induce a mucosal immune response. Objective: To compare serum and salivary IgG and IgA levels among mRNA-vaccinated individuals with or without previous SARS-CoV-2 infection. Design, Setting, and Participants: In this cohort study, SARS-CoV-2-naive participants and those with previous infection were consecutively included in the CoviCompare P and CoviCompare M mRNA vaccination trials and followed up to day 180 after vaccination with either the BNT162b2 (Pfizer-BioNTech) vaccine or the mRNA-1273 (Moderna) vaccine at the beginning of the COVID-19 vaccination campaign (from February 19 to June 8, 2021) in France. Data were analyzed from October 25, 2022, to July 13, 2023. Main Outcomes and Measures: An ultrasensitive digital enzyme-linked immunosorbent assay was used for the comparison of SARS-CoV-2 spike-specific serum and salivary IgG and IgA levels. Spike-specific secretory IgA level was also quantified at selected times. Results: A total of 427 individuals were included in 3 groups: participants with SARS-CoV-2 prior to vaccination who received 1 single dose of BNT162b2 (Pfizer-BioNTech) (n = 120) and SARS-CoV-2-naive individuals who received 2 doses of mRNA-1273 (Moderna) (n = 172) or 2 doses of BNT162b2 (Pfizer-BioNTech) (n = 135). The median age was 68 (IQR, 39-75) years, and 228 (53.4%) were men. SARS-CoV-2 spike-specific IgG saliva levels increased after 1 or 2 vaccine injections in individuals with previous infection and SARS-CoV-2-naive individuals. After vaccination, SARS-CoV-2-specific saliva IgA levels, normalized with respect to total IgA levels, were significantly higher in participants with previous infection, as compared with the most responsive mRNA-1273 (Moderna) recipients (median normalized levels, 155 × 10-5 vs 37 × 10-5 at day 29; 107 × 10-5 vs 54 × 10-5 at day 57; and 104 × 10-5 vs 70 × 10-5 at day 180 [P < .001]). In contrast, compared with day 1, spike-specific IgA levels in the BNT162b2-vaccinated SARS-CoV-2-naive group increased only at day 57 (36 × 10-5 vs 49 × 10-5 [P = .01]). Bona fide multimeric secretory IgA levels were significantly higher in individuals with previous infection compared with SARS-CoV-2-naive individuals after 2 antigenic stimulations (median optical density, 0.36 [IQR, 0.16-0.63] vs 0.16 [IQR, 0.10-0.22]; P < .001). Conclusions and Relevance: The findings of this cohort study suggest that mRNA vaccination was associated with mucosal immunity in individuals without prior SARS-CoV-2 infection, but at much lower levels than in previously infected individuals. Further studies are needed to determine the association between specific saliva IgA levels and prevention of infection or transmission.

Telomere crisis in kidney epithelial cells promotes the acquisition of a microRNA signature retrieved in aggressive renal cell carcinomas.

Telomere shortening is a major source of chromosome instability (CIN) at early stages during carcinogenesis. However, the mechanisms through which telomere-driven CIN (T-CIN) contributes to the acquisition of tumor phenotypes remain uncharacterized. We discovered that human epithelial kidney cells undergoing T-CIN display massive microRNA (miR) expression changes that are not related to local losses or gains. This widespread miR deregulation encompasses a miR-200-dependent epithelial-to-mesenchymal transition (EMT) that confers to immortalized pre-tumoral cells phenotypic traits of metastatic potential. Remarkably, a miR signature of these cells, comprising a downregulation of miRs with conserved expression in kidney, was retrieved in poorly differentiated aggressive renal cell carcinomas. Our results reveal an unanticipated connection between telomere crisis and the activation of the EMT program that occurs at pre-invasive stages of epithelial cancers, through mechanisms that involve miR deregulation. Thus, this study provides a new rational into how telomere instability contributes to the acquisition of the malignant phenotype.

Multiple-checkpoint inhibition of thymic stromal lymphopoietin-induced TH2 response by TH17-related cytokines.

BACKGROUND: The interplay between allergy and autoimmunity has been a matter of long debate. Epidemiologic studies point to a decreased frequency of allergy in patients with autoimmune diseases. However, recent studies suggest that IL-17 and related cytokines, which play a central role in autoimmunity, might also promote allergy. OBJECTIVE: To address this controversy, we systematically studied the interactions between T(H)17-related cytokines and the thymic stromal lymphopoietin (TSLP)-mediated proallergic pathway. METHODS: We used human primary dendritic cells (DCs), T cells, and skin explants. A novel geometric representation and multivariate ANOVA were used to analyze the T(H) cytokine profile. RESULTS: We show that IL-17A specifically inhibits TSLP production but increases proinflammatory IL-8 production in human skin explants exposed to TNF-α and IL-4. This inhibitory activity was confirmed in cultured skin explants of atopic dermatitis lesions. At the T-cell level, T(H)17-polarizing cytokines (IL-1ß, IL-6, TGF-ß, and IL-23) inhibited T(H)2 differentiation induced by TSLP-activated DCs. This led to a global dominance of a T(H)17-polarizing environment over TSLP-activated DCs, as revealed by clustering and computational analysis. CONCLUSIONS: Our data indicate that T(H)17-related cytokines are negative regulators of the TSLP immune pathway. This might explain the decreased frequency of allergy in patients with autoimmunity and suggests new means of manipulating proallergic responses.

Patient Derived Xenografts (PDX) Models as an Avatar to Assess Personalized Therapy Options in Uveal Melanoma: A Feasibility Study.

Uveal melanoma is the most common primary intraocular malignancy in adults. Up to 50% of UM patients develop metastatic disease, usually in the liver. When metastatic, the prognosis is poor, and few treatment options exist. Here, we investigated the feasibility of establishing patient-derived xenografts (PDXs) from a patient's tumor in order to screen for therapies that the patient could benefit from. Samples obtained from 29 primary tumors and liver metastases of uveal melanoma were grafted into SCID mice. PDX models were successfully established for 35% of primary patient tumors and 67% of liver metastases. The tumor take rate was proportional to the risk of metastases. PDXs showed the same morphology, the same GNAQ/11, BAP1, and SF3B1 mutations, and the same chromosome 3 and 8q status as the corresponding patient samples. Six PDX models were challenged with two compounds for 4 weeks. We show that, for 31% of patients with high or intermediate risk of metastasis, the timing to obtain efficacy results on PDX models derived from their primary tumors was compatible with the selection of the therapy to treat the patient after relapse. PDXs could thus be a valid tool ("avatar") to select the best personalized therapy for one third of patients that are most at risk of relapse.

Extrachromosomal amplification mechanisms in a glioma with amplified sequences from multiple chromosome loci.

Accumulation of extrachromosomal DNA molecules (double minute) is often responsible for gene amplification in cancers, but the mechanisms leading to their formation are still largely unknown. By using quantitative PCR, chromosome walking, in situ hybridization on metaphase chromosomes and whole genome analysis, we studied a glioma containing four extrachromosomally amplified loci (7p11, 1q32.1, 5p15 and 9p2). Complex extrachromosomal DNA molecules were formed by the fusion of several syntenic or non-syntenic DNA fragments from 7p11, 5p15 to 9p2. Fragments ranged from a few base pairs to megabase pairs. Scars of the amplification process remained at the original locus in the form of deletions or chromosome rearrangements. Chromosome fragmentation, due to replication stress, could explain this complex situation. In contrast, at 1q32.1, the initial extrachromosomal DNA molecule resulted from the circularization of a single fragment associated with an intrachromosomal deletion including, but larger than, the amplified sequence. The nature of the sequences involved in these rearrangements suggests that a V(D)J-like illegitimate recombination contributes to its formation.

SMETHILLIUM: spatial normalization method for Illumina infinium HumanMethylation BeadChip.

SUMMARY: DNA methylation is a major epigenetic modification in human cells. Illumina HumanMethylation27 BeadChip makes it possible to quantify the methylation state of 27 578 loci spanning 14 495 genes. We developed a non-parametric normalization method to correct the spatial background noise in order to improve the signal-to-noise ratio. The prediction performance of the proposed method was assessed on three fully methylated samples and three fully unmethylated DNA samples. We demonstrate that the spatial normalization outperforms BeadStudio to predict the methylation state of a given locus. AVAILABILITY AND IMPLEMENTATION: A R script and the data are available at the following address: http://bioinfo.curie.fr/projects/smethillium.

Influence of FLG loss-of-function mutations in host-microbe interactions during atopic skin inflammation.

BACKGROUND: Loss-of-function mutations in the filaggrin (FLG) gene directly alter skin barrier function and critically influence atopic inflammation. While skin barrier dysfunction, Th2-associated inflammation and bacterial dysbiosis are well-known characteristics of atopic dermatitis (AD), the mechanisms interconnecting genotype, transcriptome and microbiome remain largely elusive. OBJECTIVE: In-depth analysis of FLG genotype-associated skin gene expression alterations and host-microbe interactions in AD. METHODS: Multi-omics characterization of a cohort of AD patients carrying heterozygous loss-of-function mutations in the FLG gene (ADMut) (n = 15), along with matched wild-type (ADWt) patients and healthy controls. Detailed clinical characterization, microarray gene expression and 16 S rRNA-based microbial marker gene data were generated and analyzed. RESULTS: In the context of filaggrin dysfunction, the transcriptome was characterized by dysregulation of barrier function and water homeostasis, while the lesional skin of ADWt demonstrated the specific upregulation of pro-inflammatory cytokines and T-cell proliferation. S. aureus dominated the microbiome in both patient groups, however, shifting microbial communities could be observed when comparing healthy with non-lesional ADWt or ADMut skin, offering the opportunity to identify microbe-associated transcriptomic signatures. Moreover, an AD core signature with 28 genes, including CCL13, CCL18, BTC, SCIN, RAB31 and PCLO was identified. CONCLUSIONS: Our integrative approach provides molecular insights for the concept that FLG loss-of-function mutations are a genetic shortcut to atopic inflammation and unravels the complex interplay between genotype, transcriptome and microbiome in the human holobiont.

Multiparametric analysis of cytokine-driven human Th17 differentiation reveals a differential regulation of IL-17 and IL-22 production.

T helper 17 (Th17) cells produce IL-17 but can also make tumor necrosis factor, interleukin (IL)-6, IL-10, IL-21, and IL-22. These cytokines collectively contribute to the functional outcome of the Th response. IL-22 plays a critical role in some Th17-associated diseases, such as psoriasis, but its relationship to IL-17 remains controversial. Here, we used a systematic multiparametric analysis of Th-17-associated cytokines, which revealed the unexpected finding that the regulation pattern of IL-22 was most closely related to interferon-gamma, the prototypical Th1 cytokine, and not to IL-17. To explain this observation, we systematically tested the role of Th1- and Th17-inducing cytokines. We could show that IL-12 and IL-23 induced high levels of IL-22 but no IL-17. Conversely, transforming growth factor-beta inhibited IL-22 production but promoted IL-17. Thus, IL-17 and IL-22 are differentially regulated during cytokine-induced Th cell differentiation. This has important implications for the understanding and pharmacologic manipulation of Th17-associated pathologies.

Frequent genomic structural alterations at HPV insertion sites in cervical carcinoma.

To investigate whether integration of HPV DNA in cervical carcinoma is responsible for structural alterations of the host genome at the insertion site, a series of 34 primary cervical carcinomas and eight cervical cancer-derived cell lines were analysed. DNA copy number profiles were assessed using the Affymetrix GeneChip Human Mapping 250K Sty array. HPV 16, 18 or 45 integration sites were determined using the DIPS-PCR technique. The genome status at integration sites was classified as follows: no change, amplification, transition normal/gain, normal/loss or gain/LOH. A single HPV integration site was found in 34 cases; two sites were found in seven cases; and three sites in one case (51 sites). Comparison between integration sites and DNA copy number profiles showed that the genome status was altered at 17/51 (33%) integration sites, corresponding to 16/42 cases (38%). Alterations detected were amplification in nine cases, transition normal/loss in four cases, normal/gain in three cases, and gain/LOH in one case. A highly significant association was found between genomic rearrangement and integration of HPV DNA (p < 10(-10)). Activation of the replication origin located in viral integrated sequences in a cell line derived from one of the primary cervical carcinomas induced an increase of the amplification level of both viral and cellular DNA sequences flanking the integration locus. This mechanism may be implicated in the triggering of genome amplification at the HPV integration site in cervical carcinoma. Structural alterations of the host genome are frequently observed at the integration site of HPV DNA in cervical cancer and may act in oncogenesis.

Geniac: Automatic Configuration GENerator and Installer for nextflow pipelines.

With the advent of high-throughput biotechnological platforms and their ever-growing capacity, life science has turned into a digitized, computational and data-intensive discipline. As a consequence, standard analysis with a bioinformatics pipeline in the context of routine production has become a challenge such that the data can be processed in real-time and delivered to the end-users as fast as possible. The usage of workflow management systems along with packaging systems and containerization technologies offer an opportunity to tackle this challenge. While very powerful, they can be used and combined in many multiple ways which may differ from one developer to another. Therefore, promoting the homogeneity of the workflow implementation requires guidelines and protocols which detail how the source code of the bioinformatics pipeline should be written and organized to ensure its usability, maintainability, interoperability, sustainability, portability, reproducibility, scalability and efficiency. Capitalizing on Nextflow, Conda, Docker, Singularity and the nf-core initiative, we propose a set of best practices along the development life cycle of the bioinformatics pipeline and deployment for production operations which target different expert communities including i) the bioinformaticians and statisticians ii) the software engineers and iii) the data managers and core facility engineers. We implemented Geniac (Automatic Configuration GENerator and Installer for nextflow pipelines) which consists of a toolbox with three components: i) a technical documentation available at https://geniac.readthedocs.io to detail coding guidelines for the bioinformatics pipeline with Nextflow, ii) a command line interface with a linter to check that the code respects the guidelines, and iii) an add-on to generate configuration files, build the containers and deploy the pipeline. The Geniac toolbox aims at the harmonization of development practices across developers and automation of the generation of configuration files and containers by parsing the source code of the Nextflow pipeline.

biogitflow: development workflow protocols for bioinformatics pipelines with git and GitLab.

The use of a bioinformatics pipeline as a tool to support diagnostic and theranostic decisions in the healthcare process requires the definition of detailed development workflow guidelines. Therefore, we implemented protocols that describe step-by-step all the command lines and actions that the developers have to follow. Our protocols capitalized on two powerful and widely used tools: git and GitLab. They address two use cases: a nominal mode to develop a new feature in the bioinformatics pipeline and a hotfix mode to correct a bug that occurred in the production environment. The protocols are available as a comprehensive documentation at https://biogitflow.readthedocs.io and the main concepts, steps and principles are presented in this report.