Preclinical workup using long-read amplicon sequencing provides families with de novo pathogenic variants access to universal preimplantation genetic testing.

STUDY QUESTION: Can long-read amplicon sequencing be beneficial for preclinical preimplantation genetic testing (PGT) workup in couples with a de novo pathogenic variant in one of the prospective parents? SUMMARY ANSWER: Long-read amplicon sequencing represents a simple, rapid and cost-effective preclinical PGT workup strategy that provides couples with de novo pathogenic variants access to universal genome-wide haplotyping-based PGT programs. WHAT IS KNOWN ALREADY: Universal PGT combines genome-wide haplotyping and copy number profiling to select embryos devoid of both familial pathogenic variants and aneuploidies. However, it cannot be directly applied in couples with a de novo pathogenic variant in one of the partners due to the absence of affected family members required for phasing the disease-associated haplotype. STUDY DESIGN, SIZE, DURATION: This is a prospective study, which includes 32 families that were enrolled in the universal PGT program at the University Hospital of Leuven between 2018 and 2022. We implemented long-read amplicon sequencing during the preclinical PGT workup to deduce the parental origin of the disease-associated allele in the affected partner, which can then be traced in embryos during clinical universal PGT cycles. PARTICIPANTS/MATERIALS, SETTING, METHODS: To identify the parental origin of the disease-associated allele, genomic DNA from the carrier of the de novo pathogenic variant and his/her parent(s) was used for preclinical PGT workup. Primers flanking the de novo variant upstream and downstream were designed for each family. Following long-range PCR, amplicons that ranged 5-10 kb in size, were sequenced using Pacific Bioscience and/or Oxford Nanopore platforms. Next, targeted variant calling and haplotyping were performed to identify parental informative single-nucleotide variants (iSNVs) linked to the de novo mutation. Following the preclinical PGT workup, universal PGT via genome-wide haplotyping was performed for couples who proceeded with clinical PGT cycle. In parallel, 13 trophectoderm (TE) biopsies from three families that were analyzed by universal PGT, were also used for long-read amplicon sequencing to explore this approach for embryo direct mutation detection coupled with targeted long-read haplotyping. MAIN RESULTS AND THE ROLE OF CHANCE: The parental origin of the mutant allele was identified in 24/32 affected individuals during the preclinical PGT workup stage, resulting in a 75% success rate. On average, 5.95 iSNVs (SD = 4.5) were detected per locus of interest, and the average distance of closest iSNV to the de novo variant was ∼1750 bp. In 75% of those cases (18/24), the de novo mutation occurred on the paternal allele. In the remaining eight families, the risk haplotype could not be established due to the absence of iSNVs linked to the mutation or inability to successfully target the region of interest. During the time of the study, 12/24 successfully analyzed couples entered the universal PGT program, and three disease-free children have been born. In parallel to universal PGT analysis, long-read amplicon sequencing of 13 TE biopsies was also performed, confirming the segregation of parental alleles in the embryo and the results of the universal PGT. LIMITATIONS, REASONS FOR CAUTION: The main limitation of this approach is that it remains targeted with the need to design locus-specific primers. Because of the restricted size of target amplicons, the region of interest may also remain non-informative in the absence of iSNVs. WIDER IMPLICATIONS OF THE FINDINGS: Targeted haplotyping via long-read amplicon sequencing, particularly using Oxford Nanopore Technologies, provides a valuable alternative for couples with de novo pathogenic variants that allows access to universal PGT. Moreover, the same approach can be used for direct mutation analysis in embryos, as a second line confirmation of the preclinical PGT result or as a potential alternative PGT procedure in couples, where additional family members are not available. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by KU Leuven funding (no. C1/018 to J.R.V.) and Fonds Wetenschappelijk Onderzoek (1241121N to O.T.). J.R.V. is co-inventor of a patent ZL910050-PCT/EP2011/060211-WO/2011/157846 'Methods for haplotyping single-cells' and ZL913096-PCT/EP2014/068315-WO/2015/028576 'Haplotyping and copy number typing using polymorphic variant allelic frequencies' licensed to Agilent Technologies. All other authors have no conflict of interest to declare. TRIAL REGISTRATION NUMBER: N/A.

