A Faster Circular Binary Segmentation Algorithm for the Analysis of Array CGH Data

Motivation: Array CGH technologies enable the simultaneous measurement of DNA copy number for thousands of sites on a genome. We developed the circular binary segmentation (CBS) algorithm to divide the genome into regions of equal copy number (Olshen {\it et~al}, 2004). The algorithm tests for change-points using a maximal $t$-statistic with a permutation reference distribution to obtain the corresponding $p$-value. The number of computations required for the maximal test statistic is $O(N^2),$ where $N$ is the number of markers. This makes the full permutation approach computationally prohibitive for the newer arrays that contain tens of thousands markers and highlights the need for a faster. algorithm. Results: We present a hybrid approach to obtain the $p$-value of the test statistic in linear time. We also introduce a rule for stopping early when there is strong evidence for the presence of a change. We show through simulations that the hybrid approach provides a substantial gain in speed with only a negligible loss in accuracy and that the stopping rule further increases speed. We also present the analysis of array CGH data from a breast cancer cell line to show the impact of the new approaches on the analysis of real data. Availability: An R (R Development Core Team, 2006) version of the CBS algorithm has been implemented in the ``DNAcopy'' package of the Bioconductor project (Gentleman {\it et~al}, 2004). The proposed hybrid method for the $p$-value is available in version 1.2.1 or higher and the stopping rule for declaring a change early is available in version 1.5.1 or higher.

Bayesian Hidden Markov Modeling of Array CGH Data

Genomic alterations have been linked to the development and progression of cancer. The technique of Comparative Genomic Hybridization (CGH) yields data consisting of fluorescence intensity ratios of test and reference DNA samples. The intensity ratios provide information about the number of copies in DNA. Practical issues such as the contamination of tumor cells in tissue specimens and normalization errors necessitate the use of statistics for learning about the genomic alterations from array-CGH data. As increasing amounts of array CGH data become available, there is a growing need for automated algorithms for characterizing genomic profiles. Specifically, there is a need for algorithms that can identify gains and losses in the number of copies based on statistical considerations, rather than merely detect trends in the data. We adopt a Bayesian approach, relying on the hidden Markov model to account for the inherent dependence in the intensity ratios. Posterior inferences are made about gains and losses in copy number. Localized amplifications (associated with oncogene mutations) and deletions (associated with mutations of tumor suppressors) are identified using posterior probabilities. Global trends such as extended regions of altered copy number are detected. Since the posterior distribution is analytically intractable, we implement a Metropolis-within-Gibbs algorithm for efficient simulation-based inference. Publicly available data on pancreatic adenocarcinoma, glioblastoma multiforme and breast cancer are analyzed, and comparisons are made with some widely-used algorithms to illustrate the reliability and success of the technique.

A Pseudolikelihood Approach for Simultaneous Analysis of Array Comparative Genomic Hybridizations (aCGH)

DNA sequence copy number has been shown to be associated with cancer development and progression. Array-based Comparative Genomic Hybridization (aCGH) is a recent development that seeks to identify the copy number ratio at large numbers of markers across the genome. Due to experimental and biological variations across chromosomes and across hybridizations, current methods are limited to analyses of single chromosomes. We propose a more powerful approach that borrows strength across chromosomes and across hybridizations. We assume a Gaussian mixture model, with a hidden Markov dependence structure, and with random effects to allow for intertumoral variation, as well as intratumoral clonal variation. For ease of computation, we base estimation on a pseudolikelihood function. The method produces quantitative assessments of the likelihood of genetic alterations at each clone, along with a graphical display for simple visual interpretation. We assess the characteristics of the method through simulation studies and through analysis of a brain tumor aCGH data set. We show that the pseudolikelihood approach is superior to existing methods both in detecting small regions of copy number alteration and in accurately classifying regions of change when intratumoral clonal variation is present.

Gene duplication: a drive for phenotypic diversity and cause of human disease

Gene duplication is one of the key factors driving genetic innovation, i.e., producing novel genetic variants. Although the contribution of whole-genome and segmental duplications to phenotypic diversity across species is widely appreciated, the phenotypic spectrum and potential pathogenicity of small-scale duplications in individual genomes are less well explored. This review discusses the nature of small-scale duplications and the phenotypes produced by such duplications. Phenotypic variation and disease phenotypes induced by duplications are more diverse and widespread than previously anticipated, and duplications are a major class of disease-related genomic variation. Pathogenic duplications particularly involve dosage-sensitive genes with both similar and dissimilar over- and underexpression phenotypes, and genes encoding proteins with a propensity to aggregate. Phenotypes related to human-specific copy number variation in genes regulating environmental responses and immunity are increasingly recognized. Small genomic duplications containing defense-related genes also contribute to complex common phenotypes.

