Molecular analysis of human leukocyte antigen class I and class II allele frequencies and haplotype distribution in Pakistani population

AIM:

Distribution of HLA class I and II alleles and haplotype was studied in Pakistani population and compared with the data reported for Caucasoid, Africans, Orientals and Arab populations.

MATERIALS AND METHODS:

HLA class I and II polymorphisms in 1000 unrelated Pakistani individuals was studied using sequence-specific primers and polymerase chain reaction and assay.

RESULTS:

The most frequent class I alleles observed were A*02, B*35 and CW*07, with frequencies of 19.2, 13.7 and 20%, respectively. Fifteen distinct HLA-DRB1 alleles and eight HLA-DQB1 alleles were recognized. The most frequently observed DRB1 alleles which represented more than 60% of the subjects were DRB1 *03, *07, *11 and *15. The rare DRB1 alleles detected in this study were HLADRB1 *08 and *09, having frequencies of 0.9 and 1.7%, respectively. In addition, at DRB1-DQB1 loci there were 179 different haplotypes and 285 unique genotypes and the most common haplotype was DRB1*15-DQB1*06 which represented 17% of the total DRB1-DQB1 haplotypes. In our population, haplotype A*33-B*58-Cw*03 comprised 2.8% of the total class I haplotypes observed. This haplotype was seen only in the oriental populations and has not been reported in the African or European Caucasoid.

CONCLUSION:

Our study showed a close similarity of HLA class I and II alleles with that of European Caucasoid and Orientals. In Pakistani population, two rare loci and three haplotypes were identified, whereas haplotypes characteristic of Caucasians, Africans and Orientals were also found, suggesting an admixture of different races due to migration to and from this region.     

Hydraulic limitation on maximum height of Pinus strobus trees in northern Minnesota,USA

Key message

Using comparisons within and between trees, the authors show evidence for hydraulic limitation of tree height in a humid-climate species that is far from the global maximum tree height.

Abstract

We measured water status and two indicators of drought stress as a function of height within the canopies of four tall (32–35 m) eastern white pines (Pinus strobus) at an old-growth site in northern Minnesota, USA. Pre-dawn and midday xylem pressure potential measured on terminal shoots (Ψ _shoot), needle length, and foliar carbon isotope discrimination (δ ¹³C) all showed within-canopy gradients consistent with increasing drought stress with height. Midday Ψ _shoot near tree tops was ?1.8 MPa, close to values associated with stomatal closure for other temperate conifers. Pre-dawn Ψ _shoot decreased with height at >2× the gradient in gravitational potential. δ ¹³C was strongly correlated with height and weakly correlated with light. Needles were 15–25 % shorter at canopy top compared to the bottom of the canopy. Midday Ψ _shoot and needle length showed significant differences in regression model coefficients from tree to tree. The patterns are consistent with hydraulic constraints on height growth of white pine at this site.     

Haptoglobin gene promoter polymorphism and haplotypes are unique in different populations

We investigated the distribution of haptoglobin (HP) alleles and haplotypes among Africans (Ghanaians), Europeans (from South Africa), and Chinese. HP*1F was present only in Africans and Europeans, whereas HP*del was unique to Chinese. Six base substitutions at the promoter region were population specific. Only 3 out of 18 haplotypes were shared among the populations. A probable application of HP in human population genetics appears legitimate.     

Efficient haplotype block recognition of very long and dense genetic sequences

Background

The new sequencing technologies enable to scan very long and dense genetic sequences, obtaining datasets of genetic markers that are an order of magnitude larger than previously available. Such genetic sequences are characterized by common alleles interspersed with multiple rarer alleles. This situation has renewed the interest for the identification of haplotypes carrying the rare risk alleles. However, large scale explorations of the linkage-disequilibrium (LD) pattern to identify haplotype blocks are not easy to perform, because traditional algorithms have at least Θ(n ²) time and memory complexity.

Results

