HLA-G UTR Haplotype Conservation in the Malian Population: Association with Soluble HLA-G

The HLA-G molecule plays an important role in immunomodulation. In a previous study carried out on a southern French population our team showed that HLA-G haplotypes, defined by SNPs in the coding region and specific SNPs located in 5′URR and 3′UTR regulatory regions, are associated with differential soluble HLA-G expression (sHLA-G). Furthermore, the structure of these HLA-G haplotypes appears to be conserved in geographically distant populations.The aim of our study is to confirm these expectations in a sub-Saharan African population and to explore additional factors, such as HLA-A alleles, that might influence sHLA-G expression.DNA and plasma samples were collected from 229 Malians; HLA-G and HLA-A genotyping were respectively performed by the Snap Shot® method and by Luminex™ technology. sHLA-G dosage was performed using an ELISA kit. HLA-G and HLA-A allelic and haplotypic frequencies were estimated using an EM algorithm from the Gene[Rate] program. Associations between genetic and non genetic parameters with sHLA-G were performed using a non-parametric test with GRAPH PAD Prism 5.Our results reveal a good conservation of the HLA-G UTR haplotype structure in populations with different origins and demographic histories. These UTR haplotypes appear to be involved in different sHLA-G expression patterns. Specifically, the UTR-2 haplotype was associated with low sHLA-G levels, displaying a dominant negative effect. Furthermore, an allelic effect of both HLA-G and HLA-A, as well as non genetic parameters, such as age and gender possibly linked to osteogenesis and sexual hormones, also seem to be involved in the modulation of sHLA-G.These data suggest that further investigation in larger cohorts and in populations from various ethnical backgrounds is necessary not only to detect new functional polymorphism in HLA-G regulatory regions, but also to reveal the extent of biological phenomena that influence sHLA-G secretion and this might therefore have an impact on transplantation practice.     

Polymorphic Sites at the 3’ Untranslated Region of the HLA-G Gene Are Associated with Differential hla-g Soluble Levels in the Brazilian and French Population

HLA-G molecule has well-recognized tolerogenic properties, and the encoding gene shows lower frequency of polymorphism at the coding region but higher variability at regulatory 5’ and 3’ untranslated (3’UTR) regions. At least three 3’UTR polymorphic sites have been associated with HLA-G mRNA regulation, including the 14 base pair (14bp) Insertion/Deletion, +3142C-G and +3187A-G. We studied the association of polymorphic sites at 3’UTR (sequencing analysis, encompassing the 14bp Ins-Del/+3003T-C/+3010C-G/+3027C-A/+3035C-T/+3142C-G/+3187A-G/+3196C-G polymorphic sites) with plasma soluble HLA-G levels (sHLA-G, detected by ELISA) in 187 French and 153 Brazilian healthy individuals. Allele and genotype frequencies were closely similar in both populations; however, Brazilians showed a higher HLA-G 3’UTR haplotype diversity. Considering sHLA-G levels in both populations altogether, individuals presenting 14bp Del/Del showed higher levels compared to 14bpIns/Ins genotype (P <0.05); those presenting +3010C/G showed higher levels compared to the +3010C-C genotype (P< 0.05); those presenting +3027C-C showed higher levels than the +3027A-A genotype (P< 0.05); and those bearing +3035C-C showed higher levels compared to the +3035C-T (P < 0.01) and +3035T-T (P < 0.05) genotypes. The analyses of 3’UTR haplotypes showed that UTR-1 (DelTGCCCGC) was associated with higher expression of sHLA-G, whereas UTR-5 (InsTCCTGAC) and UTR-7 (InsTCATGAC) with lower expression and other UTRs (UTR-2/3/4/6) exhibited intermediate levels. Since the differential expression of HLA-G may be beneficial or harmful depending on the underlying condition, the identification of individuals genetically programmed to differentially express HLA-G may help on defining novel strategies to control the immune response against the underlying disorder.     

DNA sequence variation at a duplicated gene: excess of replacement polymorphism and extensive haplotype structure in the Drosophila melanogaster bicoid region

