The Interplay between Natural Selection and Susceptibility to Melanoma on Allele 374F of SLC45A2 Gene in a South European Population

We aimed to study the selective pressures interacting on SLC45A2 to investigate the interplay between selection and susceptibility to disease. Thus, we enrolled 500 volunteers from a geographically limited population (Basques from the North of Spain) and by resequencing the whole coding region and intron 5 of the 34 most and the 34 least pigmented individuals according to the reflectance distribution, we observed that the polymorphism Leu374Phe (L374F, rs16891982) was statistically associated with skin color variability within this sample. In particular, allele 374F was significantly more frequent among the individuals with lighter skin. Further genotyping an independent set of 558 individuals of a geographically wider population with known ancestry in the Spanish population also revealed that the frequency of L374F was significantly correlated with the incident UV radiation intensity. Selection tests suggest that allele 374F is being positively selected in South Europeans, thus indicating that depigmentation is an adaptive process. Interestingly, by genotyping 119 melanoma samples, we show that this variant is also associated with an increased susceptibility to melanoma in our populations. The ultimate driving force for this adaptation is unknown, but it is compatible with the vitamin D hypothesis. This shows that molecular evolution analysis can be used as a useful technology to predict phenotypic and biomedical consequences in humans.

Avaliação da presença de mutações de resistência no gene da NS5B e do prognóstico da infecção pelo HCV através da IL-28B em pacientes monoinfectados com HCV no RJ

Estima-se que a prevalência global da população mundial com hepatite C é de 3%. Pouco se sabe sobre a resposta ao tratamento com respeito à resistência viral. Algumas mutações no fragmento de 109 aminoácidos da NS5B são associadas com resistência ao interferon (IFN) e ribavirina (RBV). Estudos moleculares e clínicos identificaram fatores associados com o hospedeiro e vírus relacionados associada com a resposta ao tratamento, tal como o gene que codifica a IL-28B. Este estudo foi dividido em duas fases, cujos objetivos foram caracterizar a frequência de mutações que conferem resistência ao HCV e avaliar a relevância das mutações em pacientes Respondedores (R) ou Não Respondedores (NR) ao tratamento e caracterizar geneticamente as populações sobre polimorfismos genéticos nos SNPs da IL-28B em relação ao prognóstico da resposta ao tratamento. As amostras dos pacientes foram submetidas a testes de genotipagem e carga viral. As sequências geradas foram comparadas no BLAST e no banco de dados Los Alamos HCV. Realizamos o alinhamento das sequências homólogas e as mutações identificadas. Com base no genótipo e carga viral determinamos a classificação dos pacientes de acordo com a resposta à terapia. O DNA genômico foi isolado a partir de sangue periférico para a realização da tipagem de SNPs de IL-28B. A metodologia utilizada foi de PCR em tempo real utilizando sondas TaqMan SNP específico. A análise dos dados foi realizada utilizando GraphPad Prism com qui-quadrado, risco relativo (RR), Odds Ratio (OR) e intervalo de confiança de 95%, com um nível de significância de P <0,05. Foi encontrado na primeira fase deste estudo uma taxa significativa mutações associadas ao tratamento nas amostras estudadas. A prevalência de mutações associadas à resistência ao IFN e RBV bem como a novos medicamentos antivirais localizados no fragmento de 109 aminoácidos da NS5B foi examinado em 69 indivíduos infectados naïve no Rio de Janeiro, Brasil. Na segunda fase, as mutações foram clinicamente relevantes. Desde então, procuramos observar as diferenças entre melhor ou pior prognóstico de acordo com a imunogenética que mostrou diferenciação entre os grupos R e NR ao tratamento em relação ao prognóstico da resposta terapêutica. Quando as diferenças entre as sequências da NS5B e a resposta ao tratamento foram consideradas verificou-se que associada a mutação R254K, estava a C316N que poderia conduzir a uma não resposta à terapia no genótipo 1b. Os nossos dados também suportaram forte associação de IL-28B rs12979860, com elevada probabilidade de resposta à terapia de IFN + RBV. Nossos dados evidenciam a presença de pacientes virgens de tratamento que abrigam mutações de resistência previamente descritas na literatura. A análise dos fatores preditores de resposta virológica mostrou que a predição de boa resposta ou não ao tratamento e ainda da progressão da doença é dependente de uma importante interação entre a genética viral e a do hospedeiro. Fato este importante para que no momento de avaliação de diagnóstico e conduta terapêutica, o médico possa tomar medidas apropriadas para o tratamento de cada paciente individualmente independentemente do genótipo do HCV em questão.

