Pascal Belleau

Reproducibility, bioinformatic analysis and power of the SAGE method to evaluate changes in transcriptome

The serial analysis of gene expression (SAGE) method is used to study global gene expression in cells or tissues in various experimental conditions. However, its reproducibility has not yet been definitively assessed. In this study, we have evaluated the reproducibility of the SAGE method and identified the factors that affect it. The determination coefficient (R2) for the reproducibility of SAGE is 0.96. However, there are some factors that can affect the reproducibility of SAGE, such as the replication of concatemers and ditags, the number of sequenced tags and double PCR amplification of ditags. Thus, corrections for these factors must be made to ensure the reproducibility and accuracy of SAGE results. A bioinformatic analysis of SAGE data is also presented in order to eliminate these artifacts. Finally, the current study shows that increasing the number of sequenced tags improves the power of the method to detect transcripts and their regulation by experimental conditions.

Analysis of microsatellite markers and single nucleotide polymorphisms in candidate genes for susceptibility to bipolar affective disorder in the chromosome 12Q24.31 region.

Previous results from our genetic analyses using pedigrees from a French Canadian population suggested that the interval delimited by markers D12S86 and D12S378 on chromosome 12 was the most probable genomic region to contain a susceptibility gene for affective disorders. Here we present a more detailed genetic analysis of a 7.7 Mb genomic region located on 12q24.31. This region was saturated with 20 microsatellite markers to refine the candidate region and linkage analysis performed in 41 families from the Saguenay‐Lac‐St‐Jean (SLSJ) region of Quebec. The results of two point parametric analysis using MFLINK supported the presence of a susceptibility locus on chromosome 12q24.31. Association studies with microsatellite markers using a case/control sample from the same population (n = 401) and analyzed with CLUMP revealed significant allelic associations between the bipolar phenotype and markers NBG6 (P = 0.008) and NBG12 (P < 10−3). According to these results, we investigated candidate genes in the NBG12 area. We analyzed 32 genes for the presence of polymorphisms in coding sequences and intron/exon junctions and genotyped 22 non‐synonymous SNPs in the SLSJ case/control sample. Two uncommon polymorphisms (minor allele frequency ≤ 0.03) found in KIAA1595 and FLJ22471 genes, gave P‐values below 0.05 with the T1 statistic. Moreover, using haplotype analysis, a nearly significant haplotypic association was observed at the HM74 gene. These results do not give strong support for a role in the susceptibility to bipolar disorder of any of these genes analyzed. However, the significance of rare polymorphisms should be explored by further analyses.

PITX2 Is Involved in Stress Response in Cultured Human Trabecular Meshwork Cells through Regulation of SLC13A3

PURPOSE Mutations of the PITX2 gene cause Axenfeld-Rieger syndrome (ARS) and glaucoma. In this study, the authors investigated genes directly regulated by the PITX2 transcription factor to gain insight into the mechanisms underlying these disorders. METHODS RNA from nonpigmented ciliary epithelium cells transfected with hormone-inducible PITX2 and activated by mifepristone was subjected to microarray analyses. Data were analyzed using dCHIP algorithms to detect significant differences in expression. Genes with significantly altered expression in multiple microarray experiments in the presence of activated PITX2 were subjected to in silico and biochemical analyses to validate them as direct regulatory targets. One target gene was further characterized by studying the effect of its knockdown in a cell model of oxidative stress, and its expression in zebrafish embryos was analyzed by in situ hybridization. RESULTS Solute carrier family 13 sodium-dependent dicarboxylate transporter member 3 (SLC13A3) was identified as 1 of 47 potential PITX2 target genes in ocular cells. PITX2 directly regulates SLC13A3 expression, as demonstrated by luciferase reporter and chromatin immunoprecipitation assays. Reduction of PITX2 or SLC13A3 levels by small interfering RNA (siRNA)-mediated knockdown augmented the death of transformed human trabecular meshwork cells exposed to hydrogen peroxide. Zebrafish slc13a3 is expressed in anterior ocular regions in a pattern similar to that of pitx2. CONCLUSIONS The results indicate that SLC13A3 is a direct downstream target of PITX2 transcriptional regulation and that levels of PITX2 and SLC13A3 modulate responses to oxidative stress in ocular cells.

Genetic variants and haplotype analyses of the ZBRK1/ZNF350 gene in high-risk non BRCA1/2 French Canadian breast and ovarian cancer families.

