Association of IRGM Gene Mutations with Inflammatory Bowel Disease in the Indian Population

Background

Mutations in the IRGM gene have been associated with Crohn''s disease in several populations but have not been explored in Indian patients with this disease. This study examined the association of IRGM mutations with ulcerative colitis and Crohn''s disease in Indian patients with inflammatory bowel disease.

Methods

The IRGM gene was amplified in four segments and Sanger-sequenced in 101 participants (42 Crohn''s disease, 39 ulcerative colitis, and 20 healthy controls). Ten single nucleotide polymorphisms (SNP) were genotyped in 1200 participants (352 Crohn''s disease, 400 ulcerative colitis, and 448 healthy controls) using Sequenom MassARRAY iPLEX. Disease associations were evaluated for each of the ten SNPs.

Results

Thirty one mutations were identified in the IRGM gene, of which two had not hitherto been reported (150226250- ss947429272 & 150227858- ss947429273). Ten SNPs (6 from the above and 4 from the literature) were evaluated. Significant associations with Crohn''s disease were noted with the T allele of rs1000113 (OR 1.46, 95% CI 1.12–1.90), T allele of rs9637876 (OR 1.25, 95% CI 1.005–1.561) and C allele of rs 13361189 (OR 1.33, 95% CI 1.07–1.669). Two SNPs – rs11747270 and rs180802994 – did not exhibit Hardy-Weinberg equilibrium but were associated with both Crohn''s disease and ulcerative colitis in this population. The remaining SNPs did not show significant associations with either Crohn''s disease or ulcerative colitis.

Conclusions

Association of IRGM gene SNPs with Crohn''s disease is reported for the first time in Indian patients. We also report, for the first time, an association of rs 9637876 in the IRGM gene with Crohn''s disease.     

Autophagy Gene Variant IRGM ?261T Contributes to Protection from Tuberculosis Caused by Mycobacterium tuberculosis but Not by M. africanum Strains

The human immunity-related GTPase M (IRGM) has been shown to be critically involved in regulating autophagy as a means of disposing cytosolic cellular structures and of reducing the growth of intracellular pathogens in vitro. This includes Mycobacterium tuberculosis, which is in agreement with findings indicating that M. tuberculosis translocates from the phagolysosome into the cytosol of infected cells, where it becomes exposed to autophagy. To test whether IRGM plays a role in human infection, we studied IRGM gene variants in 2010 patients with pulmonary tuberculosis (TB) and 2346 unaffected controls. Mycobacterial clades were classified by spoligotyping, IS6110 fingerprinting and genotyping of the pks1/15 deletion. The IRGM genotype −261TT was negatively associated with TB caused by M. tuberculosis (OR 0.66, CI 0.52–0.84, P_nominal 0.0009, P_corrected 0.0045) and not with TB caused by M. africanum or M. bovis (OR 0.95, CI 0.70–1.30. P 0.8). Further stratification for mycobacterial clades revealed that the protective effect applied only to M. tuberculosis strains with a damaged pks1/15 gene which is characteristic for the Euro-American (EUAM) subgroup of M. tuberculosis (OR 0.63, CI 0.49–0.81, P_nominal 0.0004, P_corrected 0.0019). Our results, including those of luciferase reporter gene assays with the IRGM variants −261C and −261T, suggest a role for IRGM and autophagy in protection of humans against natural infection with M. tuberculosis EUAM clades. Moreover, they support in vitro findings indicating that TB lineages capable of producing a distinct mycobacterial phenolic glycolipid that occurs exclusively in strains with an intact pks1/15 gene inhibit innate immune responses in which IRGM contributes to the control of autophagy. Finally, they raise the possibility that the increased frequency of the IRGM −261TT genotype may have contributed to the establishment of M. africanum as a pathogen in the West African population.     

Primate immunoglobulin heavy chain constant gamma genes: an hypothesis of their evolution

Four Immunoglobulin gamma heavy chain isotypes are present in humans; the true phylogenetic relationship between the genes are not known because of the complex concerted evolution of the Ig multigene locus. Here we present data obtained from Southern blot analysis of the gamma genes in several primate species including prosimians (Lemur catta), New World Monkeys (Saguinus oedipus) and Old World Monkeys (Cercopitecus aethiops andMacaca fascicularis). Our data show the presence of a single IgG gene inLemur and probably inS.l oedipus, and of multiple genes in the two Cercopithecinae. These findings suggest that a single IgG gene was present in the ancestor of primates: we suppose that this IgG gene was later duplicated several times during primate evolution.     

Single-cell genomics of a rare environmental alphaproteobacterium provides unique insights into Rickettsiaceae evolution

