Differences In Gene Copy Number Research Articles

An Improved Quantitative Real-Time PCR Assay for the Enumeration of Heterosigma akashiwo (Raphidophyceae) Cysts Using a DNA Debris Removal Method and a Cyst-Based Standard Curve.

The identification and quantification of Heterosigma akashiwo cysts in sediments by light microscopy can be difficult due to the small size and morphology of the cysts, which are often indistinguishable from those of other types of algae. Quantitative real-time PCR (qPCR) based assays represent a potentially efficient method for quantifying the abundance of H. akashiwo cysts, although standard curves must be based on cyst DNA rather than on vegetative cell DNA due to differences in gene copy number and DNA extraction yield between these two cell types. Furthermore, qPCR on sediment samples can be complicated by the presence of extracellular DNA debris. To solve these problems, we constructed a cyst-based standard curve and developed a simple method for removing DNA debris from sediment samples. This cyst-based standard curve was compared with a standard curve based on vegetative cells, as vegetative cells may have twice the gene copy number of cysts. To remove DNA debris from the sediment, we developed a simple method involving dilution with distilled water and heating at 75°C. A total of 18 sediment samples were used to evaluate this method. Cyst abundance determined using the qPCR assay without DNA debris removal yielded results up to 51-fold greater than with direct counting. By contrast, a highly significant correlation was observed between cyst abundance determined by direct counting and the qPCR assay in conjunction with DNA debris removal (r2 = 0.72, slope = 1.07, p < 0.001). Therefore, this improved qPCR method should be a powerful tool for the accurate quantification of H. akashiwo cysts in sediment samples.

Chromosome position effects on gene expression in Escherichia coli K-12.

In eukaryotes, the location of a gene on the chromosome is known to affect its expression, but such position effects are poorly understood in bacteria. Here, using Escherichia coli K-12, we demonstrate that expression of a reporter gene cassette, comprised of the model E. coli lac promoter driving expression of gfp, varies by ∼300-fold depending on its precise position on the chromosome. At some positions, expression was more than 3-fold higher than at the natural lac promoter locus, whereas at several other locations, the reporter cassette was completely silenced: effectively overriding local lac promoter control. These effects were not due to differences in gene copy number, caused by partially replicated genomes. Rather, the differences in gene expression occur predominantly at the level of transcription and are mediated by several different features that are involved in chromosome organization. Taken together, our findings identify a tier of gene regulation above local promoter control and highlight the importance of chromosome position effects on gene expression profiles in bacteria.

Evaluating the bottlenecks of recombinant IgM production in mammalian cells

Despite the fact, that monoclonal antibodies are the fastest growing group of biopharmaceuticals in development, this is not true for the IgM class, which remains as enigmatic as ever. While more examples of usefulness of IgMs for medical applications are emerging, their recombinant production is still not common. In our study, stable monoclonal IgM producing CHO DG44 and HEK 293 cell lines, expressing two model IgM molecules (IgM-617 and IgM-012) were established. Recombinant cell lines were compared in regard of specific productivity, specific growth rate, maximal achieved antibody titer, gene copy numbers and transcription levels of transgene. IgM-617 cell lines were identified as high while IgM-012 clones were low producers. Although differences in gene copy numbers as well as in transcription levels were observed, they did not seem to be a limitation. Levels of relevant endoplasmic reticulum-stress related proteins were analyzed and no indications of unfolded protein response were detected. This could indicate that the difference in the intrinsic protein stability of our model proteins (as was previously observed on purified samples) might cause lower yields of IgM-012. Transcriptomics and/or proteomics follow up studies might be necessary for identification of potential bottlenecks in IgM producing cell lines.

Diversity of Class I Major Histocompatibility Genes From Divergent Haplotypes in Zebrafish

