Polyploidization as a retraction force in plant genome evolution: sequence rearrangements in triticale

Background

Polyploidization is a major evolutionary process in plants where hybridization and chromosome doubling induce enormous genomic stress and can generate genetic and epigenetic modifications. However, proper evaluation of DNA sequence restructuring events and the precise characterization of sequences involved are still sparse.

Methodology/Principal Findings

Inter Retrotransposons Amplified Polymorphism (IRAP), Retrotransposons Microsatellite Amplified Polymorphism (REMAP) and Inter Simple Sequence Repeat (ISSR) largely confirmed the absence of any intraspecific variation in wheat, rye and triticale. The comparative analysis of banding profiles between wheat and rye inbred lines revealed 34% of monomorphic (common to both parental species) bands for the ten different primer combinations used. The analysis of triticale plants uncovered nearly 51% of rearranged bands in the polyploid, being the majority of these modifications, due to the loss of rye bands (83%). Sequence analysis of rye fragments absent in triticale revealed for instance homology with hydroxyproline-rich glycoproteins (HRGP), a protein that belongs to a major family of inducible defence response proteins. Conversely, a wheat-specific band absent in triticale comprises a nested structure of copia-like retrotransposons elements, namely Claudia and Barbara. Sequencing of a polyploid-specific band (absent in both parents) revealed a microsatellite related sequence. Cytological studies using Fluorescent In Situ Hybridization (FISH) with REMAP products revealed a widespread distribution of retrotransposon and/or microsatellite flanking sequences on rye chromosomes, with a preferential accumulation in heterochromatic sub-telomeric domains.

Conclusions/Significance

Here, we used PCR-based molecular marker techniques involving retrotransposons and microsatellites to uncover polyploidization induced genetic restructuring in triticale. Sequence analysis of rearranged genomic fragments either from rye or wheat origin showed these to be retrotransposon-related as well as coding sequences. Further FISH analysis revealed possible chromosome hotspots for sequence rearrangements. The role of chromatin condensation on the origin of genomic rearrangements mediated by polyploidization in triticale is also discussed.     

Comparative cytogenetic analysis of the genomes of the model grass Brachypodium distachyon and its close relatives

Background and Aims

Brachypodium is a small genus of temperate grasses that comprises 12–15 species. Brachypodium distachyon is now well established as a model species for temperate cereals and forage grasses. In contrast to B. distachyon, other members of the genus have been poorly investigated at the chromosome level or not at all.

Methods

Twenty accessions comprising six species and two subspecies of Brachypodium were analysed cytogenetically. Measurements of nuclear genome size were made by flow cytometry. Chromosomal localization of 18–5·8–25S rDNA and 5S rDNA loci was performed by dual-colour fluorescence in situ hybridization (FISH) on enzymatically digested root-tip meristematic cells. For comparative phylogenetic analyses genomic in situ hybridization (GISH) applied to somatic chromosome preparations was used.

Key Results

All Brachypodium species examined have rather small genomes and chromosomes. Their chromosome numbers and genome sizes vary from 2n = 10 and 0·631 pg/2C in B. distachyon to 2n = 38 and 2·57 pg/2C in B. retusum, respectively. Genotypes with 18 and 28 chromosomes were found among B. pinnatum accessions. GISH analysis revealed that B. pinnatum with 28 chromosomes is most likely an interspecific hybrid between B. distachyon (2n = 10) and B. pinnatum (2n = 18). Two other species, B. phoenicoides and B. retusum, are also allopolyploids and B. distachyon or a close relative seems to be one of their putative ancestral species. In chromosomes of all species examined the 45S rDNA loci are distally distributed whereas loci for 5S rDNA are pericentromeric.

Conclusions

The increasing significance of B. distachyon as a model grass emphasizes the need to understand the evolutionary relationships in the genus Brachypodium and to ensure consistency in the biological nomenclature of its species. Modern molecular cytogenetic techniques such as FISH and GISH are suitable for comparative phylogenetic analyses and may provide informative chromosome- and/or genome-specific landmarks.     

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon

Background

The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon.

Results

B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner.

Conclusions

B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.     