We derived three incremental optimizations of the widely used haplotype block recognition algorithm proposed by Gabriel et al. in 2002. Our most efficient solution, called MIG ⁺⁺, has only Θ(n) memory complexity and, on a genome-wide scale, it omits >80% of the calculations, which makes it an order of magnitude faster than the original algorithm. Differently from the existing software, the MIG ⁺⁺ analyzes the LD between SNPs at any distance, avoiding restrictions on the maximal block length. The haplotype block partition of the entire HapMap II CEPH dataset was obtained in 457 hours. By replacing the standard likelihood-based D ^′ variance estimator with an approximated estimator, the runtime was further improved. While producing a coarser partition, the approximate method allowed to obtain the full-genome haplotype block partition of the entire 1000 Genomes Project CEPH dataset in 44 hours, with no restrictions on allele frequency or long-range correlations. These experiments showed that LD-based haplotype blocks can span more than one million base-pairs in both HapMap II and 1000 Genomes datasets. An application to the North American Rheumatoid Arthritis Consortium (NARAC) dataset shows how the MIG ⁺⁺ can support genome-wide haplotype association studies.

Conclusions

The MIG ⁺⁺ enables to perform LD-based haplotype block recognition on genetic sequences of any length and density. In the new generation sequencing era, this can help identify haplotypes that carry rare variants of interest. The low computational requirements open the possibility to include the haplotype block structure into genome-wide association scans, downstream analyses, and visual interfaces for online genome browsers.     

Targeted Routine Antenatal Anti-D Prophylaxis in the Prevention of RhD Immunisation - Outcome of a New Antenatal Screening and Prevention Program

Objective

To estimate the incidence of RhD immunisation after implementation of first trimester non-invasive fetal RHD screening to select only RhD negative women carrying RHD positive fetuses for routine antenatal anti-D prophylaxis (RAADP).

Materials and Methods

We present a population-based prospective observational cohort study with historic controls including all maternity care centres and delivery hospitals in the Stockholm region, Sweden. All RhD negative pregnant women were screened for fetal RHD genotype in the first trimester of pregnancy. Anti-D immunoglobulin (250–300 µg) was administered intramuscularly in gestational week 28–30 to participants with RHD positive fetuses. Main outcome measure was the incidence of RhD immunisation developing during or after pregnancy.

Results

During the study period 9380 RhD negative women gave birth in Stockholm. Non-invasive fetal RHD genotyping using cell-free fetal DNA in maternal plasma was performed in 8374 pregnancies of which 5104 (61%) were RHD positive and 3270 (39%) RHD negative. In 4590 pregnancies with an RHD positive test the women received antenatal anti-D prophylaxis. The incidence of RhD immunisation in the study cohort was 0.26 percent (24/9380) (95% CI 0.15–0.36%) compared to 0.46 percent (86/18546) (95% CI 0.37 to 0.56%) in the reference cohort. The risk ratio (RR) for sensitisation was 0.55 (95% CI 0.35 to 0.87) and the risk reduction was statistically significant (p = 0.009). The absolute risk difference was 0.20 percent, corresponding to a number needed to treat (NNT) of 500.

Conclusions

Using first trimester non-invasive antenatal screening for fetal RHD to target routine antenatal anti-D prophylaxis selectively to RhD negative women with RHD positive fetuses significantly reduces the incidence of new RhD immunisation. The risk reduction is comparable to that reported in studies evaluating the outcome of non selective RAADP to all RhD negative women. The cost-effectiveness of this targeted approach remains to be studied.     

Effects of NO 3 − -N on the growth and salinity tolerance of Tamarix laxa Willd

