Wide-Scale Analysis of Human Functional Transcription Factor Binding Reveals a Strong Bias towards the Transcription Start Site

We introduce a novel method to screen the promoters of a set of genes with shared biological function, against a precompiled library of motifs, and find those motifs which are statistically over-represented in the gene set. The gene sets were obtained from the functional Gene Ontology (GO) classification; for each set and motif we optimized the sequence similarity score threshold, independently for every location window (measured with respect to the TSS), taking into account the location dependent nucleotide heterogeneity along the promoters of the target genes. We performed a high throughput analysis, searching the promoters (from 200bp downstream to 1000bp upstream the TSS), of more than 8000 human and 23,000 mouse genes, for 134 functional Gene Ontology classes and for 412 known DNA motifs. When combined with binding site and location conservation between human and mouse, the method identifies with high probability functional binding sites that regulate groups of biologically related genes. We found many location-sensitive functional binding events and showed that they clustered close to the TSS. Our method and findings were put to several experimental tests. By allowing a "flexible" threshold and combining our functional class and location specific search method with conservation between human and mouse, we are able to identify reliably functional TF binding sites. This is an essential step towards constructing regulatory networks and elucidating the design principles that govern transcriptional regulation of expression. The promoter region proximal to the TSS appears to be of central importance for regulation of transcription in human and mouse, just as it is in bacteria and yeast.Comment: 31 pages, including Supplementary Information and figure

Characterization of Inhibitory Anti-Duffy Binding Protein II Immunity: Approach to Plasmodium vivax Vaccine Development in Thailand

Plasmodium vivax Duffy binding protein region II (DBPII) is an important vaccine candidate for antibody-mediated immunity against vivax malaria. A significant challenge for vaccine development of DBPII is its highly polymorphic nature that alters sensitivity to neutralizing antibody responses. Here, we aim to characterize naturally-acquired neutralizing antibodies against DBPII in individual Thai residents to give insight into P. vivax vaccine development in Thailand. Anti-DBPII IgG significantly increased in acute vivax infections compared to uninfected residents and naive controls. Antibody titers and functional anti-DBPII inhibition varied widely and there was no association between titer and inhibition activity. Most high titer plasmas had only a moderate to no functional inhibitory effect on DBP binding to erythrocytes, indicating the protective immunity against DBPII binding is strain specific. Only 5 of 54 samples were highly inhibitory against DBP erythrocyte-binding function. Previously identified target epitopes of inhibitory anti-DBPPII IgG (H1, H2 and H3) were localized to the dimer interface that forms the DARC binding pocket. Amino acid polymorphisms (monomorphic or dimorphic) in H1 and H3 protective epitopes change sensitivity of immune inhibition by alteration of neutralizing antibody recognition. The present study indicates Thai variant H1.T1 (R308S), H3.T1 (D384G) and H3.T3 (K386N) are the most important variants for a DBPII candidate vaccine needed to protect P. vivax in Thai residents

Early Evolution of Conserved Regulatory Sequences Associated with Development in Vertebrates

Comparisons between diverse vertebrate genomes have uncovered thousands of highly conserved non-coding sequences, an increasing number of which have been shown to function as enhancers during early development. Despite their extreme conservation over 500 million years from humans to cartilaginous fish, these elements appear to be largely absent in invertebrates, and, to date, there has been little understanding of their mode of action or the evolutionary processes that have modelled them. We have now exploited emerging genomic sequence data for the sea lamprey, Petromyzon marinus, to explore the depth of conservation of this type of element in the earliest diverging extant vertebrate lineage, the jawless fish (agnathans). We searched for conserved non-coding elements (CNEs) at 13 human gene loci and identified lamprey elements associated with all but two of these gene regions. Although markedly shorter and less well conserved than within jawed vertebrates, identified lamprey CNEs are able to drive specific patterns of expression in zebrafish embryos, which are almost identical to those driven by the equivalent human elements. These CNEs are therefore a unique and defining characteristic of all vertebrates. Furthermore, alignment of lamprey and other vertebrate CNEs should permit the identification of persistent sequence signatures that are responsible for common patterns of expression and contribute to the elucidation of the regulatory language in CNEs. Identifying the core regulatory code for development, common to all vertebrates, provides a foundation upon which regulatory networks can be constructed and might also illuminate how large conserved regulatory sequence blocks evolve and become fixed in genomic DNA