Abstract 1058The zebrafish is an important animal model for stem cell biology, cancer, and immunology research. Histocompatibility represents a key intersection of these disciplines, particularly in the context of transplantation experiments that distinguish between autologous and allogeneic tissues. Major histocompatibility (MH) genes are considered the most polymorphic genes in vertebrates, yet this immense variation occurs while maintaining conserved roles in antigen presentation. Histocompatibility in zebrafish remains poorly understood, requiring the identification of the classical MH genes as well as an evaluation of the variation between haplotypes. Although at least 11 putative zebrafish Class I MH U lineage genes have been isolated from cDNA libraries, their genomic organization and haplotype assignments remain uncharacterized for the majority.We focused our study on a set of diverse zebrafish Class I MH genes that segregate with specific haplotypes at chromosome 19, and for which donor-recipient matching was previously shown to be associated with improved engraftment after transplantation. Inside the conserved psmb8 and tpsn flanking gene regions on chromosome 19, Class I MH haplotypes can differ markedly among zebrafish strains including Tubingen and AB. Interestingly, the distinct haplotypes at chromosome 19 appear to maintain non-overlapping sets of genes, and also have gene copy number differences. For example, haplotype A contains the genes mhc1uda, mhc1uea, and mhc1ufa, with tap2 genes located in between each set of MH genes. In contrast, haplotype B contains mhc1uba and mhc1uca genes that are separated instead by a tapbp gene.Among these two highly divergent haplotypes, mhc1uea and mhc1uda from haplotype A appear to be more closely related to one another than to either of the MH genes on haplotype B. Similarly, the mhc1uba and mhc1uca genes from haplotype B share higher levels of sequence identity than in comparison with genes from haplotype A. These findings differ markedly from the MH genes in mammals where the highest similarity is between alleles of the same Class I gene (for example, HLA-A in humans), instead of the other genes on the same haplotype (eg. HLA-B and HLA-C). In addition, the zebrafish MH gene sequences are very closely related to MH genes from other fish species. As an example, mhc1ufa is more similar to a Class I MH gene from gibuna carp than to other zebrafish MH genes. These data in zebrafish indicate that selected Class I MH gene polymorphisms have been preserved within the species since before the species diverged from other teleost fishes, suggesting a selective advantage for this unique form of MH diversity. To determine mRNA expression levels of the Class I MH genes, quantitative PCR was performed on liver, spleen, kidney marrow, testis, gill, intestine and heart from zebrafish that were homozygous for either haplotype A or haplotype B. Within haplotype A, expression of mhc1uda, mhc1uea and mhc1ufa is relatively equivalent, with mhc1ufa expressed at the highest levels overall. For haplotype B, mhc1uba has the highest levels in all tissues examined, while mhc1uca expression is 10–100 fold lower and somewhat variable between tissues. Putative peptide anchor residues are highly conserved between species. Those anchor residues that are conserved in all species including mammals, as well as two additional residues that are conserved in teleosts, are also all conserved in mhc1uda, mhc1uba, and mhc1uca. In contrast, mhc1uea and mhce1ufa each have single amino acid substitutions at critical residues, K146N and Y59F, respectively. These substitutions do not necessarily disqualify these molecules from consideration as classical MH Class I genes.In conclusion, we propose that mhc1uba and mhc1uda function as classical MH Class I genes within divergent haplotypes, based on mRNA expression levels and tissue distribution, as well as sequence properties including conservation of putative peptide anchor residues. Using these same criteria, additional zebrafish Class I MH genes may be identified to serve similar roles. These predictions require verification via additional functional transplantation experiments that are currently underway in our lab. Defining the functional zebrafish Class I MH genes will provide an important foundation for future studies in immunology and transplantation. Disclosures:No relevant conflicts of interest to declare.

Deciphering the genomic structure, function and evolution of carotenogenesis related phytoene synthases in grasses

