Survival Analysis with Large Dimensional Covariates: An Application in Microarray Studies

Use of microarray technology often leads to high-dimensional and low- sample size data settings. Over the past several years, a variety of novel approaches have been proposed for variable selection in this context. However, only a small number of these have been adapted for time-to-event data where censoring is present. Among standard variable selection methods shown both to have good predictive accuracy and to be computationally efficient is the elastic net penalization approach. In this paper, adaptation of the elastic net approach is presented for variable selection both under the Cox proportional hazards model and under an accelerated failure time (AFT) model. Assessment of the two methods is conducted through simulation studies and through analysis of microarray data obtained from a set of patients with diffuse large B-cell lymphoma where time to survival is of interest. The approaches are shown to match or exceed the predictive performance of a Cox-based and an AFT-based variable selection method. The methods are moreover shown to be much more computationally efficient than their respective Cox- and AFT- based counterparts.

A radiation hybrid map of sheep chromosome 23 based on ovine BAC-end sequences

More than 375,000 BAC-end sequences (BES) of the CHORI-243 ovine BAC library have been deposited in public databases. blastn searches with these BES against HSA18 revealed 1806 unique and significant hits. We used blastn-anchored BES for an in silico prediction of gene content and chromosome assignment of comparatively mapped ovine BAC clones. Ovine BES were selected at approximately 1.3-Mb intervals of HSA18 and incorporated into a human-sheep comparative map. An ovine 5000-rad whole-genome radiation hybrid panel (USUoRH5000) was typed with 70 markers, all of which mapped to OAR23. The resulting OAR23 RH map included 43 markers derived from BES with high and unique BLAST hits to the sequence of the orthologous HSA18, nine EST-derived markers, 16 microsatellite markers taken from the ovine linkage map and two bovine microsatellite markers. Six new microsatellite markers derived from the 43 mapped BES and the two bovine microsatellite markers were linkage-mapped using the International Mapping Flock (IMF). Thirteen additional microsatellite markers were derived from other ovine BES with high and unique BLAST hits to the sequence of the orthologous HSA18 and also positioned on the ovine linkage map but not incorporated into the OAR23 RH map. This resulted in 24 markers in common and in the same order between the RH and linkage maps. Eight of the BES-derived markers were mapped using fluorescent in situ hybridization (FISH), to thereby align the RH and cytogenetic maps. Comparison of the ovine chromosome 23 RH map with the HSA18 map identified and localized three major breakpoints between HSA18 and OAR23. The positions of these breakpoints were equivalent to those previously shown for syntenic BTA24 and HSA18. This study presents evidence for the usefulness of ovine BES when constructing a high-resolution comprehensive map for a single sheep chromosome. The comparative analysis confirms and refines knowledge about chromosomal conservation and rearrangements between sheep, cattle and human. The constructed RH map demonstrates the resolution and utility of the newly constructed ovine RH panel.

The horse genome project - sequence based insights into male reproductive mechanisms

The growing knowledge on physiology, cell biology and biochemistry of the reproductive organs has provided many insights into molecular mechanisms that are required for successful reproduction. Research directed at the investigation of reproduction physiology in domestic animals was hampered in the past by a lack of species-specific genomic information. The genome sequences of dog, cattle and horse have become publicly available in 2005, 2006 and 2007 respectively. Although the gene content of mammalian genomes is generally very similar, genes involved in reproduction tend to be less conserved than the average mammalian gene. The availability of genome sequences provides a valuable resource to check whether any protein that may be known from human or mouse research is present in cattle and/or horse as well. Currently there are more than 200 genes known that are involved in the production of fertile sperm cells. Great progress has been made in the understanding of genetic aberrations that lead to male infertility. Additionally, the first genetic mechanisms are being discovered that contribute to the quantitative variation of fertility traits in fertile male animals. Here, I will review some selected aspects of genetic research in male fertility and offer some perspectives for the use of genomic sequence information.

Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Differential response of porcine osteoblasts and chondrocytes in cell or tissue culture after 5-aminolevulinic acid-based photodynamic therapy