Alien introgression in the grasses Lolium perenne (perennial ryegrass) and Festuca pratensis (meadow fescue): the development of seven monosomic substitution lines and their molecular and cytological characterization

Background and Aims

To address the issues associated with food security, environmental change and bioenergy in the context of crop plants, the production, identification and evaluation of novel plant phenotypes is fundamental. One of the major routes to this end will be wide hybridization and introgression breeding. The transfer of chromosomes and chromosome segments between related species (chromosome engineering or alien introgression) also provides an important resource for determining the genetic control of target traits. However, the realization of the full potential of chromosome engineering has previously been hampered by the inability to identify and characterize interspecific introgressions accurately.

Methods

Seven monosomic substitution lines have been generated comprising Festuca pratensis as the donor species and Lolium perenne as the recipient. Each of the seven lines has a different L. perenne chromosome replaced by the homoeologous F. pratensis chromosome (13 L. perenne + 1 F. pratensis chromosome). Molecular markers and genomic in situ hybridization (GISH) were used to assign the F. pratensis chromosomes introgressed in each of the monosomic substitutions to a specific linkage group. Cytological observations were also carried out on metaphase I of meiosis in each of the substitution lines.

Results

A significant level of synteny was found at the macro-level between L. perenne and F. pratensis. The observations at metaphase I revealed the presence of a low level of interspecific chromosomal translocations between these species.

Discussion

The isolation of the seven monosomic substitution lines provides a resource for dissecting the genetic control of important traits and for gene isolation. Parallels between the L. perenne/F. pratensis system and the Pooideae cereals such as wheat, barley, rye, oats and the model grass Brachypodium distachyon present opportunities for a comparison across the species in terms of genotype and phenotype.     

Characterization of the caleosin gene family in the Triticeae

Background

The caleosin genes encode proteins with a single conserved EF hand calcium-binding domain and comprise small gene families found in a wide range of plant species. Some members of the gene family have been shown to be upregulated by environmental stresses including low water availability and high salinity. Caleosin 3 from wheat has been shown to interact with the α-subunit of the heterotrimeric G proteins, and to act as a GTPase activating protein (GAP). This study characterizes the size and diversity of the gene family in wheat and related species and characterizes the differential tissue-specific expression of members of the gene family.

Results

A total of 34 gene family members that belong to eleven paralogous groups of caleosins were identified in the hexaploid bread wheat, T. aestivum. Each group was represented by three homeologous copies of the gene located on corresponding homeologous chromosomes, except the caleosin 10, which has four gene copies. Ten gene family members were identified in diploid barley, Hordeum vulgare, and in rye, Secale cereale, seven in Brachypodium distachyon, and six in rice, Oryza sativa. The analysis of gene expression was assayed in triticale and rye by RNA-Seq analysis of 454 sequence sets and members of the gene family were found to have diverse patterns of gene expression in the different tissues that were sampled in rye and in triticale, the hybrid hexaploid species derived from wheat and rye. Expression of the gene family in wheat and barley was also previously determined by microarray analysis, and changes in expression during development and in response to environmental stresses are presented.

Conclusions

The caleosin gene family had a greater degree of expansion in the Triticeae than in the other monocot species, Brachypodium and rice. The prior implication of one member of the gene family in the stress response and heterotrimeric G protein signaling, points to the potential importance of the caleosin gene family. The complexity of the family and differential expression in various tissues and under conditions of abiotic stress suggests the possibility that caleosin family members may play diverse roles in signaling and development that warrants further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-239) contains supplementary material, which is available to authorized users.     

An efficient method for transient gene expression in monocots applied to modify the Brachypodium distachyon cell wall

Background

Agrobacterium-mediated transformation is widely used to produce insertions into plant genomes. There are a number of well-developed Agrobacterium-mediated transformation methods for dicotyledonous plants, but there are few for monocotyledonous plants.

Methods

Three hydrolase genes were transiently expressed in Brachypodium distachyon plants using specially designed vectors that express the gene product of interest and target it to the plant cell wall. Expression of functional hydrolases in genotyped plants was confirmed using western blotting, activity assays, cell wall compositional analysis and digestibility tests.

Key Results