Variantes no gene MBL2 codificador da Lectina Ligante de Manose (MBL): implicações na leishmaniose visceral humana

A leishmaniose visceral (LV) ou calazar é uma doença endêmica, crônica, grave e de alta letalidade se não tratada. Os estudos apontam a proteína Lectina Ligante de Manose (MBL), codificada pelo gene MBL2, como uma peça-chave na imunidade inata, dada a sua função no reconhecimento microbiano, na eliminação, inflamação e morte celular. Neste trabalho realizamos um estudo do tipo caso-controle que teve como objetivo investigar a associação entre variantes no gene MBL2 e a suscetibilidade à LV em indivíduos residentes em áreas endêmicas da Ilha de São Luís-MA. A amostra foi constituída por 322 indivíduos, sendo 161 casos com LV, não aparentados, de ambos os sexos, residentes em áreas endêmicas da doença na Ilha de São Luís e 161 controles saudáveis, não infectados e não aparentados da mesma região. A identificação dos casos de LV se deu por meio do contato constante com os principais hospitais e ambulatórios de referência para a doença na cidade. Também foram feitas buscas de pacientes com LV em ambiente domiciliar, a partir de registros da FUNASA-MA. A análise molecular consistiu na genotipagem de 6 variantes localizadas na região promotora [posições -550 (C>G), -221(G>C), +4(C>T)] e codificadora [códons 52 (C>T), 54 (G>A) e 57 (G>A)] do gene MBL2, através da reação em cadeia da polimerase e sequenciamento automático. A dosagem da proteína MBL no soro foi realizada pelo teste de ELISA. Verificamos que os fenótipos MBL dependem do conjunto de alelos presentes no gene MBL2, sendo nítido o efeito que as variantes defectivas causam nos níveis da proteína. Não encontramos diferença significativa entre casos e controles em relação à distribuição dos genótipos MBL2 e dos níveis séricos de MBL. As frequências alélicas das variantes exônicas na amostra total mostram que o alelo A é o mais comum (74,8%) e que os alelos defectivos (B, C e D) se encontram principalmente em heterozigose (36,6%), o que reforça a ideia de que alelos MBL2 defectivos são mantidos na população por conferirem vantagem seletiva aos heterozigotos. Em relação aos 3 principais polimorfismos existentes na região promotora, verificamos ser a variante -221G (Y) a mais frequente (88%) seguida de +4C (P) (73%) e de -550C (L) (67%). Identificamos oito haplótipos em MBL2 num total de 644 cromossomos avaliados, em 30 combinações diferentes, sendo HYPA e LYQA os mais frequentes e HYPD e HYPB os mais raros. Todos os portadores de combinações de haplótipos homozigotos para alelos defectivos apresentaram níveis séricos de MBL indetectáveis. Os genótipos LYQA/LYQA e HYPA/HYPA apresentaram as maiores concentrações médias de MBL no soro. A combinação entre SNPs no éxon 1 e na região promotora do gene MBL2 resulta em grande variação nas concentrações de MBL em indivíduos saudáveis. Consideramos que o conjunto de dados gerados é uma contribuição valiosa que poderá ser expandida para outros cenários.

The testis-specific apoptosis related gene TTL.6 underwent adaptive evolution in the lineage leading to humans

The TTL.6 gene is a member of the tubulin-tyrosine ligase (TTL) family involved in apoptosis and preferentially expressed in the testis. We sequenced the coding region and part of the introns of TTL.6 in world wide human populations and five representativ

Sequence evolution of the CCR5 chemokine receptor gene in primates

