A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Postural changes in women with chronic pelvic pain: a case control study

Background: Chronic pelvic pain (CPP) is a lower abdominal pain lasting at least 6 months, occurring continuously or intermittently and not associated exclusively with menstruation or intercourse. Although the musculoskeletal system has been found to be involved in CPP, few studies have assessed the contribution of posture in women with CPP. We aimed to determine if the frequency of postural changes was higher in women with CPP than healthy subjects. Methods: A case-control study included 108 women with CPP of more than 6 months' duration (CPP group) who consecutively attended at the Hospital of the University of Sao Paulo and 48 healthy female volunteers (control group). Postural assessment was noninvasive and performed in the standing position, with the reference points of Kendall used as normal parameters. Factors associated with CPP were assessed by logistic regression analysis. Results: Logistic regression showed that the independent factors associated with CPP were postural changes in the cervical spine (OR 4.1; 95% CI 1.6-10.7; p < 0.01) and scapulae (OR 2.9; 95% CI 1.1-7.6; p < 0.05). Conclusion: Musculoskeletal changes were associated with CPP in 34% of women. These findings suggest that a more detailed assessment of women with CPP is necessary for better diagnosis and for more effective treatment.

Cost-effectiveness analysis of PCR for the rapid diagnosis of pulmonary tuberculosis

Background: Tuberculosis is one of the most prominent health problems in the world, causing 1.75 million deaths each year. Rapid clinical diagnosis is important in patients who have comorbidities such as Human Immunodeficiency Virus (HIV) infection. Direct microscopy has low sensitivity and culture takes 3 to 6 weeks [1-3]. Therefore, new tools for TB diagnosis are necessary, especially in health settings with a high prevalence of HIV/TB co-infection. Methods: In a public reference TB/HIV hospital in Brazil, we compared the cost-effectiveness of diagnostic strategies for diagnosis of pulmonary TB: Acid fast bacilli smear microscopy by Ziehl-Neelsen staining (AFB smear) plus culture and AFB smear plus colorimetric test (PCR dot-blot). From May 2003 to May 2004, sputum was collected consecutively from PTB suspects attending the Parthenon Reference Hospital. Sputum samples were examined by AFB smear, culture, and PCR dot-blot. The gold standard was a positive culture combined with the definition of clinical PTB. Cost analysis included health services and patient costs. Results: The AFB smear plus PCR dot-blot require the lowest laboratory investment for equipment (US$ 20,000). The total screening costs are 3.8 times for AFB smear plus culture versus for AFB smear plus PCR dot blot costs (US$ 5,635,760 versus US$ 1,498, 660). Costs per correctly diagnosed case were US$ 50,773 and US$ 13,749 for AFB smear plus culture and AFB smear plus PCR dot-blot, respectively. AFB smear plus PCR dot-blot was more cost-effective than AFB smear plus culture, when the cost of treating all correctly diagnosed cases was considered. The cost of returning patients, which are not treated due to a negative result, to the health service, was higher in AFB smear plus culture than for AFB smear plus PCR dot-blot, US$ 374,778,045 and US$ 110,849,055, respectively. Conclusion: AFB smear associated with PCR dot-blot associated has the potential to be a cost-effective tool in the fight against PTB for patients attended in the TB/HIV reference hospital.

Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR

Background: Considering the broad variation in the expression of housekeeping genes among tissues and experimental situations, studies using quantitative RT-PCR require strict definition of adequate endogenous controls. For glioblastoma, the most common type of tumor in the central nervous system, there was no previous report regarding this issue. Results: Here we show that amongst seven frequently used housekeeping genes TBP and HPRT1 are adequate references for glioblastoma gene expression analysis. Evaluation of the expression levels of 12 target genes utilizing different endogenous controls revealed that the normalization method applied might introduce errors in the estimation of relative quantities. Genes presenting expression levels which do not significantly differ between tumor and normal tissues can be considered either increased or decreased if unsuitable reference genes are applied. Most importantly, genes showing significant differences in expression levels between tumor and normal tissues can be missed. We also demonstrated that the Holliday Junction Recognizing Protein, a novel DNA repair protein over expressed in lung cancer, is extremely over-expressed in glioblastoma, with a median change of about 134 fold. Conclusion: Altogether, our data show the relevance of previous validation of candidate control genes for each experimental model and indicate TBP plus HPRT1 as suitable references for studies on glioblastoma gene expression.

