Is Rapid Adaptation to New Environments Fueled by Old Mutations? A Case Study of Copper Tolerance on Mimulus guttatus

Rapid adaptation and tolerance is a phenomenon experienced by a variety of organisms typically because of new and harsh environments. Mimulus guttatus, a plant commonly seen on the west coast of the United States, is a prime example as it has rapidly evolved to soil contamination by copper due to mining in California in the last 150 years. There have been two hypotheses posed by researchers as to the genetic basis of how organisms have evolved so quickly which I set out to study: 1) There is a low frequency of tolerant genotypes in the ancestral population otherwise known as standing variation or 2) new mutations occurred once exposed to a new environment. In the past, researchers found it difficult to distinguish between the two because they lacked the technology we have today for DNA analysis. I used four different populations of M. guttatus from varying locations in order to address which hypothesis was valid. I conducted both survival assays of these populations and DNA analysis of known tolerant and non-tolerant lines using a copper oxidase gene. I found that there was at least some degree of tolerance in all populations in the survival assays, supporting the hypothesis of standing variation. I also found patterns within DNA analysis suggesting the copper oxidase gene would be useful for further study to verify the standing variation hypothesis. The results from this experiment helps in understanding rapid evolution not just in the context of soil contamination by metals but also ties back to why an alarming number of species are not able to adapt to our constantly changing world.

PCR-Based Analysis of Mitochondrial DNA Copy Number, Mitochondrial DNA Damage, and Nuclear DNA Damage.

Because of the role that DNA damage and depletion play in human disease, it is important to develop and improve tools to assess these endpoints. This unit describes PCR-based methods to measure nuclear and mitochondrial DNA damage and copy number. Long amplicon quantitative polymerase chain reaction (LA-QPCR) is used to detect DNA damage by measuring the number of polymerase-inhibiting lesions present based on the amount of PCR amplification; real-time PCR (RT-PCR) is used to calculate genome content. In this unit, we provide step-by-step instructions to perform these assays in Homo sapiens, Mus musculus, Rattus norvegicus, Caenorhabditis elegans, Drosophila melanogaster, Danio rerio, Oryzias latipes, Fundulus grandis, and Fundulus heteroclitus, and discuss the advantages and disadvantages of these assays.

In situ synthesis of DNA microarray on functionalized cyclic olefin copolymer substrate.

Thermoplastic materials such as cyclic-olefin copolymers (COC) provide a versatile and cost-effective alternative to the traditional glass or silicon substrate for rapid prototyping and industrial scale fabrication of microdevices. To extend the utility of COC as an effective microarray substrate, we developed a new method that enabled for the first time in situ synthesis of DNA oligonucleotide microarrays on the COC substrate. To achieve high-quality DNA synthesis, a SiO(2) thin film array was prepatterned on the inert and hydrophobic COC surface using RF sputtering technique. The subsequent in situ DNA synthesis was confined to the surface of the prepatterned hydrophilic SiO(2) thin film features by precision delivery of the phosphoramidite chemistry using an inkjet DNA synthesizer. The in situ SiO(2)-COC DNA microarray demonstrated superior quality and stability in hybridization assays and thermal cycling reactions. Furthermore, we demonstrate that pools of high-quality mixed-oligos could be cleaved off the SiO(2)-COC microarrays and used directly for construction of DNA origami nanostructures. It is believed that this method will not only enable synthesis of high-quality and low-cost COC DNA microarrays but also provide a basis for further development of integrated microfluidics microarrays for a broad range of bioanalytical and biofabrication applications.

Separating DNA with different topologies by atomic force microscopy in comparison with gel electrophoresis.