Myotonic dystrophy type 1 embryonic stem cells show decreased myogenic potential, increased CpG methylation at the DMPK locus and RNA mis-splicing.

Skeletal muscle tissue is severely affected in myotonic dystrophy type 1 (DM1) patients, characterised by muscle weakness, myotonia and muscle immaturity in the most severe congenital form of the disease. Previously, it was not known at what stage during myogenesis the DM1 phenotype appears. In this study we differentiated healthy and DM1 human embryonic stem cells to myoblasts and myotubes and compared their differentiation potential using a comprehensive multi-omics approach. We found myogenesis in DM1 cells to be abnormal with altered myotube generation compared to healthy cells. We did not find differentially expressed genes between DM1 and non-DM1 cell lines within the same developmental stage. However, during differentiation we observed an aberrant inflammatory response and increased CpG methylation upstream of the CTG repeat at the myoblast level and RNA mis-splicing at the myotube stage. We show that early myogenesis modelled in hESC reiterates the early developmental manifestation of DM1.

MSH2 knock-down shows CTG repeat stability and concomitant upstream demethylation at the DMPK locus in myotonic dystrophy type 1 human embryonic stem cells.

Myotonic dystrophy type 1 (DM1) is caused by expansion of a CTG repeat in the DMPK gene, where expansion size and somatic mosaicism correlates with disease severity and age of onset. While it is known that the mismatch repair protein MSH2 contributes to the unstable nature of the repeat, its role on other disease-related features, such as CpG methylation upstream of the repeat, is unknown. In this study, we investigated the effect of an MSH2 knock-down (MSH2KD) on both CTG repeat dynamics and CpG methylation pattern in human embryonic stem cells (hESC) carrying the DM1 mutation. Repeat size in MSH2 wild-type (MSH2WT) and MSH2KD DM1 hESC was determined by PacBio sequencing and CpG methylation by bisulfite massive parallel sequencing. We found stabilization of the CTG repeat concurrent with a gradual loss of methylation upstream of the repeat in MSH2KD cells, while the repeat continued to expand and upstream methylation remained unchanged in MSH2WT control lines. Repeat instability was re-established and biased towards expansions upon MSH2 transgenic re-expression in MSH2KD lines while upstream methylation was not consistently re-established. We hypothesize that the hypermethylation at the mutant DM1 locus is promoted by the MMR machinery and sustained by a constant DNA repair response, establishing a potential mechanistic link between CTG repeat instability and upstream CpG methylation. Our work represents a first step towards understanding how epigenetic alterations and repair pathways connect and contribute to the DM1 pathology.

Multi-centre evaluation of a comprehensive preimplantation genetic test through haplotyping-by-sequencing.