Micronutrient status in lactating mothers before and after introduction of fortified flour: cross-sectional surveys in Maela refugee camp

Background Deficiency of micronutrients is common in refugee populations. Objectives Identify deficiencies and whether provided supplements and wheat flour fortified with 10 micronutrients impacts upon status among breast-feeding women from Maela refugee camp. Methods Two sequential cross-sectional studies were conducted in different groups of lactating mothers at 12 weeks postpartum. The first survey was before and the second 4-5 months after micronutrient fortified flour (MFF) had been provided to the camp (in addition to the regular food basket). Iron status and micronutrients were measured in serum, whole blood, and in breast milk samples. Results Iron and zinc deficiency and anemia were highly prevalent while low serum retinol and thiamine deficiency were rarely detected. Iron and zinc deficiency were associated with anemia, and their proportions were significantly lower after the introduction of MFF (21 vs. 35% with soluble transferrin receptor (sTfR)>8.5 mg/L, P = 0.042, and 50 vs. 73% with serum zinc<0.66 mg/L, P = 0.001). Serum sTfR, whole-blood thiamine diphosphate (TDP) and serum β-carotene were significant predictors (P<0.001) of milk iron, thiamine and β-carotene, respectively. Lower prevalence of iron deficiency in the MFF group was associated with significantly higher iron and thiamine in breast milk. Conclusions High whole-blood TDP and breast milk thiamine reflected good compliance to provided thiamine; high prevalence of iron deficiency suggested insufficient dietary iron and low acceptance to ferrous sulfate supplements. MFF as an additional food ration in Maela refugee camp seemed to have an effect in reducing both iron and zinc deficiency postpartum. © Springer-Verlag 2012

Functional annotations of diabetes nephropathy susceptibility loci through analysis of genome-wide renal gene expression in rat models of diabetes mellitus

<p>Abstract</p> <p>Background</p> <p>Hyperglycaemia in diabetes mellitus (DM) alters gene expression regulation in various organs and contributes to long term vascular and renal complications. We aimed to generate novel renal genome-wide gene transcription data in rat models of diabetes in order to test the responsiveness to hyperglycaemia and renal structural changes of positional candidate genes at selected diabetic nephropathy (DN) susceptibility loci.</p> <p>Methods</p> <p>Both Affymetrix and Illumina technologies were used to identify significant quantitative changes in the abundance of over 15,000 transcripts in kidney of models of spontaneous (genetically determined) mild hyperglycaemia and insulin resistance (Goto-Kakizaki-GK) and experimentally induced severe hyperglycaemia (Wistar-Kyoto-WKY rats injected with streptozotocin [STZ]).</p> <p>Results</p> <p>Different patterns of transcription regulation in the two rat models of diabetes likely underlie the roles of genetic variants and hyperglycaemia severity. The impact of prolonged hyperglycaemia on gene expression changes was more profound in STZ-WKY rats than in GK rats and involved largely different sets of genes. These included genes already tested in genetic studies of DN and a large number of protein coding sequences of unknown function which can be considered as functional and, when they map to DN loci, positional candidates for DN. Further expression analysis of rat orthologs of human DN positional candidate genes provided functional annotations of known and novel genes that are responsive to hyperglycaemia and may contribute to renal functional and/or structural alterations.</p> <p>Conclusion</p> <p>Combining transcriptomics in animal models and comparative genomics provides important information to improve functional annotations of disease susceptibility loci in humans and experimental support for testing candidate genes in human genetics.</p

Gene expression of NMDA receptor subunits in the cerebellum of elderly patients with schizophrenia

