Estimating the net contribution of interleukin-28B variation to spontaneous hepatitis C virus clearance

The identification of associations between interleukin-28B (IL-28B) variants and the spontaneous clearance of hepatitis C virus (HCV) raises the issues of causality and the net contribution of host genetics to the trait. To estimate more precisely the net effect of IL-28B genetic variation on HCV clearance, we optimized genotyping and compared the host contributions in multiple- and single-source cohorts to control for viral and demographic effects. The analysis included individuals with chronic or spontaneously cleared HCV infections from a multiple-source cohort (n = 389) and a single-source cohort (n = 71). We performed detailed genotyping in the coding region of IL-28B and searched for copy number variations to identify the genetic variant or haplotype carrying the strongest association with viral clearance. This analysis was used to compare the effects of IL-28B variation in the two cohorts. Haplotypes characterized by carriage of the major alleles at IL-28B single-nucleotide polymorphisms (SNPs) were highly overrepresented in individuals with spontaneous clearance versus those with chronic HCV infections (66.1% versus 38.6%, P = 6 × 10(-9) ). The odds ratios for clearance were 2.1 [95% confidence interval (CI) = 1.6-3.0] and 3.9 (95% CI = 1.5-10.2) in the multiple- and single-source cohorts, respectively. Protective haplotypes were in perfect linkage (r(2) = 1.0) with a nonsynonymous coding variant (rs8103142). Copy number variants were not detected. CONCLUSION: We identified IL-28B haplotypes highly predictive of spontaneous HCV clearance. The high linkage disequilibrium between IL-28B SNPs indicates that association studies need to be complemented by functional experiments to identify single causal variants. The point estimate for the genetic effect was higher in the single-source cohort, which was used to effectively control for viral diversity, sex, and coinfections and, therefore, offered a precise estimate of the net host genetic contribution.

qPCR-based mitochondrial DNA quantification: Influence of template DNA fragmentation on accuracy

Real-time PCR (qPCR) is the method of choice for quantification of mitochondrial DNA (mtDNA) by relative comparison of a nuclear to a mitochondrial locus. Quantitative abnormal mtDNA content is indicative of mitochondrial disorders and mostly confines in a tissue-specific manner. Thus handling of degradation-prone bioptic material is inevitable. We established a serial qPCR assay based on increasing amplicon size to measure degradation status of any DNA sample. Using this approach we can exclude erroneous mtDNA quantification due to degraded samples (e.g. long post-exicision time, autolytic processus, freeze-thaw cycles) and ensure abnormal DNA content measurements (e.g. depletion) in non-degraded patient material. By preparation of degraded DNA under controlled conditions using sonification and DNaseI digestion we show that erroneous quantification is due to the different preservation qualities of the nuclear and the mitochondrial genome. This disparate degradation of the two genomes results in over- or underestimation of mtDNA copy number in degraded samples. Moreover, as analysis of defined archival tissue would allow to precise the molecular pathomechanism of mitochondrial disorders presenting with abnormal mtDNA content, we compared fresh frozen (FF) with formalin-fixed paraffin-embedded (FFPE) skeletal muscle tissue of the same sample. By extrapolation of measured decay constants for nuclear DNA (λnDNA) and mtDNA (λmtDNA) we present an approach to possibly correct measurements in degraded samples in the future. To our knowledge this is the first time different degradation impact of the two genomes is demonstrated and which evaluates systematically the impact of DNA degradation on quantification of mtDNA copy number.

Role of KCNMA1 in breast cancer