The bicoid (bcd) gene of Drosophila has played an important role in understanding the system of developmental regulatory genes that controls segmentation in the fruit fly. Several studies in Drosophila and closely related insects suggest that bcd may be the result of a gene duplication in the Dipteran lineage. In addition, the presence of a large, conserved secondary structure in the 3' untranslated region (UTR) makes the bcd gene a good candidate for studying compensatory evolution and the relationship between RNA secondary structure and patterns of standing variation in natural populations. Despite these interesting aspects, a population-level analysis has until now not been performed on bcd. In this study, DNA sequence variation was examined for a 4-kb region of the bcd gene, including a portion of the 5' UTR, the entire coding region, and the 3' UTR, for 25 Drosophila melanogaster isofemale lines from Zimbabwe and one allele from D. simulans. Statistical tests revealed a significant excess of replacement polymorphisms in the D. melanogaster lineage that are clustered in two putative linker regions of the Bicoid protein. This result is consistent with a relaxation of selective constraints in these regions. In addition, we found a distinct haplotype structure and a significantly smaller number of haplotypes than predicted by the standard neutral model. It is unlikely that the haplotype structure is maintained by epistatic selection acting on the secondary structure in the 3' UTR or by the association of the bcd gene with polymorphic inversions. Instead, our two main observations, namely the occurrence of a haplotype structure and the excess of replacement polymorphisms, may indicate that the selective history of this gene is rather complex, involving both the relaxation of purifying selection in some parts of the protein and the action of positive selection in other parts of the gene region.     

Molecular comparison of the negative-acting nitrogen control gene,nmr, inNeurospora crassa and otherNeurospora and fungal species

In Neurospora crassa, the expression of unlinked structural genes which encode nitrogen catabolic enzymes is subject to genetic and metabolic regulation. The negative-acting nmr regulatory gene appears to play a role in nitrogen catabolite repression. Using the N. crassa nmr gene as a probe, homologous sequences were identified in a variety of other filamentous fungi. The polymerase chain reaction was used to isolate the nmr-like gene from the exotic Mauriceville strain of N. crassa and from the two related species, N. intermedia and N. sitophila. Sequence comparisons were carried out with a 1.7-kb DNA segment which includes the entire coding region of nmr plus 5' and 3' noncoding sequences. The size of the nmr coding region was identical in all three Neurospora species. Approximately 30 nucleotide base substitutions were found in the coding region of the nmr gene of each of the sister species when compared to the standard N. crassa sequence. However, most of the base changes occurred in third codon positions and were silent. The NMR proteins of N. sitophila and of N. intermedia display only three and four amino acid substitutions, respectively, from the N. crassa protein. Two regions of high variability, which include deletions and insertions of bases, were found in the 5' and 3' noncoding regions of the gene.     

Regulatory elements of copia retrotransposons, controlling the level of its expression in Drosophila melanogaster testes

Expression of the lacZ reporter gene under the control of five deletion derivatives of the copia regulatory region including the 5' long terminal repeat (LTR) and the 5' untranslated region (UTR) was assayed in the testes of transgenic Drosophila melanogaster males (larvae and imago). The full-length copia regulatory region (LTR + UTR) ensured expression of the reporter gene in testes of both larvae and adult males. Deletion of UTR or 3' end of LTR increased lacZ expression in the testes, whereas deletion of the 5' end of LTR increased it. This indicated that a positive regulator of copia expression is at the 5' end of LTR and that negative regulators are at the 3' end of LTR and in UTR. The effects of the fragments of the copia regulatory region on reporter gene expression in the testes in vivo did not completely coincide with the effects observed earlier in cultured cells. We suggest that this difference is due to different regulation of expression of the fusion constructs integrated into chromatin as compared to their transient expression.     

Identification and in silico analysis of functional SNPs of the BRCA1 gene