To determine if NMDA receptor alterations are present in the cerebellum in schizophrenia, we measured NMDA receptor binding and gene expression of the NMDA receptor subunits in a post-mortem study of elderly patients with schizophrenia and non-affected subjects. Furthermore, we assessed influence of genetic variation in the candidate gene neuregulin-1 (NRG1) on the expression of the NMDA receptor in an exploratory study. Post-mortem samples from the cerebellar cortex of ten schizophrenic patients were compared with nine normal subjects. We investigated NMDA receptor binding by receptor autoradiography and gene expression of the NMDA receptor subunits NR1, NR2A, NR2B, NR2C and NR2D by in situ hybridization. For the genetic study, we genotyped the NRG1 polymorphism rs35753505 (SNP8NRG221533). Additionally, we treated rats with the antipsychotics haloperidol or clozapine and assessed cerebellar NMDA receptor binding and gene expression of subunits to examine the effects of antipsychotic treatment. Gene expression of the NR2D subunit was increased in the right cerebellum of schizophrenic patients compared to controls. Individuals carrying at least one C allele of rs35753505 (SNP8NRG221533) showed decreased expression of the NR2C subunit in the right cerebellum, compared to individuals homozygous for the T allele. Correlation with medication parameters and the animal model revealed no treatment effects. In conclusion, increased NR2D expression results in a hyperexcitable NMDA receptor suggesting an adaptive effect due to receptor hypofunction. The decreased NR2C expression in NRG1 risk variant may cause a deficit in NMDA receptor function. This supports the hypothesis of an abnormal glutamatergic neurotransmission in the right cerebellum in the pathophysiology of schizophrenia

Genetic variation in the odorant receptors family 13 and the mhc loci influence mate selection in a multiple sclerosis dataset

<p>Abstract</p> <p>Background</p> <p>When selecting mates, many vertebrate species seek partners with major histocompatibility complex (MHC) genes different from their own, presumably in response to selective pressure against inbreeding and towards MHC diversity. Attempts at replication of these genetic results in human studies, however, have reached conflicting conclusions.</p> <p>Results</p> <p>Using a multi-analytical strategy, we report validated genome-wide relationships between genetic identity and human mate choice in 930 couples of European ancestry. We found significant similarity between spouses in the MHC at class I region in chromosome 6p21, and at the odorant receptor family 13 locus in chromosome 9. Conversely, there was significant dissimilarity in the MHC class II region, near the <it>HLA-DQA1 </it>and -<it>DQB1 </it>genes. We also found that genomic regions with significant similarity between spouses show excessive homozygosity in the general population (assessed in the HapMap CEU dataset). Conversely, loci that were significantly dissimilar among spouses were more likely to show excessive heterozygosity in the general population.</p> <p>Conclusions</p> <p>This study highlights complex patterns of genomic identity among partners in unrelated couples, consistent with a multi-faceted role for genetic factors in mate choice behavior in human populations.</p

Measuring the Evolutionary Rewiring of Biological Networks

We have accumulated a large amount of biological network data and expect even more to come. Soon, we anticipate being able to compare many different biological networks as we commonly do for molecular sequences. It has long been believed that many of these networks change, or “rewire”, at different rates. It is therefore important to develop a framework to quantify the differences between networks in a unified fashion. We developed such a formalism based on analogy to simple models of sequence evolution, and used it to conduct a systematic study of network rewiring on all the currently available biological networks. We found that, similar to sequences, biological networks show a decreased rate of change at large time divergences, because of saturation in potential substitutions. However, different types of biological networks consistently rewire at different rates. Using comparative genomics and proteomics data, we found a consistent ordering of the rewiring rates: transcription regulatory, phosphorylation regulatory, genetic interaction, miRNA regulatory, protein interaction, and metabolic pathway network, from fast to slow. This ordering was found in all comparisons we did of matched networks between organisms. To gain further intuition on network rewiring, we compared our observed rewirings with those obtained from simulation. We also investigated how readily our formalism could be mapped to other network contexts; in particular, we showed how it could be applied to analyze changes in a range of “commonplace” networks such as family trees, co-authorships and linux-kernel function dependencies

Achlya mitochondrial DNA: gene localization and analysis of inverted repeats

Mitochondrial DNA from four strains of the oomycete Achlya has been compared and nine gene loci mapped, including that of the ribosomal protein gene, var1 . Examination of the restriction enzyme site maps showed the presence of four insertions relative to a map common to all four strains. All the insertions were found in close proximity to genic regions. The four strains also cotained the inverted repeat first observed in A. ambisexualis (Hudspeth et al. 1983), allowing an examination by analysis of retained restriction sites of the evolutionary stability of repeated DNA sequences relative to single copy sequences. Although the inverted repeat is significantly more stable than single copy sequences, more detailed analysis indicated that this stability is limited to the portion encoding the ribosomal RNA genes. Thus, the apparent evolutionary stability of the repeat does not appear to derive from the inverted repeat structure per se.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/47563/1/438_2004_Article_BF00330510.pd