STUDY QUESTION: Can reduced representation genome sequencing offer an alternative to single nucleotide polymorphism (SNP) arrays as a generic and genome-wide approach for comprehensive preimplantation genetic testing for monogenic disorders (PGT-M), aneuploidy (PGT-A) and structural rearrangements (PGT-SR) in human embryo biopsy samples? SUMMARY ANSWER: Reduced representation genome sequencing, with OnePGT, offers a generic, next-generation sequencing-based approach for automated haplotyping and copy-number assessment, both combined or independently, in human single blastomere and trophectoderm samples. WHAT IS KNOWN ALREADY: Genome-wide haplotyping strategies, such as karyomapping and haplarithmisis, have paved the way for comprehensive PGT, i.e. leveraging PGT-M, PGT-A and PGT-SR in a single workflow. These methods are based upon SNP array technology. STUDY DESIGN, SIZE, DURATION: This multi-centre verification study evaluated the concordance of PGT results for a total of 225 embryos, including 189 originally tested for a monogenic disorder and 36 tested for a translocation. Concordance for whole chromosome aneuploidies was also evaluated where whole genome copy-number reference data were available. Data analysts were kept blind to the results from the reference PGT method. PARTICIPANTS/MATERIALS, SETTING, METHODS: Leftover blastomere/trophectoderm whole genome amplified (WGA) material was used, or secondary trophectoderm biopsies were WGA. A reduced representation library from WGA DNA together with bulk DNA from phasing references was processed across two study sites with the Agilent OnePGT solution. Libraries were sequenced on an Illumina NextSeq500 system, and data were analysed with Agilent Alissa OnePGT software. The embedded PGT-M pipeline utilises the principles of haplarithmisis to deduce haplotype inheritance whereas both the PGT-A and PGT-SR pipelines are based upon read-count analysis in order to evaluate embryonic ploidy. Concordance analysis was performed for both analysis strategies against the reference PGT method. MAIN RESULTS AND THE ROLE OF CHANCE: PGT-M analysis was performed on 189 samples. For nine samples, the data quality was too poor to analyse further, and for 20 samples, no result could be obtained mainly due to biological limitations of the haplotyping approach, such as co-localisation of meiotic crossover events and nullisomy for the chromosome of interest. For the remaining 160 samples, 100% concordance was obtained between OnePGT and the reference PGT-M method. Equally for PGT-SR, 100% concordance for all 36 embryos tested was demonstrated. Moreover, with embryos originally analysed for PGT-M or PGT-SR for which genome-wide copy-number reference data were available, 100% concordance was shown for whole chromosome copy-number calls (PGT-A). LIMITATIONS, REASONS FOR CAUTION: Inherent to haplotyping methodologies, processing of additional family members is still required. Biological limitations caused inconclusive results in 10% of cases. WIDER IMPLICATIONS OF THE FINDINGS: Employment of OnePGT for PGT-M, PGT-SR, PGT-A or combined as comprehensive PGT offers a scalable platform, which is inherently generic and thereby, eliminates the need for family-specific design and optimisation. It can be considered as both an improvement and complement to the current methodologies for PGT. STUDY FUNDING/COMPETING INTEREST(S): Agilent Technologies, the KU Leuven (C1/018 to J.R.V. and T.V.) and the Horizon 2020 WIDENLIFE (692065 to J.R.V. and T.V). H.M. is supported by the Research Foundation Flanders (FWO, 11A7119N). M.Z.E, J.R.V. and T.V. are co-inventors on patent applications: ZL910050-PCT/EP2011/060211- WO/2011/157846 'Methods for haplotyping single cells' and ZL913096-PCT/EP2014/068315 'Haplotyping and copy-number typing using polymorphic variant allelic frequencies'. T.V. and J.R.V. are co-inventors on patent application: ZL912076-PCT/EP2013/070858 'High-throughput genotyping by sequencing'. Haplarithmisis ('Haplotyping and copy-number typing using polymorphic variant allelic frequencies') has been licensed to Agilent Technologies. The following patents are pending for OnePGT: US2016275239, AU2014345516, CA2928013, CN105874081, EP3066213 and WO2015067796. OnePGT is a registered trademark. D.L., J.T. and R.L.R. report personal fees during the conduct of the study and outside the submitted work from Agilent Technologies. S.H. and K.O.F. report personal fees and other during the conduct of the study and outside the submitted work from Agilent Technologies. J.A. reports personal fees and other during the conduct of the study from Agilent Technologies and personal fees from Agilent Technologies and UZ Leuven outside the submitted work. B.D. reports grants from IWT/VLAIO, personal fees during the conduct of the study from Agilent Technologies and personal fees and other outside the submitted work from Agilent Technologies. In addition, B.D. has a patent 20160275239 - Genetic Analysis Method pending. The remaining authors have no conflicts of interest.

