Recommendations for locus-specific databases and their curation.

Expert curation and complete collection of mutations in genes that affect human health is essential for proper genetic healthcare and research. Expert curation is given by the curators of gene-specific mutation databases or locus-specific databases (LSDBs). While there are over 700 such databases, they vary in their content, completeness, time available for curation, and the expertise of the curator. Curation and LSDBs have been discussed, written about, and protocols have been provided for over 10 years, but there have been no formal recommendations for the ideal form of these entities. This work initiates a discussion on this topic to assist future efforts in human genetics. Further discussion is welcome.

Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

Examining ERBB2 as a candidate gene for susceptibility to leprosy (Hansen’s disease) in Brazil

Leprosy remains prevalent in Brazil. ErbB2 is a receptor for leprosy bacilli entering Schwann cells, which mediates Mycobacterium leprae-induced demyelination and the ERBB2 gene lies within a leprosy susceptibility locus on chromosome 17q11-q21. To determine whether polymorphisms at the ERBB2 locus contribute to this linkage peak, three haplotype tagging single nucleotide polymorphisms (tag-SNPs) (rs2517956, rs2952156, rs1058808) were genotyped in 72 families (208 cases; 372 individuals) from the state of Pará (PA). All three tag-SNPs were associated with leprosy per se [best SNP rs2517959 odds ratio (OR) = 2.22; 95% confidence interval (CI) 1.37-3.59; p = 0.001]. Lepromatous (LL) (OR = 3.25; 95% CI 1.37-7.70; p = 0.007) and tuberculoid (TT) (OR = 1.79; 95% CI 1.04-3.05; p = 0.034) leprosy both contributed to the association, which is consistent with the previous linkage to chromosome 17q11-q21 in the population from PA and supports the functional role of ErbB2 in disease pathogenesis. To attempt to replicate these findings, six SNPs (rs2517955, rs2517956, rs1810132, rs2952156, rs1801200, rs1058808) were genotyped in a population-based sample of 570 leprosy cases and 370 controls from the state of Rio Grande do Norte (RN) and the results were analysed using logistic regression analysis. However, none of the associations were replicated in the RN sample, whether analysed for leprosy per se, LL leprosy, TT leprosy, erythema nodosum leprosum or reversal reaction conditions. The role of polymorphisms at ERBB2 in controlling susceptibility to leprosy in Brazil therefore remains unclear.

Injection of partially purified estrogen receptor protein from Xenopus liver nuclei into oocytes activates the silent vitellogenin locus.

Injection of extracts from Xenopus liver nuclei that are enriched 2000 times in estradiol receptor into Xenopus oocytes induces transcription of the silent vitellogenin locus, which is activated in liver by estradiol, but not of the albumin locus, which is active in liver but suppressed by high levels of estradiol. Transcription initiates within the 5'-end region of the gene we have studied and probably continues into the 3' third. The activation seems to be very efficient, but most of the primary transcripts are probably rapidly and inaccurately processed. New proteins are also made and secreted by the oocytes.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.

Effects of epigenetic regulatory elements on telomeric gene expression