BackgroundCarotenoids are isoprenoid pigments, essential for photosynthesis and photoprotection in plants. The enzyme phytoene synthase (PSY) plays an essential role in mediating condensation of two geranylgeranyl diphosphate molecules, the first committed step in carotenogenesis. PSY are nuclear enzymes encoded by a small gene family consisting of three paralogous genes (PSY1-3) that have been widely characterized in rice, maize and sorghum.ResultsIn wheat, for which yellow pigment content is extremely important for flour colour, only PSY1 has been extensively studied because of its association with QTLs reported for yellow pigment whereas PSY2 has been partially characterized. Here, we report the isolation of bread wheat PSY3 genes from a Renan BAC library using Brachypodium as a model genome for the Triticeae to develop Conserved Orthologous Set markers prior to gene cloning and sequencing. Wheat PSY3 homoeologous genes were sequenced and annotated, unravelling their novel structure associated with intron-loss events and consequent exonic fusions. A wheat PSY3 promoter region was also investigated for the presence of cis-acting elements involved in the response to abscisic acid (ABA), since carotenoids also play an important role as precursors of signalling molecules devoted to plant development and biotic/abiotic stress responses. Expression of wheat PSYs in leaves and roots was investigated during ABA treatment to confirm the up-regulation of PSY3 during abiotic stress.ConclusionsWe investigated the structural and functional determinisms of PSY genes in wheat. More generally, among eudicots and monocots, the PSY gene family was found to be associated with differences in gene copy numbers, allowing us to propose an evolutionary model for the entire PSY gene family in Grasses.

Genomic diversity of the human intestinal parasite Entamoeba histolytica

BackgroundEntamoeba histolytica is a significant cause of disease worldwide. However, little is known about the genetic diversity of the parasite. We re-sequenced the genomes of ten laboratory cultured lines of the eukaryotic pathogen Entamoeba histolytica in order to develop a picture of genetic diversity across the genome.ResultsThe extreme nucleotide composition bias and repetitiveness of the E. histolytica genome provide a challenge for short-read mapping, yet we were able to define putative single nucleotide polymorphisms in a large portion of the genome. The results suggest a rather low level of single nucleotide diversity, although genes and gene families with putative roles in virulence are among the more polymorphic genes. We did observe large differences in coverage depth among genes, indicating differences in gene copy number between genomes. We found evidence indicating that recombination has occurred in the history of the sequenced genomes, suggesting that E. histolytica may reproduce sexually.ConclusionsE. histolytica displays a relatively low level of nucleotide diversity across its genome. However, large differences in gene family content and gene copy number are seen among the sequenced genomes. The pattern of polymorphism indicates that E. histolytica reproduces sexually, or has done so in the past, which has previously been suggested but not proven.

Compensation for differences in gene copy number among yeast ribosomal proteins is encoded within their promoters.

Coordinate regulation of ribosomal protein (RP) genes is key for controlling cell growth. In yeast, it is unclear how this regulation achieves the required equimolar amounts of the different RP components, given that some RP genes exist in duplicate copies, while others have only one copy. Here, we tested whether the solution to this challenge is partly encoded within the DNA sequence of the RP promoters, by fusing 110 different RP promoters to a fluorescent gene reporter, allowing us to robustly detect differences in their promoter activities that are as small as ~10%. We found that single-copy RP promoters have significantly higher activities, suggesting that proper RP stoichiometry is indeed partly encoded within the RP promoters. Notably, we also partially uncovered how this regulation is encoded by finding that RP promoters with higher activity have more nucleosome-disfavoring sequences and characteristic spatial organizations of these sequences and of binding sites for key RP regulators. Mutations in these elements result in a significant decrease of RP promoter activity. Thus, our results suggest that intrinsic (DNA-dependent) nucleosome organization may be a key mechanism by which genomes encode biologically meaningful promoter activities. Our approach can readily be applied to uncover how transcriptional programs of other promoters are encoded.

The draft genome and transcriptome of Cannabis sativa

BackgroundCannabis sativa has been cultivated throughout human history as a source of fiber, oil and food, and for its medicinal and intoxicating properties. Selective breeding has produced cannabis plants for specific uses, including high-potency marijuana strains and hemp cultivars for fiber and seed production. The molecular biology underlying cannabinoid biosynthesis and other traits of interest is largely unexplored.ResultsWe sequenced genomic DNA and RNA from the marijuana strain Purple Kush using shortread approaches. We report a draft haploid genome sequence of 534 Mb and a transcriptome of 30,000 genes. Comparison of the transcriptome of Purple Kush with that of the hemp cultivar 'Finola' revealed that many genes encoding proteins involved in cannabinoid and precursor pathways are more highly expressed in Purple Kush than in 'Finola'. The exclusive occurrence of Δ9-tetrahydrocannabinolic acid synthase in the Purple Kush transcriptome, and its replacement by cannabidiolic acid synthase in 'Finola', may explain why the psychoactive cannabinoid Δ9-tetrahydrocannabinol (THC) is produced in marijuana but not in hemp. Resequencing the hemp cultivars 'Finola' and 'USO-31' showed little difference in gene copy numbers of cannabinoid pathway enzymes. However, single nucleotide variant analysis uncovered a relatively high level of variation among four cannabis types, and supported a separation of marijuana and hemp.ConclusionsThe availability of the Cannabis sativa genome enables the study of a multifunctional plant that occupies a unique role in human culture. Its availability will aid the development of therapeutic marijuana strains with tailored cannabinoid profiles and provide a basis for the breeding of hemp with improved agronomic characteristics.

