On the molecular mass of the extracellular hemoglobin of Glossoscolex paulistus: Analytical ultracentrifugation reexamination

The giant extracellular hemoglobin of Glossoscolex paulistus (HbGp) is constituted by subunits containing heme groups with molecular masses (M) in the range of 15 to 19 kDa, monomers of 16 kDa (d), and trimers of 51 to 52 kDa (abc) linked by nonheme structures named linkers of 24 to 32 kDa (L). HbGp is homologous to Lumbricus terrestris hemoglobin (HbLt). Several reports propose M of HbLt in the range of 3.6 to 4.4 MDa. Based on subunits M determined by mass spectrometry and assuming HbGp stoichiometry of 12(abcd)₃L₃ (Vinogradov model) plus 144 heme groups, a value of M for HbGp oligomer of 3560 kDa can be predicted. This value is nearly 500 kDa higher than the unique HbGp M value reported in the literature. In the current work, sedimentation velocity analytical ultracentrifugation (AUC) experiments were performed to obtain M for HbGp in oxy and cyano-met forms. s⁰_20,w values of 58.1 ± 0.2 S and 59.6 ± 0.2 S, respectively, for the two oxidation forms were obtained. The ratio between sedimentation and diffusion coefficients supplied values for M of approximately 3600 ± 100 and 3700 ± 100 kDa for oxy and cyano-met HbGp forms, respectively. An independent determination of the partial specific volume, V_bar, for HbGp was performed based on density measurements, providing a value of 0.764 ± 0.008, in excellent agreement with the estimates from SEDFIT software. Our results show total consistency between M obtained by AUC and recent partial characterization by mass spectrometry. Therefore, HbGp possesses M very close to that of HbLt, suggesting an oligomeric assembly in agreement with the Vinogradov model.     

Towards a resolution of the long-standing controversy regarding the molecular mass of extracellular erythrocruorin of the earthworm Lumbricus terrestris

The published molecular mass of erythrocruorin of Lumbricus terrestris and related earthworm species covers a bewildering range of 3.23-4.5 MDa. A critical reexamination reveals that some mass determinations were underestimated and the results do cluster, not at one, but at two values of the molecular mass. One cluster corresponds to approximately 3.6 MDa, as predicted for a stoichiometry of 144 globin and 36 linker chains-the Vinogradov model for the hexagonal bilayer (HBL) assembly of Lumbricus erythrocruorin-and as estimated from the crystal structure of HBL at 5.5 A resolution [Proc. Natl. Acad. Sci. U. S. A. 97 (2000) 7107]. The other cluster corresponds to approximately 4.4 MDa. In addition, a molecular mass of 4.1 MDa, determined by multiangle laser light scattering (MALLS), stands apart of the two clusters, separated from the masses obtained by other methods of molecular mass determination.We propose a stoichiometry of 192 globin and 36 linker chains for the 4.4-MDa molecule. The 36 linkers and 144 out of 192 globin chains are identified with the HBL and the remaining 48 globins are allotted equally to the two halves of the axial cavity above and below the central torus of the structure. The proposed model is supported by the occurrence in some annelid species of erythrocruorin with centrally placed subunits [Biochim. Biophys. Acta 359 (1974) 210], and by the oxidation-dependent shedding of subunits in Lumbricus erythrocruorin. We propose further that the 4.1 MDa determination represents the weight average molecular mass of a population of molecules resulting from a partial dissociation of 4.4-MDa erythrocruorin. This interpretation seems reasonable on the background of the very low protein concentrations ( approximately 100 microg/ml and lower) prevailing at the MALLS experiment.     

Mass distributions of a macromolecular assembly based on electrospray ionization mass spectrometric masses of the constituent subunits

Macromolecular assemblies containing multiple protein subunits and having masses in the megadalton (MDa) range are involved in most of the functions of a living cell. Because of variation in the number and masses of subunits, macromolecular assemblies do not have a unique mass, but rather a mass distribution. The giant extracelular erythrocruorins (Ers), ∼ 3.5 MDa, comprized of at least 180 polypeptide chains, are one of the best characterized assemblies. Three-dimensional reconstructions from cryoelectron microscopic images show them to be hexagonal bilayer complexes of 12 subassemblies, each comprised of 12 globin chains, anchored to a subassembly of 36 nonglobin linker chains. We have calculated the most probable mass distributions forLumbricus andRiftia assemblies and their globin and linker subassemblies, based on theLumbricus Er stoichiometry and using accurate subunit masses obtained by electrospray ionization mass spectrometry. The expected masses ofLumbricus andRiftia Ers are 3.517 MDa and 3.284 MDa, respectively, with a possible variation of ∼ 9% due to the breadth of the mass distributions. TheLumbricus Er mass is in astonishingly good agreement with the mean of 23 known masses, 3.524 ± 0.481 MDa.     