Prognostic relevance of molecular subtypes and master regulators in pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic cancer is poorly characterized at genetic and non-genetic levels. The current study evaluates in a large cohort of patients the prognostic relevance of molecular subtypes and key transcription factors in pancreatic ductal adenocarcinoma (PDAC). METHODS: We performed gene expression analysis of whole-tumor tissue obtained from 118 surgically resected PDAC and 13 histologically normal pancreatic tissue samples. Cox regression models were used to study the effect on survival of molecular subtypes and 16 clinicopathological prognostic factors. In order to better understand the biology of PDAC we used iRegulon to identify transcription factors (TFs) as master regulators of PDAC and its subtypes. RESULTS: We confirmed the PDAssign gene signature as classifier of PDAC in molecular subtypes with prognostic relevance. We found molecular subtypes, but not clinicopathological factors, as independent predictors of survival. Regulatory network analysis predicted that HNF1A/B are among thousand TFs the top enriched master regulators of the genes expressed in the normal pancreatic tissue compared to the PDAC regulatory network. On immunohistochemistry staining of PDAC samples, we observed low expression of HNF1B in well differentiated towards no expression in poorly differentiated PDAC samples. We predicted IRF/STAT, AP-1, and ETS-family members as key transcription factors in gene signatures downstream of mutated KRAS. CONCLUSIONS: PDAC can be classified in molecular subtypes that independently predict survival. HNF1A/B seem to be good candidates as master regulators of pancreatic differentiation, which at the protein level loses its expression in malignant ductal cells of the pancreas, suggesting its putative role as tumor suppressor in pancreatic cancer. TRIAL REGISTRATION: The study was registered at ClinicalTrials.gov under the number NCT01116791 (May 3, 2010).

Genetic predisposition for beta cell fragility underlies type 1 and type 2 diabetes.

Type 1 (T1D) and type 2 (T2D) diabetes share pathophysiological characteristics, yet mechanistic links have remained elusive. T1D results from autoimmune destruction of pancreatic beta cells, whereas beta cell failure in T2D is delayed and progressive. Here we find a new genetic component of diabetes susceptibility in T1D non-obese diabetic (NOD) mice, identifying immune-independent beta cell fragility. Genetic variation in Xrcc4 and Glis3 alters the response of NOD beta cells to unfolded protein stress, enhancing the apoptotic and senescent fates. The same transcriptional relationships were observed in human islets, demonstrating the role of beta cell fragility in genetic predisposition to diabetes.

Human pancreatic cancer contains a side population expressing cancer stem cell-associated and prognostic genes.

In many types of cancers, a side population (SP) has been identified based on high efflux capacity, thereby enriching for chemoresistant cells as well as for candidate cancer stem cells (CSC). Here, we explored whether human pancreatic ductal adenocarcinoma (PDAC) contains a SP, and whether its gene expression profile is associated with chemoresistance, CSC and prognosis. After dispersion into single cells and incubation with Hoechst dye, we analyzed human PDAC resections specimens using flow cytometry (FACS). We identified a SP and main population (MP) in all human PDAC resection specimens (n = 52) analyzed, but detected immune (CD45(+)) and endothelial (CD31(+)) cells in this fraction together with tumor cells. The SP and MP cells, or more purified fractions depleted from CD31(+)/CD45(+) cells (pSP and pMP), were sorted by FACS and subjected to whole-genome expression analysis. This revealed upregulation of genes associated with therapy resistance and of markers identified before in putative pancreatic CSC. pSP gene signatures of 32 or 10 up- or downregulated genes were developed and tested for discriminatory competence between pSP and pMP in different sets of PDAC samples. The prognostic value of the pSP genes was validated in a large independent series of PDAC patients (n = 78) using nCounter analysis of expression (in tumor versus surrounding pancreatic tissue) and Cox regression for disease-free and overall survival. Of these genes, expression levels of ABCB1 and CXCR4 were correlated with worse patient survival. Thus, our study for the first time demonstrates that human PDAC contains a SP. This tumor subpopulation may represent a valuable therapeutic target given its chemoresistance- and CSC-associated gene expression characteristics with potential prognostic value.

The placental growth factor as a target against hepatocellular carcinoma in a diethylnitrosamine-induced mouse model.

BACKGROUND & AIMS: The placental growth factor (PlGF) is a member of the vascular endothelial growth factor (VEGF) family known to stimulate endothelial cell growth, migration and survival, attract angiocompetent macrophages, and determine the metastatic niche. Unlike VEGF, genetic studies have shown that PlGF is specifically involved in pathologic angiogenesis, thus its inhibition would not affect healthy blood vessels, providing an attractive drug candidate with a good safety profile. METHODS: We assess whether inhibition of PlGF could be used as a potential therapy against hepatocellular carcinoma (HCC), by using PlGF knockout mice and monoclonal anti-PlGF antibodies in a mouse model for HCC. In addition, the effect of PlGF antibodies is compared to that of sorafenib, as well as the combination of both therapies. RESULTS: We have found that both in a transgenic knockout model and in a treatment model, targeting PlGF significantly decreases tumor burden. This was achieved not only by inhibiting neovascularisation, but also by decreasing hepatic macrophage recruitment and by normalising the remaining blood vessels, thereby decreasing hypoxia and reducing the prometastatic potential of HCC. CONCLUSIONS: Considering the favourable safety profile and its pleiotropic effect on vascularisation, metastasis and inflammation, PlGF inhibition could become a valuable therapeutic strategy against HCC.