Detection of hemi/homozygotes through heteroduplex formation in high-resolution melting analysis

Heteroduplex formation, required for the complete detection of hemi/homozygotes using high-resolution melting analysis, can be induced either by pre-PCR mixing of genomic DNAs or by post-PCR mixing of PCR products from unknown and reference samples. This study investigates the effects of both methods using two single nucleotide polymorphisms in X-linked DMD gene. The results show that both methods resulted in the same effect when mixing samples with the same gene copy number. Mixing samples with different gene copy numbers has not been previously explored and we show that post-PCR mixing is insensitive to gene copy number differences as compared to pre-PCR mixing.

Auxin and monocot development.

Monocots are known to respond differently to auxinic herbicides; hence, certain herbicides kill broadleaf (i.e., dicot) weeds while leaving lawns (i.e., monocot grasses) intact. In addition, the characters that distinguish monocots from dicots involve structures whose development is controlled by auxin. However, the molecular mechanisms controlling auxin biosynthesis, homeostasis, transport, and signal transduction appear, so far, to be conserved between monocots and dicots, although there are differences in gene copy number and expression leading to diversification in function. This article provides an update on the conservation and diversification of the roles of genes controlling auxin biosynthesis, transport, and signal transduction in root, shoot, and reproductive development in rice and maize.

Dosage compensation and the global re-balancing of aneuploid genomes

Diploid genomes are exquisitely balanced systems of gene expression. The dosage-compensation systems that evolved along with monosomic sex chromosomes exemplify the intricacies of compensating for differences in gene copy number by transcriptional regulation.

Auditioning of CHO host cell lines using the artificial chromosome expression (ACE) technology

In order to maximize recombinant protein expression in mammalian cells many factors need to be considered such as transfection method, vector construction, screening techniques and culture conditions. In addition, the host cell line can have a profound effect on the protein expression. However, auditioning or directly comparing host cell lines for optimal protein expression may be difficult since most transfection methods are based on random integration of the gene of interest into the host cell genome. Thus it is not possible to determine whether differences in expression between various host cell lines are due to the phenotype of the host cell itself or genetic factors such as gene copy number or gene location. To improve cell line generation, the ACE System was developed based on pre-engineered artificial chromosomes with multiple recombination acceptor sites. This system allows for targeted transfection and has been effectively used to rapidly generate stable CHO cell lines expressing high levels of monoclonal antibody. A key feature of the ACE System is the ability to isolate and purify ACEs containing the gene(s) of interest and transfect the same ACEs into different host cell lines. This feature allows the direct auditioning of host cells since the host cells have been transfected with ACEs that contain the same number of gene copies in the same genetic environment. To investigate this audition feature, three CHO host cell lines (CHOK1SV, CHO-S and DG44) were transfected with the same ACE containing gene copies of a human monoclonal IgG1 antibody. Clonal cell lines were generated allowing a direct comparison of antibody expression and stability between the CHO host cells. Results showed that the CHOK1SV host cell line expressed antibody at levels of more than two to five times that for DG44 and CHO-S host cell lines, respectively. To confirm that the ACE itself was not responsible for the low antibody expression seen in the CHO-S based clones, the ACE was isolated and purified from these cells and transfected back into fresh CHOK1SV cells. The resulting expression of the antibody from the ACE newly transfected into CHOK1SV increased fivefold compared to its expression in CHO-S and confirmed that the differences in expression between the different CHO host cells was due to the cell phenotype rather than differences in gene copy number and/or location. These results demonstrate the utility of the ACE System in providing a rapid and direct technique for auditioning host cell lines for optimal recombinant protein expression.