An efficient, new, Agrobacterium-mediated approach was developed for transient gene expression in the grass B. distachyon, using co-cultivation of mature seeds with bacterial cells. This method allows transformed tissues to be obtained rapidly, within 3–4 weeks after co-cultivation. Also, the plants carried transgenic tissue and maintained transgenic protein expression throughout plant maturation. The efficiency of transformation was estimated at around 5 % of initially co-cultivated seeds. Application of this approach to express three Aspergillus nidulans hydrolases in the Brachypodium cell wall successfully confirmed its utility and resulted in the expected expression of active microbial proteins and alterations of cell wall composition. Cell wall modifications caused by expression of A. nidulans α-arabinofuranosidase and α-galactosidase increased the biodegradability of plant biomass.

Conclusions

This newly developed approach is a quick and efficient technique for expressing genes of interest in Brachypodium plants, which express the gene product throughout development. In the future, this could be used for broad functional genomics studies of monocots and for biotechnological applications, such as plant biomass modification for biofuel production.     

Next-generation sequencing of flow-sorted wheat chromosome 5D reveals lineage-specific translocations and widespread gene duplications

Background

The ~17 Gb hexaploid bread wheat genome is a high priority and a major technical challenge for genomic studies. In particular, the D sub-genome is relatively lacking in genetic diversity, making it both difficult to map genetically, and a target for introgression of agriculturally useful traits. Elucidating its sequence and structure will therefore facilitate wheat breeding and crop improvement.

Results

We generated shotgun sequences from each arm of flow-sorted Triticum aestivum chromosome 5D using 454 FLX Titanium technology, giving 1.34× and 1.61× coverage of the short (5DS) and long (5DL) arms of the chromosome respectively. By a combination of sequence similarity and assembly-based methods, ~74% of the sequence reads were classified as repetitive elements, and coding sequence models of 1314 (5DS) and 2975 (5DL) genes were generated. The order of conserved genes in syntenic regions of previously sequenced grass genomes were integrated with physical and genetic map positions of 518 wheat markers to establish a virtual gene order for chromosome 5D.

Conclusions

The virtual gene order revealed a large-scale chromosomal rearrangement in the peri-centromeric region of 5DL, and a concentration of non-syntenic genes in the telomeric region of 5DS. Although our data support the large-scale conservation of Triticeae chromosome structure, they also suggest that some regions are evolving rapidly through frequent gene duplications and translocations.

Sequence accessions

EBI European Nucleotide Archive, Study no. ERP002330

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1080) contains supplementary material, which is available to authorized users.     

Deciphering heterogeneity in pig genome assembly Sscrofa9 by isochore and isochore-like region analyses

Background

The isochore, a large DNA sequence with relatively small GC variance, is one of the most important structures in eukaryotic genomes. Although the isochore has been widely studied in humans and other species, little is known about its distribution in pigs.

Principal Findings

In this paper, we construct a map of long homogeneous genome regions (LHGRs), i.e., isochores and isochore-like regions, in pigs to provide an intuitive version of GC heterogeneity in each chromosome. The LHGR pattern study not only quantifies heterogeneities, but also reveals some primary characteristics of the chromatin organization, including the followings: (1) the majority of LHGRs belong to GC-poor families and are in long length; (2) a high gene density tends to occur with the appearance of GC-rich LHGRs; and (3) the density of LINE repeats decreases with an increase in the GC content of LHGRs. Furthermore, a portion of LHGRs with particular GC ranges (50%–51% and 54%–55%) tend to have abnormally high gene densities, suggesting that biased gene conversion (BGC), as well as time- and energy-saving principles, could be of importance to the formation of genome organization.

Conclusion

This study significantly improves our knowledge of chromatin organization in the pig genome. Correlations between the different biological features (e.g., gene density and repeat density) and GC content of LHGRs provide a unique glimpse of in silico gene and repeats prediction.     

Genome-wide exploration of the molecular evolution and regulatory network of mitogen-activated protein kinase cascades upon multiple stresses in Brachypodium distachyon

Background

Brachypodium distachyon is emerging as a widely recognized model plant that has very close relations with several economically important Poaceae species. MAPK cascade is known to be an evolutionarily conserved signaling module involved in multiple stresses. Although the gene sequences of MAPK and MAPKK family have been fully identified in B. distachyon, the information related to the upstream MAPKKK gene family especially the regulatory network among MAPKs, MAPKKs and MAPKKKs upon multiple stresses remains to be understood.