STUDIES ON INVERTEBRATE HEMOGLOBINS (ERYTHROCRUORINS)

1. Two high molecular invertebrate hemoglobins (the erythrocruorins of Lumbricus terrestris and of Nereis virens) as well as the low molecular erythrocruorin of Glycera dibranchiata Ehlers were studied. Their physical chemical properties were compared with those of vertebrate hemoglobin. 2. The hemin of the blood pigment of Glycera dibranchiata Ehlers was shown to be identical with that of vertebrate hemoglobin. 3. The dissociation rates of Glycera and human oxyhemoglobin were measured in the reaction meter of DuBois and t ₅₀ (half time of the reaction) was found to be identical (0.027 second) for the two pigments. The t ₅₀ value for the high molecular Lumbricus erythrocruorin was 0.070 second. 4. The chemical constitution and physical chemical properties of erythrocruorins were compared with those of vertebrate hemoglobin and of hemocyanin.     

Seasonal acclimation in resting metabolism of the skink, Mabuya brevicollis (Reptilia: Scincidae) from southwestern Saudi Arabia

The resting metabolic rate (RMR) of seasonally-acclimated Mabuya brevicollis of various body masses was determined at 20, 25, 30, 35 and 40 °C, using open-flow respirometry. RMR (ml g⁻¹ h⁻¹) decreased with increasing mass at each temperature. RMRs increaProd. Type: FTPsed as temperature increased. The highest and lowest Q₁₀ values were obtained for the temperature ranges 20–25 °C and 30–35 °C for the summer-acclimated lizards. The exponent of mass “b” in the metabolism-body mass relation ranged from 0.41 to 0.61. b values were lower in the autumn and winter-acclimated lizards than in spring and summer-acclimated lizards. Seasonal acclimation effects were evident at all temperatures (20–40 °C) for M. brevicollis. Winter-acclimated skinks had the lowest metabolic rates at different temperatures. The pattern of acclimation exhibited by M. brevicollis may represent a useful adaptation for lizards inhabiting subtropical deserts to promote activity during their active seasons.     

pH-induced quaternary assembly of Vitreoscilla hemoglobin: The monomer exhibits better peroxidase activity

pH-dependent (pH 6.0–8.0) quaternary structural changes of ferric Vitreoscilla hemoglobin (VHb) have been investigated using dynamic light scattering. The VHb exhibits a monomeric state under neutral conditions at pH 7.0, while the protein forms distinct homodimeric species at pH 6.0 and 8.0, respectively. The dissociation constant obtained using the Bio-Layer Interferometry technology indicates that, at pH 7.0, the monomer–monomer dissociation of VHb is about 6-fold or 5-fold higher (K_D = 6.34 μM) compared with that at slightly acidic pH (K_D = 1.05 μM) or slightly alkaline pH (K_D = 1.22 μM). The pH-dependent absorption spectra demonstrate that the heme microenvironment of VHb is sensitive to the changes of pH value. The maximum absorption band of heme group of VHb shifts from 402 nm to 407 nm when pH changes from 6.0 to 8.0. In addition, the fluorescence emission spectra of VHb, taken at excitation wavelength of 295 nm, suggest that the single Trp122 fluorescence quantum yields in VHb are decreased due to the formation of the homodimeric species. However, the circular dichroism spectra data display that the secondary structures of VHb are little affected by pH transitions. The pH-dependent peroxidase activity of VHb was also investigated in this study. The optimum pH for VHb using 2,2′-azino-bis-(3-ethylbenzothiazoline-6-sulfonic acid) as substrate is 7.0, which implies that the monomer state of VHb would exhibit better peroxidase activity than the homodimeric species of VHb at pH 6.0 and 8.0.     

Invasive Bombus terrestris (Hymenoptera: Apidae) parasitized by a flagellate (Euglenozoa: Kinetoplastea) and a neogregarine (Apicomplexa: Neogregarinorida)

