Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa.

PURPOSE: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. METHODS: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. RESULTS: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. CONCLUSIONS: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Addressing trypsin bias in large scale (phospho)proteome analysis by size exclusion chromatography and secondary digestion of large post-trypsin peptides.

In the vast majority of bottom-up proteomics studies, protein digestion is performed using only mammalian trypsin. Although it is clearly the best enzyme available, the sole use of trypsin rarely leads to complete sequence coverage, even for abundant proteins. It is commonly assumed that this is because many tryptic peptides are either too short or too long to be identified by RPLC-MS/MS. We show through in silico analysis that 20-30% of the total sequence of three proteomes (Schizosaccharomyces pombe, Saccharomyces cerevisiae, and Homo sapiens) is expected to be covered by Large post-Trypsin Peptides (LpTPs) with M(r) above 3000 Da. We then established size exclusion chromatography to fractionate complex yeast tryptic digests into pools of peptides based on size. We found that secondary digestion of LpTPs followed by LC-MS/MS analysis leads to a significant increase in identified proteins and a 32-50% relative increase in average sequence coverage compared to trypsin digestion alone. Application of the developed strategy to analyze the phosphoproteomes of S. pombe and of a human cell line identified a significant fraction of novel phosphosites. Overall our data indicate that specific targeting of LpTPs can complement standard bottom-up workflows to reveal a largely neglected portion of the proteome.

An integrated encyclopedia of DNA elements in the human genome.

The human genome encodes the blueprint of life, but the function of the vast majority of its nearly three billion bases is unknown. The Encyclopedia of DNA Elements (ENCODE) project has systematically mapped regions of transcription, transcription factor association, chromatin structure and histone modification. These data enabled us to assign biochemical functions for 80% of the genome, in particular outside of the well-studied protein-coding regions. Many discovered candidate regulatory elements are physically associated with one another and with expressed genes, providing new insights into the mechanisms of gene regulation. The newly identified elements also show a statistical correspondence to sequence variants linked to human disease, and can thereby guide interpretation of this variation. Overall, the project provides new insights into the organization and regulation of our genes and genome, and is an expansive resource of functional annotations for biomedical research.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution.

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Sequence homologies in the region preceding the transcription initiation site of the liver estrogen-responsive vitellogenin and apo-VLDLII genes.

In the liver of oviparous vertebrates vitellogenin gene expression is controlled by estrogen. The nucleotide sequence of the 5' flanking region of the Xenopus laevis vitellogenin genes A1, A2, B1 and B2 has been determined. These sequences have been compared to each other and to the equivalent region of the chicken vitellogenin II and apo-VLDLII genes which are also expressed in the liver in response to estrogen. The homology between the 5' flanking region of the Xenopus genes B1 and B2 is higher than between the corresponding regions of the other closely related genes A1 and A2. Four short blocks of sequence homology which are present at equivalent positions in the vitellogenin genes of both Xenopus laevis and chicken are characterized. A short sequence with two-fold rotational symmetry (GGTCANNNTGACC) was found at similar positions upstream of the five vitellogenin genes and is also present in two copies close to the 5' end of the chicken apo-VLDLII gene. The possible functional significance of this sequence, common to liver estrogen-responsive genes, is discussed.

Low Y chromosome variation in Saudi-Arabian hamadryas baboons (Papio hamadryas hamadryas).

It is important to characterise the amount of variation on the mammalian Y chromosome in order to assess its potential for use in evolutionary studies. We report very low levels of polymorphism on the Y chromosome of Saudi-Arabian hamadryas baboons, Papio hamadryas hamadryas. We found no segregating sites on the Y, despite sequence analysis of 3 kb noncontiguous intron sequence in 16 males with divergent autosomal microsatellite genotypes, and a further analysis of 1.1 kb intron sequence in 97 males from four populations by SSCP. In addition, we tested seven human-derived Y-linked microsatellites in baboons. Only four of these loci were male-specific and only one was polymorphic in our 97 male sample set. Polymorphism on the Y chromosome of Arabian hamadryas appears to be low compared to other primate species for which data are available (eg humans, chimpanzees and bonobos). Low effective population size (Ne) of paternal genes due to polygyny and female-biased adult sex ratio is a potential reason for low Y chromosome variation in this species. However, low Ne for the Y should be counterbalanced to some extent by the species' atypical pattern of male philopatry and female-biased dispersal. Allelic richness averaged over seven loci was not significantly different between an African and an Arabian population, suggesting that loss of variation during the colonisation of Arabia does not explain low Y variation. Finally, in the absence of nucleotide polymorphism, it is unclear to what extent selection could be responsible for low Y variation in this species.

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data.

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.

Relationship between estrogen receptor alpha location and gene induction reveals the importance of downstream sites and cofactors.

