Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

Effect of single nucleotide polymorphisms in cytochrome P450 isoenzyme and N-acetyltransferase 2 genes on the metabolism of artemisinin-based combination therapies in malaria patients from Cambodia and Tanzania.

The pharmacogenetics of antimalarial agents are poorly known, although the application of pharmacogenetics might be critical in optimizing treatment. This population pharmacokinetic-pharmacogenetic study aimed at assessing the effects of single nucleotide polymorphisms (SNPs) in cytochrome P450 isoenzyme genes (CYP, namely, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, and CYP3A5) and the N-acetyltransferase 2 gene (NAT2) on the pharmacokinetics of artemisinin-based combination therapies in 150 Tanzanian patients treated with artemether-lumefantrine, 64 Cambodian patients treated with artesunate-mefloquine, and 61 Cambodian patients treated with dihydroartemisinin-piperaquine. The frequency of SNPs varied with the enzyme and the population. Higher frequencies of mutant alleles were found in Cambodians than Tanzanians for CYP2C9*3, CYP2D6*10 (100C → T), CYP3A5*3, NAT2*6, and NAT2*7. In contrast, higher frequencies of mutant alleles were found in Tanzanians for CYP2D6*17 (1023C → T and 2850C → T), CYP3A4*1B, NAT2*5, and NAT2*14. For 8 SNPs, no significant differences in frequencies were observed. In the genetic-based population pharmacokinetic analyses, none of the SNPs improved model fit. This suggests that pharmacogenetic data need not be included in appropriate first-line treatments with the current artemisinin derivatives and quinolines for uncomplicated malaria in specific populations. However, it cannot be ruled out that our results represent isolated findings, and therefore more studies in different populations, ideally with the same artemisinin-based combination therapies, are needed to evaluate the influence of pharmacogenetic factors on the clearance of antimalarials.

Rapid evolution of cancer/testis genes on the X chromosome.

BACKGROUND: Cancer/testis (CT) genes are normally expressed only in germ cells, but can be activated in the cancer state. This unusual property, together with the finding that many CT proteins elicit an antigenic response in cancer patients, has established a role for this class of genes as targets in immunotherapy regimes. Many families of CT genes have been identified in the human genome, but their biological function for the most part remains unclear. While it has been shown that some CT genes are under diversifying selection, this question has not been addressed before for the class as a whole. RESULTS: To shed more light on this interesting group of genes, we exploited the generation of a draft chimpanzee (Pan troglodytes) genomic sequence to examine CT genes in an organism that is closely related to human, and generated a high-quality, manually curated set of human:chimpanzee CT gene alignments. We find that the chimpanzee genome contains homologues to most of the human CT families, and that the genes are located on the same chromosome and at a similar copy number to those in human. Comparison of putative human:chimpanzee orthologues indicates that CT genes located on chromosome X are diverging faster and are undergoing stronger diversifying selection than those on the autosomes or than a set of control genes on either chromosome X or autosomes. CONCLUSION: Given their high level of diversifying selection, we suggest that CT genes are primarily responsible for the observed rapid evolution of protein-coding genes on the X chromosome.

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes.

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

Transcriptional response to cardiac injury in the zebrafish: systematic identification of genes with highly concordant activity across in vivo models.