Adaptation of Lactobacillus plantarum IMDO 130201, a wheat sourdough isolate, to growth in wheat sourdough simulation medium at different pH values through differential gene expression.

Sourdough is a very competitive and challenging environment for microorganisms. Usually, a stable microbiota composed of lactic acid bacteria (LAB) and yeasts dominates this ecosystem. Although sourdough is rich in carbohydrates, thus providing an ideal environment for microorganisms to grow, its low pH presents a particular challenge. The nature of the adaptation to this low pH was investigated for Lactobacillus plantarum IMDO 130201, an isolate from a laboratory wheat sourdough fermentation. Batch fermentations were carried out in wheat sourdough simulation medium, and total RNA was isolated from mid-exponential-growth-phase cultures, followed by differential gene expression analysis using a LAB functional gene microarray. At low pH values, an increased expression of genes involved in peptide and amino acid metabolism was found as well as that of genes involved in plantaricin production and lipoteichoic acid biosynthesis. The results highlight cellular mechanisms that allow L. plantarum to function at a low environmental pH.

Metatranscriptome analysis for insight into whole-ecosystem gene expression during spontaneous wheat and spelt sourdough fermentations.

Lactic acid bacteria (LAB) are of industrial importance in the production of fermented foods, including sourdough-derived products. Despite their limited metabolic capacity, LAB contribute considerably to important characteristics of fermented foods, such as extended shelf-life, microbial safety, improved texture, and enhanced organoleptic properties. Triggered by the considerable amount of LAB genomic information that became available during the last decade, transcriptome and, by extension, metatranscriptome studies have become one of the most appropriate research approaches to study whole-ecosystem gene expression in more detail. In this study, microarray analyses were performed using RNA sampled during four 10-day spontaneous sourdough fermentations carried out in the laboratory with an in-house-developed LAB functional gene microarray. For data analysis, a new algorithm was developed to calculate a net expression profile for each of the represented genes, allowing use of the microarray analysis beyond the species level. In addition, metabolite target analyses were performed on the sourdough samples to relate gene expression with metabolite production. The results revealed the activation of different key metabolic pathways, the ability to use carbohydrates other than glucose (e.g., starch and maltose), and the conversion of amino acids as a contribution to redox equilibrium and flavor compound generation in LAB during sourdough fermentation.

A seven-gene set associated with chronic hypoxia of prognostic importance in hepatocellular carcinoma.

PURPOSE: Hepatocellular carcinomas (HCC) have an unpredictable clinical course, and molecular classification could provide better insights into prognosis and patient-directed therapy. We hypothesized that in HCC, certain microenvironmental regions exist with a characteristic gene expression related to chronic hypoxia which would induce aggressive behavior. EXPERIMENTAL DESIGN: We determined the gene expression pattern for human HepG2 liver cells under chronic hypoxia by microarray analysis. Differentially expressed genes were selected and their clinical values were assessed. In our hypothesis-driven analysis, we included available independent microarray studies of patients with HCC in one single analysis. Three microarray studies encompassing 272 patients were used as training sets to determine a minimal prognostic gene set, and one recent study of 91 patients was used for validation. RESULTS: Using computational methods, we identified seven genes (out of 3,592 differentially expressed under chronic hypoxia) that showed correlation with poor prognostic indicators in all three training sets (65/139/73 patients) and this was validated in a fourth data set (91 patients). Retrospectively, the seven-gene set was associated with poor survival (hazard ratio, 1.39; P = 0.007) and early recurrence (hazard ratio, 2.92; P = 0.007) in 135 patients. Moreover, using a hypoxia score based on this seven-gene set, we found that patients with a score of >0.35 (n = 42) had a median survival of 307 days, whereas patients with a score of < or =0.35 (n = 93) had a median survival of 1,602 days (P = 0.005). CONCLUSIONS: We identified a unique, liver-specific, seven-gene signature associated with chronic hypoxia that correlates with poor prognosis in HCCs.