The bacterial family Rickettsiaceae includes a group of well-known etiological agents of many human and vertebrate diseases, including epidemic typhus-causing pathogen Rickettsia prowazekii. Owing to their medical relevance, rickettsiae have attracted a great deal of attention and their host-pathogen interactions have been thoroughly investigated. All known members display obligate intracellular lifestyles, and the best-studied genera, Rickettsia and Orientia, include species that are hosted by terrestrial arthropods. Their obligate intracellular lifestyle and host adaptation is reflected in the small size of their genomes, a general feature shared with all other families of the Rickettsiales. Yet, despite that the Rickettsiaceae and other Rickettsiales families have been extensively studied for decades, many details of the origin and evolution of their obligate host-association remain elusive. Here we report the discovery and single-cell sequencing of ‘Candidatus Arcanobacter lacustris'', a rare environmental alphaproteobacterium that was sampled from Damariscotta Lake that represents a deeply rooting sister lineage of the Rickettsiaceae. Intriguingly, phylogenomic and comparative analysis of the partial ‘Candidatus Arcanobacter lacustris'' genome revealed the presence chemotaxis genes and vertically inherited flagellar genes, a novelty in sequenced Rickettsiaceae, as well as several host-associated features. This finding suggests that the ancestor of the Rickettsiaceae might have had a facultative intracellular lifestyle. Our study underlines the efficacy of single-cell genomics for studying microbial diversity and evolution in general, and for rare microbial cells in particular.     

Predicting the Proteins of Angomonas deanei,Strigomonas culicis and Their Respective Endosymbionts Reveals New Aspects of the Trypanosomatidae Family

Endosymbiont-bearing trypanosomatids have been considered excellent models for the study of cell evolution because the host protozoan co-evolves with an intracellular bacterium in a mutualistic relationship. Such protozoa inhabit a single invertebrate host during their entire life cycle and exhibit special characteristics that group them in a particular phylogenetic cluster of the Trypanosomatidae family, thus classified as monoxenics. In an effort to better understand such symbiotic association, we used DNA pyrosequencing and a reference-guided assembly to generate reads that predicted 16,960 and 12,162 open reading frames (ORFs) in two symbiont-bearing trypanosomatids, Angomonas deanei (previously named as Crithidia deanei) and Strigomonas culicis (first known as Blastocrithidia culicis), respectively. Identification of each ORF was based primarily on TriTrypDB using tblastn, and each ORF was confirmed by employing getorf from EMBOSS and Newbler 2.6 when necessary. The monoxenic organisms revealed conserved housekeeping functions when compared to other trypanosomatids, especially compared with Leishmania major. However, major differences were found in ORFs corresponding to the cytoskeleton, the kinetoplast, and the paraflagellar structure. The monoxenic organisms also contain a large number of genes for cytosolic calpain-like and surface gp63 metalloproteases and a reduced number of compartmentalized cysteine proteases in comparison to other TriTryp organisms, reflecting adaptations to the presence of the symbiont. The assembled bacterial endosymbiont sequences exhibit a high A+T content with a total of 787 and 769 ORFs for the Angomonas deanei and Strigomonas culicis endosymbionts, respectively, and indicate that these organisms hold a common ancestor related to the Alcaligenaceae family. Importantly, both symbionts contain enzymes that complement essential host cell biosynthetic pathways, such as those for amino acid, lipid and purine/pyrimidine metabolism. These findings increase our understanding of the intricate symbiotic relationship between the bacterium and the trypanosomatid host and provide clues to better understand eukaryotic cell evolution.     

Bovine carboxylesterases: Evidence for two CES1 and five families of CES genes on chromosome 18

Predicted bovine carboxylesterase (CES) protein and gene sequences were derived from bovine (Bos taurus) genomic sequence data. Two bovine CES1 genes (CES1.1 and CES1.2) were located on chromosome 18 encoding amino acid sequences that were 81% identical. Two forms of CES1.2 were also observed apparently caused by an indel polymorphism encoded at the C-terminus end. Two CES gene clusters were observed on chromosome 18: CES5–CES1.1–CES1.2 and CES2–CES3–CES6. Bovine CES1, CES2, CES3, CES5 and CES6 shared 39–45% identity with each other, but showed 71–76% identity with each of the five corresponding human CES family members. Phylogeny studies indicated that bovine CES genes originated from five ancestral gene duplication events which predated the eutherian mammalian common ancestor. In addition, a subsequent CES1 gene duplication event is proposed during mammalian evolution prior to the appearance of the Bovidae common ancestor ~ 20 MY ago.     