BACKGROUND: Zebrafish is a clinically-relevant model of heart regeneration. Unlike mammals, it has a remarkable heart repair capacity after injury, and promises novel translational applications. Amputation and cryoinjury models are key research tools for understanding injury response and regeneration in vivo. An understanding of the transcriptional responses following injury is needed to identify key players of heart tissue repair, as well as potential targets for boosting this property in humans. RESULTS: We investigated amputation and cryoinjury in vivo models of heart damage in the zebrafish through unbiased, integrative analyses of independent molecular datasets. To detect genes with potential biological roles, we derived computational prediction models with microarray data from heart amputation experiments. We focused on a top-ranked set of genes highly activated in the early post-injury stage, whose activity was further verified in independent microarray datasets. Next, we performed independent validations of expression responses with qPCR in a cryoinjury model. Across in vivo models, the top candidates showed highly concordant responses at 1 and 3 days post-injury, which highlights the predictive power of our analysis strategies and the possible biological relevance of these genes. Top candidates are significantly involved in cell fate specification and differentiation, and include heart failure markers such as periostin, as well as potential new targets for heart regeneration. For example, ptgis and ca2 were overexpressed, while usp2a, a regulator of the p53 pathway, was down-regulated in our in vivo models. Interestingly, a high activity of ptgis and ca2 has been previously observed in failing hearts from rats and humans. CONCLUSIONS: We identified genes with potential critical roles in the response to cardiac damage in the zebrafish. Their transcriptional activities are reproducible in different in vivo models of cardiac injury.

No interaction between alcohol consumption and HDL-related genes on HDL cholesterol levels.

OBJECTIVE: To assess the relationships and possible interactions between polymorphisms related to HDL levels and alcohol consumption. METHODS: Cross-sectional population-based study including 2863 women and 2546 men aged 35-75 years (CoLaus study). Alcohol intake was assessed by the reported alcohol consumption of the last 7 days. Nineteen candidate genes known to influence HDL levels were studied. RESULTS: Alcohol consumption increased HDL cholesterol levels in both genders. After multivariate adjustment for gender, age, body mass index, smoking, hypolipidaemic drug treatment, physical activity and alcohol consumption, APOA5, CETP, LIPC and LPL gene polymorphisms were significantly (10(-5) threshold) related with HDL cholesterol levels, while no genexalcohol intake interaction was found for all SNPs studied. ABCA1 polymorphisms were related to HDL cholesterol levels on bivariate analysis but the relationship was no longer significant after multivariate analysis. CONCLUSION: Our data confirm the association of alcohol consumption and of APOA5, CETP, LIPC and LPL gene polymorphisms with HDL cholesterol levels. Conversely, no genexalcohol consumption interactions were found, suggesting that the effect of alcohol consumption on HDL cholesterol levels is not mediated via a modulation of HDL related genes.

Male dominance linked to size and age, but not to 'good genes' in brown trout (Salmo trutta).

BACKGROUND: Males that are successful in intra-sexual competition are often assumed to be of superior quality. In the mating system of most salmonid species, intensive dominance fights are common and the winners monopolise most mates and sire most offspring. We drew a random sample of mature male brown trout (Salmo trutta) from two wild populations and determined their dominance hierarchy or traits linked to dominance. The fish were then stripped and their sperm was used for in vitro fertilisations in two full-factorial breeding designs. We recorded embryo viability until hatching in both experiments, and juvenile survival during 20 months after release into a natural streamlet in the second experiment. Since offspring of brown trout get only genes from their fathers, we used offspring survival as a quality measure to test (i) whether males differ in their genetic quality, and if so, (ii) whether dominance or traits linked to dominance reveal 'good genes'. RESULTS: We found significant additive genetic variance on embryo survival, i.e. males differed in their genetic quality. Older, heavier and larger males were more successful in intra-sexual selection. However, neither dominance nor dominance indicators like body length, weight or age were significantly linked to genetic quality measured as embryo or juvenile survival. CONCLUSION: We found no evidence that females can improve their offspring's genetic viability by mating with large and dominant males. If there still were advantages of mating with dominant males, they may be linked to non-genetic benefits or to genetic advantages that are context dependent and therefore possibly not revealed under our experimental conditions - even if we found significant additive genetic variation for embryo viability under such conditions.

Reversal of the silencing of tetracycline-controlled genes requires the coordinate action of distinctly acting transcription factors.