Community dynamics of bacteria in sourdough fermentations as revealed by their metatranscriptome.

The lactic acid bacterial (LAB) community dynamics of two wheat and two spelt sourdough fermentations that were daily back-slopped were monitored during a period of 10 days by hybridizing time-related RNA samples, representing the metatranscriptome, to an LAB functional gene microarray. To indicate the species present in each hybridized sample, annotation information for the 2,269 oligonucleotides on the microarray was used. The overall hybridization data revealed that after a transition phase of 5 days, in which atypical sourdough LAB species, including Enterococcus species, were found, a stabilized ecosystem was established with Lactobacillus plantarum and Lactobacillus fermentum as the dominating LAB species. Compared with the combined outcome of culture-dependent and culture-independent identification techniques, the microarray data revealed a functional role for Lactococcus lactis in the early stage ecosystem and the dominance of Pediococcus pentosaceus in most of the fermentations, besides L. plantarum and L. fermentum. Consequently, metatranscriptome hybridization data obtained using an LAB functional gene microarray was shown to be an interesting alternative to microbiological analysis of the community dynamics of complex food ecosystems.

An experimental loop design for the detection of constitutional chromosomal aberrations by array CGH.

BACKGROUND: Comparative genomic hybridization microarrays for the detection of constitutional chromosomal aberrations is the application of microarray technology coming fastest into routine clinical application. Through genotype-phenotype association, it is also an important technique towards the discovery of disease causing genes and genomewide functional annotation in human. When using a two-channel microarray of genomic DNA probes for array CGH, the basic setup consists in hybridizing a patient against a normal reference sample. Two major disadvantages of this setup are (1) the use of half of the resources to measure a (little informative) reference sample and (2) the possibility that deviating signals are caused by benign copy number variation in the "normal" reference instead of a patient aberration. Instead, we apply an experimental loop design that compares three patients in three hybridizations. RESULTS: We develop and compare two statistical methods (linear models of log ratios and mixed models of absolute measurements). In an analysis of 27 patients seen at our genetics center, we observed that the linear models of the log ratios are advantageous over the mixed models of the absolute intensities. CONCLUSION: The loop design and the performance of the statistical analysis contribute to the quick adoption of array CGH as a routine diagnostic tool. They lower the detection limit of mosaicisms and improve the assignment of copy number variation for genetic association studies.

Development and validation of a species-independent functional gene microarray that targets lactic acid bacteria.

During the last few years, genome-related information has become available for many microorganisms, including important food-related bacteria. Lactic acid bacteria (LAB) are important industrially in the production of fermented foods such as dairy products, sausages, sourdoughs, and vegetables. Despite their limited metabolic capacity, LAB contribute considerably to important characteristics of fermented foods, such as flavor and texture. In the present study, a species-independent functional gene microarray was developed that targets 406 genes that play key roles in the production of sugar catabolites, bacteriocins, exopolysaccharides, and aromas, in probiotic and biosafety characteristics, and in the stress response. Also, genes linked to negative traits, such as antibiotic resistance and virulence, are represented. As LAB ecosystems contain a variety of species, there was a more global focus on these specific functional properties. Thus, an algorithm was used to design gene-specific oligonucleotides that preferably hybridize with multiple LAB species, thereby allowing controlled cross-hybridization. For proof of concept, the microarray composed of 2,269 30-mer oligonucleotides focused on LAB species that are prevalent in sourdough ecosystems. Validation hybridizations using DNA and RNA from 18 LAB strains, covering 86% of all the oligonucleotides, showed that there were wide ranges in intensity and high reproducibility between microarrays.

Thermodynamic scaling behavior in genechips.

