A Revised Evolutionary History of the CYP1A Subfamily: Gene Duplication, Gene Conversion, and Positive Selection

Members of cytochrome P450 subfamily 1A (CYP1As) are involved in detoxification and bioactivation of common environmental pollutants. Understanding the functional evolution of these genes is essential to predicting and interpreting species differences in sensitivity to toxicity caused by such chemicals. The CYP1A gene subfamily comprises a single ancestral representative in most fish species and two paralogs in higher vertebrates, including birds and mammals. Phylogenetic analysis of complete coding sequences suggests that mammalian and bird paralog pairs (CYP1A1/2 and CYP1A4/5, respectively) are the result of independent gene duplication events. However, comparison of vertebrate genome sequences revealed that CYP1A genes lie within an extended region of conserved fine-scale synteny, suggesting that avian and mammalian CYP1A paralogs share a common genomic history. Algorithms designed to detect recombination between nucleotide sequences indicate that gene conversion has homogenized most of the length of the chicken CYP1A genes, as well as the 5′ end of mammalian CYP1As. Together, these data indicate that avian and mammalian CYP1A paralog pairs resulted from a single gene duplication event and that extensive gene conversion is responsible for the exceptionally high degree of sequence similarity between CYP1A4 and CYP1A5. Elevated nonsynonymous/synonymous substitution ratios within a putatively unconverted stretch of ∼250 bp suggests that positive selection may have reduced the effective rate of gene conversion in this region, which contains two substrate recognition sites. This work significantly alters our understanding of functional evolution in the CYP1A subfamily, suggesting that gene conversion and positive selection have been the dominant processes of sequence evolution. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Yves Van de Peer]     

Common Genetic Variants of the Human Steroid 21-Hydroxylase Gene (CYP21A2) Are Related to Differences in Circulating Hormone Levels

Purpose

Systematic evaluation of the potential relationship between the common genetic variants of CYP21A2 and hormone levels.

Methods

The relationships of CYP21A2 intron 2 polymorphisms and haplotypes with diverse baseline and stimulated blood hormone levels were studied in 106 subjects with non-functioning adrenal incidentaloma (NFAI). The rationale for using NFAI subjects is dual: i) their baseline hormone profiles do not differ from those of healthy subjects and ii) hormone levels after stimulation tests are available.

Results

The carriers (N = 27) of a well-defined CYP21A2 haplotype cluster (c5) had significantly elevated levels of cortisol (p = 0.0110), and 17-hydroxyprogesterone (p = 0.0001) after ACTH stimulation, and 11-deoxycortisol after metyrapone administration (p = 0.0017), but the hormone values were in normal ranges. In addition, the carriers (N = 33) of the C allele of the rs6462 polymorphism had a higher baseline aldosterone level (p = 0.0006). The prevalence of these genetic variants of CYP21A2 did not differ between NFAI and healthy subjects.

Conclusions

The common CYP21A2 variants presumably exert the same effect on hormone levels in the healthy and disease-affected populations. Therefore, they may contribute to complex diseases such as some cardiovascular diseases, and may influence the genotype-phenotype correlation in patients with congenital adrenal hyperplasia (CAH) including the individual need for hormone substitution.     

A General Model of Negative Frequency Dependent Selection Explains Global Patterns of Human ABO Polymorphism

The ABO locus in humans is characterized by elevated heterozygosity and very similar allele frequencies among populations scattered across the globe. Using knowledge of ABO protein function, we generated a simple model of asymmetric negative frequency dependent selection and genetic drift to explain the maintenance of ABO polymorphism and its loss in human populations. In our models, regardless of the strength of selection, models with large effective population sizes result in ABO allele frequencies that closely match those observed in most continental populations. Populations must be moderately small to fall out of equilibrium and lose either the A or B allele (N_e ≤ 50) and much smaller (N_e ≤ 25) for the complete loss of diversity, which nearly always involved the fixation of the O allele. A pattern of low heterozygosity at the ABO locus where loss of polymorphism occurs in our model is consistent with small populations, such as Native American populations. This study provides a general evolutionary model to explain the observed global patterns of polymorphism at the ABO locus and the pattern of allele loss in small populations. Moreover, these results inform the range of population sizes associated with the recent human colonization of the Americas.     