The influence of NO ₃ ^? -N on growth and osmotic adjustment was studied in Tamarix laxa Willd., a halophyte with salt glands on its twigs. Seedlings of T. laxa Willd. were exposed to 1 mM (control) or 300 mM NaCl, with 0.05, 1 or 10 mM NO ₃ ^? -N for 24 days. The relative growth rate of seedlings at 300 mM NaCl was lower than that of control plants at all NO ₃ ^? -N levels, but the concentrations of organic N and total N in the twigs did not differ between the two NaCl treatments. Increasing NO ₃ ^? supply under 300 mM NaCl improved the growth of T. laxa, indicating that NO ₃ ^? played positive roles in improving salt resistance of the plant. The twigs of T. laxa Willd. accumulated mainly inorganic ions, especially Na⁺ and Cl^?, to lower osmotic potential (Ψs): the contributions of Na⁺ and Cl^? to Ψs were estimated at 31% and 27% respectively, at the highest levels of supply of both NaCl and NO ₃ ^? -N. The estimated contribution of NO ₃ ^? -N to Ψs was as high as 20% in the twigs in these conditions, indicating that NO ₃ ^? was also involved in osmotic adjustment in the twigs. Furthermore, increases in tissue NO ₃ ^? were accompanied by decreases in tissue Cl^? and proline under 300 mM NaCl. The estimated contribution of proline to Ψs declined as with NO ₃ ^? -N supply increased from 1 to 10 mM, while the contributions of nitrate to Ψs were enhanced under 300 mM NaCl. This suggested that higher accumulation of nitrate in the vacuole alleviated the effects of salinity stress on the plant by balancing the osmotic potential. In conclusion, NO ₃ ^? -N played both nutritional and osmotic roles in T. laxa Willd. in saline conditions.     

Investigation of Caucasian rheumatoid arthritis susceptibility loci in African patients with the same disease

Introduction

The largest genetic risk to develop rheumatoid arthritis (RA) arises from a group of alleles of the HLA DRB1 locus ('shared epitope', SE). Over 30 non-HLA single nucleotide polymorphisms (SNPs) predisposing to disease have been identified in Caucasians, but they have never been investigated in West/Central Africa. We previously reported a lower prevalence of the SE in RA patients in Cameroon compared to European patients and aimed in the present study to investigate the contribution of Caucasian non-HLA RA SNPs to disease susceptibility in Black Africans.

Methods

RA cases and controls from Cameroon were genotyped for Caucasian RA susceptibility SNPs using Sequenom MassArray technology. Genotype data were also available for 5024 UK cases and 4281 UK controls and for 119 Yoruba individuals in Ibadan, Nigeria (YRI, HapMap). A Caucasian aggregate genetic-risk score (GRS) was calculated as the sum of the weighted risk-allele counts.

Results

After genotyping quality control procedures were performed, data on 28 Caucasian non-HLA susceptibility SNPs were available in 43 Cameroonian RA cases and 44 controls. The minor allele frequencies (MAF) were tightly correlated between Cameroonian controls and YRI individuals (correlation coefficient 93.8%, p = 1.7E-13), and they were pooled together. There was no correlation between MAF of UK and African controls; 13 markers differed by more than 20%. The MAF for markers at PTPN22, IL2RA, FCGR2A and IL2/IL21 was below 2% in Africans. The GRS showed a strong association with RA in the UK. However, the GRS did not predict RA in Africans (OR = 0.71, 95% CI 0.29 - 1.74, p = 0.456). Random sampling from the UK cohort showed that this difference in association is unlikely to be explained by small sample size or chance, but is statistically significant with p<0.001.

Conclusions

The MAFs of non-HLA Caucasian RA susceptibility SNPs are different between Caucasians and Africans, and several polymorphisms are barely detectable in West/Central Africa. The genetic risk of developing RA conferred by a set of 28 Caucasian susceptibility SNPs is significantly different between the UK and Africa with p<0.001. Taken together, these observations strengthen the hypothesis that the genetic architecture of RA susceptibility is different in different ethnic backgrounds.     

RHD positive haplotypes in D negative Europeans

Background  