BACKGROUND: Affymetrix Genechips are characterized by probe pairs, a perfect match (PM) and a mismatch (MM) probe differing by a single nucleotide. Most of the data preprocessing algorithms neglect MM signals, as it was shown that MMs cannot be used as estimators of the non-specific hybridization as originally proposed by Affymetrix. The aim of this paper is to study in detail on a large number of experiments the behavior of the average PM/MM ratio. This is taken as an indicator of the quality of the hybridization and, when compared between different chip series, of the quality of the chip design. RESULTS: About 250 different GeneChip hybridizations performed at the VIB Microarray Facility for Homo sapiens, Drosophila melanogaster, and Arabidopsis thaliana were analyzed. The investigation of such a large set of data from the same source minimizes systematic experimental variations that may arise from differences in protocols or from different laboratories. The PM/MM ratios are derived theoretically from thermodynamic laws and a link is made with the sequence of PM and MM probe, more specifically with their central nucleotide triplets. CONCLUSION: The PM/MM ratios subdivided according to the different central nucleotides triplets follow qualitatively those deduced from the hybridization free energies in solution. It is shown also that the PM and MM histograms are related by a simple scale transformation, in agreement with what is to be expected from hybridization thermodynamics. Different quantitative behavior is observed on the different chip organisms analyzed, suggesting that some organism chips have superior probe design compared to others.

The transcriptional repressor ARR1-SRDX suppresses pleiotropic cytokinin activities in Arabidopsis.

The signal transduction of the phytohormone cytokinin is mediated by a multistep histidine-to-aspartate phosphorelay system. One component of this system are B-type response regulators, transcription factors mediating at least part of the response to cytokinin. In planta functional analysis of this family is hampered by the high level of functional redundancy of its 11 members. We generated a dominant repressor version of the Arabidopsis (Arabidopsis thaliana) response regulator ARR1 (ARR1-SRDX) using chimeric repressor silencing technology in order to study the extent of the contribution of B-type response regulators to cytokinin activities. In a protoplast test system, ARR1-SRDX suppressed ARR6:beta-glucuronidase reporter gene activation by different B-type ARRs. 35S:ARR1-SRDX transgenic Arabidopsis plants showed phenotypic changes reminiscent of plants with a reduced cytokinin status, such as a strongly reduced leaf size, an enhanced root system, and larger seeds. Several bioassays showed that 35S:ARR1-SRDX plants have an increased resistance toward cytokinin. The rapid induction of a large part of the cytokinin response genes was dampened. The transcript levels of more than 500 genes were more than 2.5-fold reduced in 35S:ARR1-SRDX transgenic seedlings, suggesting a broad function of B-type ARRs. Collectively, the suppression of pleiotropic cytokinin activities by a dominant repressor version of a B-type ARR indicates that this protein family is involved in mediating most, if not all, of the cytokinin activities in Arabidopsis. In addition, a role for B-type ARRs in mediating cross talk with other pathways is supported by the resistance of 35S:ARR1-SRDX seeds to phytochrome B-mediated inhibition of germination by far-red light. This study demonstrates the usefulness of chimeric repressor silencing technology to overcome redundancy in transcription factor families for functional studies.

CATMA, a comprehensive genome-scale resource for silencing and transcript profiling of Arabidopsis genes.