The flagellate Crithidia bombi and the neogregarine Apicystis bombi have been found in individuals of Bombus terrestris, a Palaearctic species of bumble bee commercially reared and shipped worldwide for pollination services. B. terrestris has recently entered into the northwestern Patagonia region of Argentina from Chile, where it was introduced in 1998. Prevalence was 21.6% for C. bombi and 3.6% for A. bombi (n = 111). The pathogens were not detected in 441 bumble bees belonging to five of the eight known Argentine native species (Bombus atratus, Bombus morio, Bombus bellicosus, Bombus opifex, Bombus tucumanus) collected elsewhere in the country. Although the absence of natural occurrence of C. bombi and A. bombi in Argentine native bumble bees cannot be ascertained at present due to the limited surveys performed, it is important to report their detection in invasive B. terrestris. The invasion event is relatively recent and the accompanying pathogens are not species specific within the genus Bombus.     

EmrE dimerization depends on membrane environment

The small multi-drug resistant (SMR) transporter EmrE functions as a homodimer. Although the small size of EmrE would seem to make it an ideal model system, it can also make it challenging to work with. As a result, a great deal of controversy has surrounded even such basic questions as the oligomeric state. Here we show that the purified protein is a homodimer in isotropic bicelles with a monomer–dimer equilibrium constant (K_MD^2D) of 0.002–0.009 mol% for both the substrate-free and substrate-bound states. Thus, the dimer is stabilized in bicelles relative to detergent micelles where the K_MD^2D is only 0.8–0.95 mol% (Butler et al. 2004). In dilauroylphosphatidylcholine (DLPC) liposomes K_MD^2D is 0.0005–0.0008 mol% based on Förster resonance energy transfer (FRET) measurements, slightly tighter than bicelles. These results emphasize the importance of the lipid membrane in influencing dimer affinity.     

Gene Structure and Molecular Phylogeny of the Linker Chains from the Giant Annelid Hexagonal Bilayer Hemoglobins

Giant extracellular hexagonal bilayer hemoglobin (HBL-Hb), found only in annelids, is an ∼3500-kDa heteropolymeric structure involved in oxygen transport. The HBL-Hbs are comprised of globin and linker chains, the latter being required for the assembly of the quaternary structure. The linker chains, varying in size from 225 to 283 amino acids, have a conserved cysteine-rich domain within their N-terminal moiety that is homologous to the cysteine-rich modules constituting the ligand binding domain of the low-density lipoprotein receptor (LDLR) protein family found in many metazoans. We have investigated the gene structure of linkers from Arenicola marina, Alvinella pompejana, Nereis diversicolor, Lumbricus terrestris, and Riftia pachyptila. We found, contrary to the results obtained earlier with linker genes from N. diversicolor and L. terrestris, that in all of the foregoing cases, the linker LDL-A module is flanked by two phase 1 introns, as in the human LDLR gene, with two more introns in the 3′ side whose positions varied with the species. In addition, we obtained 13 linker cDNAs that have been determined experimentally or found in the EST database LumbriBASE. A molecular phylogenetic analysis of the linker primary sequences demonstrated that they cluster into two distinct families of linker proteins. We propose that the common gene ancestor to annelid linker genes exhibited a four-intron and five-exon structure and gave rise to the two families subsequent to a duplication event.     

Effects of UV radiation and nitrate limitation on the production of biogenic sulfur compounds by marine phytoplankton

We tested the effects of UV radiation (UVR) and nitrate limitation on the production of dimethylsulfide (DMS), particulate dimethylsulfoniopropionate (DMSPp), and particulate dimethylsulfoxide (DMSOp) in natural seawater from the Gulf of Mexico and in phytoplankton cultures. DMS/Chl a ratios in PAR-only and PAR + UV-exposed seawater were 0.44–2.0 and 0.46–1.9 nmol DMS μg⁻¹ Chl a, respectively, whereas the ratios in cultures of Amphidinium carterae were 1.0–2.2 nmol μg⁻¹ in PAR-exposed samples and 0.91–2.1 nmol μg⁻¹ in PAR + UV-exposed samples. These results suggested that UVR did not substantially affect DMS/Chl a ratios in seawater and A. carterae culture samples. Similarly, UVR had no significant effect on DMSOp/Chl a in seawater samples (0.83–1.6 nmol DMSO μg⁻¹ Chl a for PAR + UV vs. 0.70–1.5 nmol μg⁻¹ for PAR-exposed seawater samples, respectively) or in A. carterae cultures (0.20–1.3 and 0.19–0.88 nmol DMSO μg⁻¹ Chl a in PAR + UV- and PAR-exposed cultures, respectively). In an experiment with the diatom, Thalassiosira oceanica, the culture was grown in high nitrate (30 μM) or low nitrate (6 μM) media and exposed to PAR-only or PAR + UV. The low nitrate, PAR-only samples showed an increase of intracellular dimethylsulfoniopropionate (DMSP) concentration from 2.1 to 15 mmol L⁻¹ in 60 h, but the increase occurred only after cultures reached the stationary phase. Cultures of T. oceanica grown under UVR had lower growth rates than those under PAR-only (μ′ = 0.17 and 0.32 d⁻¹, respectively) and perhaps did not experience severe nitrate limitation even in the low nitrate treatment. These results suggest that the elevated UVR in low nitrate environments could result in reduction of DMSP in some species, whereas DMSP concentrations would not be affected in eutrophic areas.     