Prevalence of non-communicable diseases in Brazilian children: follow-up at school age of two Brazilian birth cohorts of the 1990's

Background: Few cohort studies have been conducted in low and middle-income countries to investigate non-communicable diseases among school-aged children. This article aims to describe the methodology of two birth cohorts, started in 1994 in Ribeirao Preto (RP), a more developed city, and in 1997/98 in Sao Luis (SL), a less developed town. Methods: Prevalences of some non-communicable diseases during the first follow-up of these cohorts were estimated and compared. Data on singleton live births were obtained at birth (2858 in RP and 2443 in SL). The follow-up at school age was conducted in RP in 2004/05, when the children were 9-11 years old and in SL in 2005/06, when the children were 7-9 years old. Follow-up rates were 68.7% in RP (790 included) and 72.7% in SL (673 participants). The groups of low (<2500 g) and high (>= 4250 g) birthweight were oversampled and estimates were corrected by weighting. Results: In the more developed city there was a higher percentage of non-nutritive sucking habits (69.1% vs 47.9%), lifetime bottle use (89.6% vs 68.3%), higher prevalence of primary headache in the last 15 days (27.9% vs 13.0%), higher positive skin tests for allergens (44.3% vs 25.3%) and higher prevalence of overweight (18.2% vs 3.6%), obesity (9.5% vs 1.8%) and hypertension (10.9% vs 4.6%). In the less developed city there was a larger percentage of children with below average cognitive function (28.9% vs 12.2%), mental health problems (47.4% vs 38.4%), depression (21.6% vs 6.0%) and underweight (5.8% vs 3.6%). There was no difference in the prevalence of bruxism, recurrent abdominal pain, asthma and bronchial hyperresponsiveness between cities. Conclusions: Some non-communicable diseases were highly prevalent, especially in the more developed city. Some high rates suggest that the burden of non-communicable diseases will be high in the future, especially mental health problems.

Allopatric speciation in ticks: genetic and reproductive divergence between geographic strains of Rhipicephalus (Boophilus) microplus

Background: The cattle tick, Rhipicephalus (Boophilus) microplus, economically impact cattle industry in tropical and subtropical regions of the world. The morphological and genetic differences among R. microplus strains have been documented in the literature, suggesting that biogeographical and ecological separation may have resulted in boophilid ticks from America/Africa and those from Australia being different species. To test the hypothesis of the presence of different boophilid species, herein we performed a series of experiments to characterize the reproductive performance of crosses between R. microplus from Australia, Africa and America and the genetic diversity of strains from Australia, Asia, Africa and America. Results: The results showed that the crosses between Australian and Argentinean or Mozambican strains of boophilid ticks are infertile while crosses between Argentinean and Mozambican strains are fertile. These results showed that tick strains from Africa (Mozambique) and America (Argentina) are the same species, while ticks from Australia may actually represent a separate species. The genetic analysis of mitochondrial 12S and 16S rDNA and microsatellite loci were not conclusive when taken separately, but provided evidence that Australian tick strains were genetically different from Asian, African and American strains. Conclusion: The results reported herein support the hypothesis that at least two different species share the name R. microplus. These species could be redefined as R. microplus (Canestrini, 1887) (for American and African strains) and probably the old R. australis Fuller, 1899 (for Australian strains), which needs to be redescribed. However, experiments with a larger number of tick strains from different geographic locations are needed to corroborate these results.

T cells, adhesion molecules and modulation of apoptosis in visceral leishmaniasis glomerulonephritis