Large genomic fibrillin-1 (FBN1) gene deletions provide evidence for true haploinsufficiency in Marfan syndrome

Mutations in the FBN1 gene are the major cause of Marfan syndrome (MFS), an autosomal dominant connective tissue disorder, which displays variable manifestations in the cardiovascular, ocular, and skeletal systems. Current molecular genetic testing of FBN1 may miss mutations in the promoter region or in other noncoding sequences as well as partial or complete gene deletions and duplications. In this study, we tested for copy number variations by successively applying multiplex ligation-dependent probe amplification (MLPA) and the Affymetrix Human Mapping 500 K Array Set, which contains probes for approximately 500,000 single-nucleotide polymorphisms (SNPs) across the genome. By analyzing genomic DNA of 101 unrelated individuals with MFS or related phenotypes in whom standard genetic testing detected no mutation, we identified FBN1 deletions in two patients with MFS. Our high-resolution approach narrowed down the deletion breakpoints. Subsequent sequencing of the junctional fragments revealed the deletion sizes of 26,887 and 302,580 bp, respectively. Surprisingly, both deletions affect the putative regulatory and promoter region of the FBN1 gene, strongly indicating that they abolish transcription of the deleted allele. This expectation of complete loss of function of one allele, i.e. true haploinsufficiency, was confirmed by transcript analyses. Our findings not only emphasize the importance of screening for large genomic rearrangements in comprehensive genetic testing of FBN1 but, importantly, also extend the molecular etiology of MFS by providing hitherto unreported evidence that true haploinsufficiency is sufficient to cause MFS.

Fatal adenovirus hepatitis during maintenance therapy for childhood acute lymphoblastic leukemia

Disseminated adenoviral infection with hepatitis is rare in children undergoing standard chemotherapy. We report on a 3(1/2)-year-old male with fatal adenovirus hepatitis receiving maintenance chemotherapy for acute lymphoblastic leukemia (ALL). Adenoviral hepatitis was proven by histology, viral culture, and PCR in a liver biopsy. Quantitative real-time PCR in the peripheral blood showed adenoviral DNA copy number >10(9)/ml. Despite aggressive supportive care and antiviral treatment with cidofovir, the patient died rapidly due to fulminant liver failure. Diagnostic and treatment options for adenovirus infection remain unsatisfactory for these patients. We propose suggestions for diagnosis and therapy.

Enhancement of heterologous expression of alkaline phytase in Pichia pastors

Phytic acid is the major storage form of phosphorus and inositol in seeds and legumes. It forms insoluble phytate salts by chelating with positively charged mineral ions. Non-ruminant animals are not able to digest phytate due to the lack of phytases in their GI tracks, thus the undigested phytate is excreted leading to environmental contamination. Supplementation with phytases in animal feed has proven to be an effective strategy to alleviate nutritional and environmental issues. The unique catalytic and thermal stability properties of alkaline phytase from lily pollen (LlALP) suggest that it has the potential to be useful as a feed supplement. Our goal is to develop a method for the production of substantial amounts of rLlALP for animal feed and structural studies. rLlALP2 has been successfully expressed in the yeast, Pichia pastoris. However, expression yield was modest (8-10 mg/L). Gene copy number has been identified as an important parameter in enhancing protein yields. Multicopy clones were selected using Zeocin-resistance-based vectors and challenging transformants to high Zeocin levels under different conditions. Data indicate that increasing selection pressure led to the generation of clones with amplification of both rLlAlp2 and Zeor genes and the two genes were not equally amplified. Additionally, clones generated by step-wise methods led to clones with greater amplification. The effects of transgene copy number and gene sequence optimization on expression levels of rLlALP2 were examined. The data indicate that increasing the copy number of rLlAlp2 in transformed clones was detrimental to expression level. The use of a sequence-optimized rLlAlp2 (op-rLlAlp2) increased expression yield of the active enzyme by 25-50%, suggesting that transcription and translation efficiency are not major bottlenecks in the production of rLlALP2. Lowering induction temperature to 20 oC led to an increase in enzyme activity of 1.2 to 20-fold, suggesting that protein folding or post-translational processes may be limiting factors for rLlALP2 production. Cumulatively, optimization of copy number, gene sequence optimization and reduced temperature led to increase of rLlALP2 enzyme activity by three-fold (25-30 mg/L). In an effort to simplify the purification process of rLlALP2, extracellular expression of phytase was investigated. Extracellular expression is dependent on the presence of an appropriate secretion signal upstream of the transgene native signal peptide(s) present in the transgene may also influence secretion efficiency. The data suggest that deletion of both N- and C-terminal signal peptides of rLlALP2 enhanced α-mating factor (α-MF)-driven secretion of LlALP2 by four-fold. The secretion signal peptide of chicken egg white lysozyme was ineffective in secretion rLlALP2 in P. pastoris. To enhance rLlALP2 secretion, effectiveness of the strong inducible promoter (PAOX1) was compared with the constitutive promoter (PGAP). The intracellular yield of rLlALP2 was about four-fold greater under the control of PGAP compared to PAOX1 and extracellular expression level of rLlALP2 was around eight-fold (75-100 mg/L) greater. The successful production of active rLlALP2 in P. pastoris will allow us to conduct the animal feed supplementation studies and structural studies.

