Ecologically and evolutionarily important SNPs identified in natural populations

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.     

Genotypic and Allelic Variability in CYP19A1 among Populations of African and European Ancestry

CYP19A1 facilitates the bioconversion of estrogens from androgens. CYP19A1 intron single nucleotide polymorphisms (SNPs) may alter mRNA splicing, resulting in altered CYP19A1 activity, and potentially influencing disease susceptibility. Genetic studies of CYP19A1 SNPs have been well documented in populations of European ancestry; however, studies in populations of African ancestry are limited. In the present study, ten ‘candidate’ intronic SNPs in CYP19A1 from 125 African Americans (AA) and 277 European Americans (EA) were genotyped and their frequencies compared. Allele frequencies were also compared with HapMap and ASW 1000 Genomes populations. We observed significant differences in the minor allele frequencies between AA and EA in six of the ten SNPs including rs10459592 (p<0.0001), rs12908960 (p<0.0001), rs1902584 (p = 0.016), rs2470144 (p<0.0001), rs1961177 (p<0.0001), and rs6493497 (p = 0.003). While there were no significant differences in allele frequencies between EA and CEU in the HapMap population, a 1.2- to 19-fold difference in allele frequency for rs10459592 (p = 0.004), rs12908960 (p = 0.0006), rs1902584 (p<0.0001), rs2470144 (p = 0.0006), rs1961177 (p<0.0001), and rs6493497 (p = 0.0092) was observed between AA and the Yoruba (YRI) population. Linkage disequilibrium (LD) blocks and haplotype clusters that is unique to the EA population but not AA was also observed. In summary, we demonstrate that differences in the allele frequencies of CYP19A1 intron SNPs are not consistent between populations of African and European ancestry. Thus, investigations into whether CYP19A1 intron SNPs contribute to variations in cancer incidence, outcomes and pharmacological response seen in populations of different ancestry may prove beneficial.     

A new look towards BAC-based array CGH through a comprehensive comparison with oligo-based array CGH

Background

Natural populations of the teleost fish Fundulus heteroclitus tolerate a broad range of environmental conditions including temperature, salinity, hypoxia and chemical pollutants. Strikingly, populations of Fundulus inhabit and have adapted to highly polluted Superfund sites that are contaminated with persistent toxic chemicals. These natural populations provide a foundation to discover critical gene pathways that have evolved in a complex natural environment in response to environmental stressors.

Results

We used Fundulus cDNA arrays to compare metabolic gene expression patterns in the brains of individuals among nine populations: three independent, polluted Superfund populations and two genetically similar, reference populations for each Superfund population. We found that up to 17% of metabolic genes have evolved adaptive changes in gene expression in these Superfund populations. Among these genes, two (1.2%) show a conserved response among three polluted populations, suggesting common, independently evolved mechanisms for adaptation to environmental pollution in these natural populations.

Conclusion

Significant differences among individuals between polluted and reference populations, statistical analyses indicating shared adaptive changes among the Superfund populations, and lack of reduction in gene expression variation suggest that common mechanisms of adaptive resistance to anthropogenic pollutants have evolved independently in multiple Fundulus populations. Among three independent, Superfund populations, two genes have a common response indicating that high selective pressures may favor specific responses.     

Exploring genomic variation associated with drought stress in Picea mariana populations

Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype‐by‐environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.     

Profiling Deleterious Non-synonymous SNPs of Smoker's Gene CYP1A1

CYP1A1 gene belongs to the cytochrome P450 family and is known better as smokers’ gene due to its hyperactivation as a consequence of long term smoking. The expression of CYP1A1 induces polycyclic aromatic hydrocarbon production in the lungs, which when over expressed, is known to cause smoking related diseases, such as cardiovascular pathologies, cancer, and diabetes. Single nucleotide polymorphisms (SNPs) are the simplest form of genetic variations that occur at a higher frequency, and are denoted as synonymous and non-synonymous SNPs on the basis of their effects on the amino acids. This study adopts a systematic in silico approach to predict the deleterious SNPs that are associated with disease conditions. It is inferred that four SNPs are highly deleterious, among which the SNP with rs17861094 is commonly predicted to be harmful by all tools. Hydrophobic (isoleucine) to hydrophilic (serine) amino acid variation was observed in the candidate gene. Hence, this investigation aims to characterize a candidate gene from 159 SNPs of CYP1A1.     