BACKGROUND: Regulation of genes transferred to eukaryotic organisms is often limited by the lack of consistent expression levels in all transduced cells, which may result in part from epigenetic gene silencing effects. This reduces the efficacy of ligand-controlled gene switches designed for somatic gene transfers such as gene therapy. METHODS: A doxycycline-controlled transgene was stably introduced in human cells, and clones were screened for epigenetic silencing of the transgene. Various regulatory proteins were targeted to the silent transgene, to identify those that would mediate regulation by doxycycline. RESULTS: A doxycycline-controlled minimal promoter was found to be prone to gene silencing, which prevents activation by a fusion of the bacterial TetR DNA-binding domain with the VP16 activator. DNA modification studies indicated that the silenced transgene adopts a poorly accessible chromatin structure. Several cellular transcriptional activators were found to restore an accessible DNA structure when targeted to the silent transgene, and they cooperated with Tet-VP16 to mediate regulation by doxycycline. CONCLUSIONS: Reversal of the silencing of a tetracycline-regulated minimal promoter requires a chromatin-remodeling activity for subsequent promoter activation by the Tet-VP16 fusion protein. Thus, distinct regulatory elements may be combined to obtain long-term regulation and persistent expression of exogenous genes in eukaryotic cells.

Efficient Subtractive Cloning of Genes Activated by Lipopolysaccharide and Interferon γ in Primary-Cultured Cortical Cells of Newborn Mice.

Innate immune responses play a central role in neuroprotection and neurotoxicity during inflammatory processes that are triggered by pathogen-associated molecular pattern-exhibiting agents such as bacterial lipopolysaccharide (LPS) and that are modulated by inflammatory cytokines such as interferon γ (IFNγ). Recent findings describing the unexpected complexity of mammalian genomes and transcriptomes have stimulated further identification of novel transcripts involved in specific physiological and pathological processes, such as the neural innate immune response that alters the expression of many genes. We developed a system for efficient subtractive cloning that employs both sense and antisense cRNA drivers, and coupled it with in-house cDNA microarray analysis. This system enabled effective direct cloning of differentially expressed transcripts, from a small amount (0.5 µg) of total RNA. We applied this system to isolation of genes activated by LPS and IFNγ in primary-cultured cortical cells that were derived from newborn mice, to investigate the mechanisms involved in neuroprotection and neurotoxicity in maternal/perinatal infections that cause various brain injuries including periventricular leukomalacia. A number of genes involved in the immune and inflammatory response were identified, showing that neonatal neuronal/glial cells are highly responsive to LPS and IFNγ. Subsequent RNA blot analysis revealed that the identified genes were activated by LPS and IFNγ in a cooperative or distinctive manner, thereby supporting the notion that these bacterial and cellular inflammatory mediators can affect the brain through direct but complicated pathways. We also identified several novel clones of apparently non-coding RNAs that potentially harbor various regulatory functions. Characterization of the presently identified genes will give insights into mechanisms and interventions not only for perinatal infection-induced brain damage, but also for many other innate immunity-related brain disorders.

A role for clock genes in sleep homeostasis.

The timing and quality of both sleep and wakefulness are thought to be regulated by the interaction of two processes. One of these two processes keeps track of the prior sleep-wake history and controls the homeostatic need for sleep while the other sets the time-of-day that sleep preferably occurs. The molecular pathways underlying the latter, circadian process have been studied in detail and their key role in physiological time-keeping has been well established. Analyses of sleep in mice and flies lacking core circadian clock gene proteins have demonstrated, however, that besides disrupting circadian rhythms, also sleep homeostatic processes were affected. Subsequent studies revealed that sleep loss alters both the mRNA levels and the specific DNA-binding of the key circadian transcriptional regulators to their target sequences in the mouse brain. The fact that sleep loss impinges on the very core of the molecular circadian circuitry might explain why both inadequate sleep and disrupted circadian rhythms can similarly lead to metabolic pathology. The evidence for a role for clock genes in sleep homeostasis will be reviewed here.

Predator-prey chemical warfare determines the expression of biocontrol genes by rhizosphere-associated Pseudomonas fluorescens.

