Study of gene expression in thyrotropin-stimulated thyroid cells by cDNA expression array: ID3 transcription modulating factor as an early response protein and tumor marker in thyroid carcinomas.

Induction of cell proliferation by mitogen or growth factor stimulation leads to the specific induction or repression of a large number of genes. To identify genes differentially regulated by the cAMP-dependent transduction pathway, which is poorly characterized so far, we used the cDNA expression array technology. Hybridizations of Atlas human cDNA expression arrays with (32)P-labeled cDNA probes derived from control or thyrotropin (TSH)-stimulated dog thyrocytes in primary culture generated expression profiles of hundreds of genes simultaneously. Among the genes that displayed modified expression, we selected the transcription factor ID3, whose expression was increased by a cAMP-dependent stimulus. ID3 overexpression after TSH stimulation was first verified by Northern blotting analysis, and its mRNA regulation was then investigated in response to a variety of agents acting on thyrocyte proliferation and/or differentiation. We show that: (1) ID3 mRNA induction was stronger after stimulation of the cAMP cascade, but was not restricted to this signaling pathway, as phorbol myristate ester (TPA) and insulin also stimulated mRNA accumulation; (2) in contrast, powerful mitogens for thyroid cells, epidermal growth factor and hepatocyte growth factor, did not significantly modify ID3 mRNA levels; (3) ID3 protein levels closely parallelled mRNA levels, as revealed by immunofluorescence experiments showing a nuclear signal regulated by TSH; (4) in papillary thyroid carcinomas, ID3 mRNA was downregulated. Our results suggest that ID3 expression might be more related to the differentiating process induced by TSH than to the proliferative action of this hormone.

Molecular diagnosis for heterogeneous genetic diseases with high-throughput DNA sequencing applied to retinitis pigmentosa

BACKGROUND:
The genetic heterogeneity of many Mendelian disorders, such as retinitis pigmentosa which results from mutations in over 40 genes, is a major obstacle to obtaining a molecular diagnosis in clinical practice. Targeted high-throughput DNA sequencing offers a potential solution and was used to develop a molecular diagnostic screen for patients with retinitis pigmentosa.
METHODS:
A custom sequence capture array was designed to target the coding regions of all known retinitis pigmentosa genes and used to enrich these sequences from DNA samples of five patients. Enriched DNA was subjected to high-throughput sequencing singly or in pools, and sequence variants were identified by alignment of up to 10 million reads per sample to the normal reference sequence. Potential pathogenicity was assessed by functional predictions and frequency in controls.
RESULTS AND CONCLUSIONS:
Known homozygous PDE6B and compound heterozygous CRB1 mutations were detected in two patients. A novel homozygous missense mutation (c.2957A?T; p.N986I) in the cyclic nucleotide gated channel ß1 (CNGB1) gene predicted to have a deleterious effect and absent in 720 control chromosomes was detected in one case in which conventional genetic screening had failed to detect mutations. The detection of known and novel retinitis pigmentosa mutations in this study establishes high-throughput DNA sequencing with DNA pooling as an effective diagnostic tool for heterogeneous genetic diseases.

Distinct methylation patterns in genes that affect mitochondrial function are associated with kidney disease in blood-derived DNA from individuals with Type 1 diabetes

AIMS: Epigenetic modifications, such as DNA methylation, can influence the risk of developing kidney disease. We studied methylation profiles in genes related to mitochondrial function to assess whether differences in these epigenetic features were associated with diabetic kidney disease in people with Type 1 diabetes.

METHODS: A case-control association study was undertaken (n = 196 individuals with diabetic kidney disease vs. n = 246 individuals without renal disease). Participants were White and diagnosed with Type 1 diabetes before 31 years of age. Genes that encode mitochondrial proteins (n = 780) were downloaded from mitoproteome. org. DNA methylation profiles from blood-derived DNA were generated using the Illumina Infinium HumanMethylation450 (262 samples) and Illumina Infinium HumanMethylation27 (192 samples) arrays. Beta values (β) were calculated and quality control was conducted, including evaluating blind duplicate DNA samples.