Differential modes of selection on the rhodopsin gene in coastal Baltic and North Sea populations of the sand goby,Pomatoschistus minutus

An excellent model to elucidate the mechanisms and importance of evolution in the marine environment is the spectral tuning mechanism of the visual pigment in vertebrates. In the sand goby Pomatoschistus minutus (Teleostei; Gobiidae), a distribution‐wide study showed that spatial variation at the rhodopsin gene (RH1) matches the characteristics of specific light environments. This match suggests that populations are locally adapted to selective light regimes targeting the RH1 gene. If so, then the direction of selection should depend on the regional spatial and temporal stability of the light conditions. We tested this prediction by comparing goby populations from two regions: the Baltic Sea, characterized by divergent, but temporally stable light conditions, and the North Sea, characterized by locally heterogeneous and temporally variable light conditions. RH1 sequences of 491 Pomatoschistus minutus individuals from 15 locations were analysed. We found that variation at the RH1 gene in the Baltic populations showed signatures of diversifying selection, whereas the RH1 gene in the North Sea showed signatures of stabilizing selection. These different modes of selection are consistent with the regional light conditions and hence support our predictions, but may also be influenced by migration between the open sea and more turbid estuarine environments. An interesting observation is that within one gene, synonymous and non‐synonymous SNPs show a totally different pattern between populations. Population differentiation based on non‐synonymous SNPs of the RH1 gene correlated with spectral variation of the local environment of the sand goby populations. In contrast, the differentiation based on synonymous SNPs of RH1 reflects more the neutral historical pattern of the species.     

Identification of functional DNA variants in the constitutive promoter region of MDM2

Although mutations in the oncoprotein murine double minute 2 (MDM2) are rare, MDM2 gene overexpression has been observed in several human tumors. Given that even modest changes in MDM2 levels might influence the p53 tumor suppressor signaling pathway, we postulated that sequence variation in the promoter region of MDM2 could lead to disregulated expression and variation in gene dosage. Two promoters have been reported for MDM2; an internal promoter (P2), which is located near the end of intron 1 and is p53-responsive, and an upstream constitutive promoter (P1), which is p53-independent. Both promoter regions contain DNA variants that could influence the expression levels of MDM2, including the well-studied single nucleotide polymorphism (SNP) SNP309, which is located in the promoter P2; i.e., upstream of exon 2. In this report, we screened the promoter P1 for DNA variants and assessed the functional impact of the corresponding SNPs. Using the dbSNP database and genotyping validation in individuals of European descent, we identified three common SNPs (?1494?G?>?A; indel 40?bp; and ?182?C?>?G). Three major promoter haplotypes were inferred by using these three promoter SNPs together with rs2279744 (SNP309). Following subcloning into a gene reporter system, we found that two of the haplotypes significantly influenced MDM2 promoter activity in a haplotype-specific manner. Site-directed mutagenesis experiments indicated that the 40?bp insertion/deletion variation is causing the observed allelic promoter activity. This study suggests that part of the variability in the MDM2 expression levels could be explained by allelic p53-independent P1 promoter activity.     

Differential vitamin D 24-hydroxylase/CYP24A1 gene promoter methylation in endothelium from benign and malignant human prostate