Background: Immune complex deposition is the accepted mechanism of pathogenesis of VL glomerulopathy however other immune elements may participate. Further in the present study, no difference was seen between immunoglobulin and C3b deposit intensity in glomeruli between infected and non-infected dogs thus T cells, adhesion molecules and parameters of proliferation and apoptosis were analysed in dogs with naturally acquired VL from an endemic area. The dog is the most important domestic reservoir of the protozoa Leishmania (L.) chagasi that causes visceral leishmaniasis (VL). The similarity of VL manifestation in humans and dogs renders the study of canine VL nephropathy of interest with regard to human pathology. Methods: From 55 dogs with VL and 8 control non-infected dogs from an endemic area, kidney samples were analyzed by immunohistochemistry for immunoglobulin and C3b deposits, staining for CD4+ and CD8+ T cells, ICAM-1, P-selectin and quantified using morphometry. Besides proliferation marker Ki-67, apoptosis markers M30 and TUNEL staining, and related cytokines TNF-alpha, IL-1 alpha were searched and quantified. Results: We observed similar IgG, IgM and IgA and C3b deposit intensity in dogs with VL and non-infected control dogs. However we detected the Leishmania antigen in cells in glomeruli in 54, CD4+ T cells in the glomeruli of 44, and CD8+ T cells in 17 of a total of 55 dogs with VL. Leishmania antigen was absent and T cells were absent/scarse in eight non-infected control dogs. CD 4+ T cells predominate in proliferative patterns of glomerulonephritis, however the presence of CD4+ and CD8+ T cells were not different in intensity in different patterns of glomerulonephritis. The expression of ICAM-1 and P-selectin was significantly greater in the glomeruli of infected dogs than in control dogs. In all patterns of glomerulonephritis the expression of ICAM-1 ranged from minimum to moderately severe and P-selectin from absent to severe. In the control animals the expression of these molecules ranged from absent to medium intensity. It was not observed any correlation between severity of the disease and these markers. There was a correlation between the number of Leishmania antigen positive cells and CD4+ T cells, and between the number of CD4+ T cells and CD8+ T cells. In dogs presenting different histopathological patterns of glomerulonephritis, parameters of proliferation and apoptosis were studied. Ki-67, a proliferative marker, was not detected locally, but fewer apoptotic cells and lower TNF-alpha expression were seen in infected animals than in non-infected controls. Conclusion: Immunopathogenic mechanisms of VL glomerulonephritis are complex and data in the present study suggest no clear participation of immunoglobulin and C3b deposits in these dogs but the possible migration of CD4+ T cells into the glomeruli, participation of adhesion molecules, and diminished apoptosis of cells contributing to determine the proliferative pattern of glomerulonephritis in VL.

Can type of school be used as an alternative indicator of socioeconomic status in dental caries studies? A cross-sectional study

Background: Despite the importance of collecting individual data of socioeconomic status (SES) in epidemiological oral health surveys with children, this procedure relies on the parents as respondents. Therefore, type of school (public or private schools) could be used as an alternative indicator of SES, instead of collecting data individually. The aim of this study was to evaluate the use of the variable type of school as an indicator of socioeconomic status as a substitute of individual data in an epidemiological survey about dental caries in Brazilian preschool children. Methods: This study followed a cross-sectional design, with a random sample of 411 preschool children aged 1 to 5 years, representative of Catalao, Brazil. A calibrated examiner evaluated the prevalence of dental caries and parents or guardians provided information about several individual socioeconomic indicators by means of a semi-structured questionnaire. A multilevel approach was used to investigate the association among individual socioeconomic variables, as well as the type of school, and the outcome. Results: When all significant variables in the univariate analysis were used in the multiple model, only mother's schooling and household income (individual socioeconomic variables) presented significant associations with presence of dental caries, and the type of school was not significantly associated. However, when the type of school was used alone, children of public school presented significantly higher prevalence of dental caries than those enrolled in private schools. Conclusions: The type of school used as an alternative indicator for socioeconomic status is a feasible predictor for caries experience in epidemiological dental caries studies involving preschool children in Brazilian context.

RNA interference-mediated knockdown of CD49e (alpha 5 integrin chain) in human thymic epithelial cells modulates the expression of multiple genes and decreases thymocyte adhesion