Analysis of CYP21A1P and the duplicated CYP21A2 genes

The RCCX module on chromosome 6p21.3 has 3 possible forms: monomodular, bimodular, and trimodular. Chromosomes with 4 RCCX modules are very rare. In the monomodule, most of the CYP21A1P genes do not exist. However, haplotypes of the RCCX module with more than one CYP21A2 gene were observed. Obviously, the gene located downstream of the XA gene can possibly include the CYP21A2 as well as the CYP21A1P gene.     

Analysis of the CYP21A1P pseudogene: indication of mutational diversity and CYP21A2-like and duplicated CYP21A2 genes

The CYP21A1P gene downstream of the XA gene, carrying 15 deteriorated mutations, is a nonfunctional pseudogene that shares 98% nucleotide sequence homology with CYP21A2 located on chromosome 6p21.3. However, these mutations in the CYP21A1P gene are not totally involved in each individual. From our analysis of 100 healthy ethnic Chinese (i.e., Taiwanese) (n = 200 chromosomes) using the polymerase chain reaction (PCR) products combined with an amplification-created restriction site (ACRS) method and DNA sequencing, we found that approximately 10% of CYP21A1P alleles (n = 195 chromosomes) presented the CYP21A2 sequence; frequencies of P30, V281, Q318, and R356 in that locus were approximately 24%, 21%, 11%, and 34%, respectively, and approximately 90% of the CYP21A1P alleles had 15 mutated loci. In addition, approximately 2.5% (n = 5 chromosomes) showed four haplotypes of the 3.7-kb TaqI-produced fragment of the CYP21A2-like gene and one duplicated CYP21A2 gene. We conclude that the pseudogene of the CYP21A1P mutation presents diverse variants. Moreover, the existence of the CYP21A2-like gene is more abundant than that of the duplicated CYP21A2 gene downstream of the XA gene and could not be distinguished from the CYP21A2–TNXB gene; thus, it may be misdiagnosed by previously established methods for congenital adrenal hyperplasia caused by a 21-hydroxylase deficiency.     

A Genetic Polymorphism in Coumarin 7-Hydroxylation: Sequence of the Human CYP2A Gnes and Identification of Variant CYP2A6 Alleles

A group of human cytochrome P450 genes encompassing the CYP2A, CYP2B, and CYP2F subfamilies were cloned and assembled into a 350-kb contig localized on the long arm of chromosome 19. Three complete CYP2A genes—CYP2A6, CYP2A7, and CYP2A13—plus two pseudogenes truncated after exon 5, were identified and sequenced. A variant CYP2A6 allele that differed from the corresponding CYP2A6 and CYP2A7 cDNAs previously sequenced was found and was designated CYP2A6ν2. Sequence differences in the CYP2A6ν2 gene are restricted to regions encompassing exons 3, 6, and 8, which bear sequence relatedness with the corresponding exons of the CYP2A7 gene, located downstream and centromeric of CYP2A6ν2, suggesting recent gene-conversion events. The sequencing of all the CYP2A genes allowed the design of a PCR diagnostic test for the normal CYP2A6 allele, the CYP2A6ν2 allele, and a variant—designated CYP2A6ν1—that encodes an enzyme with a single inactivating amino acid change. These variant alleles were found in individuals who were deficient in their ability to metabolize the CYP2A6 probe drug coumarin. The allelic frequencies of CYP2A6ν1 and CYP2A6ν2 differed significantly between Caucasian, Asian, and African-American populations. These studies establish the existence of a new cytochrome P450 genetic polymorphism.     

Differential Gene Expression Patterns of EBV Infected EBNA-3A Positive and Negative Human B Lymphocytes