Epigenetic alterations occur in tumor-associated vessels in the tumor microenvironment. Methylation of the CYP24A1 gene promoter differs in endothelial cells isolated from tumors and non-tumor microenvironments in mice. The epigenetic makeup of endothelial cells of human tumor-associated vasculature is unknown due to difficulty of isolating endothelial cells populations from a heterogeneous tissue microenvironment. To ascertain CYP24A1 promoter methylation in tumor-associated endothelium, we utilized laser microdissection guided by CD31 immunohistochemistry to procure endothelial cells from human prostate tumor specimens. Prostate tissues were obtained following robotic radical prostatectomy from men with clinically localized prostate cancer. Adjacent histologically benign prostate tissues were used to compare endothelium from benign versus tumor microenvironments. Sodium bisulfite sequencing of CYP24A1 promoter region showed that the average CYP24A1 promoter methylation in the endothelium was 20% from the tumor microenvironment compared with 8.2% in the benign microenvironment (p&lt;0.05). A 2-fold to 17-fold increase in CYP24A1 promoter methylation was observed in the prostate tumor endothelium compared with the matched benign prostate endothelium in four patient samples, while CYP24A1 remained unchanged in two patient sample. In addition, there is no correlation of the level of CYP24A1 promoter methylation in prostate tumor-associated endothelium with that of epithelium/stroma. This study demonstrates that the CYP24A1 promoter is methylated in tumor-associated endothelium, indicating that epigenetic alterations in CYP24A1 may play a role in determining the phenotype of tumor-associated vasculature in the prostate tumor microenvironment.     

Identification of genetic polymorphisms of CYP2S1 in a Finnish Caucasian population

CYP2S1 is a recently discovered member of the cytochrome P450 (CYP) gene superfamily. Interestingly, even though the DNA sequence identifies it as the sole member of the new CYP2S family, CYP2S1 exhibits many features typical to CYP1 family members, e.g. dioxin-inducibility mediated by the aryl hydrocarbon receptor (AHR) and the aryl hydrocarbon receptor nuclear translocator (ARNT). In addition, CYP2S1 metabolises some aromatic hydrocarbons as well as cellular substances. These characteristics, together with a wide extrahepatic tissue distribution, suggest that CYP2S1 may have an important role in both exogenous and endogenous metabolism. This is the first study characterising CYP2S1 alleles and naming them with the recommended CYP allele nomenclature. We used denaturing gradient gel electrophoresis (DGGE) and direct sequencing to investigate genetic variation of CYP2S1 in 100 male Finnish Caucasians. Those exons in which variation was found were examined in subsequent 100 subjects. The coding region of all of the nine exons, as well as a 449 bp fragment of the proximal promoter region, was analysed. This systematic investigation revealed eight single nucleotide polymorphisms (SNPs), which comprise nine different variant alleles (haplotypes), in addition to the wild-type allele. Seven of the SNPs occurred in the protein-coding areas and one in the proximal 3' untranslated region (3'UTR). Two of these sequence variations (10347C > T and 13106C > T) result in non-conservative amino acid substitutions, i.e. Arg380Cys and Pro466Leu, respectively. The respective allelic variants, CYP2S1*2 ([10347C > T]) and CYP2S1*3 (13106C > T; 13255A > G]), occurred in our study population at frequencies of 0.50 and 3.75%, respectively. The most common of the variant alleles was CYP2S1*1H (23.8%), harbouring a 13255A > G substitution located in the 3'UTR.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

Development and characterization of a chimaeric tissue-specific promoter in wheat and rice endosperm