OBJECTIVE: Outcome in osteochondral allografting is limited by the immunological incompatibility of the grafted tissue. Based on a resistance of chondrocytes to photodynamic therapy in cell culture it is proposed that 5-aminolevulinic acid-based photodynamic therapy (5-ALA-PDT) might be used to inactivate bone while maintaining viability of chondrocytes and thus immunomodulate bone selectively. METHODS: Chondrocytes and osteoblasts from porcine humeral heads were either isolated (cell culture) or treated in situ (tissue culture). To quantify cytotoxic effects of 5-ALA-PDT (0-20J/cm(2), 100mW/cm(2)) an (3-(4,5-dimethylthiazol-2-yl)-2,5-di-phenyltetrazolium bromide) (MTT)-assay was used in cell culture and in situ hybridization in tissue culture to assess metabolic active cells (functional osteoblasts: colalpha(1)(I) mRNA, functional chondrocytes: colalpha(1)(II) mRNA). RESULTS: In cell culture, survival after 5-ALA-PDT was significantly higher for chondrocytes (5J/cm(2): 87+/-12% compared to untreated cells) than for osteoblasts (5J/cm(2): 12+/-11%). In tissue culture, the percentage of functional chondrocytes in cartilage showed a decrease after 5-ALA-PDT (direct fixation: 92+/-2%, 20J/cm(2): 35+/-15%; P<0.0001). A significant decrease in the percentage of bone surfaces covered by functional osteoblasts was observed in freshly harvested (31+/-3%) compared to untreated tissues maintained in culture (11+/-4%, P<0.0001), with no further decrease after 5-ALA-PDT. CONCLUSION: Chondrocytes were more resistant to 5-ALA-PDT than osteoblasts in cell culture, while in tissue culture a loss of functional chondrocytes was observed after 5-ALA-PDT. Since osteoblasts - but not chondrocytes - were sensitive to the tissue culture conditions, devitalized bone with functional cartilage might already be achieved by applying specific tissue culture conditions even without 5-ALA-PDT.

Hypermethylation of the DPYD promoter region is not a major predictor of severe toxicity in 5-fluorouracil based chemotherapy

BACKGROUND: The activity of dihydropyrimidine dehydrogenase (DPD), the key enzyme of pyrimidine catabolism, is thought to be an important determinant for the occurrence of severe toxic reactions to 5-fluorouracil (5-FU), which is one of the most commonly prescribed chemotherapeutic agents for the treatment of solid cancers. Genetic variation in the DPD gene (DPYD) has been proposed as a main factor for variation in DPD activity in the population. However, only a small proportion of severe toxicities in 5-FU based chemotherapy can be explained with such rare deleterious DPYD mutations resulting in severe enzyme deficiencies. Recently, hypermethylation of the DPYD promoter region has been proposed as an alternative mechanism for DPD deficiency and thus as a major cause of severe 5-FU toxicity. METHODS: Here, the prognostic significance of this epigenetic marker with respect to severe 5-FU toxicity was assessed in 27 cancer patients receiving 5-FU based chemotherapy, including 17 patients experiencing severe toxic side effects following drug administration, none of which were carriers of a known deleterious DPYD mutation, and ten control patients. The methylation status of the DPYD promoter region in peripheral blood mononuclear cells was evaluated by analysing for each patient between 19 and 30 different clones of a PCR-amplified 209 base pair fragment of the bisulfite-modified DPYD promoter region. The fragments were sequenced to detect bisulfite-induced, methylation-dependent sequence differences. RESULTS: No evidence of DPYD promoter methylation was observed in any of the investigated patient samples, whereas in a control experiment, as little as 10% methylated genomic DNA could be detected. CONCLUSION: Our results indicate that DYPD promoter hypermethylation is not of major importance as a prognostic factor for severe toxicity in 5-FU based chemotherapy.

Systemic analysis of microarray data sets towards identifying genes regulating root development

Nitrogen and water are essential for plant growth and development. In this study, we designed experiments to produce gene expression data of poplar roots under nitrogen starvation and water deprivation conditions. We found low concentration of nitrogen led first to increased root elongation followed by lateral root proliferation and eventually increased root biomass. To identify genes regulating root growth and development under nitrogen starvation and water deprivation, we designed a series of data analysis procedures, through which, we have successfully identified biologically important genes. Differentially Expressed Genes (DEGs) analysis identified the genes that are differentially expressed under nitrogen starvation or drought. Protein domain enrichment analysis identified enriched themes (in same domains) that are highly interactive during the treatment. Gene Ontology (GO) enrichment analysis allowed us to identify biological process changed during nitrogen starvation. Based on the above analyses, we examined the local Gene Regulatory Network (GRN) and identified a number of transcription factors. After testing, one of them is a high hierarchically ranked transcription factor that affects root growth under nitrogen starvation. It is very tedious and time-consuming to analyze gene expression data. To avoid doing analysis manually, we attempt to automate a computational pipeline that now can be used for identification of DEGs and protein domain analysis in a single run. It is implemented in scripts of Perl and R.

Analysis of independent microarray datasets of renal biopsies identifies a robust transcript signature of acute allograft rejection