The chemokine receptor CCR5 can serve as a coreceptor for M-tropic HIV-1 infection and both M-tropic and T-tropic SIV infection. We sequenced the entire CCR5 gene from 10 nonhuman primates: Pongo pygmaeus, Hylobates leucogenys, Trachypithecus francoisi, Trachypithecus phayrei, Pygathrix nemaeus, Rhinopithecus roxellanae, Rhinopithecus bieti, Rhinopithecus avunculus, Macaca assamensis, and Macaca arctoides. When compared with CCR5 sequences from humans and other primates, our results demonstrate that:(1) nucleotide and amino acid sequences of CCR5 among primates are highly homologous, with variations slightly concentrated on the amino and carboxyl termini; and (2) site Asp13, which is critical for CD4-independent binding of SIV gp120 to Macaca mulatta CCR5, was also present in all other nonhuman primates tested here, suggesting that those nonhuman primate CCR5s might also bind SIV gp120 without the presence of CD4. The topologies of CCR5 gene trees constructed here conflict with the putative opinion that the snub-nosed langurs compose a monophyletic group, suggesting that the CCR5 gene may not be a good genetic marker for low-level phylogenetic analysis. The evolutionary rate of CCR5 was calculated, and our results suggest a slowdown in primates after they diverged from rodents. The synonymous mutation rate of CCR5 in primates is constant, about 1.1 x 10(-9) synonymous mutations per site per year. Comparisons of K-a and K-s suggest that the CCR5 genes have undergone negative or purifying selection. K-a/K-s ratios from cercopithecines and colobines are significantly different, implying that selective pressures have played different roles in the two lineages.

Evolution of the DRD2 gene haplotype and its association with alcoholism in Mexican Americans

The human D2 dopamine receptor gene (DRD2) plays a central role in the neuromodulation of appetitive behaviors and is implicated in having a possible role in susceptibility to alcoholism. We genotyped an SNP in DRD2 Exon 8 in 251 nonalcoholic, unrelated, healthy controls and 200 alcoholic Mexican Americans. The DRD2 haplotypes were analyzed using the Exon 8 genotype in combination with five other SNP genotypes, which were obtained from our previous study. The ancestral origins of the DRD2 polymorphisms have been determined by sequencing the homologous region in other higher primates. Twenty DRD2 haplotypes, defined as H1 to H20 based on their frequency from high to low, were obtained in this major minority population. The ancestral haplotype "I-132-G-C-G-A1" and two one-step mutation haplotypes were absent in our study population. The haplotype H1, "I-B1-T-C-A-A1", with the highest frequency in the population, is a three-step mutation from the ancestral form. The first five or eight major haplotypes make up 87% or 95% of the entire population, respectively. The prevalence of the haplotype H1+ (H1/H1 and H1/Hn genotypes) is significantly higher in alcoholics and alcoholic subgroups, including early onset drinkers and benders, than in their respective control groups. The Promoter -141C allele is in linkage disequilibrium (LD) with five other loci in the nonalcoholic group, but not in the alcoholic group. All of the other five loci are in LD in both the alcoholic and control groups. The DRD2 TaqI B allele is in complete LD with the allele located in intron 6. Five SNPs, Promoter -141C, TaqI B (or Intron 6), Exon 7, Exon 8, and TaqI A, are sufficient to define the DRD2 haplotypes in Mexican Americans. Our data indicate that the DRD2 haplotypes are associated with alcoholism in Mexican Americans. (c) 2005 Elsevier Inc. All rights reserved.

Evolution of cd59 gene in mammals

The CD59-coding sequences were obtained from 5 mammals by PCR and BLAST, and combined with the available sequences in GenBank, the nucleotide substitution rates of mammalian cd59 were calculated. Results of synonymous and nonsynonymous substitution rates

Motilin and ghrelin gene experienced episodic evolution during primitive placental mammal evolution

Motilin and ghrelin, members of a structure-function-related hormone family, play important roles in gastrointestinal function, regulation of energy homeostasis and growth hormone secretion. We observed episodic evolution in both of their prehormone gene sequences during primitive placental mammal evolution, during which most of the nonsynonymous changes result in radical substitution. Of note, a functional obestatin hormone might have only originated after this episodic evolution event. Early in placental mammal evolution, a series of biology complexities evolved. At the same time the motilin and ghrelin prehormone genes, which play important roles in several of these processes, experienced episodic evolution with dramatic changes in their coding sequences. These observations suggest that some of the lineage-specific physiological adaptations are due to episodic evolution of the motilin and ghrelin genes.

The role of MC1R gene in buffalo coat color