Results

In this study, we have identified MAPKKKs which belong to the biggest gene family of MAPK cascade kinases. We have systematically investigated the evolution of whole MAPK cascade kinase gene family in terms of gene structures, protein structural organization, chromosomal localization, orthologs construction and gene duplication analysis. Our results showed that most BdMAPK cascade kinases were located at the low-CpG-density region, and the clustered members in each group shared similar structures of the genes and proteins. Synteny analysis showed that 62 or 21 pairs of duplicated orthologs were present between B. distachyon and Oryza sativa, or between B. distachyon and Arabidopsis thaliana respectively. Gene expression data revealed that BdMAPK cascade kinases were rapidly regulated by stresses and phytohormones. Importantly, we have constructed a regulation network based on co-expression patterns of the expression profiles upon multiple stresses performed in this study.

Conclusions

BdMAPK cascade kinases were involved in the signaling pathways of multiple stresses in B. distachyon. The network of co-expression regulation showed the most of duplicated BdMAPK cascade kinase gene orthologs demonstrated their convergent function, whereas few of them developed divergent function in the evolutionary process. The molecular evolution analysis of identified MAPK family genes and the constructed MAPK cascade regulation network under multiple stresses provide valuable information for further investigation of the functions of BdMAPK cascade kinase genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1452-1) contains supplementary material, which is available to authorized users.     

Barcoding bugs: DNA-based identification of the true bugs (Insecta: Hemiptera: Heteroptera)

Background

DNA barcoding, the analysis of sequence variation in the 5′ region of the mitochondrial cytochrome c oxidase I (COI) gene, has been shown to provide an efficient method for the identification of species in a wide range of animal taxa. In order to assess the effectiveness of barcodes in the discrimination of Heteroptera, we examined 344 species belonging to 178 genera, drawn from specimens in the Canadian National Collection of Insects.

Methodology/Principal Findings

Analysis of the COI gene revealed less than 2% intra-specific divergence in 90% of the taxa examined, while minimum interspecific distances exceeded 3% in 77% of congeneric species pairs. Instances where barcodes fail to distinguish species represented clusters of morphologically similar species, except one case of barcode identity between species in different genera. Several instances of deep intraspecific divergence were detected suggesting possible cryptic species.

Conclusions/Significance

Although this analysis encompasses 0.8% of the described global fauna, our results indicate that DNA barcodes will aid the identification of Heteroptera. This advance will be useful in pest management, regulatory and environmental applications and will also reveal species that require further taxonomic research.     

Hospital and community ampicillin-resistant Enterococcus faecium are evolutionarily closely linked but have diversified through niche adaptation

Background

Ampicillin-resistant Enterococcus faecium (ARE) has emerged as a nosocomial pathogen. Here, we quantified ARE carriage in different community sources and determined genetic relatedness with hospital ARE.

Methods and Results

ARE was recovered from rectal swabs of 24 of 79 (30%) dogs, 11 of 85 (13%) cats and 0 of 42 horses and from 3 of 40 (8%) faecal samples of non-hospitalized humans receiving amoxicillin. Multi-locus Sequence Typing revealed 21 sequence types (STs), including 5 STs frequently associated with hospital-acquired infections. Genes previously found to be enriched in hospital ARE, such as IS16, orf903, orf905, orf907, were highly prevalent in community ARE (≥79%), while genes with a proposed role in pathogenesis, such as esp, hyl and ecbA, were found rarely (≤5%) in community isolates. Comparative genome analysis of 2 representative dog isolates revealed that the dog strain of ST192 was evolutionarily closely linked to two previously sequenced hospital ARE, but had, based on gene content, more genes in common with the other, evolutionarily more distantly related, dog strain (ST266).

Conclusion

ARE were detected in dogs, cats and sporadically in healthy humans, with evolutionary linkage to hospital ARE. Yet, their accessory genome has diversified, probably as a result of niche adaptation.     

Muscle hypertrophy in prepubescent tennis players: a segmentation MRI study