Blood group genotyping is increasingly utilized for prenatal diagnosis and after recent transfusions, but still lacks the specificity of serology. In whites, the presence of antigen D is predicted, if two or more properly selected RHD-specific polymorphism are detected. This prediction must fail, if an antigen D negative RHD positive allele is encountered. Excluding RHDψ and Cde ^S frequent only in individuals of African descent, most of these alleles are unknown and the population frequency of any such allele has not been determined.     

Bias in effect size of systemic lupus erythematosus susceptibility loci across Europe: a case-control study

Introduction

We aimed to investigate whether the effect size of the systemic lupus erythematosus (SLE) risk alleles varies across European subpopulations.

Methods

European SLE patients (n = 1,742) and ethnically matched healthy controls (n = 2,101) were recruited at 17 centres from 10 different countries. Only individuals with self-reported ancestry from the country of origin were included. In addition, participants were genotyped for top ancestry informative markers and for 25 SLE associated SNPs. The results were used to compare effect sizes between the Central Eureopan and Southern European subgroups.

Results

Twenty of the 25 SNPs showed independent association with SLE, These SNPs showed a significant bias to larger effect sizes in the Southern subgroup, with 15/20 showing this trend (P = 0.019) and a larger mean odds ratio of the 20 SNPs (1.46 vs. 1.34, P = 0.02) as well as a larger difference in the number of risk alleles (2.06 vs. 1.63, P = 0.027) between SLE patients and controls than for Central Europeans. This bias was reflected in a very significant difference in the cumulative genetic risk score (4.31 vs. 3.48, P = 1.8 × 10^-32). Effect size bias was accompanied by a lower number of SLE risk alleles in the Southern subjects, both patients and controls, the difference being more marked between the controls (P = 1.1 × 10^-8) than between the Southern and Central European patients (P = 0.016). Seven of these SNPs showed significant allele frequency clines.

Conclusion

Our findings showed a bias to larger effect sizes of SLE loci in the Southern Europeans relative to the Central Europeans together with clines of SLE risk allele frequencies. These results indicate the need to study risk allele clines and the implications of the polygenic model of inheritance in SLE.     

Patterns of linkage disequilibrium and haplotype distribution in disease candidate genes

Background

The adequacy of association studies for complex diseases depends critically on the existence of linkage disequilibrium (LD) between functional alleles and surrounding SNP markers.

Results

We examined the patterns of LD and haplotype distribution in eight candidate genes for osteoporosis and/or obesity using 31 SNPs in 1,873 subjects. These eight genes are apolipoprotein E (APOE), type I collagen α1 (COL1A1), estrogen receptor-α (ER-α), leptin receptor (LEPR), parathyroid hormone (PTH)/PTH-related peptide receptor type 1 (PTHR1), transforming growth factor-β1 (TGF-β1), uncoupling protein 3 (UCP3), and vitamin D (1,25-dihydroxyvitamin D₃) receptor (VDR). Yin yang haplotypes, two high-frequency haplotypes composed of completely mismatching SNP alleles, were examined. To quantify LD patterns, two common measures of LD, D' and r², were calculated for the SNPs within the genes. The haplotype distribution varied in the different genes. Yin yang haplotypes were observed only in PTHR1 and UCP3. D' ranged from 0.020 to 1.000 with the average of 0.475, whereas the average r² was 0.158 (ranging from 0.000 to 0.883). A decay of LD was observed as the intermarker distance increased, however, there was a great difference in LD characteristics of different genes or even in different regions within gene.

Conclusion

The differences in haplotype distributions and LD patterns among the genes underscore the importance of characterizing genomic regions of interest prior to association studies.     

Deployment of ACT antimalarials for treatment of malaria: challenges and opportunities

Background