BACKGROUND: The Complete Arabidopsis Transcript MicroArray (CATMA) initiative combines the efforts of laboratories in eight European countries 1 to deliver gene-specific sequence tags (GSTs) for the Arabidopsis research community. The CATMA initiative offers the power and flexibility to regularly update the GST collection according to evolving knowledge about the gene repertoire. These GST amplicons can easily be reamplified and shared, subsets can be picked at will to print dedicated arrays, and the GSTs can be cloned and used for other functional studies. This ongoing initiative has already produced approximately 24,000 GSTs that have been made publicly available for spotted microarray printing and RNA interference. RESULTS: GSTs from the CATMA version 2 repertoire (CATMAv2, created in 2002) were mapped onto the gene models from two independent Arabidopsis nuclear genome annotation efforts, TIGR5 and PSB-EuGène, to consolidate a list of genes that were targeted by previously designed CATMA tags. A total of 9,027 gene models were not tagged by any amplified CATMAv2 GST, and 2,533 amplified GSTs were no longer predicted to tag an updated gene model. To validate the efficacy of GST mapping criteria and design rules, the predicted and experimentally observed hybridization characteristics associated to GST features were correlated in transcript profiling datasets obtained with the CATMAv2 microarray, confirming the reliability of this platform. To complete the CATMA repertoire, all 9,027 gene models for which no GST had yet been designed were processed with an adjusted version of the Specific Primer and Amplicon Design Software (SPADS). A total of 5,756 novel GSTs were designed and amplified by PCR from genomic DNA. Together with the pre-existing GST collection, this new addition constitutes the CATMAv3 repertoire. It comprises 30,343 unique amplified sequences that tag 24,202 and 23,009 protein-encoding nuclear gene models in the TAIR6 and EuGène genome annotations, respectively. To cover the remaining untagged genes, we identified 543 additional GSTs using less stringent design criteria and designed 990 sequence tags matching multiple members of gene families (Gene Family Tags or GFTs) to cover any remaining untagged genes. These latter 1,533 features constitute the CATMAv4 addition. CONCLUSION: To update the CATMA GST repertoire, we designed 7,289 additional sequence tags, bringing the total number of tagged TAIR6-annotated Arabidopsis nuclear protein-coding genes to 26,173. This resource is used both for the production of spotted microarrays and the large-scale cloning of hairpin RNA silencing vectors. All information about the resulting updated CATMA repertoire is available through the CATMA database http://www.catma.org.

The effect of a sertoli cell-selective knockout of the androgen receptor on testicular gene expression in prepubertal mice.

To unravel the molecular mechanisms mediating the effects of androgens on spermatogenesis, testicular gene expression was compared in mice with Sertoli cell-selective androgen receptor knockout (SCARKO) and littermate controls on postnatal d 10. Microarray analysis identified 692 genes with significant differences in expression. Of these, 28 appeared to be down-regulated and 12 up-regulated at least 2-fold in SCARKOs compared with controls. For nine of the more than 2-fold down-regulated genes, androgen regulation was confirmed by treatment of wild-type mice with an antiandrogen (flutamide). Some of them were previously described to be androgen regulated or essential for spermatogenesis. Serine-type protease inhibitors were markedly overrepresented in this down-regulated subgroup. A time study (d 8-20), followed by cluster analysis, allowed identification of distinct expression patterns of differentially expressed genes. Three genes with a pattern closely resembling that of Pem, a prototypical androgen-regulated gene expressed in Sertoli cells, were selected for confirmation by quantitative RT-PCR and additional analysis. The data confirm that the SCARKO model allows identification of novel androgen-regulated genes in the testis. Moreover, they suggest that protease inhibitors and other proteins related to tubular restructuring and cell junction dynamics may be controlled in part by androgens.

Benchmarking the CATMA microarray. A novel tool for Arabidopsis transcriptome analysis.

Transcript profiling is crucial to study biological systems, and various platforms have been implemented to survey mRNAs at the genome scale. We have assessed the performance of the CATMA microarray designed for Arabidopsis (Arabidopsis thaliana) transcriptome analysis and compared it with the Agilent and Affymetrix commercial platforms. The CATMA array consists of gene-specific sequence tags of 150 to 500 bp, the Agilent (Arabidopsis 2) array of 60mer oligonucleotides, and the Affymetrix gene chip (ATH1) of 25mer oligonucleotide sets. We have matched each probe repertoire with the Arabidopsis genome annotation (The Institute for Genomic Research release 5.0) and determined the correspondence between them. Array performance was analyzed by hybridization with labeled targets derived from eight RNA samples made of shoot total RNA spiked with a calibrated series of 14 control transcripts. CATMA arrays showed the largest dynamic range extending over three to four logs. Agilent and Affymetrix arrays displayed a narrower range, presumably because signal saturation occurred for transcripts at concentrations beyond 1,000 copies per cell. Sensitivity was comparable for all three platforms. For Affymetrix GeneChip data, the RMA software package outperformed Microarray Suite 5.0 for all investigated criteria, confirming that the information provided by the mismatch oligonucleotides has no added value. In addition, taking advantage of replicates in our dataset, we conducted a robust statistical analysis of the platform propensity to yield false positive and false negative differentially expressed genes, and all gave satisfactory results. The results establish the CATMA array as a mature alternative to the Affymetrix and Agilent platforms.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.