PHYLOGENETIC ANALYSES, GENE DOSAGE EFFECT AND SELECTION IN MEMBERS OF THE LATERAL BOUNDARY DOMAIN (LBD) GENES TRANSCRIPTION FACTOR FAMILY

Wood formation is an economically and environmentally important process and has played a significant role in the evolution of terrestrial plants. Despite its significance, the molecular underpinnings of the process are still poorly understood. We have previously shown that four Lateral Boundary Domain (LBD) transcription factors have important roles in the regulation of wood formation with two (LBD1 and LBD4) involved in secondary phloem and ray cell development and two (LBD15 and LBD18) in secondary xylem formation. Here, we used comparative phylogenetic analyses to test potential roles of the four LBD genes in the evolution of woodiness. We studied the copy number and variation in DNA and amino acid sequences of the four LBDs in a wide range of woody and herbaceous plant taxa with fully sequenced and annotated genomes. LBD1 showed the highest gene copy number across the studied species, and LBD1 gene copy number was strongly and significantly correlated with the level of ray seriation. The lianas, cucumber and grape, with multiseriate ray cells showed the highest gene copy number (12 and 11, respectively). Because lianas’ growth habit requires significant twisting and bending, the less lignified ray parenchyma cells likely facilitate stem flexibility and maintenance of xylem conductivity. We further demonstrate conservation of amino acids in the LBD18 protein sequences that are specific to woody taxa. Neutrality tests showed evidence for strong purifying selection on these gene regions across various orders, indicating adaptive convergent evolution of LBD18. Structural modeling demonstrates that the conserved amino acids have a significant impact on the tertiary protein structure and thus are likely of significant functional importance.

Regulation of phosphoinositide 3-kinase expression in health and disease

Both the biology and the therapeutic potential of the phosphoinositide 3-kinase (PI3K) signalling axis have been the subject of intense investigation; however, little is known about the regulation of PI3K expression. Emerging evidence indicates that PI3K levels change in response to cellular stimulation with insulin and nuclear receptor ligands, and during various physiological and pathological processes including differentiation, regeneration, hypertension and cancer. Recently identified mechanisms that control PI3K production include increased gene copy number in cancer, and transcriptional regulation of the p110alpha PI3K gene by FOXO3a, NF-kappaB and p53, and of the PI3K regulatory subunits by STAT3, EBNA-2 and SREBP. In most instances, however, the impact of alterations in PI3K expression on PI3K signalling and disease remains to be established.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Clustering subspecies of Aeromonas salmonicida using IS630 typing.