The recently achieved significant improvement of cereal transformation protocols provides facilities to alter the protein composition of the endosperm, for example, to increase or decrease the quantity of one of its protein components or to express foreign molecules. To achieve this goal, strong endosperm-specific promoters have to be available. The aim of our work was to develop a more efficient tissue-specific promoter which is currently used. A chimaeric promoter was assembled using the 5′ UTR (1,900 bp) of the gene coding for the 1Bx17 HMW glutenin subunit protein, responsible for tissue-specific expression and the first intron of the rice actin gene (act1). The sequence around of the translation initial codon was optimized. The effect of the intron and promoter regulatory sequences, using different lengths of 1Bx17 HMW-GS promoter, were studied on the expression of uidA gene. The function of promoter elements, promoter length, and the first intron of the rice actin gene were tested by a transient expression assay in immature wheat endosperm and in stable transgenic rice plants. Results showed that insertion of the rice act1 first intron increased GUS expression by four times in transient assay. The shortest 1Bx17 HMW-GS promoter fragment (173 bp) linked to the intron and GUS reporter gene provided almost the same expression level than the intronless long 1Bx17 HMW-GS promoter. Analysis of the stable transformant plants revealed that 173 nucleotides were sufficient for endosperm-specific expression of the uidA gene, despite 13 nucleotides missing from the HMW enhancer sequence, a relevant regulatory element in the promoter region.     

Frequencies of four genetic polymorphisms in the <Emphasis Type="Italic">CYP1A2</Emphasis> gene in Turkish population

Cytochrome P450 (CYP) 1A2 gene is involved in the metabolic activation of several carcinogens and altered metabolization of some clinically used drugs. We aimed to investigate the distributions of genetic polymorphisms-3860 (G/A)(CYP1A2*1C) and-2467 (T/del)(CYP1A2*1D) in the 5′-flanking region and-739 (T/G)(CYP1A2*1E) and-163(C/A)(CYP1A2*1F) in the first intron of the CYP1A2 gene in 110 unrelated healthy Turkish volunteers by PCR-RFLP technique. The frequencies of each polymorphism in Turkish population were found as 0.04, 0.92, 0.01, 0.27 for CYP1A2*1C, CYP1A2*1D, CYP1A2*1E, CYP1A2*1F, respectively. Compared with other populations, CYP1A2*1Dhas been found to be significantly increased in Turkish population. On the other hand, in general, the frequencies of the other polymorphisms were concordant with those in the Egyptian and Caucasian populations, and were different from those in the Japanese, Chinese and Ethiopian populations. Our results suggest that due to increased frequency of CYP1A2*1D in Turkish population, unctional significance of CYP1A2*1D should be evaluated. It might be screened to determine the relationship between CYP1A2*1D and CYP1A2 related drug metabolisms in associated groups. This article was submitted by the authors in English.     

The genetic basis of population fecundity prediction across multiple field populations of Nilaparvata lugens

Identifying the molecular markers for complex quantitative traits in natural populations promises to provide novel insight into genetic mechanisms of adaptation and to aid in forecasting population dynamics. In this study, we investigated single nucleotide polymorphisms (SNPs) using candidate gene approach from high‐ and low‐fecundity populations of the brown planthopper (BPH) Nilaparvata lugens Stål (Hemiptera: Delphacidae) divergently selected for fecundity. We also tested whether the population fecundity can be predicted by a few SNPs. Seven genes (ACE, fizzy, HMGCR, LpR, Sxl, Vg and VgR) were inspected for SNPs in N. lugens, which is a serious insect pest of rice. By direct sequencing of the complementary DNA and promoter sequences of these candidate genes, 1033 SNPs were discovered within high‐ and low‐fecundity BPH populations. A panel of 121 candidate SNPs were selected and genotyped in 215 individuals from 2 laboratory populations (HFP and LFP) and 3 field populations (GZP, SGP and ZSP). Prior to association tests, population structure and linkage disequilibrium (LD) among the 3 field populations were analysed. The association results showed that 7 SNPs were significantly associated with population fecundity in BPH. These significant SNPs were used for constructing general liner models with stepwise regression. The best predictive model was composed of 2 SNPs (ACE_‐862 and VgR_‐816) with very good fitting degree. We found that 29% of the phenotypic variation in fecundity could be accounted for by only two markers. Using two laboratory populations and a complete independent field population, the predictive accuracy was 84.35–92.39%. The predictive model provides an efficient molecular method to predict BPH fecundity of field populations and provides novel insights for insect population management.