Studying copy number variations using a nanofluidic platform

Copy number variations (CNVs) in the human genome are conventionally detected using high-throughput scanning technologies, such as comparative genomic hybridization and high-density single nucleotide polymorphism (SNP) microarrays, or relatively low-throughput techniques, such as quantitative polymerase chain reaction (PCR). All these approaches are limited in resolution and can at best distinguish a twofold (or 50%) difference in copy number. We have developed a new technology to study copy numbers using a platform known as the digital array, a nanofluidic biochip capable of accurately quantitating genes of interest in DNA samples. We have evaluated the digital array's performance using a model system, to show that this technology is exquisitely sensitive, capable of differentiating as little as a 15% difference in gene copy number (or between 6 and 7 copies of a target gene). We have also analyzed commercial DNA samples for their CYP2D6 copy numbers and confirmed that our results were consistent with those obtained independently using conventional techniques. In a screening experiment with breast cancer and normal DNA samples, the ERBB2 gene was found to be amplified in about 35% of breast cancer samples. The use of the digital array enables accurate measurement of gene copy numbers and is of significant value in CNV studies.

The CGView Server: a comparative genomics tool for circular genomes

The CGView Server generates graphical maps of circular genomes that show sequence features, base composition plots, analysis results and sequence similarity plots. Sequences can be supplied in raw, FASTA, GenBank or EMBL format. Additional feature or analysis information can be submitted in the form of GFF (General Feature Format) files. The server uses BLAST to compare the primary sequence to up to three comparison genomes or sequence sets. The BLAST results and feature information are converted to a graphical map showing the entire sequence, or an expanded and more detailed view of a region of interest. Several options are included to control which types of features are displayed and how the features are drawn. The CGView Server can be used to visualize features associated with any bacterial, plasmid, chloroplast or mitochondrial genome, and can aid in the identification of conserved genome segments, instances of horizontal gene transfer, and differences in gene copy number. Because a collection of sequences can be used in place of a comparison genome, maps can also be used to visualize regions of a known genome covered by newly obtained sequence reads. The CGView Server can be accessed at http://stothard.afns.ualberta.ca/cgview_server/

Comparative Genomic Sequence and Expression Analyses ofMedicago truncatulaand Alfalfa Subspeciesfalcata COLD-ACCLIMATION-SPECIFICGenes

In Arabidopsis (Arabidopsis thaliana) the low-temperature induction of genes encoding the C-REPEAT BINDING FACTOR (CBF) transcriptional activators is a key step in cold acclimation. CBFs in turn activate a battery of downstream genes known as the CBF regulon, which collectively act to increase tolerance to low temperatures. Fundamental questions are: What determines the size and scope of the CBF regulon, and is this is a major determinant of the low-temperature tolerance capacity of individual plant species? Here we have begun to address these questions through comparative analyses of Medicago truncatula and Medicago sativa subsp. falcata. M. truncatula survived to -4 degrees C but did not cold acclimate, whereas Medicago falcata cold acclimated and survived -14 degrees C. Both species possessed low-temperature-induced CBFs but differed in the expression of the COLD-ACCLIMATION-SPECIFIC (CAS) genes, which are candidate CBF targets. M. falcata CAS30 was robustly cold-responsive whereas the MtCAS31 homolog was not. M. falcata also possessed additional CAS30 homologs in comparison to the single CAS31 gene in M. truncatula. MfCAS30 possessed multiple pairs of closely spaced C-REPEAT/DEHYDRATION RESPONSIVE ELEMENT (CRT/DRE) motifs, the cognate CBF binding site in its upstream region whereas MtCAS31 lacked one CRT/DRE partner of the two proximal partner pairs. CAS genes also shared a promoter structure comprising modules proximal and distal to the coding sequence. CAS15, highly cold-responsive in both species, harbored numerous CRT/DRE motifs, but only in the distal module. However, fusion of the MtCAS15 promoter, including the distal module, to a reporter gene did not result in low-temperature responsiveness in stably transformed Arabidopsis. In contrast, both MtCAS31 and MfCAS30 promoter fusions were low-temperature responsive, although the MfCAS31 fusion was less robust than the MfCAS30 fusion. From these studies we conclude that CAS genes harbor CRT/DRE motifs, their proximity to one another is likely key to regulatory output in Medicago, and they may be located kilobases distal to the transcriptional start site. We hypothesize that these differences in CRT/DRE copy numbers in CAS30/CAS31 upstream regions combined with differences in gene copy numbers may be a factor in determining differences in low-temperature tolerance between M. truncatula and M. falcata.