Melanocortin-1 receptor (MC1R) plays a major role in pigmentation in many species. To investigate if the MC1R gene is associated with coat color in water buffalo, the coding region of MC1R gene of 216 buffalo samples was sequenced, which included 49 black river buffalo (Murrah and Nili-Ravi), 136 swamp buffalo (Dehong, Diandongnan, Dechang, Guizhou, and Xilin) with white and gray body, and 31 hybrid offspring of river buffalo Nili-Ravi (or Murrah) and swamp buffalo. Among the three variation sites found, SNP684 was synonymous, while SNP310 and SNP384 were nonsynonymous, leading to p.S104G and p.I128M changes, respectively. Only Individuals carrying homozygote E-BR/E-BR were black. The genotype and phenotype analysis of the hybrid offspring of black river buffalo and gray swamp buffalo further revealed that the river buffalo type allele E-BR or the allele carrying the amino acid p.104S was important for the full function of MC1R. The in silico functional analysis showed that the amino acid substitutions p.G104S and p.M128I had significant impact on the function of MC1R. Above results indicate that the allele E-BR or the allele carrying the amino acid p.104S was associated with the black coat color in buffalo.

Origin and evolution of new exons in the rodent zinc finger protein 39 gene

The origin of new structures and functions is an important process in evolution. In the past decades, we have obtained some preliminary knowledge of the origin and evolution of new genes. However, as the basic unit of genes, the origin and evolution of exons remain unclear. Because young exons retain the footprints of origination, they can be good materials for studying origin and evolution of new exons. In this paper, we report two young exons in a zinc finger protein gene of rodents. Since they are unique sequences in mouse and rat genome and no homologous sequences were found in the orthologous genes of human and pig, the young exons might originate after the divergence of primates and rodents through exonization of intronic sequences. Strong positive selection was detected in the new exons between mouse and rat, suggesting that these exons have undergone significant functional divergence after the separation of the two species. On the other hand, population genetics data of mouse demonstrate that the new exons have been subject to functional constraint, indicating an important function of the new exons in mouse. Functional analyses suggest that these new exons encode a nuclear localization signal peptide, which may mediate new ways of nuclear protein transport. To our knowledge, this is the first example of the origin and evolution of young exons.

Repetitive element-mediated recombination as a mechanism for new gene origination in Drosophila

Previous studies of repetitive elements (REs) have implicated a mechanistic role in generating new chimerical genes. Such examples are consistent with the classic model for exon shuffling, which relies on non-homologous recombination. However, recent data

A de novo originated gene depresses budding yeast mating pathway and is repressed by the protein encoded by its antisense strand

Recent transcription profiling studies have revealed an unexpectedly large proportion of antisense transcripts in eukaryotic genomes. These antisense genes seem to regulate gene expression by interacting with sense genes. Previous studies have focused on the non-coding antisense genes, but the possible regulatory role of the antisense protein is poorly understood. In this study, we found that a protein encoded by the antisense gene ADF1 acts as a transcription suppressor, regulating the expression of sense gene MDF1 in Saccharomyces cerevisiae. Based on the evolutionary, genetic, cytological and biochemical evidence, we show that the protein-coding sense gene MDF1 most likely originated de novo from a previously non-coding sequence and can significantly suppress the mating efficiency of baker's yeast in rich medium by binding MAT alpha 2 and thus promote vegetative growth. These results shed new light on several important issues, including a new sense-antisense interaction mechanism, the de novo origination of a functional gene, and the regulation of yeast mating pathway.

Effects of exogenous double-stranded RNA on the basonuclin gene expression in mouse oocytes