Background: The thymus is a central lymphoid organ, in which bone marrow-derived T cell precursors undergo a complex process of maturation. Developing thymocytes interact with thymic microenvironment in a defined spatial order. A component of thymic microenvironment, the thymic epithelial cells, is crucial for the maturation of T-lymphocytes through cell-cell contact, cell matrix interactions and secretory of cytokines/chemokines. There is evidence that extracellular matrix molecules play a fundamental role in guiding differentiating thymocytes in both cortical and medullary regions of the thymic lobules. The interaction between the integrin alpha 5 beta 1 (CD49e/CD29; VLA-5) and fibronectin is relevant for thymocyte adhesion and migration within the thymic tissue. Our previous results have shown that adhesion of thymocytes to cultured TEC line is enhanced in the presence of fibronectin, and can be blocked with anti-VLA-5 antibody. Results: Herein, we studied the role of CD49e expressed by the human thymic epithelium. For this purpose we knocked down the CD49e by means of RNA interference. This procedure resulted in the modulation of more than 100 genes, some of them coding for other proteins also involved in adhesion of thymocytes; others related to signaling pathways triggered after integrin activation, or even involved in the control of F-actin stress fiber formation. Functionally, we demonstrated that disruption of VLA-5 in human TEC by CD49e-siRNA-induced gene knockdown decreased the ability of TEC to promote thymocyte adhesion. Such a decrease comprised all CD4/CD8-defined thymocyte subsets. Conclusion: Conceptually, our findings unravel the complexity of gene regulation, as regards key genes involved in the heterocellular cell adhesion between developing thymocytes and the major component of the thymic microenvironment, an interaction that is a mandatory event for proper intrathymic T cell differentiation.

Trends in socioeconomic inequalities in cancer mortality in Barcelona: 1992-2003

Background: The objective of this study was to assess trends in cancer mortality by educational level in Barcelona from 1992 to 2003. Methods: The study population comprised Barcelona inhabitants aged 20 years or older. Data on cancer deaths were supplied by the system of information on mortality. Educational level was obtained from the municipal census. Age-standardized rates by educational level were calculated. We also fitted Poisson regression models to estimate the relative index of inequality (RII) and the Slope Index of Inequalities (SII). All were calculated for each sex and period (1992-1994, 1995-1997, 1998-2000, and 2001-2003). Results: Cancer mortality was higher in men and women with lower educational level throughout the study period. Less-schooled men had higher mortality by stomach, mouth and pharynx, oesophagus, larynx and lung cancer. In women, there were educational inequalities for cervix uteri, liver and colon cancer. Inequalities of overall and specific types of cancer mortality remained stable in Barcelona; although a slight reduction was observed for some cancers. Conclusion: This study has identified those cancer types presenting the greatest inequalities between men and women in recent years and shown that in Barcelona there is a stable trend in inequalities in the burden of cancer.

Trends in dermatomyositis- and polymyositis-related mortality in the state of Sao Paulo, Brazil, 1985-2007: multiple cause-of-death analysis