Atomic force microscopy, which is normally used for DNA imaging to gain qualitative results, can also be used for quantitative DNA research, at a single-molecular level. Here, we evaluate the performance of AFM imaging specifically for quantifying supercoiled and relaxed plasmid DNA fractions within a mixture, and compare the results with the bulk material analysis method, gel electrophoresis. The advantages and shortcomings of both methods are discussed in detail. Gel electrophoresis is a quick and well-established quantification method. However, it requires a large amount of DNA, and needs to be carefully calibrated for even slightly different experimental conditions for accurate quantification. AFM imaging is accurate, in that single DNA molecules in different conformations can be seen and counted. When used carefully with necessary correction, both methods provide consistent results. Thus, AFM imaging can be used for DNA quantification, as an alternative to gel electrophoresis.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

Altered gene expression and DNA damage in peripheral blood cells from Friedreich's ataxia patients: cellular model of pathology.

The neurodegenerative disease Friedreich's ataxia (FRDA) is the most common autosomal-recessively inherited ataxia and is caused by a GAA triplet repeat expansion in the first intron of the frataxin gene. In this disease, transcription of frataxin, a mitochondrial protein involved in iron homeostasis, is impaired, resulting in a significant reduction in mRNA and protein levels. Global gene expression analysis was performed in peripheral blood samples from FRDA patients as compared to controls, which suggested altered expression patterns pertaining to genotoxic stress. We then confirmed the presence of genotoxic DNA damage by using a gene-specific quantitative PCR assay and discovered an increase in both mitochondrial and nuclear DNA damage in the blood of these patients (p<0.0001, respectively). Additionally, frataxin mRNA levels correlated with age of onset of disease and displayed unique sets of gene alterations involved in immune response, oxidative phosphorylation, and protein synthesis. Many of the key pathways observed by transcription profiling were downregulated, and we believe these data suggest that patients with prolonged frataxin deficiency undergo a systemic survival response to chronic genotoxic stress and consequent DNA damage detectable in blood. In conclusion, our results yield insight into the nature and progression of FRDA, as well as possible therapeutic approaches. Furthermore, the identification of potential biomarkers, including the DNA damage found in peripheral blood, may have predictive value in future clinical trials.

Mitochondrial DNA polymorphism A4917G is independently associated with age-related macular degeneration.

The objective of this study was to determine if MTND2*LHON4917G (4917G), a specific non-synonymous polymorphism in the mitochondrial genome previously associated with neurodegenerative phenotypes, is associated with increased risk for age-related macular degeneration (AMD). A preliminary study of 393 individuals (293 cases and 100 controls) ascertained at Vanderbilt revealed an increased occurrence of 4917G in cases compared to controls (15.4% vs.9.0%, p = 0.11). Since there was a significant age difference between cases and controls in this initial analysis, we extended the study by selecting Caucasian pairs matched at the exact age at examination. From the 1547 individuals in the Vanderbilt/Duke AMD population association study (including 157 in the preliminary study), we were able to match 560 (280 cases and 280 unaffected) on exact age at examination. This study population was genotyped for 4917G plus specific AMD-associated nuclear genome polymorphisms in CFH, LOC387715 and ApoE. Following adjustment for the listed nuclear genome polymorphisms, 4917G independently predicts the presence of AMD (OR = 2.16, 95%CI 1.20-3.91, p = 0.01). In conclusion, a specific mitochondrial polymorphism previously implicated in other neurodegenerative phenotypes (4917G) appears to convey risk for AMD independent of recently discovered nuclear DNA polymorphisms.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.

Membrane binding of plasmid DNA and endocytic pathways are involved in electrotransfection of mammalian cells.