The genome of Epstein-Barr virus (EBV) encodes 86 proteins, but only a limited set is expressed in EBV–growth transformed B cells, termed lymphoblastoid cell lines (LCLs). These cells proliferate via the concerted action of EBV nuclear antigens (EBNAs) and latent membrane proteins (LMPs), some of which are rate limiting to establish a stable homeostasis of growth promoting and anti-apoptotic activities. We show here that EBV mutants, which lack the EBNA-3A gene, are impaired but can still initiate cell cycle entry and proliferation of primary human B cells in contrast to an EBNA-2 deficient mutant virus. Surprisingly, and in contrast to previous reports, these viral mutants are attenuated in growth transformation assays but give rise to permanently growing EBNA-3A negative B cell lines which exhibit reduced proliferation rates and elevated levels of apoptosis. Expression profiles of EBNA-3A deficient LCLs are characterized by 129 down-regulated and 167 up-regulated genes, which are significantly enriched for genes involved in apoptotic processes or cell cycle progression like the tumor suppressor gene p16/INK4A, or might contribute to essential steps of the viral life cycle in the infected host. In addition, EBNA-3A cellular target genes remarkably overlap with previously identified targets of EBNA-2. This study comprises the first genome wide expression profiles of EBNA-3A target genes generated within the complex network of viral proteins of the growth transformed B cell and permits a more detailed understanding of EBNA-3A''s function and contribution to viral pathogenesis.     

CYP21和CYP21P基因启动子序列对绿色荧光蛋白基因表达的影响