Transgene expression in eukaryotic cells strongly depends on the locus of integration in the host genome and often results in limited transcription level because of unfavorable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which are believed to act on the chromatin structure and may protect transgenes from this so-called position effect. Despite being extensively used to increase transgene expression, the mechanism of action of many of these elements remains largely unknown. Here we evaluated different epigenetic regulatory DNA elements for their ability to protect transgene transcription at telomeres, a defined chromatin environment associated to low or inconsistent expression caused by the Telomere Position Effect (TPE). For the assessment of the effects of epigenetic regulators at telomeres, a novel dual reporter system had to be designed. Telomeric integration of the newly-developed dual reporter system carrying different epigenetic regulators showed that MARs (Matrix Attachment Regions), a UCOE (Ubiquitous Chromatin-Opening Element) or the chicken cHS4 insulator have strong barrier activity which prevented TPE from spreading toward the centromere, resulting in stable and in some cases increased expression of a telomeric-distal reporter gene. In addition, MARs and STAR element 40 resulted in an increase of cells expressing the telomeric-proximal reporter gene, suggesting also an anti-silencing effect. Chromatin immunoprecipitation assays revealed that at telomeres MARs actively promote the deposition of euchromatic histone marks, especially acetylation of both histone H3 and H4, which might be involved in MARs' barrier and transcriptional activator activities. Differently, the chromatin in proximity of the UCOE element was depleted of several repressive chromatin marks, such as trimethylated lysine 9 and lysine 27 on histone H3 and trimethylated lysine 20 of histone H4, possibly favoring the preservation of an open chromatin structure at the integration site. We conclude that epigenetic regulatory elements that may be used to enhance and sustain transgene expression have all a specific epigenetic signature which might be at the basis of their mechanism of action, and that a combination of different classes of epigenetic regulators might be advantageous when high levels of protein expression are required. - L'expression des transgènes dans les cellules eucaryotes est fortement influencée par leur site d'intégration dans le génome. En effet, une structure chromatinienne défavorable au niveau du site d'intégration peut fortement limiter le degré d'expression d'un transgène. Il existe toutefois des séquences d'ADN qui, en agissant sur la structure de la chromatine, permettent de limiter cet effet de position et, par conséquent, de promouvoir l'expression soutenue d'un transgène. Ces éléments génomiques, connus comme régulateurs épigénétiques, sont largement utilisés dans plusieurs domaines où une expression élevée et soutenue est requise, malgré un mode de fonctionnement parfois méconnu. Dans cette étude, j'ai évalué la capacité de différents régulateurs épigénétiques à protéger la transcription de transgènes au niveau des télomères, régions chromatiniennes bien définies qui ont été associées à un fort effet de silençage, connu comme «effet de position télomérique». Pour cela, un nouveau système à deux gènes rapporteurs a été développé. Lorsque des MARs (Matrix Attachment Regions, séquences d'ADN pouvant s'associer à la matrice nucléaire), un UCOE (Ubiquitous Chromatin-Opening Element, élément pouvant ouvrir la chromatine) ou l'isolateur génétique cHS4 (dérivé du locus de la β-globine de poulet) sont placés entre les deux gènes rapporteurs, une forte activité barrière bloquant la propagation de la chromatine répressive télomérique est observée, résultant en un plus grand nombre de cellules exprimant le gène télomérique-distal. D'autre part, une augmentation du nombre de cellules exprimant le gène télomérique-proximal, observée en présence des éléments MAR et STAR 40 (Stabilizing Anti-Repressor element 40, un élément pouvant prévenir la répression génique), suggère aussi un faible effet anti-silençage pour ces éléments. Des expériences d'immunoprécipitation de la chromatine démontrent qu'au télomère, les MARs favorisent l'assemblage de marqueurs de la chromatine active, surtout l'acétylation des histones H3 et H4, qui pourraient être à la base de l'activité barrière et de celle d'activateur transcriptionel. Différemment, la chromatine à proximité de l'élément UCOE est particulièrement pauvre en marqueurs de la chromatine silencieuse, comme la trimethylation des lysines 9 et 27 de l'histone H3, ainsi que la trimethylation de la lysine 20 de l'histone H4. Cela suggère que UCOE pourrait préserver une structure chromatinienne ouverte au site d'intégration, favorisant l'expression des gènes à sa proximité. En conclusion, les régulateurs épigénétiques analysés lors de cette étude ont tous montré une signature épigénétique spécifique qui pourrait être à la base de leurs mécanismes de fonctionnement, suggérant aussi qu'une utilisation d'éléments épigénétiques de classe différente dans un même vecteur d'expression pourrait être avantageuse lorsque de hauts et soutenus niveaux d'expression sont nécessaires.

Recessive multiple epiphyseal dysplasia (rMED) with homozygosity for C653S mutation in the DTDST gene - phenotype, molecular diagnosis and surgical treatment of habitual dislocation of multilayered patella: case report.

BACKGROUND: Multiple epiphyseal dysplasia (MED) is one of the more common generalised skeletal dysplasias. Due to its clinical heterogeneity diagnosis may be difficult. Mutations of at least six separate genes can cause MED. Joint deformities, joint pain and gait disorders are common symptoms. CASE PRESENTATION: We report on a 27-year-old male patient suffering from clinical symptoms of autosomal recessive MED with habitual dislocation of a multilayered patella on both sides, on the surgical treatment and on short-term clinical outcome. Clinical findings were: bilateral hip and knee pain, instability of femorotibial and patellofemoral joints with habitual patella dislocation on both sides, contractures of hip, elbow and second metacarpophalangeal joints. Main radiographic findings were: bilateral dislocated multilayered patella, dysplastic medial tibial plateaus, deformity of both femoral heads and osteoarthritis of the hip joints, and deformity of both radial heads. In the molecular genetic analysis, the DTDST mutation g.1984T > A (p.C653S) was found at the homozygote state. Carrier status was confirmed in the DNA of the patient's parents. The mutation could be considered to be the reason for the patient's disease. Surgical treatment of habitual patella dislocation with medialisation of the tibial tuberosity led to an excellent clinical outcome. CONCLUSIONS: The knowledge of different phenotypes of skeletal dysplasias helps to select genes for genetic analysis. Compared to other DTDST mutations, this is a rather mild phenotype. Molecular diagnosis is important for genetic counselling and for an accurate prognosis. Even in case of a multilayered patella in MED, habitual patella dislocation could be managed successfully by medialisation of the tibial tuberosity.

