Molecular evolution of 5' flanking regions of 87 candidate genes for atherosclerotic cardiovascular disease

Molecular cloning and characterization of NPR1 gene from Arachis hypogaea

The NPR1 gene was an important regulator for a plant disease resistance. The cDNA of NPR1 gene was cloned from peanut cultivar Ri Hua 1 by rapid amplification of cDNA ends-polymerase chain reaction (RACE-PCR). The full length cDNA of Arachis hypogaea NPR1 consisted of 2,078 base pairs with a 1,446 bp open-reading frame encoding 481 amino acids. The predicted NPR1 contained the highly conserved functional domains (BTB/POZ domain from M1 to D116), protein-protein interaction domains (three ankyrin repeats from K158 to L186; N187 to L217 and R221 to D250) and one NPR1-like domain (C262 to S469). The DNA sequence of the NPR1 gene was 2,332 or 2,223 bp. Both two sequences contained three introns and four exons. The NPR1 transcripts were expressed mainly in roots and leaves, while fewer signals were detected in the stems. Amount of the NPR1 transcript was significantly increased 1 h after salicylic acid challenge and was eventually 5.3 times greater than that in the control group. Both the DNA sequence and the coding sequence were obtained from eight cultivars and nine wild species of Arachis. Maximum likelihood analyses of d N/d S ratios for 25 sequences from different species showed that different selection pressures may have acted on different branches.

Controlling type-I error of the McDonald-Kreitman test in genomewide scans for selection on noncoding DNA

Departures from the assumption of homogenously interdigitated neutral and putatively selected sites in the McDonald-Kreitman test can lead to false rejections of the neutral model in the presence of intermediate levels of recombination. This problem is exacerbated by small sample sizes, nonequilibrium demography, recombination rate variation, and in comparisons involving more recently diverged species. I propose that establishing significance levels by coalescent simulation with recombination can improve the fidelity of the test in genomewide scans for selection on noncoding DNA.

Molecular population genetics of PCSK9: a signature of recent positive selection

OBJECTIVE

Proprotein convertase subtilisin-like kexin type 9 (PCSK9) is a newly discovered serine protease that plays a key role in regulating plasma low density lipoprotein (LDL) cholesterol levels. Both rare mutations and common variants in the coding regions of PCSK9 affect LDL cholesterol levels and coronary heart disease risk, as well as response to lipid-lowering therapy.

METHODS

The authors characterized the patterns of variation at the PCSK9 locus in African-Americans and European-Americans using resequenced data from the SeattleSNPs database (pga.gs.washington.edu). The authors performed a test of population differentiation and the long-range haplotype test to detect signatures of recent position selection on PCSK9.

RESULTS

A significantly high FST (a measure of population differentiation) between African-Americans and European-Americans was noted for single nucleotide polymorphism (SNP) rs505151 (FST=0.309). The long-range haplotype test was suggestive of non-neutral evolution of two SNPs, rs505151 and rs562556 that are associated with elevated LDL cholesterol levels ('gain-of-function' mutations). The modes of selection were different between African-Americans and European-Americans. The authors observed signals of recent positive selection on the ancestral allele of nonsynonymous SNP rs505151 (E670G, P=0.0227 and P=0.0001 in theoretical and empirical distribution, respectively) and the derived allele of nonsynonymous SNP rs562556 (I474V, P=0.0227 and 0.0001) in African-Americans, whereas in European-Americans the ancestral allele of SNP rs562556 (P=0.1320 and 0.0370) appeared to be under positive selection.

CONCLUSION

The authors' findings suggest that evolutionary dynamics may underlie the 'gain-of-function' mutations in PCSK9 that are associated with higher LDL cholesterol levels.

Evidence for positive selection in the C-terminal domain of the cholesterol metabolism gene PCSK9 based on phylogenetic analysis in 14 primate species

Background

Cholesterol homeostasis is maintained through finely tuned mechanisms regulating intestinal absorption, hepatic biosynthesis and secretion as well as plasma clearance. Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a secreted enzyme of the serine protease family that reduces cellular uptake of plasma low-density lipoprotein (LDL) cholesterol by promoting LDL receptor (LDL-R) degradation. Species-specific positive selection has been noted in the LDLR promoter, leading to differential expression of LDLR among primates. Whether PCSK9 experienced significant selective pressure to maintain a functional relationship with its target protein, LDL-R, is unknown.

Methodology/Principal Findings

We compiled the sequences of the coding regions of PCSK9 from 14 primate species in the clade of Hominoids, Old World monkeys and New World monkeys. To detect selective pressure at the protein level, the ratios of nonsynonymous/synonymous substitution rate (d_N/d_S) under different evolutionary models were calculated across the phylogeny of PCSK9. Maximum likelihood analyses of d_N/d_S ratios for the aligned coding region sequences among 14 primate species indicated that PCSK9 was subject to a strong functional constraint (i.e., purifying selection). However, positive selection was noted in the functional carboxyl-terminal (C-terminal) domain in many branches across the phylogeny, especially in the lineage leading to the orangutan. Furthermore, at least five positively selected amino acids were detected in this lineage using the branch-site model A. In a sliding-window analysis, several d_N/d_S peaks in the C-terminal domain in both the human and the orangutan branches were noted.