Electric field mediated gene delivery or electrotransfection is a widely used method in various studies ranging from basic cell biology research to clinical gene therapy. Yet, mechanisms of electrotransfection are still controversial. To this end, we investigated the dependence of electrotransfection efficiency (eTE) on binding of plasmid DNA (pDNA) to plasma membrane and how treatment of cells with three endocytic inhibitors (chlorpromazine, genistein, dynasore) or silencing of dynamin expression with specific, small interfering RNA (siRNA) would affect the eTE. Our data demonstrated that the presence of divalent cations (Ca(2+) and Mg(2+)) in electrotransfection buffer enhanced pDNA adsorption to cell membrane and consequently, this enhanced adsorption led to an increase in eTE, up to a certain threshold concentration for each cation. Trypsin treatment of cells at 10 min post electrotransfection stripped off membrane-bound pDNA and resulted in a significant reduction in eTE, indicating that the time period for complete cellular uptake of pDNA (between 10 and 40 min) far exceeded the lifetime of electric field-induced transient pores (∼10 msec) in the cell membrane. Furthermore, treatment of cells with the siRNA and all three pharmacological inhibitors yielded substantial and statistically significant reductions in the eTE. These findings suggest that electrotransfection depends on two mechanisms: (i) binding of pDNA to cell membrane and (ii) endocytosis of membrane-bound pDNA.

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae).

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Association between DNA damage response and repair genes and risk of invasive serous ovarian cancer.

BACKGROUND: We analyzed the association between 53 genes related to DNA repair and p53-mediated damage response and serous ovarian cancer risk using case-control data from the North Carolina Ovarian Cancer Study (NCOCS), a population-based, case-control study. METHODS/PRINCIPAL FINDINGS: The analysis was restricted to 364 invasive serous ovarian cancer cases and 761 controls of white, non-Hispanic race. Statistical analysis was two staged: a screen using marginal Bayes factors (BFs) for 484 SNPs and a modeling stage in which we calculated multivariate adjusted posterior probabilities of association for 77 SNPs that passed the screen. These probabilities were conditional on subject age at diagnosis/interview, batch, a DNA quality metric and genotypes of other SNPs and allowed for uncertainty in the genetic parameterizations of the SNPs and number of associated SNPs. Six SNPs had Bayes factors greater than 10 in favor of an association with invasive serous ovarian cancer. These included rs5762746 (median OR(odds ratio)(per allele) = 0.66; 95% credible interval (CI) = 0.44-1.00) and rs6005835 (median OR(per allele) = 0.69; 95% CI = 0.53-0.91) in CHEK2, rs2078486 (median OR(per allele) = 1.65; 95% CI = 1.21-2.25) and rs12951053 (median OR(per allele) = 1.65; 95% CI = 1.20-2.26) in TP53, rs411697 (median OR (rare homozygote) = 0.53; 95% CI = 0.35 - 0.79) in BACH1 and rs10131 (median OR( rare homozygote) = not estimable) in LIG4. The six most highly associated SNPs are either predicted to be functionally significant or are in LD with such a variant. The variants in TP53 were confirmed to be associated in a large follow-up study. CONCLUSIONS/SIGNIFICANCE: Based on our findings, further follow-up of the DNA repair and response pathways in a larger dataset is warranted to confirm these results.

No association between TGFB1 polymorphisms and late radiotherapy toxicity: a meta-analysis.

BACKGROUND: Transforming growth factor-beta 1 (TGF-β1) protein may be multifunctional and related to the development of fibrosis, induction of apoptosis, extracellular signaling and inhibition of proliferation in response to radiation-induced DNA damage. Several studies have investigated associations between single nucleotide polymorphisms (SNPs) in the TGFB1 gene and risk of late radiation-induced injury of normal tissue, but the conclusions remain controversial. METHODS: We searched three electronic databases (i.e., MEDLINE, EMBASE and EBSCO) for eligible publications and performed a meta-analysis assessing the association of three commonly studied SNPs in TGFB1 (i.e., rs1800469, rs1800470 and rs1800471) with risk of late radiation-induced injury of normal tissue. RESULTS: We finally included 28 case-only studies from 16 publications on aforementioned SNPs in TGFB1. However, we did not find statistical evidence of any significant association with overall risk of late radiotherapy toxicity in the pooled analysis or in further stratified analysis by cancer type, endpoint, ethnicity and sample size. CONCLUSIONS: This meta-analysis did not find statistical evidence for an association between SNPs in TGFB1 and risk of late radiation-induced injury of normal tissue, but this finding needs further confirmation by a single large study.