BACKGROUND: To understand cancer-related modifications to transcriptional programs requires detailed knowledge about the activation of signal-transduction pathways and gene expression programs. To investigate the mechanisms of target gene regulation by human estrogen receptor alpha (hERalpha), we combine extensive location and expression datasets with genomic sequence analysis. In particular, we study the influence of patterns of DNA occupancy by hERalpha on expression phenotypes. RESULTS: We find that strong ChIP-chip sites co-localize with strong hERalpha consensus sites and detect nucleotide bias near hERalpha sites. The localization of ChIP-chip sites relative to annotated genes shows that weak sites are enriched near transcription start sites, while stronger sites show no positional bias. Assessing the relationship between binding configurations and expression phenotypes, we find binding sites downstream of the transcription start site (TSS) to be equally good or better predictors of hERalpha-mediated expression as upstream sites. The study of FOX and SP1 cofactor sites near hERalpha ChIP sites shows that induced genes frequently have FOX or SP1 sites. Finally we integrate these multiple datasets to define a high confidence set of primary hERalpha target genes. CONCLUSION: Our results support the model of long-range interactions of hERalpha with the promoter-bound cofactor SP1 residing at the promoter of hERalpha target genes. FOX motifs co-occur with hERalpha motifs along responsive genes. Importantly we show that the spatial arrangement of sites near the start sites and within the full transcript is important in determining response to estrogen signaling.

Genotyping microarray for CSNB-associated genes.

PURPOSE: Congenital stationary night blindness (CSNB) is a clinically and genetically heterogeneous retinal disease. Although electroretinographic (ERG) measurements can discriminate clinical subgroups, the identification of the underlying genetic defects has been complicated for CSNB because of genetic heterogeneity, the uncertainty about the mode of inheritance, and time-consuming and costly mutation scanning and direct sequencing approaches. METHODS: To overcome these challenges and to generate a time- and cost-efficient mutation screening tool, the authors developed a CSNB genotyping microarray with arrayed primer extension (APEX) technology. To cover as many mutations as possible, a comprehensive literature search was performed, and DNA samples from a cohort of patients with CSNB were first sequenced directly in known CSNB genes. Subsequently, oligonucleotides were designed representing 126 sequence variations in RHO, CABP4, CACNA1F, CACNA2D4, GNAT1, GRM6, NYX, PDE6B, and SAG and spotted on the chip. RESULTS: Direct sequencing of genes known to be associated with CSNB in the study cohort revealed 21 mutations (12 novel and 9 previously reported). The resultant microarray containing oligonucleotides, which allow to detect 126 known and novel mutations, was 100% effective in determining the expected sequence changes in all known samples assessed. In addition, investigation of 34 patients with CSNB who were previously not genotyped revealed sequence variants in 18%, of which 15% are thought to be disease-causing mutations. CONCLUSIONS: This relatively inexpensive first-pass genetic testing device for patients with a diagnosis of CSNB will improve molecular diagnostics and genetic counseling of patients and their families and gives the opportunity to analyze whether, for example, more progressive disorders such as cone or cone-rod dystrophies underlie the same gene defects.

GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Detection of three Allexivirus species infecting garlic in Brazil

Garlic viruses often occur in mixed infections under field conditions. In this study, garlic samples collected in three geographical areas of Brazil were tested by Dot-ELISA for the detection of allexiviruses using monoclonal specific antibodies to detect Garlic virus A (GarV-A), Garlic virus B (GarV-B), Garlic virus C (GarV-C) and a polyclonal antiserum able to detect the three virus species mentioned plus Garlic virus D (GarV-D). The detected viruses were biologically isolated by successive passages through Chenopodium quinoa. Reverse Transcriptase Polimerase Chain Reaction (RT-PCR) was performed using primers designed from specific regions of the coat protein genes of Japanese allexiviruses available in the Genetic Bank of National Center of Biotechnology Information (NCBI). By these procedures, individual garlic virus genomes were isolated and sequenced. The nucleotide and amino acid sequence analysis and the one with serological data revealed the presence of three distinct allexiviruses GarV-C, GarV-D and a recently described allexivirus, named Garlic mite-borne filamentous virus (GarMbFV), in Brazil.

The evolution of the thyroid hormone distributor protein transthyretin in the order insectivora, class mammalia.