In late 2002, the health authorities of Mozambique implemented sulphadoxine-pyrimethamine (SP)/amodiaquine (AQ) as first-line treatment against uncomplicated falciparum malaria. In 2004, this has been altered to SP/artesunate in line with WHO recommendations of using Artemisinin Combination Therapies (ACTs), despite the fact that all the neighbouring countries have abandoned SP-drug combinations due to high levels of SP drug resistance. In the study area, one year prior to the change to SP/AQ, SP alone was used to treat uncomplicated malaria cases. The study described here investigated the immediate impact of the change to SP on the frequency of SP and CQ resistance-related haplotypes in the Plasmodium falciparum genes Pfdhfr, Pfdhps and Pfcrt before and a year after the introduction of SP.

Methods

Samples were collected during two cross sectional surveys in early 2002 and 2003 involving 796 and 692 children one year or older and adults randomly selected living in Maciana, an area located in Manhiça district, Southern Mozambique. Out of these, 171 and 173 P. falciparum positive samples were randomly selected to measure the frequency of resistance- related haplotypes in Pfdhfr, Pfdhps and Pfcrt based on results obtained by nested PCR followed by sequence-specific oligonucleotide probe (SSOP)-ELISA.

Results

The frequency of the SP-resistance associated Pfdhps double mutant (SGEAA) haplotype increased significantly from 14% to 35% (P < 0.001), while the triple mutant Pfdhfr haplotype (CIRNI) remained high and only changed marginally from 46% to 53% (P = 0.405) after one year with SP as first-line treatment in the study area. Conversely, the combined Pfdhfr/Pfdhps quintuple mutant haplotype increased from 8% to 26% (P = 0.005). The frequency of the chloroquine resistance associated Pfcrt -CVIET haplotype was above 90% in both years.

Conclusion

These retrospective findings add to the general concern on the lifespan of the combination of SP/artesunate in Mozambique. The high frequency of quintuple Pfdhfr/Pfdhps haplotypes found here as early as 2002 most likely cause ideal conditions for the development of artesunate tolerance in the P. falciparum populations and may even endanger the sensitivity to the second-line drug, Coartem^®.     

Haplotypic analysis of the TNF locus by association efficiency and entropy

Background

To understand the causal basis of TNF associations with disease, it is necessary to understand the haplotypic structure of this locus. We genotyped 12 single-nucleotide polymorphisms (SNPs) distributed over 4.3 kilobases in 296 healthy, unrelated Gambian and Malawian adults. We generated 592 high-quality haplotypes by integrating family- and population-based reconstruction methods.

Results

We found 32 different haplotypes, of which 13 were shared between the two populations. Both populations were haplotypically diverse (gene diversity = 0.80, Gambia; 0.85, Malawi) and significantly differentiated (p < 10^-5 by exact test). More than a quarter of marker pairs showed evidence of intragenic recombination (29% Gambia; 27% Malawi). We applied two new methods of analyzing haplotypic data: association efficiency analysis (AEA), which describes the ability of each SNP to detect every other SNP in a case-control scenario; and the entropy maximization method (EMM), which selects the subset of SNPs that most effectively dissects the underlying haplotypic structure. AEA revealed that many SNPs in TNF are poor markers of each other. The EMM showed that 8 of 12 SNPs (Gambia) and 7 of 12 SNPs (Malawi) are required to describe 95% of the haplotypic diversity.

Conclusions

The TNF locus in the Gambian and Malawi sample is haplotypically diverse and has a rich history of intragenic recombination. As a consequence, a large proportion of TNF SNPs must be typed to detect a disease-modifying SNP at this locus. The most informative subset of SNPs to genotype differs between the two populations.     

Ancient pathogen-driven adaptation triggers increased susceptibility to non-celiac wheat sensitivity in present-day European populations

Background

Non-celiac wheat sensitivity is an emerging wheat-related syndrome showing peak prevalence in Western populations. Recent studies hypothesize that new gliadin alleles introduced in the human diet by replacement of ancient wheat with modern varieties can prompt immune responses mediated by the CXCR3-chemokine axis potentially underlying such pathogenic inflammation. This cultural shift may also explain disease epidemiology, having turned European-specific adaptive alleles previously targeted by natural selection into disadvantageous ones.