Purpose

To asses if tennis at prepubertal age elicits the hypertrophy of dominant arm muscles.

Methods

The volume of the muscles of both arms was determined using magnetic resonance imaging (MRI) in 7 male prepubertal tennis players (TP) and 7 non-active control subjects (CG) (mean age 11.0±0.8 years, Tanner 1–2).

Results

TP had 13% greater total muscle volume in the dominant than in the contralateral arm. The magnitude of inter-arm asymmetry was greater in TP than in CG (13 vs 3%, P<0.001). The dominant arm of TP was 16% greater than the dominant arm of CG (P<0.01), whilst non-dominant arms had similar total muscle volumes in both groups (P = 0.25), after accounting for height as covariate. In TP, dominant deltoid (11%), forearm supinator (55%) and forearm flexors (21%) and extensors (25%) were hypertrophied compared to the contralateral arm (P<0.05). In CG, the dominant supinator muscle was bigger than its contralateral homonimous (63%, P<0.05).

Conclusions

Tennis at prepubertal age is associated with marked hypertrophy of the dominant arm, leading to a marked level of asymmetry (+13%), much greater than observed in non-active controls (+3%). Therefore, tennis particpation at prepubertal age is associated with increased muscle volumes in dominant compared to the non-dominant arm, likely due to selectively hypertrophy of the loaded muscles.     

The cervical microbiome over 7 years and a comparison of methodologies for its characterization

Background

The rapidly expanding field of microbiome studies offers investigators a large choice of methods for each step in the process of determining the microorganisms in a sample. The human cervicovaginal microbiome affects female reproductive health, susceptibility to and natural history of many sexually transmitted infections, including human papillomavirus (HPV). At present, long-term behavior of the cervical microbiome in early sexual life is poorly understood.

Methods

The V6 and V6–V9 regions of the 16S ribosomal RNA gene were amplified from DNA isolated from exfoliated cervical cells. Specimens from 10 women participating in the Natural History Study of HPV in Guanacaste, Costa Rica were sampled successively over a period of 5–7 years. We sequenced amplicons using 3 different platforms (Sanger, Roche 454, and Illumina HiSeq 2000) and analyzed sequences using pipelines based on 3 different classification algorithms (usearch, RDP Classifier, and pplacer).

Results

Usearch and pplacer provided consistent microbiome classifications for all sequencing methods, whereas RDP Classifier deviated significantly when characterizing Illumina reads. Comparing across sequencing platforms indicated 7%–41% of the reads were reclassified, while comparing across software pipelines reclassified up to 32% of the reads. Variability in classification was shown not to be due to a difference in read lengths. Six cervical microbiome community types were observed and are characterized by a predominance of either G. vaginalis or Lactobacillus spp. Over the 5–7 year period, subjects displayed fluctuation between community types. A PERMANOVA analysis on pairwise Kantorovich-Rubinstein distances between the microbiota of all samples yielded an F-test ratio of 2.86 (p<0.01), indicating a significant difference comparing within and between subjects’ microbiota.

Conclusions

Amplification and sequencing methods affected the characterization of the microbiome more than classification algorithms. Pplacer and usearch performed consistently with all sequencing methods. The analyses identified 6 community types consistent with those previously reported. The long-term behavior of the cervical microbiome indicated that fluctuations were subject dependent.     

Evaluation of the Use of the Polyubiquitin Genes,Ubi4 and Ubi10 as Reference Genes for Expression Studies in Brachypodium distachyon

Background

Brachypodium distachyon is emerging as the model plant for temperate grass research and the genome of the community line Bd21 has been sequenced. Additionally, techniques have been developed for Agrobacterium-mediated transformation for the generation of T-DNA insertional lines. Recently, it was reported that expression of the polyubiquitin genes, Ubi4 and Ubi10 are stable in different tissues and growth hormone-treated plant samples, leading to the conclusion that both Ubi4 and Ubi10 are good reference genes for normalization of gene expression data using real-time, quantitative PCR (qPCR).

Principal Findings