Our current understanding of breast cancer susceptibility involves mutations in the 2 major genes BRCA1 and BRCA2, found in about 25% of high‐risk families, as well as few other low penetrance genes such as ATM and CHEK2. Approximately two‐thirds of the multiple cases families remain to be explained by mutations in still unknown genes. In a candidate gene approach to identify new genes potentially involved in breast cancer susceptibility, we analyzed genomic variants in the ZBRK1 gene, a co‐repressor implicated in BRCA1‐mediated repression of GADD45. Direct sequencing of ZBRK1 entire coding region in affected breast cancer individuals from 97 high‐risk French Canadian breast/ovarian cancer families and 94 healthy controls led to the identification of 18 genomic variants. Haplotype analyses, using PHASE, COCAPHASE and HaploStats programs, put in evidence 3 specific haplotypes which could potentially modulate breast cancer risk, and among which 2 that are associated with a potential protective effect (p = 0.01135 and p = 0.00268), while another haplotype is over‐represented in the case group (p = 0.00143). Further analyses of these haplotypes indicated that a strong component of the observed difference between both groups emerge from the first 5 variants (out of 12 used for haplotype determination). The present study also permitted to determine a set of tagging SNPs that could be useful for subsequent analyses in large scale association studies. Additional studies in large cohorts and other populations will however be needed to further evaluate if common and/or rare ZBRK1 sequence variants and haplotypes could be associated with a modest/intermediate breast cancer risk.

Polymorphic variations in the FANCA gene in high‐risk non‐BRCA1/2 breast cancer individuals from the French Canadian population

The majority of genes associated with breast cancer susceptibility, including BRCA1 and BRCA2 genes, are involved in DNA repair mechanisms. Moreover, among the genes recently associated with an increased susceptibility to breast cancer, four are Fanconi Anemia (FA) genes: FANCD1/BRCA2, FANCJ/BACH1/BRIP1, FANCN/PALB2 and FANCO/RAD51C. FANCA is implicated in DNA repair and has been shown to interact directly with BRCA1. It has been proposed that the formation of FANCA/G (dependent upon the phosphorylation of FANCA) and FANCB/L sub‐complexes altogether with FANCM, represent the initial step for DNA repair activation and subsequent formation of other sub‐complexes leading to ubiquitination of FANCD2 and FANCI. As only approximately 25% of inherited breast cancers are attributable to BRCA1/2 mutations, FANCA therefore becomes an attractive candidate for breast cancer susceptibility. We thus analyzed FANCA gene in 97 high‐risk French Canadian non‐BRCA1/2 breast cancer individuals by direct sequencing as well as in 95 healthy control individuals from the same population. Among a total of 85 sequence variants found in either or both series, 28 are coding variants and 19 of them are missense variations leading to amino acid change. Three of the amino acid changes, namely Thr561Met, Cys625Ser and particularly Ser1088Phe, which has been previously reported to be associated with FA, are predicted to be damaging by the SIFT and PolyPhen softwares. cDNA amplification revealed significant expression of 4 alternative splicing events (insertion of an intronic portion of intron 10, and the skipping of exons 11, 30 and 31). In silico analyzes of relevant genomic variants have been performed in order to identify potential variations involved in the expression of these spliced transcripts. Sequence variants in FANCA could therefore be potential spoilers of the Fanconi‐BRCA pathway and as a result, they could in turn have an impact in non‐BRCA1/2 breast cancer families.

Organoid profiling identifies common responders to chemotherapy in pancreatic cancer

Pancreatic cancer is the most lethal common solid malignancy. Systemic therapies are often ineffective, and predictive biomarkers to guide treatment are urgently needed. We generated a pancreatic cancer patient-derived organoid (PDO) library that recapitulates the mutational spectrum and transcriptional subtypes of primary pancreatic cancer. New driver oncogenes were nominated and transcriptomic analyses revealed unique clusters. PDOs exhibited heterogeneous responses to standard-of-care chemotherapeutics and investigational agents. In a case study manner, we found that PDO therapeutic profiles paralleled patient outcomes and that PDOs enabled longitudinal assessment of chemosensitivity and evaluation of synchronous metastases. We derived organoid-based gene expression signatures of chemosensitivity that predicted improved responses for many patients to chemotherapy in both the adjuvant and advanced disease settings. Finally, we nominated alternative treatment strategies for chemorefractory PDOs using targeted agent therapeutic profiling. We propose that combined molecular and therapeutic profiling of PDOs may predict clinical response and enable prospective therapeutic selection.Significance: New approaches to prioritize treatment strategies are urgently needed to improve survival and quality of life for patients with pancreatic cancer. Combined genomic, transcriptomic, and therapeutic profiling of PDOs can identify molecular and functional subtypes of pancreatic cancer, predict therapeutic responses, and facilitate precision medicine for patients with pancreatic cancer. Cancer Discov; 8(9); 1112-29. ©2018 AACR.See related commentary by Collisson, p. 1062This article is highlighted in the In This Issue feature, p. 1047.

Germline copy number variation in the YTHDC2 gene: does it have a role in finding a novel potential molecular target involved in pancreatic adenocarcinoma susceptibility?