特异性扩增CYP21基因和CYP21P基因启动子区域－770bp--1bp片段,去除pEGFP-N1载体中的CMV启动子,构建含CYP21基因启动子的pEGFP-N1载体（pCYP21)和CYP21P基因启动子的pEGFP-N1载体（pCYP21P),分别将上述两种构建载体、野生型pEGFP-N1(阳性对照）质粒及阴性对照转染入肾上腺皮质来源的Y1细胞系中,用倒置荧光显微镜,以及激光共聚焦显微镜等方法观测绿色荧光蛋白的表达。转染后,在荧光倒置显微镜下首次发现Y1细胞中出现绿色荧光蛋白的时间阳性对照为3小时,pCYP21为7小时,pCYP21P与阳性对照（未转染任何载体的Y1细胞）始终未观测到绿色荧光蛋白。阳性对照和pCYP21的绿色荧光蛋白表达强于pCYP21,pCYP21P与阴性对照始终未观测到绿色荧光蛋白。阳性对照和pCYP21的绿色荧光蛋白在胞核中的荧光强度高于胞浆。上述结果进一步表明,含有CYP21和CYP21P两种基因启动子的GFP质粒在Y1细胞中表达存在显著差异。     

T细胞阳性选择和阴性选择的发生比较

胸腺细胞的阳性选择与阴性选择在T细胞的发育成熟过程中起着关键作用。自骨髓进入胸腺的胸腺细胞在增殖后大多数(约95％)在阳性选择或阴性选择中凋亡,其中大约90％在阳性选择中因胸腺细胞的T细胞受体(TCR)不能识别自身MHC而凋亡,约5％在阴性选择中因胸腺细胞的TCR能够识别自身肽而凋亡。只有5％的胸腺细胞存活下来,成为具有自身限制性和自身耐受特性的成熟T细胞离开胸腺。过去传统的观点认为,胸腺细胞的阳性选择与阴性选择是发生于不同解剖位置、不同发育时期的两个截然分开的过程。近年来,随着转基因动物技术、体外器官培养等一些新技术的应用,对阳性选择与阴性选择的发生机制有了新的认识。     

人MIF基因启动子多态性与结核病易感性的关系

为了研究巨噬细胞移动抑制因子(MIF)与结核病易感性的关系,分析了肺结核病人和正常人MIF基因启动子多态性。DNA测序分析的结果显示,正常人MIF基因启动子-794区CATT序列重复次数存在明显的差异,5、6和7个拷贝的几率分别是22.2%、44.4%和33.3%。7个肺结核病人中6人有6个拷贝,1人有7个拷贝,出现的频率分别为85.6%和14.3%。对比分析提示,肺结核病人MIF基因启动子区CATT低拷贝数(低于或等于6个)出现的频率明显高于正常人(85.6%比66.6%),可能同结核病的致病相关。     

Functional Polymorphisms in the CYP2C19 Gene Contribute to Digestive System Cancer Risk: Evidence from 11,042 Subjects

Background

CYP2C19 belongs to the cytochrome P450 superfamily of enzymes involved in activating and detoxifying many carcinogens and endogenous compounds, which has attracted considerable attention as a candidate gene for digestive system cancer. CYP2C19 has two main point mutation sites (CYP2C19*2, CYP2C19*3) leading to poor metabolizer (PM) phenotype. In the past decade, the relationship between CYP2C19 polymorphism and digestive system cancer has been reported in various ethnic groups; however, these studies have yielded contradictory results.

Methods

To clarify this inconsistency, we performed this meta-analysis. Databases including Pubmed, EMBASE, Web of Science and China National Knowledge Infrastructure (CNKI) were searched to find relevant studies. Odds ratios (ORs) with 95% confidence intervals (CIs) were used to assess the strength of association.

Results

In total, 18 studies with 4,414 cases and 6,628 controls were included. Overall, significantly elevated digestive system cancer risk was associated CYP2C19 PM with OR of 1.66 (95%CI: 1.31–2.10, P<10⁻⁵) when all studies were pooled into the meta-analysis. There was strong evidence of heterogeneity (P = 0.006), which largely disappeared after stratification by cancer type. In the stratified analyses according to cancer type, ethnicity, control source and sample size, significantly increased risks were found.

Conclusions

In summary, our meta-analysis suggested that the PM phenotype caused by the variation on CYP2C19 gene is associated with increased risk of digestive system cancer, especially in East Asians.     

Involvement of the Modifier Gene of a Human Mendelian Disorder in a Negative Selection Process

Background

Identification of modifier genes and characterization of their effects represent major challenges in human genetics. SAA1 is one of the few modifiers identified in humans: this gene influences the risk of renal amyloidosis (RA) in patients with familial Mediterranean fever (FMF), a Mendelian autoinflammatory disorder associated with mutations in MEFV. Indeed, the SAA1 α homozygous genotype and the p.Met694Val homozygous genotype at the MEFV locus are two main risk factors for RA.

Methodology/Principal Findings

Here, we investigated Armenian FMF patients and controls from two neighboring countries: Armenia, where RA is frequent (24%), and Karabakh, where RA is rare (2.5%). Sequencing of MEFV revealed similar frequencies of p.Met694Val homozygotes in the two groups of patients. However, a major deficit of SAA1 α homozygotes was found among Karabakhian patients (4%) as compared to Armenian patients (24%) (p = 5.10⁻⁵). Most importantly, we observed deviations from Hardy-Weinberg equilibrium (HWE) in the two groups of patients, and unexpectedly, in opposite directions, whereas, in the two control populations, genotype distributions at this locus were similar and complied with (HWE).

Conclusions/Significance

The excess of SAA1α homozygotes among Armenian patients could be explained by the recruitment of patients with severe phenotypes. In contrast, a population-based study revealed that the deficit of α/α among Karabakhian patients would result from a negative selection against carriers of this genotype. This study, which provides new insights into the role of SAA1 in the pathophysiology of FMF, represents the first example of deviations from HWE and selection involving the modifier gene of a Mendelian disorder.     

Bioengineered Nisin A Derivatives with Enhanced Activity against Both Gram Positive and Gram Negative Pathogens

Nisin is a bacteriocin widely utilized in more than 50 countries as a safe and natural antibacterial food preservative. It is the most extensively studied bacteriocin, having undergone decades of bioengineering with a view to improving function and physicochemical properties. The discovery of novel nisin variants with enhanced activity against clinical and foodborne pathogens has recently been described. We screened a randomized bank of nisin A producers and identified a variant with a serine to glycine change at position 29 (S29G), with enhanced efficacy against S. aureus SA113. Using a site-saturation mutagenesis approach we generated three more derivatives (S29A, S29D and S29E) with enhanced activity against a range of Gram positive drug resistant clinical, veterinary and food pathogens. In addition, a number of the nisin S29 derivatives displayed superior antimicrobial activity to nisin A when assessed against a range of Gram negative food-associated pathogens, including E. coli, Salmonella enterica serovar Typhimurium and Cronobacter sakazakii. This is the first report of derivatives of nisin, or indeed any lantibiotic, with enhanced antimicrobial activity against both Gram positive and Gram negative bacteria.