Expression profiles reveal parallel evolution of epistatic interactions involving the CRP regulon in Escherichia coli

The extent and nature of epistatic interactions between mutations are issues of fundamental importance in evolutionary biology. However, they are difficult to study and their influence on adaptation remains poorly understood. Here, we use a systems-level approach to examine epistatic interactions that arose during the evolution of Escherichia coli in a defined environment. We used expression arrays to compare the effect on global patterns of gene expression of deleting a central regulatory gene, crp. Effects were measured in two lineages that had independently evolved for 20,000 generations and in their common ancestor. We found that deleting crp had a much more dramatic effect on the expression profile of the two evolved lines than on the ancestor. Because the sequence of the crp gene was unchanged during evolution, these differences indicate epistatic interactions between crp and mutations at other loci that accumulated during evolution. Moreover, a striking degree of parallelism was observed between the two independently evolved lines; 115 genes that were not crp-dependent in the ancestor became dependent on crp in both evolved lines. An analysis of changes in crp dependence of well-characterized regulons identified a number of regulatory genes as candidates for harboring beneficial mutations that could account for these parallel expression changes. Mutations within three of these genes have previously been found and shown to contribute to fitness. Overall, these findings indicate that epistasis has been important in the adaptive evolution of these lines, and they provide new insight into the types of genetic changes through which epistasis can evolve. More generally, we demonstrate that expression profiles can be profitably used to investigate epistatic interactions.     

The evolution of ADAM gene family in eukaryotes

The ADAM (A Disintegrin And Metalloprotease) gene family encodes proteins with adhesion and proteolytic functions. ADAM proteins are associated with diseases like cancers. Twenty ADAM genes have been identified in humans. However, little is known about the evolution of the family. We analyzed the repertoire of ADAM genes in a vast number of eukaryotic genomes to clarify the main gene copy number expansions. For the first time, we provide compelling evidence that early-branching green algae (Mamiellophyceae) have ADAM genes, suggesting that they originated in the last common ancestor of eukaryotes, before the split of plants, fungi and animals. The ADAM family expanded in early metazoans, with the most significative gene expansion happening during the first steps of vertebrate evolution. We concluded that most of mammal ADAM diversity can be explained by gene duplications in early bone fish. Our data suggest that ADAM genes were lost early in green plant evolution.     

Identification and characterization of complement factor H in Branchiostoma belcheri

Complement factor H (CFH) is an essential regulator of the complement system and plays very important roles in animal innate immunity. Although the complement system of amphioxus has been extensively studied, the expression in amphioxus and evolution of CFH gene remain unknown. In this study, we identified and characterized an amphioxus (Branchiostoma belcheri) CFH gene (designated as AmphiCFH). Our results showed that the full-length cDNA of AmphiCFH gene consists of 1295 bp nucleotides containing an 855 bp open reading frame (ORF) that was predicted to encode a 284 amino acid protein. The putative AmphiCFH protein possessed the characteristic of the CFH protein family, including typical CCP (complement control protein) domain. Real-time PCR analysis showed that the AmphiCFH was ubiquitously and differentially expressed in five investigated tissues (intestine, gills, notochord, muscles, and hepatic cecum). The expression level of the AmphiCFH gene was induced upon lipopolysaccharide stimulation, indicating that the AmphiCFH gene might be involved in innate immunity. In addition, phylogenetic analysis showed that the AmphiCFH gene was located between that of invertebrates and vertebrates, suggesting that the AmphiCFH gene is a member of the CFH gene family. In conclusion, our findings provided an insight into animal innate immunity and evolution of the CFH gene family.     

Seven New Complete Plastome Sequences Reveal Rampant Independent Loss of the ndh Gene Family across Orchids and Associated Instability of the Inverted Repeat/Small Single-Copy Region Boundaries

Earlier research has revealed that the ndh loci have been pseudogenized, truncated, or deleted from most orchid plastomes sequenced to date, including in all available plastomes of the two most species-rich subfamilies, Orchidoideae and Epidendroideae. This study sought to resolve deeper-level phylogenetic relationships among major orchid groups and to refine the history of gene loss in the ndh loci across orchids. The complete plastomes of seven orchids, Oncidium sphacelatum (Epidendroideae), Masdevallia coccinea (Epidendroideae), Sobralia callosa (Epidendroideae), Sobralia aff. bouchei (Epidendroideae), Elleanthus sodiroi (Epidendroideae), Paphiopedilum armeniacum (Cypripedioideae), and Phragmipedium longifolium (Cypripedioideae) were sequenced and analyzed in conjunction with all other available orchid and monocot plastomes. Most ndh loci were found to be pseudogenized or lost in Oncidium, Paphiopedilum and Phragmipedium, but surprisingly, all ndh loci were found to retain full, intact reading frames in Sobralia, Elleanthus and Masdevallia. Character mapping suggests that the ndh genes were present in the common ancestor of orchids but have experienced independent, significant losses at least eight times across four subfamilies. In addition, ndhF gene loss was correlated with shifts in the position of the junction of the inverted repeat (IR) and small single-copy (SSC) regions. The Orchidaceae have unprecedented levels of homoplasy in ndh gene presence/absence, which may be correlated in part with the unusual life history of orchids. These results also suggest that ndhF plays a role in IR/SSC junction stability.     