Background: Dermatomyositis (DM) and polymyositis (PM) are rare systemic autoimmune rheumatic diseases with high fatality rates. There have been few population-based mortality studies of dermatomyositis and polymyositis in the world, and none have been conducted in Brazil. The objective of the present study was to employ multiple-cause of-death methodology in the analysis of trends in mortality related to dermatomyositis and polymyositis in the state of Sao Paulo, Brazil, between 1985 and 2007. Methods: We analyzed mortality data from the Sao Paulo State Data Analysis System, selecting all death certificates on which DM or PM was listed as a cause of death. The variables sex, age and underlying, associated or total mentions of causes of death were studied using mortality rates, proportions and historical trends. Statistical analysis were performed by chi-square and H Kruskal-Wallis tests, variance analysis and linear regression. A p value less than 0.05 was regarded as significant. Results: Over a 23-year period, there were 318 DM-related deaths and 316 PM-related deaths. Overall, DM/PM was designated as an underlying cause in 55.2% and as an associated cause in 44.8%; among 634 total deaths females accounted for 71.5%. During the study period, age-and gender-adjusted DM mortality rates did not change significantly, although PM as an underlying cause and total mentions of PM trended lower (p < 0.05). The mean ages at death were 47.76 +/- 20.81 years for DM and 54.24 +/- 17.94 years for PM (p = 0.0003). For DM/PM, respectively, as underlying causes, the principal associated causes of death were as follows: pneumonia (in 43.8%/33.5%); respiratory failure (in 34.4%/32.3%); interstitial pulmonary diseases and other pulmonary conditions (in 28.9%/17.6%); and septicemia (in 22.8%/15.9%). For DM/PM, respectively, as associated causes, the following were the principal underlying causes of death: respiratory disorders (in 28.3%/26.0%); circulatory disorders (in 17.4%/20.5%); neoplasms (in 16.7%/13.7%); infectious and parasitic diseases (in 11.6%/9.6%); and gastrointestinal disorders (in 8.0%/4.8%). Of the 318 DM-related deaths, 36 involved neoplasms, compared with 20 of the 316 PM-related deaths (p = 0.03). Conclusions: Our study using multiple cause of deaths found that DM/PM were identified as the underlying cause of death in only 55.2% of the deaths, indicating that both diseases were underestimated in the primary mortality statistics. We observed a predominance of deaths in women and in older individuals, as well as a trend toward stability in the mortality rates. We have confirmed that the risk of death is greater when either disease is accompanied by neoplasm, albeit to lesser degree in individuals with PM. The investigation of the underlying and associated causes of death related to DM/PM broaden the knowledge of the natural history of both diseases and could help integrate mortality data for use in the evaluation of control measures for DM/PM.

Genetic parameters for postweaning traits in Braunvieh cattle

Genetic parameters for traits related to postweaning growth in Braunvieh cattle, reared under tropical and sub-tropical conditions in Brazil, were studied. Weight traits were weight at 365 days of age (W365, N = 4055), at 450 days (W450, N = 3453), and at 550 days (W550, N = 1946), while weight gains were gain from weaning to 365 days of age (WGW365, N = 3060), from weaning to 450 days (WGW450, N = 2764), from weaning to 550 days (WGW550, N = 1531), from 365 to 550 days of age (WG365550, N = 1528), from 365 to 450 days (WG365450, N = 2401), and from 450 to 550 days (WG450550, N = 1563). A full animal model was used for estimating the variance components, using the MTDFREML software. The dataset contained 18,688 animals with phenotypic measures and 35,188 animals in the relationship matrix. Heritability estimates for postweaning weights decreased with age. For W365, W450 and W550, respectively, the direct heritability estimates were 0.29 +/- 0.061, 0.25 +/- 0.057, 0.16 +/- 0.060, maternal heritability was 0.20 +/- 0.035, 0.18 +/- 0.035, 0.13 +/- 0.052, and total heritability was 0.30, 0.35, 0.26. In this breed, maternal influence was found to be important up to 550 days of age. The greater genetic correlations between weights were observed for weights measured at shorter intervals. A large environmental effect was observed for weight gain between weaning and 550 days; this effect was greater for the gains between 365 and 550 days.

Characterization of new IS elements and studies of their dispersion in two subspecies of Leifsonia xyli