Single-nucleotide polymorphisms (SNPs) play a major role in the understanding of the genetic basis of many complex human diseases. Also, the genetics of human phenotype variation could be understood by knowing the functions of these SNPs. It is still a major challenge to identify the functional SNPs in a disease-related gene. In this work, we have analyzed the genetic variation that can alter the expression and the function of the BRCA1 gene using computational methods. Of the total 477 SNPs, 65 were found to be nonsynonymous (ns) SNPs. Among the 14 SNPs in the untranslated region, 4 were found in the 5' and 10 were found in the 3' untranslated region (UTR). It was found that 16.9% of the nsSNPs were damaging, by both the SIFT and the PolyPhen servers. The UTR Resource tool suggested that 2 of 4 SNPs in the 5' UTR and 3 of 10 SNPs in the 3' UTR might change the protein expression levels. We identified major mutations from proline to serine at positions 1776 and 1812 of the native protein of the BRCA1 gene. From a comparison of the stabilizing residues of the native and mutant proteins, we propose that an nsSNP (rs1800751) could be an important candidate for the breast cancer caused by the BRCA1 gene.     

Interplay between hnRNP A1 and a cis-acting element in the 3' UTR of CYP2A5 mRNA is central for high expression of the gene

Our previous evidence suggests that heterogeneous nuclear ribonucleoprotein (hnRNP) A1 plays a part in the regulation of the Cyp2a5 gene by interacting with the 3' untranslated region (UTR) of the CYP2A5 mRNA. However, the exact role of this interaction is not clear. The aim of the present work was to gain further insight into the regulation process of Cyp2a5. For this purpose the 3' UTR of CYP2A5 was fused to the coding region of luciferase mRNA. Luciferase recombinants containing either the full length 3' UTR, or the 3' UTR lacking a previously described 71 nucleotide (nt) region (the hnRNP A1 primary binding site), were transiently expressed in cells expressing or lacking hnRNP A1. The expression of the luciferase recombinants was examined both at mRNA and enzyme activity levels. The results disclosed that the presence of hnRNP A1 was required for the high expression of the recombinant carrying the full length 3' UTR of CYP2A5. Deletion of the hnRNP A1 primary binding site dramatically modified the expression pattern: the mRNA levels and luciferase activities of the deletion mutant were independent from hnRNP A1. These results conclusively demonstrate that the 71 nt region in the 3' UTR of CYP2A5 mRNA can confer hnRNP A1-dependent regulation to a gene. In addition, comparison of RNA levels and luciferase activities suggested that regions flanking the hnRNP A1 binding site could regulate translation of the CYP2A5 mRNA. These results are consistent with a model in which the binding of hnRNP A1 to the 71 nt putative hairpin-loop region in the CYP2A5 mRNA 3' UTR upregulates mRNA levels possibly by protecting the mRNA from degradation.     

Mining the LIPG allelic spectrum reveals the contribution of rare and common regulatory variants to HDL cholesterol

Genome-wide association studies (GWAS) have successfully identified loci associated with quantitative traits, such as blood lipids. Deep resequencing studies are being utilized to catalogue the allelic spectrum at GWAS loci. The goal of these studies is to identify causative variants and missing heritability, including heritability due to low frequency and rare alleles with large phenotypic impact. Whereas rare variant efforts have primarily focused on nonsynonymous coding variants, we hypothesized that noncoding variants in these loci are also functionally important. Using the HDL-C gene LIPG as an example, we explored the effect of regulatory variants identified through resequencing of subjects at HDL-C extremes on gene expression, protein levels, and phenotype. Resequencing a portion of the LIPG promoter and 5' UTR in human subjects with extreme HDL-C, we identified several rare variants in individuals from both extremes. Luciferase reporter assays were used to measure the effect of these rare variants on LIPG expression. Variants conferring opposing effects on gene expression were enriched in opposite extremes of the phenotypic distribution. Minor alleles of a common regulatory haplotype and noncoding GWAS SNPs were associated with reduced plasma levels of the LIPG gene product endothelial lipase (EL), consistent with its role in HDL-C catabolism. Additionally, we found that a common nonfunctional coding variant associated with HDL-C (rs2000813) is in linkage disequilibrium with a 5' UTR variant (rs34474737) that decreases LIPG promoter activity. We attribute the gene regulatory role of rs34474737 to the observed association of the coding variant with plasma EL levels and HDL-C. Taken together, the findings show that both rare and common noncoding regulatory variants are important contributors to the allelic spectrum in complex trait loci.