RESULTS: Fifty-four Cytosine-phosphate-Guanine probes across 51 unique genes were significantly associated (P ≤ 10(-8) ) with diabetic kidney disease across both the 450K and the 27K methylation arrays. A subanalysis, employing the 450K array, identified 755 Cytosine-phosphate-Guanine probes in 374 genes that were significantly associated (P ≤ 10(-8) ) with end-stage renal disease. Forty-six of the top-ranked variants for diabetic kidney disease were also identified as being differentially methylated in individuals with end-stage renal disease. The largest change in methylation (Δβ = 0.2) was observed for cg03169527 in the TAMM41 gene, chromosome 3p25.2. Three genes, PMPCB, TSFM and AUH, were observed with differential methylation at multiple Cytosine-phosphate-Guanine sites each (P < 10(-12) ).

CONCLUSIONS: Differential methylation in genes that influence mitochondrial function are associated with kidney disease in individuals with Type 1 diabetes. 

AURKA overexpression accompanies dysregulation of DNA-damage response genes in invasive urothelial cell carcinoma

Invasive urothelial cell carcinoma (UCC) is characterized by increased chromosomal instability and follows an aggressive clinical course in contrast to non-invasive disease. To identify molecular processes that confer and maintain an aggressive malignant phenotype, we used a high-throughput genome-wide approach to interrogate a cohort of high and low clinical risk UCC tumors. Differential expression analyses highlighted cohesive dysregulation of critical genes involved in the G(2)/M checkpoint in aggressive UCC. Hierarchical clustering based on DNA Damage Response (DDR) genes separated tumors according to a pre-defined clinical risk phenotype. Using array-comparative genomic hybridization, we confirmed that the DDR was disrupted in tumors displaying high genomic instability. We identified DNA copy number gains at 20q13.2-q13.3 (AURKA locus) and determined that overexpression of AURKA accompanied dysregulation of DDR genes in high risk tumors. We postulated that DDR-deficient UCC tumors are advantaged by a selective pressure for AURKA associated override of M phase barriers and confirmed this in an independent tissue microarray series. This mechanism that enables cancer cells to maintain an aggressive phenotype forms a rationale for targeting AURKA as a therapeutic strategy in advanced stage UCC.

A new look towards BAC-based array CGH through a comprehensive comparison with oligo-based array CGH

BACKGROUND: Currently, two main technologies are used for screening of DNA copy number; the BAC (Bacterial Artificial Chromosome) and the recently developed oligonucleotide-based CGH (Chromosomal Comparative Genomic Hybridization) arrays which are capable of detecting small genomic regions with amplification or deletion. The correlation as well as the discriminative power of these platforms has never been compared statistically on a significant set of human patient samples.

RESULTS: In this paper, we present an exhaustive comparison between the two CGH platforms, undertaken at two independent sites using the same batch of DNA from 19 advanced prostate cancers. The comparison was performed directly on the raw data and a significant correlation was found between the two platforms. The correlation was greatly improved when the data were averaged over large chromosomic regions using a segmentation algorithm. In addition, this analysis has enabled the development of a statistical model to discriminate BAC outliers that might indicate microevents. These microevents were validated by the oligo platform results.

CONCLUSION: This article presents a genome-wide statistical validation of the oligo array platform on a large set of patient samples and demonstrates statistically its superiority over the BAC platform for the Identification of chromosomic events. Taking advantage of a large set of human samples treated by the two technologies, a statistical model has been developed to show that the BAC platform could also detect microevents.

DNA microarray profiling of genes differentially regulated by the histone deacetylase inhibitors vorinostat and LBH589 in colon cancer cell lines

BACKGROUND: Despite the significant progress made in colon cancer chemotherapy, advanced disease remains largely incurable and novel efficacious chemotherapies are urgently needed. Histone deacetylase inhibitors (HDACi) represent a novel class of agents which have demonstrated promising preclinical activity and are undergoing clinical evaluation in colon cancer. The goal of this study was to identify genes in colon cancer cells that are differentially regulated by two clinically advanced hydroxamic acid HDACi, vorinostat and LBH589 to provide rationale for novel drug combination partners and identify a core set of HDACi-regulated genes.

METHODS: HCT116 and HT29 colon cancer cells were treated with LBH589 or vorinostat and growth inhibition, acetylation status and apoptosis were analyzed in response to treatment using MTS, Western blotting and flow cytometric analyses. In addition, gene expression was analyzed using the Illumina Human-6 V2 BeadChip array and Ingenuity Pathway Analysis.