Methods

To explore this evolutionary scenario, we performed ultra-deep sequencing of genes pivotal in the CXCR3-inflammatory pathway on individuals diagnosed for non-celiac wheat sensitivity and we applied anthropological evolutionary genetics methods to sequence data from worldwide populations to investigate the genetic legacy of natural selection on these loci.

Results

Our results indicate that balancing selection has maintained two divergent CXCL10/CXCL11 haplotypes in Europeans, one responsible for boosting inflammatory reactions and another for encoding moderate chemokine expression.

Conclusions

This led to considerably higher occurrence of the former haplotype in Western people than in Africans and East Asians, suggesting that they might be more prone to side effects related to the consumption of modern wheat varieties. Accordingly, this study contributed to shed new light on some of the mechanisms potentially involved in the disease etiology and on the evolutionary bases of its present-day epidemiological patterns. Moreover, overrepresentation of disease homozygotes for the dis-adaptive haplotype plausibly accounts for their even more enhanced CXCR3-axis expression and for their further increase in disease risk, representing a promising finding to be validated by larger follow-up studies.

     

Conservation of the COP9/signalosome in budding yeast

Background

Variation at the PPARG locus may influence susceptibility to type 2 diabetes and related traits. The Pro12Ala polymorphism may modulate receptor activity and is associated with protection from type 2 diabetes. However, there have been inconsistent reports of its association with obesity. The silent C1431T polymorphism has not been as extensively studied, but the rare T allele has also been inconsistently linked to increases in weight. Both rare alleles are in linkage disequilibrium and the independent associations of these two polymorphisms have not been addressed.

Results

We have genotyped a large population with type 2 diabetes (n = 1107), two populations of non-diabetics from Glasgow (n = 186) and Dundee (n = 254) and also a healthy group undergoing physical training (n = 148) and investigated the association of genotype with body mass index. This analysis has demonstrated that the Ala12 and T1431 alleles are present together in approximately 70% of the carriers. By considering the other 30% of individuals with haplotypes that only carry one of these polymorphisms, we have demonstrated that the Ala12 allele is consistently associated with a lower BMI, whilst the T1431 allele is consistently associated with higher BMI.

Conclusion

This study has therefore revealed an opposing interaction of these polymorphisms, which may help to explain previous inconsistencies in the association of PPARG polymorphisms and body weight.     

TLR7 single-nucleotide polymorphisms in the 3' untranslated region and intron 2 independently contribute to systemic lupus erythematosus in Japanese women: a case-control association study

Introduction

The Toll-like receptor 7 (TLR7) gene, encoded on human chromosome Xp22.3, is crucial for type I interferon production. A recent multicenter study in East Asian populations, comprising Chinese, Korean and Japanese participants, identified an association of a TLR7 single-nucleotide polymorphism (SNP) located in the 3' untranslated region (3' UTR), rs3853839, with systemic lupus erythematosus (SLE), especially in males, although some difference was observed among the tested populations. To test whether additional polymorphisms contribute to SLE in Japanese, we systematically analyzed the association of TLR7 with SLE in a Japanese female population.

Methods

A case-control association study was conducted on eight tag SNPs in the TLR7 region, including rs3853839, in 344 Japanese females with SLE and 274 healthy female controls.

Results

In addition to rs3853839, two SNPs in intron 2, rs179019 and rs179010, which were in moderate linkage disequilibrium with each other (r ² = 0.53), showed an association with SLE (rs179019: P = 0.016, odds ratio (OR) 2.02, 95% confidence interval (95% CI) 1.15 to 3.54; rs179010: P = 0.018, OR 1.75, 95% CI 1.10 to 2.80 (both under the recessive model)). Conditional logistic regression analysis revealed that the association of the intronic SNPs and the 3' UTR SNP remained significant after we adjusted them for each other. When only the patients and controls carrying the risk genotypes at the 3' UTR SNPpositionwere analyzed, the risk of SLE was significantly increased when the individuals also carried the risk genotypes at both of the intronic SNPs (P = 0.0043, OR 2.45, 95% CI 1.31 to 4.60). Furthermore, the haplotype containing the intronic risk alleles in addition to the 3' UTR risk allele was associated with SLE under the recessive model (P = 0.016, OR 2.37, 95% CI 1.17 to 4.80), but other haplotypes were not associated with SLE.