Crystal structure of the single-stranded RNA binding protein HutP from Geobacillus thermodenitrificans

RNA binding proteins control gene expression by the attenuation/antitermination mechanism. HutP is an RNA binding antitermination protein. It regulates the expression of hut operon when it binds with RNA by modulating the secondary structure of single-stranded hut mRNA. HutP necessitates the presence of l-histidine and divalent metal ion to bind with RNA. Herein, we report the crystal structures of ternary complex (HutP–l-histidine–Mg²⁺) and EDTA (0.5 M) treated ternary complex (HutP–l-histidine–Mg²⁺), solved at 1.9 Å and 2.5 Å resolutions, respectively, from Geobacillus thermodenitrificans. The addition of 0.5 M EDTA does not affect the overall metal-ion mediated ternary complex structure and however, the metal ions at the non-specific binding sites are chelated, as evidenced from the results of structural features.     

Oxygen binding of the heme-containing subunits of the extracellular hemoglobin of Lumbricus, terrestris

Fully oxygenated heme-containing subunits of the extracellular hemoglobin of $\text{Lumbricus}\text{,}\text{terrestris}$ were isolated by gel filtration in 0.05 M sodium borate, pH 9.0 at 7 ± 1°C. The polypeptide chain composition of the isolated subunits was determined by SDS PAGE. Most of the subunits exhibited reversible oxygenation curves with oxygen affinities higher than the intact hemoglobin. The Hill plots were nonlinear for $\text{Lumbricus}$ hemoglobin and its subunits: the latter exhibited substantial cooperativities as evidenced by Hill constants at half-saturation in the range 2.0 to 2.8.     

The complete mitochondrial genome of Taeniogonalos taihorina (Bischoff) (Hymenoptera: Trigonalyidae) reveals a novel gene rearrangement pattern in the Hymenoptera

The family Trigonalyidae is considered to be one of the most basal lineages in the suborder Apocrita of Hymenoptera. Here, we determine the first complete mitochondrial genome of the Trigonalyidae, from the species Taeniogonalos taihorina (Bischoff, 1914). This mitochondrial genome is 15,927 bp long, with a high A + T-content of 84.60%. It contains all of the 37 typical animal mitochondrial genes and an A + T-rich region. The orders and directions of all genes are different from those of previously reported hymenopteran mitochondrial genomes. Eight tRNA genes, three protein-coding genes and the A + T-rich region were rearranged, with the dominant gene rearrangement events being translocation and local inversion. The arrangements of three tRNA clusters, trnY–trnM–trnI–trnQ, trnW–trnL2–trnC, and trnH–trnA–trnR–trnN–trnS–trnE–trnF, and the position of the cox1 gene, are novel to the Hymenoptera, even the insects. Six long intergenic spacers are present in the genome. The secondary structures of the RNA genes are normal, except for trnS2, in which the D-stem pairing is absent.     

Preparative separation of alkaloids from Nelumbo nucifera leaves by conventional and pH-zone-refining counter-current chromatography

Two modes of high-speed counter-current chromatography (HSCCC) were successfully applied to the separation of alkaloids from crude extract of Nelumbo nucifera leaves. The conventional HSCCC separations were performed with a two-phase solvent system composed of tetrachloromethane–CHCl₃–methanol–0.1 M HCl at a volume ratio of 1:3:3:2 (v/v/v/v), and 120 mg crude extract could be successfully separated. pH-Zone-refining CCC was performed with a two-phase solvent system composed of petroleum ether (60–90 °C)–ethyl acetate–methanol–water (5:5:2:8, v/v/v/v) where triethylamine (10 mM) was added to the upper organic stationary phase as a retainer and hydrochloric acid (5 mM) to the aqueous mobile phase as an eluent. From 4.0 g of the crude extract, 120 mg N-nornuciferine, 1020 mg nuciferine and 96 mg roemerine were obtained in a single run each with a purity of over 98% as determined by HPLC. The structures of the isolated compounds were identified by ESI-MS, ¹H NMR and ¹³C NMR.     