Conclusions

These results suggest that among primates, differential selective pressure has shaped evolutionary patterns in the functional domains of PCSK9, an important regulator of cholesterol homeostasis.

Comparative and evolutionary pharmacogenetics of ABCB1: complex signatures of positive selection on coding and regulatory regions

BACKGROUND

As a major mediator in the complex interplay between humans and the xenobiotic environment, the ABCBI transporter gene is an obvious candidate for comparative and evolutionary pharmacogenetic studies. It has been recently reported that common variants in its coding region, which are variously associated with drug response and disease susceptibility, may have conferred differential selective sweep in various populations. Fully profiling the alletic architecture and explicitly interrogating the natural selection at ABCBI are needed to understand its evolutionary population genetics.

METHODS AND RESULTS

Using a comprehensive single nucleotide polymorphism variants in coding and regulatory regions, as well as comparable genotype data from the Environmental Genome Project, we systematically characterized the extent and length of linkage disequilibrium throughout the ABCBI locus in three major ethnic populations (African, European, and Chinese). We observed pronounced signals of recent positive selection on the derived alleles of three common single nucleotide polymorphisms coding regions: e12/1236T, e21/2677T, and e26/3435T in the Chinese, as well as on extended haplotype homozygosity were also observed for two potentially functional common variants in the 5'f/-4489G (rs17149810) in the Chinese and 5'f/-693T (rs3213619) in the Africans, respectively, which may have shaped the phylogenetically inferred star-like haplotype structure of the 5'flanking region.

CONCLUSION

Our finding reveal complex signatures of natural selection on both coding and regulatory regions of the human ABCBI gene, point to potential functional relevance of its regulatory variants, and suggest that evolutionary dynamics and transcriptional regulation may underline the phenotypic variation in xenobiotic disposition and varying predisposition to complex in which xenobiotics play a role.

Differentially expressed genes in adipose tissues of high body weight-selected (obese) and unselected (lean) mouse lines

Recently, quantitative trait loci (QTLs) for body weight and obesity have been mapped in an intercross population between the high body weight-selected mouse line DU6i and the inbred line DBA/2. Most QTLs were highly significant, but had small effects only. Under the hypothesis that small-effect QTLs might result from changes in gene activity, our strategy to identify candidate genes for the observed effects was directed towards the identification of differentially expressed genes. Therefore, here we compare the transcription profile of about 11 000 genes in epididymal fat tissues of males of two high body weight-selected (DU6 and DU6i) and two unselected mouse lines (DUKs and DBA/2). For the hybridisation of GeneChips, we used pooled samples of 20 individual mice. By pair-wise comparisons between selected and unselected mouse lines, a set of 77 genes was identified representing genes whose level of expression differed between obese and lean mouse strains. According to the functional classification of genes, 69 differentially expressed genes were involved in regulatory and metabolic pathways, cell division, cell stability, or immune response, and thus might have an effect on body weight and fat accumulation. 14 out of these genes, occur in QTL regions for body weight or abdominal fat weight. Further analyses are necessary to discriminate between genes directly causing QTL effects and indirectly regulated differentially expressed genes.

Patterns of population differentiation of candidate genes for cardiovascular disease

BACKGROUND

The basis for ethnic differences in cardiovascular disease (CVD) susceptibility is not fully understood. We investigated patterns of population differentiation (FST) of a set of genes in etiologic pathways of CVD among 3 ethnic groups: Yoruba in Nigeria (YRI), Utah residents with European ancestry (CEU), and Han Chinese (CHB) + Japanese (JPT). We identified 37 pathways implicated in CVD based on the PANTHER classification and 416 genes in these pathways were further studied; these genes belonged to 6 biological processes (apoptosis, blood circulation and gas exchange, blood clotting, homeostasis, immune response, and lipoprotein metabolism). Genotype data were obtained from the HapMap database.

RESULTS

We calculated FST for 15,559 common SNPs (minor allele frequency > or = 0.10 in at least one population) in genes that co-segregated among the populations, as well as an average-weighted FST for each gene. SNPs were classified as putatively functional (non-synonymous and untranslated regions) or non-functional (intronic and synonymous sites). Mean FST values for common putatively functional variants were significantly higher than FST values for nonfunctional variants. A significant variation in FST was also seen based on biological processes; the processes of 'apoptosis' and 'lipoprotein metabolism' showed an excess of genes with high FST. Thus, putative functional SNPs in genes in etiologic pathways for CVD show greater population differentiation than non-functional SNPs and a significant variance of FST values was noted among pairwise population comparisons for different biological processes.

CONCLUSION

These results suggest a possible basis for varying susceptibility to CVD among ethnic groups.