In plants and less-advanced animal species, such as C.elegans, introduction of exogenous double-stranded RNA (dsRNA) into cells would trigger degradation of the mRNA with homologous sequence and interfere with the endogenous gene expression. It might represent an ancient anti-virus response which could prevent the mutation in the genome that was caused by virus infection or mobile DNA elements insertion. This phenomenon was named RNA interference, or RNAi. In this study, RNAi was used to investigate the function of basonuclin gene during oogenesis. Microinjection of dsRNA directed towards basonuclin into mouse germinal-vesicle-intact (GV) oocytes brought down the abundance of the cognate mRNA effectively in a time- and concentration-dependent manner. This reduction effect was sequence-specific and showed no negative effect on other non-homologous gene expression in oocytes, which indicated that dsRNA can recognize and cause the degradation of the transcriptional products of endogenous basonuclin gene in a sequence-specific manner. Immunofluorescence results showed that RNAi could reduce the concentration of basonuclin protein to some extent, but the effect was less efficient than the dsRNA targeting towards tPA and cMos which was also expressed in oocytes. This result might be due to the long half life of basonuclin protein in oocytes and the short reaction time which was posed by the limited life span of GV oocytes cultured in vitro. In summary, dsRNA could inhibit the expression of the cognate gene in oocytes at both mRNA and protein levels. The effect was similar to Knock-out technique which was based on homologous recombination. Furthermore, hairpin-style dsRNA targeting basonuclin gene could be produced by transcription from a recombinant plasmid and worked efficiently to deplete the cognate mRNA in oocytes. This finding offered a new way to study the function of basonuclin in the early stage of oogenesis by infection of primordial oocytes with the plasmid expressing hairpin-style basonuclin dsRNA.

Origin and evolution of the RIG-I like RNA helicase gene family

Background: The DExD/H domain containing RNA helicases such as retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5) are key cytosolic pattern recognition receptors (PRRs) for detecting nucleotide pathogen associated molecular patterns (PAMPs) of invading viruses. The RIG-I and MDA5 proteins differentially recognise conserved PAMPs in double stranded or single stranded viral RNA molecules, leading to activation of the interferon system in vertebrates. They share three core protein domains including a RNA helicase domain near the C terminus (HELICc), one or more caspase activation and recruitment domains (CARDs) and an ATP dependent DExD/H domain. The RIG-I/MDA5 directed interferon response is negatively regulated by laboratory of genetics and physiology 2 (LGP2) and is believed to be controlled by the mitochondria antiviral signalling protein (MAVS), a CARD containing protein associated with mitochondria. Results: The DExD/H containing RNA helicases including RIG-I, MDA5 and LGP2 were analysed in silico in a wide spectrum of invertebrate and vertebrate genomes. The gene synteny of MDA5 and LGP2 is well conserved among vertebrates whilst conservation of the gene synteny of RIG-I is less apparent. Invertebrate homologues had a closer phylogenetic relationship with the vertebrate RIG-Is than the MDA5/LGP2 molecules, suggesting the RIG-I homologues may have emerged earlier in evolution, possibly prior to the appearance of vertebrates. Our data suggest that the RIG-I like helicases possibly originated from three distinct genes coding for the core domains including the HELICc, CARD and ATP dependent DExD/H domains through gene fusion and gene/domain duplication. Furthermore, presence of domains similar to a prokaryotic DNA restriction enzyme III domain (Res III), and a zinc finger domain of transcription factor (TF) IIS have been detected by bioinformatic analysis. Conclusion: The RIG-I/MDA5 viral surveillance system is conserved in vertebrates. The RIG-I like helicase family appears to have evolved from a common ancestor that originated from genes encoding different core functional domains. Diversification of core functional domains might be fundamental to their functional divergence in terms of recognition of different viral PAMPs.

Molecular characterization and expression profiles in response to bacterial infection of Chinese soft-shelled turtle interleukin-8 (IL-8), the first reptilian chemokine gene

In this study, an IL-8 homologue has been cloned and identified from a reptile, Chinese soft-shelled turtle for the first time. The full-length cDNA of turtle IL-8 was 1188 bp and contained a 312 bp open reading frame (ORF) coding for a protein of 104 amino acids. The chemokine CXC domain, which contained Glu-Leu-Arg (ELR) motif and four cysteine residues, was well conserved in turtle IL-8. The 4924 bp genomic DNA of turtle IL-8 contained four exons and three introns. Phylogenetic analysis showed that the amino acid sequence of turtle IL-8 clustered together with birds. RT-PCR analysis showed that turtle IL-8 mRNA was constitutively expressed liver, spleen, kidney, heart, blood and intestine tissues of control turtles. Real-time quantitative PCR analysis further indicated that the turtle IL-8 mRNA expression was apparent in various tissues at 8 h and up-regulated significantly during 8 h-7 d after Aeromonas hydrophila infection. The present studies will help us to understand the evolution of IL-8 molecule and the inflammatory response mechanism in reptiles. (C) 2009 Elsevier Ltd. All rights reserved.