Evaporative water loss and energy metabolic in two small mammals,voles (Eothenomys miletus) and mice (Apodemus chevrieri), in Hengduan mountains region

Evaporative water loss (EWL) and energy metabolism were measured at different temperatures in Eothenomys miletus and Apodemus chevrieri in dry air. The thermal neutral zone (TNZ) of E. miletus was 22.5–30 °C and that of A. chevrieri was 20–27.5 °C. Mean body temperatures of the two species were 35.75±0.5 and 36.54±0.61 °C. Basal metabolic rates (BMR) were 1.92±0.17 and 2.7±0.5 ml O₂/g h, respectively. Average minimum thermal conductance (C_m) were 0.23±0.08 and 0.25±0.06 ml O₂/g h °C. EWL in E. miletus and A. chevrieri increased with the increase in temperature; the maximal EWL at 35 °C was 4.78±0.6 mg H₂O/g h in E. miletus, and 5.92±0.43 mg H₂O/g h in A. chevrieri. Percentage of evaporative heat loss to total heat production (EHL/HP) increased with the increase in temperature; the maximal EHL/HP was 22.45% at 30 °C in E. miletus, and in A. chevrieri it was 19.96% at 27.5 °C. The results may reflect features of small rodents in the Hengduan mountains region: both E. miletus and A. chevrieri have high levels of BMR and high levels of total thermal conductance, compared with the predicted values based on their body masses, while their body temperatures are relatively low. EWL plays an important role in temperature regulation.     

Evolutionary Origin and Host Range of Plagiotoma lumbrici (Ciliophora,Hypotrichia), an Obligate Gut Symbiont of Lumbricid Earthworms

Four common earthworm species, the anecic Lumbricus terrestris, the endogeic Octolasion tyrteum as well as the epigeic Eisenia fetida and Dendrobaena veneta, were examined for the presence of the microbial gut symbiont Plagiotoma lumbrici. The evolutionary origin of this endobiotic microbe was reconstructed, using the 18S rRNA gene, the ITS1‐5.8S‐ITS2 region, and the first two domains of the 28S rRNA gene. Plagiotoma lumbrici was exclusively detected in the anecic Lumbricus terrestris. Multigene analyses and the ITS2 secondary structure robustly determined the phylogenetic home of Plagiotoma lumbrici populations within the oxytrichid Dorsomarginalia (Spirotrichea: Hypotrichia) as a sister taxon of the free‐living Hemiurosomoida longa. This indicates that earthworms obtained their gut endosymbiont by ingesting soil/leaf litter containing oxytrichine ciliates that became adapted to the intestinal tract of earthworms. Interestingly, according to the literature data, Plagiotoma lumbrici was detected in multiple anecic and some epigeic but never in endogeic earthworms. These observations suggest that Plagiotoma lumbrici might be adapted to certain gut conditions and the lifestyle of anecic Lumbricidae, such as Lumbricus, Aporrectodea, and Scherotheca, as well as of some co‐occurring epigeic Lumbricus species.     

Null genotypes of GSTM1 and GSTT1 contribute to increased risk of diabetes mellitus: A meta-analysis

Diabetes mellitus (DM) is a common disease which results from various causes including genetic and environmental factors. Glutathione S-Transferase M1 (GSTM1) and Glutathione S-Transferase T1 (GSTT1) genes are polymorphic in human and the null genotypes lead to the absence of enzyme function. Many studies assessed the associations between GSTM1/GSTT1 null genotypes and DM risk but reported conflicting results. In order to get a more precise estimate of the associations of GSTM1/GSTT1 null genotypes with DM risk, we performed this meta-analysis. Published literature from PubMed, Embase and China Biology Medicine (CBM) databases was searched for eligible studies. Pooled odds ratios (OR) and corresponding 95% confidence intervals (95%CI) were calculated using a fixed- or random-effects model. 11 publications (a total of 2577 cases and 4572 controls) were finally included into this meta-analysis. Meta-analyses indicated that null genotypes of GSTM1/GSTT1 and dual null genotype of GSTM1–GSTT1 were all associated with increased risk of DM (GSTM1: OR _{random-effects} = 1.60, 95%CI 1.10–2.34, P_OR = 0.014; GSTT1: OR _{random-effects} = 1.47, 95%CI 1.12–1.92, P_OR = 0.005; GSTM1–GSTT1: OR _{fixed-effects} = 1.83, 95%CI 1.30–2.59, P_OR = 0.001). Subgroup by ethnicity suggested significant associations between null genotypes of GSTM1 and GSTT1 and DM risk among Asians (GSTM1: OR _{random-effects} = 1.77, 95%CI 1.24–2.53, P_OR = 0.002; GSTT1: OR _{random-effects} = 1.58, 95%CI 1.09–2.27, P_OR = 0.015). This meta-analysis suggests null genotypes of GSTM1/GSTT1 and dual null genotype of GSTM1–GSTT1 are all associated with increased risk of DM, and null genotypes of GSTM1/GSTT1 and dual null genotype of GSTM1–GSTT1 are potential biomarkers of DM.     