KCNMA1 encodes the α-subunit of the large conductance, voltage and Ca(2+)-activated (BK) potassium channel and has been reported as a target gene of genomic amplification at 10q22 in prostate cancer. To investigate the prevalence of the amplification in other human cancers, the copy number of KCNMA1 was analyzed by fluorescence-in-situ-hybridization (FISH) in 2,445 tumors across 118 different tumor types. Amplification of KCNMA1 was restricted to a small but distinct fraction of breast, ovarian and endometrial cancer with the highest prevalence in invasive ductal breast cancers and serous carcinoma of ovary and endometrium (3-7%). We performed an extensive analysis on breast cancer tissue microarrays (TMA) of 1,200 tumors linked to prognosis. KCNMA1 amplification was significantly associated with high tumor stage, high grade, high tumor cell proliferation, and poor prognosis. Immunofluorescence revealed moderate or strong KCNMA1 protein expression in 8 out of 9 human breast cancers and in the breast cancer cell line MFM223. KCNMA1-function in breast cancer cell lines was confirmed by whole-cell patch clamp recordings and proliferation assays, using siRNA-knockdown, BK channel activators such as 17ß-estradiol and the BK-channel blocker paxilline. Our findings revealed that enhanced expression of KCNMA1 correlates with and contributes to high proliferation rate and malignancy of breast cancer.

BCL-2 family member BOK is widely expressed but its loss has only minimal impact in mice

BOK/MTD was discovered as a protein that binds to the anti-apoptotic Bcl-2 family member MCL-1 and shares extensive amino-acid sequence similarity to BAX and BAK, which are essential for the effector phase of apoptosis. Therefore, and on the basis of its reported expression pattern, BOK is thought to function in a BAX/BAK-like pro-apoptotic manner in female reproductive tissues. In order to determine the function of BOK, we examined its expression in diverse tissues and investigated the consequences of its loss in Bok(-/-) mice. We confirmed that Bok mRNA is prominently expressed in the ovaries and uterus, but also observed that it is present at readily detectable levels in several other tissues such as the brain and myeloid cells. Bok(-/-) mice were produced at the expected Mendelian ratio, appeared outwardly normal and proved fertile. Histological examination revealed that major organs in Bok(-/-) mice displayed no morphological aberrations. Although several human cancers have somatically acquired copy number loss of the Bok gene and BOK is expressed in B lymphoid cells, we found that its deficiency did not accelerate lymphoma development in Eμ-Myc transgenic mice. Collectively, these results indicate that Bok may have a role that largely overlaps with that of other members of the Bcl-2 family, or may have a function restricted to specific stress stimuli and/or tissues.

SNP array profiling of childhood adrenocortical tumors reveals distinct pathways of tumorigenesis and highlights candidate driver genes

Childhood adrenocortical tumors (ACT) are rare malignancies, except in southern Brazil, where a higher incidence rate is associated to a high frequency of the founder R337H TP53 mutation. To date, copy number alterations in these tumors have only been analyzed by low-resolution comparative genomic hybridization.

Detection of exon deletions within an entire gene (CFTR) by relative quantification on the LightCycler

BACKGROUND: Cystic fibrosis (CF) is associated with at least 1 pathogen point sequence variant on each CFTR allele. Some symptomatic patients, however, have only 1 detectable pathogen sequence variant and carry, on the other allele, a large deletion that is not detected by conventional screening methods. METHODS: For relative quantitative real-time PCR detection of large deletions in the CFTR gene, we designed DNA-specific primers for each exon of the gene and primers for a reference gene (beta2-microglobulin). For PCR we used a LightCycler system (Roche) and calculated the gene-dosage ratio of CFTR to beta2-microglobulin. We tested the method by screening all 27 exons in 3 healthy individuals and 2 patients with only 1 pathogen sequence variant. We then performed specific deletion screenings in 10 CF patients with known large deletions and a blinded analysis in which we screened 24 individuals for large deletions by testing 8 of 27 exons. RESULTS: None of the ratios for control samples were false positive (for deletions or duplications); moreover, for all samples from patients with known large deletions, the calculated ratios for deleted exons were close to 0.5. In addition, the results from the blinded analysis demonstrated that our method can also be used for the screening of single individuals. CONCLUSIONS: The LightCycler assay allows reliable and rapid screening for large deletions in the CFTR gene and detects the copy number of all 27 exons.

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.

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.