RESULTS: Treatment with either vorinostat or LBH589 rapidly induced histone acetylation, cell cycle arrest and inhibited the growth of both HCT116 and HT29 cells. Bioinformatic analysis of the microarray profiling revealed significant similarity in the genes altered in expression following treatment with the two HDACi tested within each cell line. However, analysis of genes that were altered in expression in the HCT116 and HT29 cells revealed cell-line-specific responses to HDACi treatment. In addition a core cassette of 11 genes modulated by both vorinostat and LBH589 were identified in both colon cancer cell lines analyzed.

CONCLUSION: This study identified HDACi-induced alterations in critical genes involved in nucleotide metabolism, angiogenesis, mitosis and cell survival which may represent potential intervention points for novel therapeutic combinations in colon cancer. This information will assist in the identification of novel pathways and targets that are modulated by HDACi, providing much-needed information on HDACi mechanism of action and providing rationale for novel drug combination partners. We identified a core signature of 11 genes which were modulated by both vorinostat and LBH589 in a similar manner in both cell lines. These core genes will assist in the development and validation of a common gene set which may represent a molecular signature of HDAC inhibition in colon cancer.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

SCAPA (Spacially Addressable Protein Array) – a Novel Protein Array for the Differential Profiling of Ligand-receptor Induced Signalling

While protein microarray technology has been successful in demonstrating its usefulness for large scale high-throughput proteome profiling, performance of antibody/antigen microarrays has been only moderately productive. Immobilization of either the capture antibodies or the protein samples on solid supports has severe drawbacks. Denaturation of the immobilized proteins as well as inconsistent orientation of antibodies/ligands on the arrays can lead to erroneous results. This has prompted a number of studies to address these challenges by immobilizing proteins on biocompatible surfaces, which has met with limited success. Our strategy relates to a multiplexed, sensitive and high-throughput method for the screening quantification of intracellular signalling proteins from a complex mixture of proteins. Each signalling protein to be monitored has its capture moiety linked to a specific oligo ‘tag’. The array involves the oligonucleotide hybridization-directed localization and identification of different signalling proteins simultaneously, in a rapid and easy manner. Antibodies have been used as the capture moieties for specific identification of each signaling protein. The method involves covalently partnering each antibody/protein molecule with a unique DNA or DNA derivatives oligonucleotide tag that directs the antibody to a unique site on the microarray due to specific hybridization with a complementary tag-probe on the array. Particular surface modifications and optimal conditions allowed high signal to noise ratio which is essential to the success of this approach.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Array based detection of antibiotic resistance genes in Gram negative bacteria isolated from retail poultry meat in the UK and Ireland

The use of antibiotics in birds and animals intended for human consumption within the European Union (EU) and elsewhere has been subject to regulation prohibiting the use of antimicrobials as growth promoters and the use of last resort antibiotics in an attempt to reduce the spread of multi-resistant Gram negative bacteria. Given the inexorable spread of antibiotic resistance there is an increasing need for improved monitoring of our food. Using selective media, Gram negative bacteria were isolated from retail chicken of UK-Intensively reared (n = 27), Irish-Intensively reared (n = 19) and UK-Free range (n = 30) origin and subjected to an oligonucleotide based array system for the detection of 47 clinically relevant antibiotic resistance genes (ARGs) and two integrase genes. High incidences of β-lactamase genes were noted in all sample types, acc (67%), cmy (80%), fox (55%) and tem (40%) while chloramphenicol resistant determinants were detected in bacteria from the UK poultry portions and were absent in bacteria from the Irish samples. Denaturing Gradient Gel Electrophoresis (DGGE) was used to qualitatively analyse the Gram negative population in the samples and showed the expected diversity based on band stabbing and DNA sequencing. The array system proved to be a quick method for the detection of antibiotic resistance gene (ARG) burden within a mixed Gram negative bacterial population.

Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair

Recurrent submicroscopic genomic copy number changes are the result of nonallelic homologous recombination (NAHR). Nonrecurrent aberrations, however, can result from different nonexclusive recombination-repair mechanisms. We previously described small microduplications at Xq28 containing MECP2 in four male patients with a severe neurological phenotype. Here, we report on the fine-mapping and breakpoint analysis of 16 unique microduplications. The size of the overlapping copy number changes varies between 0.3 and 2.3 Mb, and FISH analysis on three patients demonstrated a tandem orientation. Although eight of the 32 breakpoint regions coincide with low-copy repeats, none of the duplications are the result of NAHR. Bioinformatics analysis of the breakpoint regions demonstrated a 2.5-fold higher frequency of Alu interspersed repeats as compared with control regions, as well as a very high GC content (53%). Unexpectedly, we obtained the junction in only one patient by long-range PCR, which revealed nonhomologous end joining as the mechanism. Breakpoint analysis in two other patients by inverse PCR and subsequent array comparative genomic hybridization analysis demonstrated the presence of a second duplicated region more telomeric at Xq28, of which one copy was inserted in between the duplicated MECP2 regions. These data suggest a two-step mechanism in which part of Xq28 is first inserted near the MECP2 locus, followed by breakage-induced replication with strand invasion of the normal sister chromatid. Our results indicate that the mechanism by which copy number changes occur in regions with a complex genomic architecture can yield complex rearrangements.

Production of DNA microarray and expression analysis of genes from Xylella fastidiosa in different culture media

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Implementation of DNA markers to produce genomically - Enhanced EPDs in nellore cattle

Background: The sequencing and publication of the cattle genome and the identification of single nucleotide polymorphism (SNP) molecular markers have provided new tools for animal genetic evaluation and genomic-enhanced selection. These new tools aim to increase the accuracy and scope of selection while decreasing generation interval. The objective of this study was to evaluate the enhancement of accuracy caused by the use of genomic information (Clarifide® - Pfizer) on genetic evaluation of Brazilian Nellore cattle. Review: The application of genome-wide association studies (GWAS) is recognized as one of the most practical approaches to modern genetic improvement. Genomic selection is perhaps most suited to the improvement of traits with low heritability in zebu cattle. The primary interest in livestock genomics has been to estimate the effects of all the markers on the chip, conduct cross-validation to determine accuracy, and apply the resulting information in GWAS either alone [9] or in combination with bull test and pedigree-based genetic evaluation data. The cost of SNP50K genotyping however limits the commercial application of GWAS based on all the SNPs on the chip. However, reasonable predictability and accuracy can be achieved in GWAS by using an assay that contains an optimally selected predictive subset of markers, as opposed to all the SNPs on the chip. The best way to integrate genomic information into genetic improvement programs is to have it included in traditional genetic evaluations. This approach combines traditional expected progeny differences based on phenotype and pedigree with the genomic breeding values based on the markers. Including the different sources of information into a multiple trait genetic evaluation model, for within breed dairy cattle selection, is working with excellent results. However, given the wide genetic diversity of zebu breeds, the high-density panel used for genomic selection in dairy cattle (Ilumina Bovine SNP50 array) appears insufficient for across-breed genomic predictions and selection in beef cattle. Today there is only one breed-specific targeted SNP panel and genomic predictions developed using animals across the entire population of the Nellore breed (www.pfizersaudeanimal.com), which enables genomically - enhanced selection. Genomic profiles are a way to enhance our current selection tools to achieve more accurate predictions for younger animals. Material and Methods: We analyzed the age at first calving (AFC), accumulated productivity (ACP), stayability (STAY) and heifer pregnancy at 30 months (HP30) in Nellore cattle fitting two different animal models; 1) a traditional single trait model, and 2) a two-trait model where the genomic breeding value or molecular value prediction (MVP) was included as a correlated trait. All mixed model analyses were performed using the statistical software ASREML 3.0. Results: Genetic correlation estimates between AFC, ACP, STAY, HP30 and respective MVPs ranged from 0.29 to 0.46. Results also showed an increase of 56%, 36%, 62% and 19% in estimated accuracy of AFC, ACP, STAY and HP30 when MVP information was included in the animal model. Conclusion: Depending upon the trait, integration of MVP information into genetic evaluation resulted in increased accuracy of 19% to 62% as compared to accuracy from traditional genetic evaluation. GE-EPD will be an effective tool to enable faster genetic improvement through more dependable selection of young animals.

Genes candidatos a marcadores tumorais na progressão do adenocarcinoma de próstata indentificados por análise de HR-CGH e CGH-ARRAY

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)