Conclusions

The TLR7 intronic SNPs rs179019 and rs179010 are associated with SLE independently of the 3' UTR SNP rs3853839 in Japanese women. Our findings support a role of TLR7 in predisposition for SLE in Asian populations.     

On the relationship between an Asian haplotype on chromosome 6 that reduces androstenone levels in boars and the differential expression of SULT2A1 in the testis

Background

Androstenone is one of the major compounds responsible for boar taint, a pronounced urine-like odor produced when cooking boar meat. Several studies have identified quantitative trait loci (QTL) for androstenone level on Sus scrofa chromosome (SSC) 6. For one of the candidate genes in the region SULT2A1, a difference in expression levels in the testis has been shown at the protein and RNA level.

Results

Haplotypes were predicted for the QTL region and their effects were estimated showing that haplotype 1 was consistently related with a lower level, and haplotype 2 with a higher level of androstenone. A recombinant haplotype allowed us to narrow down the QTL region from 3.75 Mbp to 1.94 Mbp. An RNA-seq analysis of the liver and testis revealed six genes that were differentially expressed between homozygotes of haplotypes 1 and 2. Genomic sequences of these differentially expressed genes were checked for variations within potential regulatory regions. We identified one variant located within a CpG island that could affect expression of SULT2A1 gene. An allele-specific expression analysis in the testis did not show differential expression between the alleles of SULT2A1 located on the different haplotypes in heterozygous animals. However a synonymous mutation C166T (SSC6: 49,117,861 bp in Sscrofa 10.2; C/T) was identified within the exon 2 of SULT2A1 for which the haplotype 2 only had the C allele which was higher expressed than the T allele, indicating haplotype-independent allelic-imbalanced expression between the two alleles. A phylogenetic analysis for the 1.94 Mbp region revealed that haplotype 1, associated with low androstenone level, originated from Asia.

Conclusions

Differential expression could be observed for six genes by RNA-seq analysis. No difference in the ratio of C:T expression of SULT2A1 for the haplotypes was found by the allele-specific expression analysis, however, a difference in expression between the C over T allele was found for a variation within SULT2A1, showing that the difference in androstenone levels between the haplotypes is not caused by the SNP in exon 2.     

STK39 Polymorphism Is Associated with Essential Hypertension: A Systematic Review and Meta-Analysis

Background

A recent genome-wide association study identified STK39as a candidate gene for blood pressure (BP) in Europeans. Subsequently, several studies have attempted to replicate the association across different ethnic populations. However, the results have been inconsistent.

Objective and Methods

We performed a meta-analysis to elucidate the association between the STK39 rs3754777 polymorphism (or proxy) and hypertension. Published literature from PubMed and Embase databases were retrieved and pooled odds ratio (OR) with 95% confidence interval (CI) was calculated using fixed- or random-effects model.

Results

Using appropriate inclusion/exclusion criteria, we identified 10 studies that included 21, 863 hypertensive cases and 24, 480 controls from different ethnicities. The meta-analysis showed a significant association of STK39 rs3754777 variant with hypertension (OR = 1.10, 95%CI = 1.06–1.15, p = 7.95×10⁻⁶). Further subgroup analysis by ethnicity suggested that the association was significant in Europeans (OR = 1.08, 95% CI = 1.03–1.14, p = 0.002) and in East Asians (OR = 1.16, 95%CI = 1.07–1.25, p = 4.34×10⁻⁴), but not in Africans (OR = 1.01, 95%CI 0.80–1.27, p = 0.932). We further confirmed the positive association by sensitivity analysis. No publication bias was detected (Begg’s test, p = 0.721; Egger’s test, p = 0.744).