Transcriptomics could contribute significantly to the early and specific diagnosis of rejection episodes by defining 'molecular Banff' signatures. Recently, the description of pathogenesis-based transcript sets offered a new opportunity for objective and quantitative diagnosis. Generating high-quality transcript panels is thus critical to define high-performance diagnostic classifier. In this study, a comparative analysis was performed across four different microarray datasets of heterogeneous sample collections from two published clinical datasets and two own datasets including biopsies for clinical indication, and samples from nonhuman primates. We characterized a common transcriptional profile of 70 genes, defined as acute rejection transcript set (ARTS). ARTS expression is significantly up-regulated in all AR samples as compared with stable allografts or healthy kidneys, and strongly correlates with the severity of Banff AR types. Similarly, ARTS were tested as a classifier in a large collection of 143 independent biopsies recently published by the University of Alberta. Results demonstrate that the 'in silico' approach applied in this study is able to identify a robust and reliable molecular signature for AR, supporting a specific and sensitive molecular diagnostic approach for renal transplant monitoring.

The NF-{kappa}B negative regulator TNFAIP3 (A20) is inactivated by somatic mutations and genomic deletions in marginal zone lymphomas

Unique and shared cytogenetic abnormalities have been documented for marginal zone lymphomas (MZLs) arising at different sites. Recently, homozygous deletions of the chromosomal band 6q23, involving the tumor necrosis factor alpha-induced protein 3 (TNFAIP3, A20) gene, a negative regulator of NF-kappaB, were described in ocular adnexal MZL, suggesting a role for A20 as a tumor suppressor in this disease. Here, we investigated inactivation of A20 by DNA mutations or deletions in a panel of extranodal MZL (EMZL), nodal MZL (NMZL), and splenic MZL (SMZL). Inactivating mutations encoding truncated A20 proteins were identified in 6 (19%) of 32 MZLs, including 2 (18%) of 11 EMZLs, 3 (33%) of 9 NMZLs, and 1 (8%) of 12 SMZLs. Two additional unmutated nonsplenic MZLs also showed monoallelic or biallelic A20 deletions by fluorescent in situ hybridization (FISH) and/or SNP-arrays. Thus, A20 inactivation by either somatic mutation and/or deletion represents a common genetic aberration across all MZL subtypes, which may contribute to lymphomagenesis by inducing constitutive NF-kappaB activation.

Phylogeny of the family Pasteurellaceae based on rpoB sequences

Sequences of the gene encoding the beta-subunit of the RNA polymerase (rpoB) were used to delineate the phylogeny of the family Pasteurellaceae. A total of 72 strains, including the type strains of the major described species as well as selected field isolates, were included in the study. Selection of universal rpoB-derived primers for the family allowed straightforward amplification and sequencing of a 560 bp fragment of the rpoB gene. In parallel, 16S rDNA was sequenced from all strains. The phylogenetic tree obtained with the rpoB sequences reflected the major branches of the tree obtained with the 16S rDNA, especially at the genus level. Only a few discrepancies between the trees were observed. In certain cases the rpoB phylogeny was in better agreement with DNA-DNA hybridization studies than the phylogeny derived from 16S rDNA. The rpoB gene is strongly conserved within the various species of the family of Pasteurellaceae. Hence, rpoB gene sequence analysis in conjunction with 16S rDNA sequencing is a valuable tool for phylogenetic studies of the Pasteurellaceae and may also prove useful for reorganizing the current taxonomy of this bacterial family.

Fluorescent In Situ hybridization: a new tool for the direct identification and detection of F. psychrophilum

F. psychrophilum is the causative agent of Bacterial Cold Water Disease (BCW) and Rainbow Trout Fry Syndrome (RTFS). To date, diagnosis relies mainly on direct microscopy or cultural methods. Direct microscopy is fast but not very reliable, whereas cultural methods are reliable but time-consuming and labor-intensive. So far fluorescent in situ hybridization (FISH) has not been used in the diagnosis of flavobacteriosis but it has the potential to rapidly and specifically detect F. psychrophilum in infected tissues. Outbreaks in fish farms, caused by pathogenic strains of Flavobacterium species, are increasingly frequent and there is a need for reliable and cost-effective techniques to rapidly diagnose flavobacterioses. This study is aimed at developing a FISH that could be used for the diagnosis of F. psychrophilum infections in fish. We constructed a generic probe for the genus Flavobacterium ("Pan-Flavo") and two specific probes targeting F. psychrophilum based on 16S rRNA gene sequences. We tested their specificity and sensitivity on pure cultures of different Flavobacterium and other aquatic bacterial species. After assessing their sensitivity and specificity, we established their limit of detection and tested the probes on infected fresh tissues (spleen and skin) and on paraffin-embedded tissues. The results showed high sensitivity and specificity of the probes (100% and 91% for the Pan-Flavo probe and 100% and 97% for the F. psychrophilum probe, respectively). FISH was able to detect F. psychrophilum in infected fish tissues, thus the findings from this study indicate this technique is suitable as a fast and reliable method for the detection of Flavobacterium spp. and F. psychrophilum.