A large deletion in the adRP gene PRPF31: evidence that haploinsufficiency is the cause of disease.

PURPOSE: To report a large deletion that encompasses more than 90% of PRPF31 gene and two other neighboring genes in their entirety in an adRP pedigree that appears to show only the typical clinical features of retinitis pigmentosa. METHODS: To identify PRPF31 mutation in a dominant RP family (ADRP2) previously linked to the RP11 locus, the 14 exons of PRPF31 were screened for mutations by direct sequencing. To investigate the possibility of a large deletion, microsatellite markers near PRPF31 gene were analyzed by non-denaturing PAGE. RESULTS: Initial screening of PRPF31 gene in the ADRP2 family did not reveal an obvious mutation. A large deletion was however suspected due to lack of heterozygosity for nearly all PRPF31 intragenic single nucleotide polymorphysm (SNPs). In order to estimate the size of the deletion, SNPs and microsatellite markers spanning and flanking PRPF31 were analyzed in the entire ADRP2 family. Haplotype analysis with the above markers suggested a deletion of approximately 30 kb that included the putative promoter region of a novel gene OSCAR, the entire genomic content of genes NDUFA3, TFPT and more than 90% of PRPF31 gene. Sequence analysis of the region flanking the potential deletion showed a high presence of Alu elements implicating Alu mediated recombination as the mechanism responsible for this event. CONCLUSIONS: This mutation provides evidence that haploinsufficiency rather than aberrant function of mutated proteins is the cause of disease in these adRP patients with mutations in PRPF31 gene.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Infantile spasms is associated with deletion of the MAGI2 gene on chromosome 7q11.23-q21.11.

Infantile spasms (IS) is the most severe and common form of epilepsy occurring in the first year of life. At least half of IS cases are idiopathic in origin, with others presumed to arise because of brain insult or malformation. Here, we identify a locus for IS by high-resolution mapping of 7q11.23-q21.1 interstitial deletions in patients. The breakpoints delineate a 500 kb interval within the MAGI2 gene (1.4 Mb in size) that is hemizygously disrupted in 15 of 16 participants with IS or childhood epilepsy, but remains intact in 11 of 12 participants with no seizure history. MAGI2 encodes the synaptic scaffolding protein membrane-associated guanylate kinase inverted-2 that interacts with Stargazin, a protein also associated with epilepsy in the stargazer mouse.

Regulation of the expression of Fit insect toxin locus genes in the root-associated biocontrol pseudomonad CHA0

The root-colonizing Pseudomonas fluorescens strain CHA0 is a biocontrol agent of soil-borne plant diseases caused by fungal and oomycete pathogens. Remarkably, this plant-beneficial pseudomonad is also endowed with potent insecticidal activity that depends on the production of a large protein toxin termed Fit (for P. fluorescens insecticidal toxin). In our present work, the genomic locus encoding the P. fluorescens insect toxin is subjected to a detailed molecular analysis. The Fit toxin gene fitD is flanked upstream by the fitABC genes and downstream by the fitE gene that encode the ABC transporter, membrane fusion, and outer membrane efflux components of a type I protein secretion system predicted to function in toxin export. The fitF, fitG, and fitH genes located downstream of fitE code for regulatory proteins having domain structures typical of signal transduction histidine kinases, LysR-type transcriptional regulators, and response regulators, respectively. The role of these insect toxin locus-associated control elements is being investigated with mutants defective for the regulatory genes and with GFP-based reporter fusions to putative promoter regions upstream of the transporter genes fitA and fitE, the toxin gene fitD, and the regulatory genes fitF and fitH. Our preliminary findings suggest that the three regulators interact with known global regulators of biocontrol factor expression to control Fit toxin expression and secretion.

Sense and antisense transcripts in the histone H1 (HIS-1) locus of Leishmania major.