Background: Leifsonia xyli is a xylem-inhabiting bacterial species comprised of two subspecies: L. xyli subsp. xyli (Lxx) and L. xyli subsp. cynodontis (Lxc). Lxx is the causal agent of ratoon stunting disease in sugarcane commercial fields and Lxc colonizes the xylem of several grasses causing either mild or no symptoms of disease. The completely sequenced genome of Lxx provided insights into its biology and pathogenicity. Since IS elements are largely reported as an important source of bacterial genome diversification and nothing is known about their role in chromosome architecture of L. xyli, a comparative analysis of Lxc and Lxx elements was performed. Results: Sample sequencing of Lxc genome and comparative analysis with Lxx complete DNA sequence revealed a variable number of IS transposable elements acting upon genomic diversity. A detailed characterization of Lxc IS elements and a comparative review with IS elements of Lxx are presented. Each genome showed a unique set of elements although related to same IS families when considering features such as similarity among transposases, inverted and direct repeats, and element size. Most of the Lxc and Lxx IS families assigned were reported to maintain transposition at low levels using translation regulatory mechanisms, consistent with our in silico analysis. Some of the IS elements were found associated with rearrangements and specific regions of each genome. Differences were also found in the effect of IS elements upon insertion, although none of the elements were preferentially associated with gene disruption. A survey of transposases among genomes of Actinobacteria showed no correlation between phylogenetic relatedness and distribution of IS families. By using Southern hybridization, we suggested that diversification of Lxc isolates is also mediated by insertion sequences in probably recent events. Conclusion: Collectively our data indicate that transposable elements are involved in genome diversification of Lxc and Lxx. The IS elements were probably acquired after the divergence of the two subspecies and are associated with genome organization and gene contents. In addition to enhancing understanding of IS element dynamics in general, these data will contribute to our ongoing comparative analyses aimed at understanding the biological differences of the Lxc and Lxx.

Group I introns and associated homing endonuclease genes reveals a clinal structure for Porphyra spiralis var. amplifolia (Bangiales, Rhodophyta) along the Eastern coast of South America

Background: Group I introns are found in the nuclear small subunit ribosomal RNA gene (SSU rDNA) of some species of the genus Porphyra (Bangiales, Rhodophyta). Size polymorphisms in group I introns has been interpreted as the result of the degeneration of homing endonuclease genes (HEG) inserted in peripheral loops of intron paired elements. In this study, intron size polymorphisms were characterized for different Porphyra spiralis var. amplifolia (PSA) populations on the Southern Brazilian coast, and were used to infer genetic relationships and genetic structure of these PSA populations, in addition to cox2-3 and rbcL-S regions. Introns of different sizes were tested qualitatively for in vitro self-splicing. Results: Five intron size polymorphisms within 17 haplotypes were obtained from 80 individuals representing eight localities along the distribution of PSA in the Eastern coast of South America. In order to infer genetic structure and genetic relationships of PSA, these polymorphisms and haplotypes were used as markers for pairwise Fst analyses, Mantel's test and median joining network. The five cox2-3 haplotypes and the unique rbcL-S haplotype were used as markers for summary statistics, neutrality tests Tajima's D and Fu's Fs and for median joining network analyses. An event of demographic expansion from a population with low effective number, followed by a pattern of isolation by distance was obtained for PSA populations with the three analyses. In vitro experiments have shown that introns of different lengths were able to self-splice from pre-RNA transcripts. Conclusion: The findings indicated that degenerated HEGs are reminiscent of the presence of a full-length and functional HEG, once fixed for PSA populations. The cline of HEG degeneration determined the pattern of isolation by distance. Analyses with the other markers indicated an event of demographic expansion from a population with low effective number. The different degrees of degeneration of the HEG do not refrain intron self-splicing. To our knowledge, this was the first study to address intraspecific evolutionary history of a nuclear group I intron; to use nuclear, mitochondrial and chloroplast DNA for population level analyses of Porphyra; and intron size polymorphism as a marker for population genetics.

Phylogenetic nomenclature and evolution of mannose-binding lectin (MBL2) haplotypes

Background: Polymorphisms of the mannose-binding lectin gene (MBL2) affect the concentration and functional efficiency of the protein. We recently used haplotype-specific sequencing to identify 23 MBL2 haplotypes, associated with enhanced susceptibility to several diseases. Results: In this work, we applied the same method in 288 and 470 chromosomes from Gabonese and European adults, respectively, and found three new haplotypes in the last group. We propose a phylogenetic nomenclature to standardize MBL2 studies and found two major phylogenetic branches due to six strongly linked polymorphisms associated with high MBL production. They presented high Fst values and were imbedded in regions with high nucleotide diversity and significant Tajima's D values. Compared to others using small sample sizes and unphased genotypic data, we found differences in haplotyping, frequency estimation, Fu and Li's D* and Fst results. Conclusion: Using extensive testing for selective neutrality, we confirmed that stochastic evolutionary factors have had a major role in shaping this polymorphic gene worldwide.