Mycoplasma feriruminatoris sp. nov., a fast growing Mycoplasma species isolated from wild Caprinae

Five Mycoplasma strains from wild Caprinae were analyzed: four from Alpine ibex (Capra ibex) which died at the Berlin Zoo between 1993 and 1994, one from a Rocky Mountain goat collected in the USA prior to 1987. These five strains represented a population different from the populations belonging to the 'Mycoplasma mycoides cluster' as tested using multi locus sequence typing, Matrix-assisted laser desorption/ionization time of flight mass spectrometry analysis and DNA-DNA hybridization. Analysis of the 16S rRNA gene (rrs), genomic sequence based in silico as well as laboratory DNA-DNA hybridization, and the analysis of phenotypic traits in particular their exceptionally rapid growth all confirmed that they do not belong to any Mycoplasma species described to date. We therefore suggest these strains represent a novel species, for which we propose the name Mycoplasma feriruminatoris sp. nov. The type strain is G5847(T) (=DSM 26019(T)=NCTC 1362(T)).

Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Diagnostic and prognostic implications of the PAX8-PPARγ translocation in thyroid carcinomas-a TMA-based study of 226 cases

AIMS Follicular thyroid carcinoma (FTC) has been a diagnostic challenge for decades. The PAX8-PPARγ rearrangement has been detected in FTC and classic papillary thyroid carcinomas (PTCs). The aims of this study were to assess the presence of PAX8-PPARγ by using tissue microarrays in a large cohort of different thyroid neoplasms, and to assess its diagnostic and prognostic implications. METHODS AND RESULTS Fluorescence in-situ hybridization (FISH) analysis for PAX8-PPARγ was performed on 226 thyroid tumours, comprising FTCs (n = 59), PTCs (n = 126), poorly differentiated thyroid carcinomas (PDs; n = 34), follicular thyroid adenomas (FTAs; n = 5), and follicular tumours of unknown malignant potential (FTUMPs; n = 2). PAX8-PPARγ was detected in 12% of FTCs, 1% of PTCs, 7% of PDs, and in both cases of FTUMP. There was no correlation between the extent of capsular or vascular invasion and PAX8-PPARγ, or between lymph node or haematogenous metastasis and PAX8-PPARγ. Overall survival (OS), tumour-specific survival (TSS) and relapse-free-survival (RFS) were not influenced by PAX8-PPARγ. CONCLUSIONS In this study, we demonstrate for the first time the presence of PAX8-PPARγ in PDs and FTUMPs, whereas in FTCs and PTCs the prevalence of PAX8-PPARγ is lower than previously reported. PAX8-PPARγ did not correlate with invasiveness or affect prognosis in any tumour type.

Prediction of outcome in patients with low-grade squamous intraepithelial lesions by fluorescence in situ hybridization analysis of human papillomavirus, TERC, and MYC

BACKGROUND Cytology is an excellent method with which to diagnose preinvasive lesions of the uterine cervix, but it suffers from limited specificity for clinically significant lesions. Supplementary methods might predict the natural course of the detected lesions. The objective of the current study was to test whether a multicolor fluorescence in situ hybridization (FISH) assay might help to stratify abnormal results of Papanicolaou tests. METHODS A total of 219 liquid-based cytology specimens of low-grade squamous intraepithelial lesions (LSIL), 49 atypical squamous cells of undetermined significance (ASCUS) specimens, 52 high-grade squamous intraepithelial lesion (HSIL) specimens, and 50 normal samples were assessed by FISH with probes for the human papillomavirus (HPV), MYC, and telomerase RNA component (TERC). Subtyping of HPV by polymerase chain reaction (PCR) was performed in a subset of cases (n=206). RESULTS There was a significant correlation found between HPV detection by FISH and PCR (P<.0001). In patients with LSILs, the presence of HPV detected by FISH was significantly associated with disease progression (P<.0001). An increased MYC and/or TERC gene copy number (>2 signals in>10% of cells) prevailed in 43% of ASCUS specimens and was more frequent in HSIL (85%) than in LSIL (33%) (HSIL vs LSIL: P<.0001). Increased TERC gene copy number was significantly correlated with progression of LSIL (P<.01; odds ratio, 7.44; area under the receiver operating characteristic curve, 0.73; positive predictive value, 0.30; negative predictive value, 0.94) CONCLUSIONS: The detection of HPV by FISH analysis is feasible in liquid-based cytology and is significantly correlated with HPV analysis by PCR. The analysis of TERC gene copy number may be useful for risk stratification in patients with LSIL.