Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Evidence for persistent dysfunction of wild-type aldosterone synthase gene in glucocorticoid-treated familial hyperaldosteronism type I

Background In familial hyperaldosteronism type I (FH-I), glucocorticoid treatment suppresses adrenocorticotrophic hormone-regulated hybrid gene expression and corrects hyperaldosteronism. Objective To determine whether the wild-type aldosterone synthase genes, thereby released from chronic suppression, are capable of functioning normally. Methods We compared mid-morning levels of plasma potassium, plasma aldosterone, plasma renin activity (PRA) and aldosterone : PRA ratios, measured with patients in an upright position, and responsiveness of aldosterone levels to infusion of angiotensin II (AII), for 11 patients with FH-I before and during long-term (0.8-14.3 years) treatment with 0.25-0.75 mg/day dexamethasone or 2.5-10 mg/day prednisolone. Results During glucocorticoid treatment, hypertension was corrected in all. Potassium levels, which had been low (< 3.5 mmol/l) in two patients before treatment, were normal in all during treatment (mean 4.0 +/- 0.1 mmol/l, range 3.5-4.6). Aldosterone levels during treatment [13.2 +/- 2.1 ng/100 ml (mean +/- SEM)] were lower than those before treatment (20.1 +/- 2.5 ng/100 ml, P < 0.05). PRA levels, which had been suppressed before treatment (0.5 +/- 0.2 ng/ml per h), were unsuppressed during treatment (5.1 +/- 1.5 ng/ml per h, P < 0.01) and elevated (> 4 ng/ml per h) in six patients. Aldosterone : PRA ratios, which had been elevated (> 30) before treatment (101.1 +/- 25.9), were much lower during treatment (4.1 +/- 1.0, P < 0.005) and below normal (< 5) in eight patients. Surprisingly, aldosterone level, which had not been responsive (< 50% rise) to infusion of AII for all 11 patients before treatment, remained unresponsive for 10 during treatment. Conclusions Apparently regardless of duration of glucocorticoid treatment in FH-I, aldosterone level remains poorly responsive to AII, with a higher than normal PRA and a low aldosterone : PRA ratio. This is consistent with there being a persistent defect in functioning of wild-type aldosterone synthase gene. (C) Rapid Science Publishers ISSN 0263-6352.

Identification of two novel mutations in the hydroxymethylbilane synthase gene in three patients from two unrelated families with acute intermittent porphyria

We have screened the hydroxymethylbilane synthase cDNAs of 3 patients from 2 families suffering from acute intermittent porphyria (AIP) from Scotland and South Africa using heteroduplex and chemical cleavage of mismatch analyses, Direct sequencing was used to characterise the mutations, The two novel mutations identified were a missense mutation at nucleotide position 64 in exon 3 (R22C) and a single base-pair deletion in exon 15, These mutations are predicted to affect the normal function of the enzyme and, therefore, are expected to be the primary cause of disease in these patients.

Specific PCR based detection of Phytophthora medicaginis using the intergenic spacer region of the ribosomal DNA

A technique based on the polymerase chain reaction (PCR) for the specific detection of Phytophthora medicaginis was developed using nucleotide sequence information of the ribosomal DNA (rDNA) regions. The complete IGS 2 region between the 5 S gene of one rDNA repeat and the small subunit of the adjacent repeat was sequenced for P. medicaginis and related species. The entire nucleotide sequence length of the IGS 2 of P. medicaginis was 3566 bp. A pair of oligonucleotide primers (PPED04 and PPED05), which allowed amplification of a specific fragment (364 bp) within the IGS 2 of P. medicaginis using the PCR, was designed. Specific amplification of this fragment from P. medicaginis was highly sensitive, detecting template DNA as low as 4 ng and in a host-pathogen DNA ratio of 1000000:1. Specific PCR amplification using PPED04 and PPED05 was successful in detecting P. medicaginis in lucerne stems infected under glasshouse conditions and field infected lucerne roots. The procedures developed in this work have application to improved identification and detection of a wide range of Phytophthora spp. in plants and soil.

Nucleotide sequence of the nifH gene coding for nitrogen reductase in the acetic acid bacterium Acetobacter diazotrophicus