A comparative study of genome-wide SNP, CGH microarray and protein expression analysis to explore genotypic and phenotypic mechanisms of acquired antiestrogen resistance in breast cancer

Allelic imbalance is a common feature of many malignancies. We have measured allelic imbalance in genomic DNA from the breast cancer cell lines T47D, MDA-MB-231, two antiestrogen sensitive (MCF7N and MCF7L) and two resistant MCF7 cell lines (MMU2 and LCC9) using single nucleotide polymorphism (SNP) oligonucleotide microarrays. DNA from MCF7(L) and MMU2 cells was also analysed by comparative genome hybridisation (CGH) to compare with SNP microarray data. Proteins previously determined to be involved in disease progression were quantified by Western blot and compared to array data. The SNP and CGH array both detected cytogenetic abnormalities commonly found in breast cancer: amplification of chromosomes 11q13-14.1, 17q and 20q containing cyclin D1, BCAS1 and 3 (Breast Cancer Amplified Sequence) and AIB1 (Amplified in Breast cancer) genes; losses at 6q, 9p and X chromosomes, which included ERalpha (Estrogen Receptor alpha) and p16 ( INK4A ) genes. However the SNP chip array data additionally identified regions of loss of heterozygosity (LOH) followed by duplication of the remaining allele-uniparental disomy (UPD). Good concordance between SNP arrays and CGH analyses was observed, however there was poor correlation between gene copy number and protein levels between the cell lines. There were reductions in ERalpha, cyclin D1 and p27 protein levels whilst p21 protein levels were elevated in antiestrogen resistant MCF7 cell lines. Although protein levels varied there was no difference in gene copy number. This study shows SNP and CGH array analysis are powerful tools for analysis of allelic imbalance in breast cancer. However, the antiestrogen resistant phenotype was likely to be due to changes in gene expression and protein degradation rather than in altered gene copy number.

Human SULT1A1 gene: copy number differences and functional implications.

SULT1A1, which catalyzes the sulfate conjugation of a wide variety of natural and synthetic compounds, is genetically polymorphic. Biochemical and pharmacogenetic studies have demonstrated that individual variation in the level of enzyme activity is inherited. Common single-nucleotide polymorphisms (SNPs) located in the open reading frame and in the 5'-flanking region (5'-FR) may account for a portion of this individual variation. In this study, we demonstrate the presence of SULT1A1 gene deletions and duplications, representing an additional source of variability in the metabolic activity of this enzyme. A quantitative multiplex PCR assay was used to measure the extent of copy number differences and the frequency of these events in different populations. An analysis of DNA from 362 Caucasian-American and 99 African-American showed the presence of 1 to approximately 5 copies of SULT1A1 in individual samples: 5% of Caucasian subjects contained a single copy of the gene and 26% had three or more copies, while 63% of African-American subjects had three or more copies. Analysis of the genomic region surrounding the SULT1A1 gene in three separate cases with a deletion demonstrated that the entire SULT1A1 gene was affected. Reporter assays, constructed for each of the various 5'-FR SNP haplotypes, suggest that these may also play a role in SULT1A1 activity. However, the variability in the level of enzyme activity among 23 human platelet and 267 human liver samples was best explained by gene copy number differences when all sources of genetic variability were considered (P < 0.0001). Overall, these observations have obvious implications for the effectiveness of SULT1A1 as a drug and hormone metabolizing enzyme and its potential role as a risk factor for disease.

Prognostic Value of PAI1 in Invasive Breast Cancer: Evidence That Tumor-Specific Factors Are More Important Than Genetic Variation in Regulating PAI1 Expression