Conclusions

The present meta-analysis confirms the significant association of STK39 polymorphism with susceptibility to hypertension in Europeans and East Asians. Future studies should include gene–gene and gene–environment interactions to investigate the identified association.     

The STR polymorphisms in intron 8 may provide information about the molecular evolution of RH haplotypes

We identified simple-sequence repeat polymorphisms in intron 8 of the RHD and RHCE genes, both of which contained the 5-bp repeat unit (AAAAT)n. We analyzed the polymorphisms of this short tandem repeat (STR) in 104 Japanese RhD-positive and 124 RhD-negative (87 RHD gene negative and 37 nonfunctional RHD gene positive) donors by the polymerase chain reaction (PCR) and subsequent typing by electrophoresis and silver staining. We found five alleles (10, 11, 12, 13, and 14 repeats) in the RHD gene and four (7, 8, 9, and 10 repeats) in the RHCE gene. The Rh phenotypes were closely associated with polymorphisms of the STR. The Ce allele had 12 repeats in the RHD gene and 9 repeats in the RHCE gene at high frequency. The cE allele frequently had 10–12 repeats in the RHD gene and 10 repeats in the RHCE gene. The 10 repeats in the RHCE gene were identified exclusively in the 87 RHD gene-negative donors and 9 repeats were identified only in those with the RhC antigen. These results indicate that both haplotypes of dce and dcE arose from single RHD gene deletion and recombination events, respectively. In the 37 RhD-negative donors with a nonfunctional RHD gene, 12 repeats in the RHD gene and 9 repeats in the RHCE gene were frequently observed. Thus, the RhD-negative with a nonfunctional RHD gene combination might have arisen from the DCe haplotype via a mutation that abolished RHD gene expression. These findings suggest that the STR polymorphisms might shed light upon the molecular evolution of RH haplotypes. Received: 30 November 1998 / Accepted: 8 February 1999     

Balancing Selection on a Regulatory Region Exhibiting Ancient Variation That Predates Human–Neandertal Divergence

Ancient population structure shaping contemporary genetic variation has been recently appreciated and has important implications regarding our understanding of the structure of modern human genomes. We identified a ∼36-kb DNA segment in the human genome that displays an ancient substructure. The variation at this locus exists primarily as two highly divergent haplogroups. One of these haplogroups (the NE1 haplogroup) aligns with the Neandertal haplotype and contains a 4.6-kb deletion polymorphism in perfect linkage disequilibrium with 12 single nucleotide polymorphisms (SNPs) across diverse populations. The other haplogroup, which does not contain the 4.6-kb deletion, aligns with the chimpanzee haplotype and is likely ancestral. Africans have higher overall pairwise differences with the Neandertal haplotype than Eurasians do for this NE1 locus (p<10⁻¹⁵). Moreover, the nucleotide diversity at this locus is higher in Eurasians than in Africans. These results mimic signatures of recent Neandertal admixture contributing to this locus. However, an in-depth assessment of the variation in this region across multiple populations reveals that African NE1 haplotypes, albeit rare, harbor more sequence variation than NE1 haplotypes found in Europeans, indicating an ancient African origin of this haplogroup and refuting recent Neandertal admixture. Population genetic analyses of the SNPs within each of these haplogroups, along with genome-wide comparisons revealed significant F_ST (p = 0.00003) and positive Tajima''s D (p = 0.00285) statistics, pointing to non-neutral evolution of this locus. The NE1 locus harbors no protein-coding genes, but contains transcribed sequences as well as sequences with putative regulatory function based on bioinformatic predictions and in vitro experiments. We postulate that the variation observed at this locus predates Human–Neandertal divergence and is evolving under balancing selection, especially among European populations.