BACKGROUND The insertion element IS630 found in Aeromonas salmonicida belongs to the IS630-Tc1-mariner superfamily of transposons. It is present in multiple copies and represents approximately half of the IS present in the genome of A. salmonicida subsp. salmonicida A449. RESULTS By using High Copy Number IS630 Restriction Fragment Length Polymorphism (HCN-IS630-RFLP), strains of various subspecies of Aeromonas salmonicida showed conserved or clustering patterns, thus allowing their differentiation from each other. Fingerprints of A. salmonicida subsp. salmonicida showed the highest homogeneity while 'atypical' A. salmonicida strains were more heterogeneous. IS630 typing also differentiated A. salmonicida from other Aeromonas species. The copy number of IS630 in Aeromonas salmonicida ranges from 8 to 35 and is much lower in other Aeromonas species. CONCLUSIONS HCN-IS630-RFLP is a powerful tool for subtyping of A. salmonicida. The high stability of IS630 insertions in A. salmonicida subsp. salmonicida indicates that it might have played a role in pathoadaptation of A. salmonicida which has reached an optimal configuration in the highly virulent and specific fish pathogen A. salmonicida subsp. salmonicida.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

EGFR exon-level biomarkers of the response to bevacizumab/erlotinib in non-small cell lung cancer

Activating epidermal growth factor receptor (EGFR) mutations are recognized biomarkers for patients with metastatic non-small cell lung cancer (NSCLC) treated with EGFR tyrosine kinase inhibitors (TKIs). EGFR TKIs can also have activity against NSCLC without EGFR mutations, requiring the identification of additional relevant biomarkers. Previous studies on tumor EGFR protein levels and EGFR gene copy number revealed inconsistent results. The aim of the study was to identify novel biomarkers of the response to TKIs in NSCLC by investigating whole genome expression at the exon-level. We used exon arrays and clinical samples from a previous trial (SAKK19/05) to investigate the expression variations at the exon-level of 3 genes potentially playing a key role in modulating treatment response: EGFR, V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS) and vascular endothelial growth factor (VEGFA). We identified the expression of EGFR exon 18 as a new predictive marker for patients with untreated metastatic NSCLC treated with bevacizumab and erlotinib in the first line setting. The overexpression of EGFR exon 18 in tumor was significantly associated with tumor shrinkage, independently of EGFR mutation status. A similar significant association could be found in blood samples. In conclusion, exonic EGFR expression particularly in exon 18 was found to be a relevant predictive biomarker for response to bevacizumab and erlotinib. Based on these results, we propose a new model of EGFR testing in tumor and blood.

RNA-dependent genome processing during nuclear differentiation: the model systems of stichotrichous ciliates

We introduce ciliated protozoa, and more specifically the stichotrichous ciliates Oxytricha and Stylonychia, as biological model systems for the analysis of programmed DNA-reorganization processes during nuclear differentiation. These include DNA excision, DNA elimination, reordering of gene segments and specific gene amplification. We show that small nuclear RNAs specify DNA sequences to be excised or retained, but also discuss the need for a RNA template molecule derived from the parental nucleus for these processes. This RNA template guides reordering of gene segments to become functional genes and determines gene copy number in the differentiated nucleus. Since the template is derived from the parental macronucleus, gene reordering and DNA amplification are inherited in a non-Mendelian epigenetic manner.

Estrogen receptor β expression and androgen receptor phosphorylation correlate with a poor clinical outcome in hormone-naive prostate cancer and are elevated in castration-resistant disease

Patients with advanced prostate cancer (PC) are usually treated with androgen withdrawal. While this therapy is initially effective, nearly all PCs become refractory to it. As hormone receptors play a crucial role in this process, we constructed a tissue microarray consisting of PC samples from 107 hormone-naïve (HN) and 101 castration-resistant (CR) PC patients and analyzed the androgen receptor (AR) gene copy number and the protein expression profiles of AR, Serin210-phosphorylated AR (pAR(210)), estrogen receptor (ER)β, ERα and the proliferation marker Ki67. The amplification of the AR gene was virtually restricted to CR PC and was significantly associated with increased AR protein expression (P<0.0001) and higher tumor cell proliferation (P=0.001). Strong AR expression was observed in a subgroup of HN PC patients with an adverse prognosis. In contrast, the absence of AR expression in CR PC was significantly associated with a poor overall survival. While pAR(210) was predominantly found in CR PC patients (P<0.0001), pAR(210) positivity was observed in a subgroup of HN PC patients with a poor survival (P<0.05). Epithelial ERα expression was restricted to CR PC cells (9%). ERβ protein expression was found in 38% of both HN and CR PCs, but was elevated in matched CR PC specimens. Similar to pAR(210), the presence of ERβ in HN patients was significantly associated with an adverse prognosis (P<0.005). Our results strongly suggest a major role for pAR(210) and ERβ in HN PC. The expression of these markers might be directly involved in CR tumor growth.