Identification of conserved C2H2 zinc-finger gene families in the Bilateria

Background:Identification of orthologous relationships between genes from widely divergent taxa allows partial reconstruction of the gene complement of ancestral genomes. C2H2 zinc-finger genes are one of the largest and most complex gene superfamilies in metazoan genomes, with hundreds of members in the human genome. Here we analyze C2H2 zinc-finger genes from three taxa - Drosophila, Caenorhabditis elegans and human - from which near-complete genome sequence data are available.Results:Our analyses conclusively identify 39 families of genes, of which 38 can be defined as orthology groups in that they are descended from single ancestral genes in the common ancestor of Drosophila, C. elegans and humans.Conclusions:On the basis of current metazoan phylogeny, these 39 groups represent the minimum complement of C2H2 zinc-finger genes present in the genome of the bilaterian common ancestor.     

Evolutionary and functional diversity of green fluorescent proteins in cephalochordates

Green fluorescent protein (GFP) has been widely used as a molecular marker in modern biological research. Before the recent report of one GFP gene in Branchiostoma floridae, GFP family members were cloned only from other two groups of species: Cnidaria and Copepoda. Here we describe the complete GFP gene repertoire of B. floridae which includes 13 functional genes and 2 pseudogenes, representing the largest GFP family found so far. Coupling with nine other GFP sequences from another two species of genus Branchiostoma and the sequences from Cnidaria and Copepoda, we made a deep-level phylogenetic analysis for GFP genes in cephalochordates and found: 1) GFP genes have experienced a divergent evolution in cephalochordates; 2) all amphioxus GFP genes form four main clades on the tree which had diverged before the radiation of the last common ancestor of all extant cephalochordates; 3) GFP genes in amphioxus shared a common ancestor with that in Copepoda rather than being derived from horizontal gene transfer, which indicates that our ancestor was derived from a fluorescent organism and lost this ability after its separation from Cephalochordata, and also makes GFP a rare gene which has a rather unusual evolutionary path. In addition, we also provided evidence indicating that GFP genes have evolved divergent functions by specializing their expression profile, and different fluorescent spectra by changing their emission peaks. These findings spark two interesting issues: what are GFP in vivo functions in cephalochordates and why they are lost in other examined deuterostomes?     

Evolution of the ABPA Subunit of Androgen-Binding Protein Expressed in the Submaxillary Glands in New and Old World Rodent Taxa

The salivary androgen-binding proteins (ABPs) are members of the secretoglobin gene family present in mammals. Each ABP is a heterodimer assembled as an ABPA subunit encoded by an Abpa gene and linked by disulfide bridges to an ABPBG subunit encoded by an Abpbg gene. The ABP dimers are secreted into the saliva of mice and then transferred to the pelage after grooming and subsequently to the environment allowing an animal to mark territory with a biochemical signal. The putative role of the mouse salivary ABPs is that of pheromones mediating mate selection resulting in assortative mating in the Mus musculus species complex. We focused on comparing patterns of molecular evolution between the Abpa genes expressed in the submaxillary glands of species of New World and Old World muroids. We found that in both sets of rodents the Abpa genes expressed in the submaxillary glands appear to be evolving under a similar evolutionary regime, with relatively high nonsynonymous substitution rates, suggesting that ABP might play a similar biological role in both systems. Thus, ABP could be involved with mate recognition and species isolation in New World as well as Old World muroids.     

Molecular Evolution of the Odorant and Gustatory Receptor Genes in Lepidopteran Insects: Implications for Their Adaptation and Speciation