Histone H1 in the parasitic protozoan Leishmania is a developmentally regulated protein encoded by two genes, HIS-1.1 and HIS-1.2. These genes are separated by approximately 20 kb of sequence and are located on the same DNA strand of chromosome 27. When Northern blots of parasite RNA were probed with HIS-1 strand-specific riboprobes, we detected sense and antisense transcripts that were polyadenylated and developmentally regulated. When the HIS-1.2 coding region was replaced with the coding region of the neomycin phosphotransferase gene, antisense transcription of this gene was unaffected, indicating that the regulatory elements controlling antisense transcription were located outside of the HIS-1.2 gene, and that transcription in Leishmania can occur from both DNA strands even in the presence of transcription of a selectable marker in the complementary strand. A search for other antisense transcripts within the HIS-1 locus identified an additional transcript (SC-1) within the intervening HIS-1 sequence, downstream of adenine and thymine-rich sequences. These results show that gene expression in Leishmania is not only regulated polycistronically from the sense strand of genomic DNA, but that the complementary strand of DNA also contains sequences that could drive expression of open reading frames from the antisense strand of DNA. These findings suggest that the parasite has evolved in such a way as to maximise the transcription of its genome, a mechanism that might be important for it to maintain virulence.

Genomewide association study using a high-density single nucleotide polymorphism array and case-control design identifies a novel essential hypertension susceptibility locus in the promoter region of endothelial NO synthase.

Essential hypertension is a multifactorial disorder and is the main risk factor for renal and cardiovascular complications. The research on the genetics of hypertension has been frustrated by the small predictive value of the discovered genetic variants. The HYPERGENES Project investigated associations between genetic variants and essential hypertension pursuing a 2-stage study by recruiting cases and controls from extensively characterized cohorts recruited over many years in different European regions. The discovery phase consisted of 1865 cases and 1750 controls genotyped with 1M Illumina array. Best hits were followed up in a validation panel of 1385 cases and 1246 controls that were genotyped with a custom array of 14 055 markers. We identified a new hypertension susceptibility locus (rs3918226) in the promoter region of the endothelial NO synthase gene (odds ratio: 1.54 [95% CI: 1.37-1.73]; combined P=2.58 · 10(-13)). A meta-analysis, using other in silico/de novo genotyping data for a total of 21 714 subjects, resulted in an overall odds ratio of 1.34 (95% CI: 1.25-1.44; P=1.032 · 10(-14)). The quantitative analysis on a population-based sample revealed an effect size of 1.91 (95% CI: 0.16-3.66) for systolic and 1.40 (95% CI: 0.25-2.55) for diastolic blood pressure. We identified in silico a potential binding site for ETS transcription factors directly next to rs3918226, suggesting a potential modulation of endothelial NO synthase expression. Biological evidence links endothelial NO synthase with hypertension, because it is a critical mediator of cardiovascular homeostasis and blood pressure control via vascular tone regulation. This finding supports the hypothesis that there may be a causal genetic variation at this locus.

Autosomal Recessive Posterior Microphthalmos/Nanophthalmos is Caused by Loss-of-function Mutations in LOC646960, a Novel Gene Encoding a Trypsin-like Serine Protease

Purpose: Posterior microphthalmos (MCOP)/nanophthalmos (NNO) is a developmental anomaly characterized by extreme hyperopia due to short axial length. The population of the Faroe Islands shows a high prevalence of an autosomal recessive form (arMCOP). The gene mutated in arMCOP is not yet known.Methods: Genetic mapping by linkage analysis using microsatellite and single nucleotide polymorphisms, mutation analysis by PCR and sequencing, molecular modellingResults: Having refined the position of the disease locus (MCOP6) in an interval of 250 kb in chromosome 2q37.1 in Faroese families, we detected 3 mutations in a novel gene, LOC646960: Patients of 10 different Faroese families were either homozygous (n=22) for c.926G>C (p.Trp309Ser) or compound heterozygous (n=6) for c.926G>C and c.526C>G (p.Arg176Gly), whereas a homozygous 1 bp duplication (c.1066dupC) was identified in patients with arNNO from a Tunisian family. In two unrelated patients with MCOP, no LOC646960 mutation was found. LOC646960 is expressed in the human adult retina and RPE. The expression of the mouse homologue in the eye can be first detected at E17 and is highest in adults. The predicted protein is a 603 amino acid long secreted trypsin-like serine peptidase. c.1066dupC should result in a functional null allele. Molecular modelling of the p.Trp309Ser mutant suggests that both affinity and reactivity of the enzyme towards in vivo substrates are substantially reduced.Conclusions: Postnatal growth of the eye is important for proper development of the refractive components (emmetropization), and is mainly due to elongation of the posterior segment from 10-11 mm at birth to 15-16 mm at the age of 13 years. Optical defocus leads to changes in axial length by moving the retina towards the image plane. arMCOP may theoretically be explained, in line with the expression pattern of LOC646960, by a postnatal growth retardation of the posterior segment.