The nifH gene sequence of the nitrogen-fixing bacterium Acetobacter diazotrophicus was determined with the use of the polymerase chain reaction and universal degenerate oligonucleotide primers. The gene shows highest pair-wise similarity to the nifH gene of Azospirillum brasilense. The phylogenetic relationships of the nifH gene sequences were compared with those inferred from 16S rRNA gene sequences. Knowledge of the sequence of the nifH gene contributes to the growing database of nifH gene sequences, and will allow the detection of Acet. diazotrophicus from environmental samples with nifH gene-based primers.

Concordance between dispersal and mitochondrial gene flow: isolation by distance in a tropical teleost, Lates calcarifer (Australian barramundi)

Patterns of population subdivision and the relationship between gene flow and geographical distance in the tropical estuarine fish Lares calcarifer (Centropomidae) were investigated using mtDNA control region sequences. Sixty-three putative haplotypes were resolved from a total of 270 individuals from nine localities within three geographical regions spanning the north Australian coastline. Despite a continuous estuarine distribution throughout the sampled range, no haplotypes were shared among regions. However, within regions, common haplotypes were often shared among localities. Both sequence-based (average Phi(ST)=0.328) and haplotype-based (average Phi(ST)=0.182) population subdivision analyses indicated strong geographical structuring. Depending on the method of calculation, geographical distance explained either 79 per cent (sequence-based) or 23 per cent (haplotype-based) of the variation in mitochondrial gene flow. Such relationships suggest that genetic differentiation of L. calcarifer has been generated via isolation-by-distance, possibly in a stepping-stone fashion. This pattern of genetic structure is concordant with expectations based on the life history of L. calcarifer and direct studies of its dispersal patterns. Mitochondrial DNA variation, although generally in agreement with patterns of allozyme variation, detected population subdivision at smaller spatial scales. Our analysis of mtDNA variation in L. calcarifer confirms that population genetic models can detect population structure of not only evolutionary significance but also of demographic significance. Further, it demonstrates the power of inferring such structure from hypervariable markers, which correspond to small effective population sizes.

Functional polymorphism of the human arylamine N-acetyltransferase type 1 gene caused by (CT)-T-190 and G(560)A mutations

Human N-acetyltransferase type 1 (NAT1) catalyses the N- or O-acetylation of various arylamine and heterocyclic amine substrates and is able to bioactivate several known carcinogens. Despite wide inter-individual variability in activity, historically, NAT1 was considered to be monomorphic in nature. However, recent reports of allelic variation at the NAT1 locus suggest that it may be a polymorphically expressed enzyme. In the present study, peripheral blood mononuclear cell NAT1 activity in 85 individuals was found to be bimodally distributed with approximately 8% of the population being slow acetylators. Subsequent sequencing of the individuals having slow acetylator status showed all to have either a (CT)-T-190 or G(560)A base substitution located in the protein encoding region of the NAT1 gene. The (CT)-T-190 base substitution changed a highly conserved Arg(64), which others have shown to be essential for fully functional NAT1 protein. The (CT)-T-190 mutation has not been reported previously and we have named it NAT1*17. The G(560)A mutation is associated with the base substitutions previously observed in the NAT1*10 allele and this variant (NAT1*14) encodes for a protein with reduced acetylation capacity. A novel method using linear PCR and dideoxy terminators was developed for the detection of NAT1*14 and NAT1*17. Neither of these variants was found in the rapid acetylator population. We conclude that both the (CT)-T-190 (NAT1*17) and G(560)A (NAT1*14) NAT1 structural variants are involved in a distinct NAT1 polymorphism. Because NAT1 can bioactivate several carcinogens, this polymorphism may have implications for cancer risk in individual subjects. (C) 1998 Chapman & Hall Ltd.

Molecular cloning and sequence of the ovine gastrin gene

A clone encoding ovine preprogastrin was isolated from a sheep genomic library. The deduced 104 amino acid sequence of ovine preprogastrin was 92% and 68% identical to the sequences of bovine and human preprogastrin, respectively. While the similarity was greatest in the gastrin-17 sequence, an unexpected similarity was also observed in the N-terminus of mature progastrin.

Mutational analysis of the cystathionine beta-synthase gene: A splicing mutation, two missense mutations and an insertion in patients with homocystinuria

RT-PCR and direct sequence analyses were used to define mutations in the cystathionine beta-synthase (CBS) gene in two unrelated male patients with vitamin B6 nonresponsive homocystinuria. Both patients were compound heterozygotes for CBS alleles containing point mutations. One patient had a maternally derived G-->A transition in the splice-donor site of intron 1, resulting in aberrant splicing of CBS mRNA. The other allele contained a missense mutation resulting in the previously reported E144K mutant CBS protein. The second patient had a maternally derived 4 bp insertion in exon 17, predicted to cause a CBS peptide of altered amino acid sequence. A 494G-->A transition was found in exon 4 of the other allele, predicting a C165Y substitution. Expression of recombinant CBS protein, containing the C165Y mutation, had no detectable catalytic activity. Each mutation was confirmed in genomic DNA. (C) 1998 Wiley-Liss, Inc.