Thyroid hormones are involved in the regulation of growth and metabolism in all vertebrates. Transthyretin is one of the extracellular proteins with high affinity for thyroid hormones which determine the partitioning of these hormones between extracellular compartments and intracellular lipids. During vertebrate evolution, both the tissue pattern of expression and the structure of the gene for transthyretin underwent characteristic changes. The purpose of this study was to characterize the position of Insectivora in the evolution of transthyretin in eutherians, a subclass of Mammalia. Transthyretin was identified by thyroxine binding and Western analysis in the blood of adult shrews, hedgehogs, and moles. Transthyretin is synthesized in the liver and secreted into the bloodstream, similar to the situation for other adult eutherians, birds, and diprotodont marsupials, but different from that for adult fish, amphibians, reptiles, monotremes, and Australian polyprotodont marsupials. For the characterization of the structure of the gene and the processing of mRNA for transthyretin, cDNA libraries were prepared from RNA from hedgehog and shrew livers, and full-length cDNA clones were isolated and sequenced. Sections of genomic DNA in the regions coding for the splice sites between exons 1 and 2 were synthesized by polymerase chain reaction and sequenced. The location of splicing was deduced from comparison of genomic with cDNA nucleotide sequences. Changes in the nucleotide sequence of the transthyretin gene during evolution are most pronounced in the region coding for the N-terminal region of the protein. Both the derived overall amino sequences and the N-terminal regions of the transthyretins in Insectivora were found to be very similar to those in other eutherians but differed from those found in marsupials, birds, reptiles, amphibians, and fish. Also, the pattern of transthyretin precursor mRNA splicing in Insectivora was more similar to that in other eutherians than to that in marsupials, reptiles, and birds. Thus, in contrast to the marsupials, with a different pattern of transthyretin gene expression in the evolutionarily "older" polyprotodonts compared with the evolutionarily "younger" diprotodonts, no separate lineages of transthyretin evolution could be identified in eutherians. We conclude that transthyretin gene expression in the liver of adult eutherians probably appeared before the branching of the lineages leading to modern eutherian species.

Study of simple sequence repeat markers from coffee expressed sequences associated to leaf miner resistance

The objective of this work was to identify expressed simple sequence repeats (SSR) markers associated to leaf miner resistance in coffee progenies. Identification of SSR markers was accomplished by directed searches on the Brazilian Coffee Expressed Sequence Tags (EST) database. Sequence analysis of 32 selected SSR loci showed that 65% repeats are of tetra-, 21% of tri- and 14% of dinucleotides. Also, expressed SSR are localized frequently in the 5'-UTR of gene transcript. Moreover, most of the genes containing SSR are associated with defense mechanisms. Polymorphisms were analyzed in progenies segregating for resistance to the leaf miner and corresponding to advanced generations of a Coffea arabica x Coffea racemosa hybrid. Frequency of SSR alleles was 2.1 per locus. However, no polymorphism associated with leaf miner resistance was identified. These results suggest that marker-assisted selection in coffee breeding should be performed on the initial cross, in which genetic variability is still significant.

Multicentric carpotarsal osteolysis is caused by mutations clustering in the amino-terminal transcriptional activation domain of MAFB.

Multicentric carpotarsal osteolysis (MCTO) is a rare skeletal dysplasia characterized by aggressive osteolysis, particularly affecting the carpal and tarsal bones, and is frequently associated with progressive renal failure. Using exome capture and next-generation sequencing in five unrelated simplex cases of MCTO, we identified previously unreported missense mutations clustering within a 51 base pair region of the single exon of MAFB, validated by Sanger sequencing. A further six unrelated simplex cases with MCTO were also heterozygous for previously unreported mutations within this same region, as were affected members of two families with autosomal-dominant MCTO. MAFB encodes a transcription factor that negatively regulates RANKL-induced osteoclastogenesis and is essential for normal renal development. Identification of this gene paves the way for development of novel therapeutic approaches for this crippling disease and provides insight into normal bone and kidney development.

Mutations in DNAH5 account for only 15% of a non-preselected cohort of patients with primary ciliary dyskinesia.

BACKGROUND: Primary ciliary dyskinesia (PCD) is characterised by recurrent infections of the upper respiratory airways (nose, bronchi, and frontal sinuses) and randomisation of left-right body asymmetry. To date, PCD is mainly described with autosomal recessive inheritance and mutations have been found in five genes: the dynein arm protein subunits DNAI1, DNAH5 and DNAH11, the kinase TXNDC3, and the X-linked retinitis pigmentosa GTPase regulator RPGR. METHODS: We screened 89 unrelated individuals with PCD for mutations in the coding and splice site regions of the gene DNAH5 by denaturing high performance liquid chromatography (DHPLC) and sequencing. Patients were mainly of European origin and were recruited without any phenotypic preselection. RESULTS: We identified 18 novel (nonsense, splicing, small deletion and missense) and six previously described mutations. Interestingly, these DNAH5 mutations were mainly associated with outer + inner dyneins arm ultrastructural defects (50%). CONCLUSION: Overall, mutations on both alleles of DNAH5 were identified in 15% of our clinically heterogeneous cohort of patients. Although genetic alterations remain to be identified in most patients, DNAH5 is to date the main PCD gene.