Soil bacteria are heavily consumed by protozoan predators, and many bacteria have evolved defense strategies such as the production of toxic exometabolites. However, the production of toxins is energetically costly and therefore is likely to be adjusted according to the predation risk to balance the costs and benefits of predator defense. We investigated the response of the biocontrol bacterium Pseudomonas fluorescens CHA0 to a common predator, the free-living amoeba Acanthamoeba castellanii. We monitored the effect of the exposure to predator cues or direct contact with the predators on the expression of the phlA, prnA, hcnA, and pltA genes, which are involved in the synthesis of the toxins, 2,4-diacetylphloroglucinol (DAPG), pyrrolnitrin, hydrogen cyanide, and pyoluteorin, respectively. Predator chemical cues led to 2.2-, 2.0-, and 1.2-fold increases in prnA, phlA, and hcnA expression, respectively, and to a 25% increase in bacterial toxicity. The upregulation of the tested genes was related to the antiprotozoan toxicity of the corresponding toxins. Pyrrolnitrin and DAPG had the highest toxicity, suggesting that bacteria secrete a predator-specific toxin cocktail. The response of the bacteria was elicited by supernatants of amoeba cultures, indicating that water-soluble chemical compounds were responsible for induction of the bacterial defense response. In contrast, direct contact of bacteria with living amoebae reduced the expression of the four bacterial toxin genes by up to 50%, suggesting that protozoa can repress bacterial toxicity. The results indicate that predator-prey interactions are a determinant of toxin production by rhizosphere P. fluorescens and may have an impact on its biocontrol potential.

Biogeographic origin and radiation of the Old World crocidurine shrews (Mammalia: Soricidae) inferred from mitochondrial and nuclear genes.

The crocidurine shrews include the most speciose genus of mammals, Crocidura. The origin and evolution of their radiation is, however, poorly understood because of very scant fossil records and a rather conservative external morphology between species. Here, we use an alignment of 3560 base pairs of mitochondrial and nuclear DNA to generate a phylogenetic hypothesis for the evolution of Old World shrews of the subfamily Crocidurinae. These molecular data confirm the monophyly of the speciose African and Eurasian Crocidura, which also includes the fossorial, monotypic genus Diplomesodon. The phylogenetic reconstructions give further credit to a paraphyletic position of Suncus shrews, which are placed into at least two independent clades (one in Africa and sister to Sylvisorex and one in Eurasia), at the base of the Crocidura radiation. Therefore, we recommend restricting the genus Suncus to the Palaearctic and Oriental taxa, and to consider all the African Suncus as Sylvisorex. Using molecular dating and biogeographic reconstruction analyses, we suggest a Palaearctic-Oriental origin for Crocidura dating back to the Upper Miocene (6.8 million years ago) and several subsequent colonisations of the Afrotropical region by independent lineages of Crocidura.

Genes associated with retinitis pigmentosa and allied diseases are frequently mutated in the general population.

BACKGROUND: Retinitis pigmentosa and other hereditary retinal degenerations (HRD) are rare genetic diseases leading to progressive blindness. Recessive HRD are caused by mutations in more than 100 different genes. Laws of population genetics predict that, on a purely theoretical ground, such a high number of genes should translate into an extremely elevated frequency of unaffected carriers of mutations. In this study we estimate the proportion of these individuals within the general population, via the analyses of data from whole-genome sequencing. METHODOLOGY/PRINCIPAL FINDINGS: We screened complete and high-quality genome sequences from 46 control individuals from various world populations for HRD mutations, using bioinformatic tools developed in-house. All mutations detected in silico were validated by Sanger sequencing. We identified clear-cut, null recessive HRD mutations in 10 out of the 46 unaffected individuals analyzed (∼22%). CONCLUSIONS/SIGNIFICANCE: Based on our data, approximately one in 4-5 individuals from the general population may be a carrier of null mutations that are responsible for HRD. This would be the highest mutation carrier frequency so far measured for a class of Mendelian disorders, especially considering that missenses and other forms of pathogenic changes were not included in our assessment. Among other things, our results indicate that the risk for a consanguineous couple of generating a child with a blinding disease is particularly high, compared to other genetic conditions.