Lepidoptera (comprised of butterflies and moths) is one of the largest groups of insects, including more than 160,000 described species. Chemoreception plays important roles in the adaptation of these species to a wide range of niches, e.g., plant hosts, egg-laying sites, and mates. This study investigated the molecular evolution of the lepidopteran odorant (Or) and gustatory receptor (Gr) genes using recently identified genes from Bombyx mori, Danaus plexippus, Heliconius melpomene, Plutella xylostella, Heliothis virescens, Manduca sexta, Cydia pomonella, and Spodoptera littoralis. A limited number of cases of large lineage-specific gene expansion are observed (except in the P. xylostella lineage), possibly due to selection against tandem gene duplication. There has been strong purifying selection during the evolution of both lepidopteran odorant and gustatory genes, as shown by the low ω values estimated through CodeML analysis, ranging from 0.0093 to 0.3926. However, purifying selection has been relaxed on some amino acid sites in these receptors, leading to sequence divergence, which is a precursor of positive selection on these sequences. Signatures of positive selection were detected only in a few loci from the lineage-specific analysis. Estimation of gene gains and losses suggests that the common ancestor of the Lepidoptera had fewer Or genes compared to extant species and an even more reduced number of Gr genes, particularly within the bitter receptor clade. Multiple gene gains and a few gene losses occurred during the evolution of Lepidoptera. Gene family expansion may be associated with the adaptation of lepidopteran species to plant hosts, especially after angiosperm radiation. Phylogenetic analysis of the moth sex pheromone receptor genes suggested that chromosomal translocations have occurred several times. New sex pheromone receptors have arisen through tandem gene duplication. Positive selection was detected at some amino acid sites predicted to be in the extracellular and transmembrane regions of the newly duplicated genes, which might be associated with the evolution of the new pheromone receptors.     

Expansion of a novel endogenous retrovirus throughout the pericentromeres of modern humans

BackgroundApproximately 8% of the human genome consists of sequences of retroviral origin, a result of ancestral infections of the germ line over millions of years of evolution. The most recent of these infections is attributed to members of the human endogenous retrovirus type-K (HERV-K) (HML-2) family. We recently reported that a previously undetected, large group of HERV-K (HML-2) proviruses, which are descendants of the ancestral K111 infection, are spread throughout human centromeres.ResultsStudying the genomes of certain cell lines and the DNA of healthy individuals that seemingly lack K111, we discover new HERV-K (HML-2) members hidden in pericentromeres of several human chromosomes. All are related through a common ancestor, termed K222, which is a virus that infected the germ line approximately 25 million years ago. K222 exists as a single copy in the genomes of baboons and high order primates, but not New World monkeys, suggesting that progenitor K222 infected the primate germ line after the split between New and Old World monkeys. K222 exists in modern humans at multiple loci spread across the pericentromeres of nine chromosomes, indicating it was amplified during the evolution of modern humans.ConclusionsCopying of K222 may have occurred through recombination of the pericentromeres of different chromosomes during human evolution. Evidence of recombination between K111 and K222 suggests that these retroviral sequences have been templates for frequent cross-over events during the process of centromere recombination in humans.     

The SKINT1-Like Gene Is Inactivated in Hominoids But Not in All Primate Species: Implications for the Origin of Dendritic Epidermal T Cells

Dendritic epidermal T cells, which express an invariant Vγ5Vδ1 T-cell receptor and account for 95% of all resident T cells in the mouse epidermis, play a critical role in skin immune surveillance. These γδ T cells are generated by positive selection in the fetal thymus, after which they migrate to the skin. The development of dendritic epidermal T cells is critically dependent on the Skint1 gene expressed specifically in keratinocytes and thymic epithelial cells, suggesting an indispensable role for Skint1 in the selection machinery for specific intraepithelial lymphocytes. Phylogenetically, rodents have functional SKINT1 molecules, but humans and chimpanzees have a SKINT1-like (SKINT1L) gene with multiple inactivating mutations. In the present study, we analyzed SKINT1L sequences in representative primate species and found that all hominoid species have a common inactivating mutation, but that Old World monkeys such as olive baboons, green monkeys, cynomolgus macaques and rhesus macaques have apparently functional SKINT1L sequences, indicating that SKINT1L was inactivated in a common ancestor of hominoids. Interestingly, the epidermis of cynomolgus macaques contained a population of dendritic-shaped γδ T cells expressing a semi-invariant Vγ10/Vδ1 T-cell receptor. However, this population of macaque T cells differed from rodent dendritic epidermal T cells in that their Vγ10/Vδ1 T-cell receptors displayed junctional diversity and expression of Vγ10 was not epidermis-specific. Therefore, macaques do not appear to have rodent-type dendritic epidermal T cells despite having apparently functional SKINT1L. Comprehensive bioinformatics analysis indicates that SKINT1L emerged in an ancestor of placental mammals but was inactivated or lost multiple times in mammalian evolution and that Skint1 arose by gene duplication in a rodent lineage, suggesting that authentic dendritic epidermal T cells are presumably unique to rodents.