Plasminogen activator inhibitor-1 (PAI1) can promote cancer progression, and its protein expression in tumors is an independent indicator of poor prognosis in many forms of cancer. Here, we show that high PAI1 mRNA levels also predict for shorter overall survival in two independent breast cancer data sets, highlighting the importance of its transcriptional regulation. The -675insG (4G/5G) single-nucleotide polymorphism in the PAI1 gene promoter has been shown to influence PAI1 transcription, with the 4G allele eliciting higher reporter gene expression in vitro and higher levels of circulating PAI1 in vivo. Nevertheless, its genotypic distribution in 2,539 British women with invasive breast cancer was virtually identical to that seen in 1,832 matched controls (P = 0.72), and annual mortality rates for 4G4G, 4G5G, and 5G5G cases were 2.6%, 2.8%, and 3.1% per year, respectively (P = 0.10). Thus, there was no association with breast cancer incidence or outcome, and in a separate set of breast cancers, the 4G/5G single-nucleotide polymorphism showed no association with PAI1 mRNA expression (P = 0.85). By contrast, connective tissue growth factor (CTGF), which can regulate PAI1 expression in culture, was associated with PAI1 expression in three independent cohorts (P << 0.0001). In addition, PAI1 gene copy number differences in the tumors were correlated with PAI1 mRNA expression (P = 0.0005) and seemed to affect expression independently of CTGF. Thus, local factors, such as CTGF and genomic amplification, seem to be more important than germ line genetic variation in influencing PAI1 expression and its untoward effects in breast cancer.

Gene capture and random amplification for quantitative recovery of homologous genes

The polymerase chain reaction (PCR) is instrumental in molecular analysis of microorganisms, allowing for the selective amplification of nucleic acids directly from clinical and environmental samples. However, the principles that allow for targeted amplification of DNA become a hindrance when attempting to simultaneously discriminate and quantify complex mixtures of homologous genes. Here we present a simple solution to the quantitative problem by separating the enrichment and amplification aspects of a conventional PCR reaction. In this assay, genes are enriched using a DNA oligonucleotide capture probe and subsequently amplified in a two-step random amplification protocol. In order to evaluate the quantitative aspects of the gene capture assay, we used real-time quantitative-PCR to measure initial and final concentrations of homologous genes from constructed mixtures of genomes. Upon sampling for the universal DNA-dependent RNA polymerase gene, rpoC, we were able to demonstrate quantitative recoveries from a mixed DNA sample despite differences in gene copy number ranging up to 4 orders of magnitude. This suggests that minority populations as low as 0.01% of the total community are represented as accurately as populations at higher abundance. These results offer new possibilities for accurately and quantitatively monitoring diverse mixtures of microorganisms.

Divergence in gene expression related to variation in host specificity of an ectomycorrhizal fungus

Ectomycorrhizae are formed by mutualistic interactions between fungi and the roots of woody plants. During symbiosis the two organisms exchange carbon and nutrients in a specific tissue that is formed at the contact between a compatible fungus and plant. There is considerable variation in the degree of host specificity among species and strains of ectomycorrhizal fungi. In this study, we have for the first time shown that this variation is associated with quantitative differences in gene expression, and with divergence in nucleotide sequences of symbiosis-regulated genes. Gene expression and sequence evolution were compared in different strains of the ectomycorrhizal fungus Paxillus involutus; the strains included Nau, which is not compatible with birch and poplar, and the two compatible strains Maj and ATCC200175. On a genomic level, Nau and Maj were very similar. The sequence identity was 98.9% in the 16 loci analysed, and only three out of 1075 genes analysed by microarray-based hybridizations had signals indicating differences in gene copy numbers. In contrast, 66 out of the 1075 genes were differentially expressed in Maj compared to Nau after contact with birch roots. Thirty-seven of these symbiosis-regulated genes were also differentially expressed in the ATCC strain. Comparative analysis of DNA sequences of the symbiosis-regulated genes in different strains showed that two of them have evolved at an enhanced rate in Nau. The sequence divergence can be explained by a decreased selection pressure, which in turn is determined by lower functional constraints on these proteins in Nau as compared to the compatible strains.