Objective: The vast majority of pancreatic cancers occurs sporadically. The discovery of frequent variations in germline gene copy number can significantly influence the expression levels of genes that predispose to pancreatic adenocarcinoma. We prospectively investigated whether patients with sporadic pancreatic adenocarcinoma share specific gene copy number variations (CNVs) in their germline DNA. Patients and methods: DNA samples were analyzed from peripheral leukocytes from 72 patients with a diagnosis of sporadic pancreatic adenocarcinoma and from 60 controls using Affymetrix 500K array set. Multiplex ligation-dependent probe amplification (MLPA) assay was performed using a set of self-designed MLPA probes specific for seven target sequences. Results: We identified a CNV-containing DNA region associated with pancreatic cancer risk. This region shows a deletion of 1 allele in 36 of the 72 analyzed patients but in none of the controls. This region is of particular interest since it contains the YTHDC2 gene encoding for a putative DNA/RNA helicase, such protein being frequently involved in cancer susceptibility. Interestingly, 82.6% of Sicilian patients showed germline loss of one allele. Conclusions: Our results suggest that the YTHDC2 gene could be a potential candidate for pancreatic cancer susceptibility and a useful marker for early detection as well as for the development of possible new therapeutic strategies.

Analysis of germline gene copy number variants of patients with sporadic pancreatic adenocarcinoma reveals specific variations.

Objectives: The rapid fatality of pancreatic cancer is, in large part, the result of diagnosis at an advanced stage in the majority of patients. Identification of individuals at risk of developing pancreatic adenocarcinoma would be useful to improve the prognosis of this disease. There is presently no biological or genetic indicator allowing the detection of patients at risk. Our main goal was to identify copy number variants (CNVs) common to all patients with sporadic pancreatic cancer. Methods: We analyzed gene CNVs in leukocyte DNA from 31 patients with sporadic pancreatic adenocarcinoma and from 93 matched controls. Genotyping was performed with the use of the GeneChip Human Mapping 500K Array Set (Affymetrix). Results: We identified 431 single nucleotide polymorphism (SNP) probes with abnormal hybridization signal present in the DNA of all 31 patients. Of these SNP probes, 284 corresponded to 3 or more copies and 147 corresponded to 1 or 0 copies. Several cancer-associated genes were amplified in all patients. Conversely, several genes supposed to oppose cancer development were present as single copy. Conclusions: These data suggest that a set of 431 CNVs could be associated with the disease. This set could be useful for early diagnosis.

Fine temporal analysis of DHT transcriptional modulation of the ATM/Gadd45g signaling pathways in the mouse uterus

In rodents, the uterus of a mature female undergoes changes during the uterine cycle, under the control of steroid hormones. 5α‐Dihydrotestosterone (DHT) is recognized to play an important role in the regulation of androgen action in normal endometrium. Using microarray technology, a screening analysis of genes responding to DHT in the uterus of ovariectomized mice, has allowed us to highlight multiple genes of the ATM/Gadd45g pathway that are modulated following exposure to DHT. Two phases of regulation were identified. In the early phase, the expression of genes involved in the G2/M arrest is rapidly increased, followed by the repression of genes of the G1/S checkpoint, and by the induction of transcriptional regulators. Later, i.e. from 12 to 24 hr, genes involved in G2/M transition, cytoarchitectural and lipid‐related genes are stimulated by DHT while immunity‐related genes appear to be differentially regulated by the hormone. These results show that a physiological dose of DHT induces the transcription of genes promoting the cell cycle progression in mice. Profile determination of temporal uterine gene expression at the transcriptional level enables us to suggest that the DHT modulation of genes involved in ATM/Gadd45g signaling in an ATM‐ or p53‐independent manner, could play an important role in the cyclical changes of uterine cells in the mouse uterus. Mol. Reprod. Dev. 76: 278–288, 2009.

metagene Profiles Analyses Reveal Regulatory Element's Factor-Specific Recruitment Patterns.

ChIP-Sequencing (ChIP-Seq) provides a vast amount of information regarding the localization of proteins across the genome. The aggregation of ChIP-Seq enrichment signal in a metagene plot is an approach commonly used to summarize data complexity and to obtain a high level visual representation of the general occupancy pattern of a protein. Here we present the R package metagene, the graphical interface Imetagene and the companion package similaRpeak. Together, they provide a framework to integrate, summarize and compare the ChIP-Seq enrichment signal from complex experimental designs. Those packages identify and quantify similarities or dissimilarities in patterns between large numbers of ChIP-Seq profiles. We used metagene to investigate the differential occupancy of regulatory factors at noncoding regulatory regions (promoters and enhancers) in relation to transcriptional activity in GM12878 B-lymphocytes. The relationships between occupancy patterns and transcriptional activity suggest two different mechanisms of action for transcriptional control: i) a “gradient effect” where the regulatory factor occupancy levels follow transcription and ii) a “threshold effect” where the regulatory factor occupancy levels max out prior to reaching maximal transcription. metagene, Imetagene and similaRpeak are implemented in R under the Artistic license 2.0 and are available on Bioconductor.