Microsatellite instability and loss of heterozygosity at DNA mismatch repair gene loci occurs during hepatic carcinogenesis

DNA mismatch repair is an important mechanism involved in maintaining the fidelity of genomic DNA. Defective DNA mismatch repair is implicated in a variety of gastrointestinal and other turners; however, its role in hepatocellular carcinoma (HCC) has not been assessed. Formalin-fixed, paraffin-embedded archival pathology tissues from 46 primary liver tumors were studied by microdissection and microsatellite analysis of extracted DNA to assess the degree of microsatellite instability, a marker of defective mismatch repair, and to determine the extent and timing of allelic loss of two DNA mismatch repair genes, human Mut S homologue-2 (hMSH2) and human Mut L homologue-1 (hMLH1), and the tumor suppressor genes adenomatous polyposis coli gene (APC), p53, and DPC4. Microsatellite instability was detected in 16 of the tumors (34.8%). Loss of heterozygosity at microsatellites linked to the DNA mismatch repair genes, hMSH2 and/or hMLH1, was found in 9 cases (19.6%), usually in association with microsatellite instability. Importantly, the pattern of allelic loss was uniform in 8 of these 9 tumors, suggesting that clonal loss had occurred. Moreover, loss at these loci also occurred in nonmalignant tissue adjacent to 4 of these tumors, where it was associated with marked allelic heterogeneity. There was relatively infrequent loss of APC, p53, or DPC4 loci that appeared unrelated to loss of hMSH2 or hMLH1 gene loci. Loss of heterozygosity at hMSH2 and/or hMLH1 gene loci, and the associated microsatellite instability in premalignant hepatic tissues suggests a possible causal role in hepatic carcinogenesis in a subset of hepatomas.

From sabotage to camouflage: viral evasion of cytotoxic T lymphocyte and natural killer cell-mediated immunity

The outcome of a virus infection is strongly influenced by interactions between host immune defences and virus 'anti-defence' mechanisms. For many viruses, their continued survival depends on, the speed of their attach: their capacity to replicate and transmit to uninfected hosts prior to their elimination by an effective immune response. In contrast, the success of persistent viruses lies in their capacity for immunological subterfuge: the evasion of host defence mechanisms by either mutation (covered elsewhere in this issue, by Gould and Bangham, pp. 321-328) or interference with the action of host cellular proteins that are important components of the immune response. This review will focus on the strategies employed by persistent viruses against two formidable host defences against virus infection: the CD8+ cytotoxic T lymphocyte (CTL) and natural killer (NK) cell responses.

Febrile seizures and generalized epilepsy associated with a mutation in the Na+-channel beta 1 subunit gene SCN1B

Febrile seizures affect approximately 3% of all children under six years of age and are by far the most common seizure disorder(1). A small proportion of children with febrile seizures later develop ongoing epilepsy with afebrile seizures(2). Segregation analysis suggests the majority of cases have complex inheritance(3) but rare families show apparent autosomal dominant: inheritance. Two putative loci have been mapped (FEB1 and FEB2), but specific genes have not yet been identified(4,5). We recently described a clinical subset, termed generalized epilepsy with febrile seizures plus (GEFS(+)), in which many family members have seizures with fever that may persist beyond six years of age or be associated with afebrile generalized seizures(6). We now report linkage, in another large GEFS(+) family, to chromosome region 19q13.1 and identification of a mutation in the voltage-gated sodium (Na+)-channel beta 1 subunit gene (SCN1B). The mutation changes a conserved cysteine residue disrupting a putative disulfide bridge which normally maintains an extracellular immunoglobulin-like fold. Go-expression of the mutant pr subunit with a brain Na+-channel alpha subunit in Xenopus laevis oocytes demonstrates that the mutation interferes with the ability of the subunit to modulate channel-gating kinetics consistent with a loss-of-function allele. This observation develops the theme that idiopathic epilepsies are a family of channelopathies and raises the possibility of involvement of other Na+-channel subunit genes in febrile seizures and generalized epilepsies with complex inheritance patterns.