Biochemical characterization of a low molecular weight aspartic protease inhibitor from thermo-tolerant Bacillus licheniformis: Kinetic interactions with Pepsin

The present article reports a low molecular weight aspartic protease inhibitor, API, from a newly isolated thermo-tolerant Bacillus licheniformis. The inhibitor was purified to homogeneity as shown by rp-HPLC and SDS-PAGE. API is found to be stable over a broad pH range of 2–11 and at temperature 90 °C for 2 1/2 h. It has a Mr (relative molecular mass) of 1363 Da as shown by MALDI-TOF spectra and 1358 Da as analyzed by SDS-PAGE .The amino acid analysis of the peptide shows the presence of 12 amino acid residues having Mr of 1425 Da. The secondary structure of API as analyzed by the CD spectra showed 7% α-helix, 49% β-sheet and 44% aperiodic structure. The Kinetic studies of Pepsin–API interactions reveal that API is a slow-tight binding competitive inhibitor with the IC₅₀ and K_i values 4.0 nM and (3.83 nM–5.31 nM) respectively. The overall inhibition constant K_i? value is 0.107 ± 0.015 nM. The progress curves are time-dependent and consistent with slow-tight binding inhibition: E + I ? (k₄, k₅) EI ? (k₆, k₇) EI?. Rate constant k₆ = 2.73 ± 0.32 s^− 1 reveals a fast isomerization of enzyme–inhibitor complex and very slow dissociation as proved by k₇ = 0.068 ± 0.009 s^− 1. The Rate constants from the intrinsic tryptophanyl fluorescence data is in agreement with those obtained from the kinetic analysis; therefore, the induced conformational changes were correlated to the isomerization of EI to EI?.     

Association of interleukin-13 SNP rs20541 with allergic rhinitis risk: A meta-analysis

Studies investigating the association between interleukin-13 (IL-13) single nucleotide polymorphism (SNP) rs20541 and allergic rhinitis (AR) risk have reported conflicting results. The aim of the present study was to conduct a meta-analysis assessing the possible association of IL-13 SNP rs20541 with AR risk. Eight studies were included in the present meta-analysis (2153 cases and 3931 controls). The combined results based on all studies showed that IL-13 SNP rs20541 was associated with increased AR risk (Gln versus Arg: odds ratio (OR) = 1.18, 95% confidence interval (CI) = 1.08–1.30; Gln/Gln versus Arg/Arg: OR = 1.52, 95% CI = 1.20–1.92; Arg/Gln + Gln/Gln versus Arg/Arg: OR = 1.19, 95% CI = 1.06–1.33; Gln/Gln versus Arg/Gln + Arg/Arg: OR = 1.42, 95% CI = 1.13–1.79). When stratifying for race, IL-13 SNP rs20541 exhibited increased AR risk in Asians (Gln versus Arg: OR = 1.20, 95% CI = 1.06–1.36; Gln/Gln versus Arg/Arg: OR = 1.57, 95% CI = 1.17–2.12; Arg/Gln + Gln/Gln versus Arg/Arg: OR = 1.22, 95% CI = 1.04–1.44; Gln/Gln versus Arg/Gln + Arg/Arg: OR = 1.45, 95% CI = 1.09–1.93), while no significant association was detected in Caucasians (Gln versus Arg: OR = 1.28, 95% CI = 0.93 ~ 1.78; Gln/Gln versus Arg/Arg: OR = 1.42, 95% CI = 0.96–2.11; Arg/Gln + Gln/Gln versus Arg/Arg: OR = 1.35, 95% CI = 0.89–2.05; Gln/Gln versus Arg/Gln + Arg/Arg: OR = 1.37, 95% CI = 0.93–2.02). This meta-analysis supported that IL-13 SNP rs20541 was associated with AR, particularly in Asians.