Mining of the Joint Genome Institute (JGI) 8X Brachypodium distachyon genome assembly showed that Ubi4 and Ubi10 share a high level of sequence identity (89%), and in silico analyses of the sequences of Ubi4 (Bradi3g04730) and Ubi10 (Bradi1g32860) showed that the primers used previously exhibit multiple binding sites within the coding sequences arising from the presence of tandem repeats of the coding regions. This can potentially result in over-estimation of steady-state levels of Ubi4 and Ubi10. Additionally, due to the high level of sequence identity between both genes, primers used previously for amplification of Ubi4 can bind to Ubi10 and vice versa, resulting in the formation of non-specific amplification products.

Conclusions

The results from this study indicate that the primers used previously were not sufficiently robust and specific. Additionally, their use would result in over-estimation of the steady-state expression levels of Ubi4. Our results question the validity of using the previously proposed primer sets for qPCR amplification of Ubi4 and Ubi10. We demonstrate that primers designed to target the 3′-UTRs of Ubi4 and Ubi10 are better suited for real-time normalization of steady-state expression levels in Brachypodium distachyon.     

Genotypic diversity and drug susceptibility patterns among M. tuberculosis complex isolates from South-Western Ghana

Objective

The aim of this study was to use spoligotyping and large sequence polymorphism (LSP) to study the population structure of M. tuberculosis complex (MTBC) isolates.

Methods

MTBC isolates were identified using standard biochemical procedures, IS6110 PCR, and large sequence polymorphisms. Isolates were further typed using spoligotyping, and the phenotypic drug susceptibility patterns were determined by the proportion method.

Result

One hundred and sixty-two isolates were characterised by LSP typing. Of these, 130 (80.25%) were identified as Mycobacterium tuberculosis sensu stricto (MTBss), with the Cameroon sub-lineage being dominant (N = 59/130, 45.38%). Thirty-two (19.75%) isolates were classified as Mycobacterium africanum type 1, and of these 26 (81.25%) were identified as West-Africa I, and 6 (18.75%) as West-Africa II. Spoligotyping sub-lineages identified among the MTBss included Haarlem (N = 15, 11.53%), Ghana (N = 22, 16.92%), Beijing (4, 3.08%), EAI (4, 3.08%), Uganda I (4, 3.08%), LAM (2, 1.54%), X (N = 1, 0.77%) and S (2, 1.54%). Nine isolates had SIT numbers with no identified sub-lineages while 17 had no SIT numbers. MTBss isolates were more likely to be resistant to streptomycin (p<0.008) and to any drug resistance (p<0.03) when compared to M. africanum.

Conclusion

This study demonstrated that overall 36.4% of TB in South-Western Ghana is caused by the Cameroon sub-lineage of MTBC and 20% by M. africanum type 1, including both the West-Africa 1 and West-Africa 2 lineages. The diversity of MTBC in Ghana should be considered when evaluating new TB vaccines.     

Tumor necrosis factor-alpha g308α gene polymorphism and essential hypertension: a meta-analysis involving 2244 participants

Background

The tumor necrosis factor-alpha (TNFα) G308A gene polymorphism has been implicated in susceptibility to essential hypertension (EH), but study results are still controversial.

Objective and Methods

The present meta-analysis is performed to investigate the relationship between the TNFα G308A gene polymorphism and EH. Electronic databases were searched and seven separate studies on the association of the TNF α G308A gene polymorphism with EH were analyzed. The meta-analysis involved 1092 EH patients and 1152 controls. The pooled odds ratios (ORs) and their corresponding 95% confidence interval (CI) were calculated by a fixed or random effect model.

Results

A significant relationship between the TNFα G308A gene polymorphism and EH was found in an allelic genetic model (OR: 1.45, 95% CI: 1.17 to 1.80, P = 0.0008), a recessive genetic model (OR: 3.181, 95% CI: 1.204 to 8.408, P = 0.02), and a homozygote model (OR: 3.454, 95% CI: 1.286 to 9.278, P = 0.014). No significant association between them was detected in both a dominant genetic model (OR: 1.55, 95% CI: 0.99 to 2.42, P = 0.06) or a heterozygote genetic model (OR: 1.45, 95% CI: 0.90 to 2.33, P = 0.13).

Conclusion

The TNFα G308A gene polymorphism is associated with EH susceptibility.