Search

—Biodiversity-based approaches are at the forefront of creating new varieties in soybean production. In this paper, we analyzed the natural variability in the population of cultivated soybean varieties used in breeding in the Russian Federation using genome-wide association studies. Genome regions that control the protein and oil contents in seeds, as well as productivity, have been identified, many of which have not been described previously. These results open the way to the design of new markers for marker selection of this crop.

Transposons are genomic elements that can relocate within a host genome using a ‘cut’- or ‘copy-and-paste’ mechanism. They make up a significant part of many genomes, serve as a driving force for genome evolution, and are linked with Mendelian diseases and cancers. Interactions between two specific retrotransposon types, autonomous (e.g., LINE1/L1) and nonautonomous (e.g., Alu), may lead to fluctuations in the number of these transposons in the genome over multiple cell generations. We developed and examined a simple model of retrotransposon dynamics under conditions where transposon replication machinery competed for cellular resources: namely, free ribosomes and available energy (i.e., ATP molecules). Such competition is likely to occur in stress conditions that a malfunctioning cell may experience as a result of a malignant transformation. The modeling revealed that the number of actively replicating LINE1 and Alu elements in a cell decreases with the increasing competition for resources, however, stochastic effects interfere with this simple trend. We stochastically simulated the transposon dynamics in a cell population and showed that the population splits into pools with drastically different transposon behaviors. The early extinction of active Alu elements resulted in a larger number of LINE1 copies occurring in the first pool, as there was no competition between the two types of transposons in this pool. In the other pool, the competition process remained and the number of L1 copies was kept small. As the level of available resources reached a critical value, both types of dynamics demonstrated an increase in noise levels, and both the period and the amplitude of predator–prey oscillations rose in one of the cell pools. We hypothesized that the presented dynamical effects associated with the impact of the competition for cellular resources inflicted on the dynamics of retrotransposable elements could be used as a characteristic feature to assess a cell state, or to control the transposon activity.

Modern flax cultivars are susceptible to many diseases, arguably, the most economically damaging of these is the Fusarium wilt fungal disease. Over the past decades international flax breeding initiatives resulted in the development of resistant cultivars. However, much remains to be learned about the mechanisms of resistance to Fusarium infection in flax. As a first step to uncover the genetic factors associated with resistance to Fusarium wilt disease, we performed a genome-wide association study (GWAS) using 297 accessions from the collection of the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. These genotypes were infected with a highly pathogenic Fusarium oxysporum f.sp. lini MI39 strain, the wilt symptoms were documented in the course of three successive years. Six different single-locus models implemented in GAPIT3 R package were applied to a selected subset of 72,526 SNPs. A total of 15 QTNs (Quantitative Trait Nucleotides) were detected during at least two years of observation, while eight QTNs were found during all three years of the experiment. Of these, ten QTNs occupied a region of 640 Kb at the start of chromosome 1, while the remaining QTNs mapped to chromosomes 8, 11 and 13. All stable QTNs demonstrate a statistically significant allelic effect across 3 years of the experiment. Importantly, several QTNs spanned regions that harbored genes involved in the pathogen recognition and plant immunity response, including the KIP1-like protein (Lus10025717) and NBS-LRR protein (Lus10025852). Our results provide novel insights into the genetic architecture of flax resistance to Fusarium wilt and pinpoint potential candidate genes for further in-depth studies.

A collection of flax accessions from Russian Federal Research Center for Bast Fiber Crops was characterised to evaluate its phenotypic diversity. 406 samples representing different morphotypes were selected for thorough quantitative assessment of various agronomic traits. We measured height, length of technical part of the stem, technical part weight, inflorescence length, number of bolls and seeds per plant, 1000 seed weight, the diameter of the stem, the number of internodes and finally, distance between internodes. The fiber quality was estimated by calculating stem slenderness, stem taperingness and elementary fiber length. The dataset was produced in a framework of a project focused on characterization of diversity of flax genotypes and phenotypes, as well as on identification of genomic regions associated with various traits, it is hosted on Figshare.

Fusarium wilt of flax is an aggressive disease caused by the soil-borne fungal pathogen Fusarium oxysporum f. sp. lini. It is a challenging pathogen presenting a constant threat to flax production industry worldwide. Previously, we reported chromosome-level assemblies of 5 highly pathogenic F. oxysporum f. sp. lini strains. We sought to characterize the genomic architecture of the fungus and outline evolutionary mechanisms shaping the pathogen genome. Here, we reveal the complex multi-compartmentalized genome organization and uncover its diverse evolutionary dynamics, which boosts genetic diversity and facilitates host adaptation. In addition, our results suggest that host of functions implicated in the life cycle of mobile genetic elements are main contributors to dissimilarity between proteomes of different Fusaria. Finally, our experiments demonstrate that mobile genetics elements are expressed in planta upon infection, alluding to their role in pathogenicity. On the whole, these results pave the way for further in-depth studies of evolutionary forces shaping the host–pathogen interaction.

Background: Alcohol dependence (AD) is a cluster of clinical phenomena that develop after repeated alcohol use. At least 50% of factors predisposing to AD are genetic and women affected with this disorder present with more psychiatric comorbidities, probably indicating different genetic factors involved. We aimed to run a genome-wide association study (GWAS) followed by a bioinformatic functional annotation of associated genomic regions in a cohort of patients with AD and eight additional clinical measures. Methods: Chi-square tests and the signal localization approach were used for the genome-wide screening of associations with AD. Likelihood ratio tests were used to evaluate associations with the additional clinical measures. Regulatory sequences in the associated linkage disequilibrium blocks were retrieved from bioinformatic databases. The resulting list of regulated genes was analyzed for possible interactions. Results: A genome-wide significant association of rs220677 with AD (p-value = 1.33×10-8) was discovered in female patients. Associations with six other single nucleotide polymorphisms listed in previous GWAS of psychiatric and addiction traits were differently replicated in female and male patients. The bioinformatic analysis showed that regulatory elements in the seven associated linkage disequilibrium blocks define the expression of 78 protein-coding genes. Nearly 68% of these and of 120 other previously published coding genes associated with alcohol phenotypes directly interact in a single network (i.e. without additional interactors added by the software). Conclusions: Results of this study suggest that molecular pathways implicated in the pathogenesis of AD are the same in both sexes, but different genes play more prominent roles in men and women. The discovered gene network may contain clues to the understanding of the pathogenesis and to the discovery of new treatment options of AD.

The review discusses, in a format of a timeline, the studies of different types of genetic variants, present in Homo sapiens, but absent in all other primate, mammalian, or vertebrate species, tested so far. The main characteristic of these variants is that they are found in regions of high evolutionary conservation. These sequence variations include single nucleotide substitutions (called human accelerated regions), deletions, and segmental duplications. The rationale for finding such variations in the human genome is that they could be responsible for traits, specific to our species, of which the human brain is the most remarkable. As became obvious, the vast majority of human-specific single nucleotide substitutions are found in noncoding, likely regulatory regions. A number of genes, associated with these human-specific alleles, often through novel enhancer activity, were in fact shown to be implicated in human-specific development of certain brain areas, including the prefrontal cortex. Human-specific deletions may remove regulatory sequences, such as enhancers. Segmental duplications, because of their large size, create new coding sequences, like new functional paralogs. Further functional study of these variants will shed light on evolution of our species, as well as on the etiology of neurodevelopmental disorders.

Arginine methylation widens the mechanism of control by E2F1 from a transcription factor to a regulator of alternative RNA splicing., E2F is a family of master transcription regulators involved in mediating diverse cell fates. Here, we show that residue-specific arginine methylation (meR) by PRMT5 enables E2F1 to regulate many genes at the level of alternative RNA splicing, rather than through its classical transcription-based mechanism. The p100/TSN tudor domain protein reads the meR mark on chromatin-bound E2F1, allowing snRNA components of the splicing machinery to assemble with E2F1. A large set of RNAs including spliced variants associate with E2F1 by virtue of the methyl mark. By focusing on the deSUMOylase SENP7 gene, which we identified as an E2F target gene, we establish that alternative splicing is functionally important for E2F1 activity. Our results reveal an unexpected consequence of arginine methylation, where reader-writer interplay widens the mechanism of control by E2F1, from transcription factor to regulator of alternative RNA splicing, thereby extending the genomic landscape under E2F1 control.

In the present study we conducted a genome-wide association study (GWAS) in a cohort of 505 patients with paranoid schizophrenia (SCZ), of which 95 had tardive dyskinesia (TD), and 503 healthy controls. Using data generated by the PsychENCODE Consortium (PEC) and other bioinformatic databases, we revealed a gene network, implicated in neurodevelopment and brain function, associated with both these disorders. Almost all these genes are in gene or isoform co-expression PEC network modules important for the functioning of the brain; the activity of these networks is also altered in SCZ, bipolar disorder and autism spectrum disorders. The associated PEC network modules are enriched for gene ontology terms relevant to the brain development and function (CNS development, neuron development, axon ensheathment, synapse, synaptic vesicle cycle, and signaling receptor activity) and to the immune system (inflammatory response). Results of the present study suggest that orofacial and limbtruncal types of TD seem to share the molecular network with SCZ. Paranoid SCZ and abnormal involuntary movements that indicate the orofacial type of TD are associated with the same genomic loci on chromosomes 3p22.2, 8q21.13, and 13q14.2. The limbtruncal type of TD is associated with a locus on chromosome 3p13 where the best functional candidate is FOXP1, a high-confidence SCZ gene. The results of this study shed light on common pathogenic mechanisms for SCZ and TD, and indicate that the pathogenesis of the orofacial and limbtruncal types of TD might be driven by interacting genes implicated in neurodevelopment.

The cell-surface attachment protein (Env) of the HERV-K(HML-2) lineage of endogenous retroviruses is a potentially attractive tumour-associated antigen for anti-cancer immunotherapy. The human genome contains around 100 integrated copies (called proviruses or loci) of the HERV-K(HML-2) virus and we argue that it is important for therapy development to know which and how many of these contribute to protein expression, and how this varies across tissues. We measured relative provirus expression in HERV-K(HML-2), using enriched RNA-Seq analysis with both short- and long-read sequencing, in three Mantle Cell Lymphoma cell lines (JVM2, Granta519 and REC1). We also confirmed expression of the Env protein in two of our cell lines using Western blotting, and analysed provirus expression data from all other relevant published studies. Firstly, in both our and other reanalysed studies, approximately 10% of the transcripts mapping to HERV-K(HML-2) came from Env-encoding proviruses. Secondly, in one cell line the majority of the protein expression appears to come from one provirus (12q14.1). Thirdly, we find a strong tissue-specific pattern of provirus expression. A possible dependency of Env expression on a single provirus, combined with the earlier observation that this provirus is not present in all individuals and a general pattern of tissue-specific expression among proviruses, has serious implications for future HERV-K(HML-2)-targeted immunotherapy. Further research into HERV-K(HML-2) as a possible tumour-associated antigen in blood cancers requires a more targeted, proteome-based, screening protocol that will consider these polymorphisms within HERV-K(HML-2). We include a plan (and necessary alignments) for such work.

One lineage of human endogenous retroviruses (HERVs), HERV-K(HML2), is upregulated in many cancers, some autoimmune/inflammatory diseases, and HIV-infected cells. Despite 3 decades of research, it is not known if these viruses play a causal role in disease, and there has been recent interest in whether they can be used as immunotherapy targets. Resolution of both these questions will be helped by an ability to distinguish between the effects of different integrated copies of the virus (loci). Research so far has concentrated on the 20 or so recently integrated loci that, with one exception, are in the human reference genome sequence. However, this viral lineage has been copying in the human population within the last million years, so some loci will inevitably be present in the human population but absent from the reference sequence. We therefore performed the first detailed search for such loci by mining whole-genome sequences generated by next-generation sequencing. We found a total of 17 loci, and the frequency of their presence ranged from only 2 of the 358 individuals examined to over 95% of them. On average, each individual had six loci that are not in the human reference genome sequence. Comparing the number of loci that we found to an expectation derived from a neutral population genetic model suggests that the lineage was copying until at least ∼250,000 years ago. IMPORTANCE About 5% of the human genome sequence is composed of the remains of retroviruses that over millions of years have integrated into the chromosomes of egg and/or sperm precursor cells. There are indications that protein expression of these viruses is higher in some diseases, and we need to know (i) whether these viruses have a role in causing disease and (ii) whether they can be used as immunotherapy targets in some of them. Answering both questions requires a better understanding of how individuals differ in the viruses that they carry. We carried out the first careful search for new viruses in some of the many human genome sequences that are now available thanks to advances in sequencing technology. We also compared the number that we found to a theoretical expectation to see if it is likely that these viruses are still replicating in the human population today.

Craniosynostosis, the premature fusion of the cranial sutures, is a heterogeneous disorder with a prevalence of ~1 in 2,200 (refs. 1,2). A specific genetic etiology can be identified in ~21% of cases3, including mutations of TWIST1, which encodes a class II basic helix-loop-helix (bHLH) transcription factor, and causes Saethre-Chotzen syndrome, typically associated with coronal synostosis4-6. Starting with an exome sequencing screen, we identified 38 heterozygous TCF12 mutations in 347 samples from unrelated individuals with craniosynostosis. The mutations predominantly occurred in patients with coronal synostosis and accounted for 32% and 10% of subjects with bilateral and unilateral pathology, respectively. TCF12 encodes one of three class I E-proteins that heterodimerize with class II bHLH proteins such as TWIST1. We show that TCF12 and TWIST1 act synergistically in a transactivation assay, and that mice doubly heterozygous for loss-of-function mutations in Tcf12 and Twist1 exhibit severe coronal synostosis. Hence, the dosage of TCF12/TWIST1 heterodimers is critical for coronal suture development.

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Reactome (http://www.reactome.org) is an expert-authored, peer-reviewed knowledgebase of human reactions and pathways that functions as a data mining resource and electronic textbook. Its current release includes 2975 human proteins, 2907 reactions and 4455 literature citations. A new entity-level pathway viewer and improved search and data mining tools facilitate searching and visualizing pathway data and the analysis of user-supplied high-throughput data sets. Reactome has increased its utility to the model organism communities with improved orthology prediction methods allowing pathway inference for 22 species and through collaborations to create manually curated Reactome pathway datasets for species including Arabidopsis, Oryza sativa (rice), Drosophila and Gallus gallus (chicken). Reactome's data content and software can all be freely used and redistributed under open source terms.

The genome expression pattern is strongly modified during the heat shock response (HSR) to form an adaptive state. This may be partly achieved by modulating microRNA levels that control the expression of a great number of genes that are embedded within the gene circuitry. Here, we investigated the cross-talk between two highly conserved and universal house-keeping systems, the HSR and microRNA machinery, in Drosophila melanogaster . We demonstrated that pronounced interstrain differences in the microRNA levels are alleviated after heat shock (HS) to form a uniform microRNA pattern. However, individual strains exhibit different patterns of microRNA expression during the course of recovery. Importantly, HS-regulated microRNAs may target functionally similar HS-responsive genes involved in the HSR. Despite the observed general downregulation of primary microRNA precursor expression as well as core microRNA pathway genes after HS, the levels of many mature microRNAs are upregulated. This indicates that the regulation of miRNA expression after HS occurs at transcriptional and post-transcriptional levels. It was also shown that deletion of all hsp70 genes had no significant effect on microRNA biogenesis but might influence the dynamics of microRNA expression during the HSR.

InterPro is an integrated resource for protein families, domains and functional sites, which integrates the following protein signature databases: PROSITE, PRINTS, ProDom, Pfam, SMART, TIGRFAMs, PIRSF, SUPERFAMILY, Gene3D and PANTHER. The latter two new member databases have been integrated since the last publication in this journal. There have been several new developments in InterPro, including an additional reading field, new database links, extensions to the web interface and additional match XML files. InterPro has always provided matches to UniProtKB proteins on the website and in the match XML file on the FTP site. Additional matches to proteins in UniParc (UniProt archive) are now available for download in the new match XML files only. The latest InterPro release (13.0) contains more than 13 000 entries, covering over 78% of all proteins in UniProtKB. The database is available for text- and sequence-based searches via a webserver (http://www.ebi.ac.uk/interpro), and for download by anonymous FTP (ftp://ftp.ebi.ac.uk/pub/databases/interpro). The InterProScan search tool is now also available via a web service at http://www.ebi.ac.uk/Tools/webservices/WSInterProScan.html.

InterPro, an integrated documentation resource for protein families, protein domains, and functional sites, was developed to amalgamate the individual efforts of the PROSITE, PRINTS, Pfam, and ProDom databases. InterPro can be used for the computational functional classification of newly determined amino acid sequences that lack biochemical characterization and for comparative genome analysis. InterPro contains over 3500 entries for more than 1 000 000 hits in SWISS-PROT and TrEMBL. The database is accessible for text-and sequence-based searches at http://www.ebi.ac.uk/interpro/. InterPro was used for the complete analysis of the proteome of the pathogenic microorganism Mycobacterium tuberculosis and the comparison with the predicted protein-coding sequences of the complete genomes of Bacillus subtilis and Escherichia coli. It was found that 64.8% of proteins in the proteome of M. tuberculosis matched InterPro entries and can be classified by their functions. The comparison with B. subtilis and E. coli provided information on the most common protein families and domains and on the most highly represented protein families in each organism. Thus, InterPro is a useful tool for general comparison of complete proteomes and their compositions.

Manual curation has long been held to be the “gold standard” for functional annotation of DNA sequence. Our experience with the annotation of more than 20,000 full-length cDNA sequences revealed problems with this approach, including inaccurate and inconsistent assignment of gene names, as well as many good assignments that were difficult to reproduce using only computational methods. For the FANTOM2 annotation of more than 60,000 cDNA clones, we developed a number of methods and tools to circumvent some of these problems, including an automated annotation pipeline that provides high-quality preliminary annotation for each sequence by introducing an “uninformative filter” that eliminates uninformative annotations, controlled vocabularies to accurately reflect both the functional assignments and the evidence supporting them, and a highly refined, Web-based manual annotation tool that allows users to view a wide array of sequence analyses and to assign gene names and putative functions using a consistent nomenclature. The ultimate utility of our approach is reflected in the low rate of reassignment of automated assignments by manual curation. Based on these results, we propose a new standard for large-scale annotation, in which the initial automated annotations are manually investigated and then computational methods are iteratively modified and improved based on the results of manual curation.

InterPro, an integrated documentation resource of protein families, domains and functional sites, was created in 1999 as a means of amalgamating the major protein signature databases into one comprehensive resource. PROSITE, Pfam, PRINTS, ProDom, SMART and TIGRFAMs have been manually integrated and curated and are available in InterPro for text- and sequence-based searching. The results are provided in a single format that rationalises the results that would be obtained by searching the member databases individually. The latest release of InterPro contains 5629 entries describing 4280 families, 1239 domains, 95 repeats and 15 post-translational modifications. Currently, the combined signatures in InterPro cover more than 74% of all proteins in SWISS-PROT and TrEMBL, an increase of nearly 15% since the inception of InterPro. New features of the database include improved searching capabilities and enhanced graphical user interfaces for visualisation of the data. The database is available via a webserver (http://www.ebi.ac.uk/interpro) and anonymous FTP (ftp://ftp.ebi.ac.uk/pub/databases/interpro).

The ERYTHROPOIETIN (EPO) gene is regulated by the transcription factor Hypoxia Inducible Factor-α (HIF-α). In this pathway, Prolyl Hydroxylase Domain protein 2 (PHD2) hydroxylates two prolyl residues in HIF-α, which in turn promotes HIF-α degradation by the von Hippel Lindau (VHL) protein. Evidence that HIF-2α is the important isoform for EPO regulation in humans comes from the recent observation that mutations in the HIF2A gene are associated with cases of erythrocytosis. We report here a new erythrocytosis-associated mutation, p.Asp539Glu, in the HIF2A gene. Similar to all reported cases, the affected residue is in close vicinity and C-terminal to the primary hydroxylation site in HIF-2α, Pro531. This mutation, however, is notable in producing a rather subtle amino acid substitution. Nonetheless, we find that this mutation compromises binding of HIF-2α to both PHD2 and VHL, and we propose that this mutation is the cause of erythrocytosis in this individual.

We have undertaken the first comparative pilot gene discovery analysis of approximately 25,000 random genomic and expressed sequence tags (ESTs) from three species of Plasmodium, the infectious agent that causes malaria. A total of 5482 genome survey sequences (GSSs) and 5582 ESTs were generated from mung bean nuclease (MBN) and cDNA libraries, respectively, of the ANKA line of the rodent malaria parasite Plasmodium berghei, and 10,874 GSSs generated from MBN libraries of the Salvador I and Belem lines of Plasmodium vivax, the most geographically wide-spread human malaria pathogen. These tags, together with 2438 Plasmodium falciparum sequences present in GenBank, were used to perform first-pass assembly and transcript reconstruction, and non-redundant consensus sequence datasets created. The datasets were compared against public protein databases and more than 1000 putative new Plasmodium proteins identified based on sequence similarity. Homologs of previously characterized Plasmodium genes were also identified, increasing the number of P. vivax and P. berghei sequences in public databases at least 10-fold. Comparative studies with other species of Apicomplexa identified interesting homologs of possible therapeutic or diagnostic value. A gene prediction program, Phat, was used to predict probable open reading frames for proteins in all three datasets. Predicted and non-redundant BLAST-matched proteins were submitted to InterPro, an integrated database of protein domains, signatures and families, for functional classification. Thus a partial predicted proteome was created for each species. This first comparative analysis of Plasmodium protein coding sequences represents a valuable resource for further studies on the biology of this important pathogen.

β-III spectrin is present in the brain and is known to be important in the function of the cerebellum. Heterozygous mutations in SPTBN2, the gene encoding β-III spectrin, cause Spinocerebellar Ataxia Type 5 (SCA5), an adult-onset, slowly progressive, autosomal-dominant pure cerebellar ataxia. SCA5 is sometimes known as “Lincoln ataxia,” because the largest known family is descended from relatives of the United States President Abraham Lincoln. Using targeted capture and next-generation sequencing, we identified a homozygous stop codon in SPTBN2 in a consanguineous family in which childhood developmental ataxia co-segregates with cognitive impairment. The cognitive impairment could result from mutations in a second gene, but further analysis using whole-genome sequencing combined with SNP array analysis did not reveal any evidence of other mutations. We also examined a mouse knockout of β-III spectrin in which ataxia and progressive degeneration of cerebellar Purkinje cells has been previously reported and found morphological abnormalities in neurons from prefrontal cortex and deficits in object recognition tasks, consistent with the human cognitive phenotype. These data provide the first evidence that β-III spectrin plays an important role in cortical brain development and cognition, in addition to its function in the cerebellum; and we conclude that cognitive impairment is an integral part of this novel recessive ataxic syndrome, Spectrin-associated Autosomal Recessive Cerebellar Ataxia type 1 (SPARCA1). In addition, the identification of SPARCA1 and normal heterozygous carriers of the stop codon in SPTBN2 provides insights into the mechanism of molecular dominance in SCA5 and demonstrates that the cell-specific repertoire of spectrin subunits underlies a novel group of disorders, the neuronal spectrinopathies, which includes SCA5, SPARCA1, and a form of West syndrome., Author Summary β-III spectrin is present in the brain and is known to be important in the function of the cerebellum. Mutations in β-III spectrin cause spinocerebellar ataxia type 5 (SCA5), sometimes called Lincoln ataxia because it was first described in the relatives of United States President Abraham Lincoln. This is generally an adult-onset progressive cerebellar disorder. Recessive mutations have not previously been described in any of the brain spectrins. We identified a homozygous mutation in SPTBN2, which causes a more severe disorder than SCA5, with a developmental cerebellar ataxia, which is present from childhood; in addition there is marked cognitive impairment. We call this novel condition SPARCA1 (Spectrin-associated Autosomal Recessive Cerebellar Ataxia type 1). This condition could be caused by two separate gene mutations; but we show, using a combination of genome-wide mapping, whole-genome sequencing, and detailed behavioural and neuropathological analysis of a β-III spectrin mouse knockout, that both the ataxia and cognitive impairment are caused by the recessive mutations in β-III spectrin. SPARCA1 is one of a family of neuronal spectrinopathies and illustrates the importance of spectrins in brain development and function.

We consider the problem of biological complexity via a projection of protein-coding genes of complex organisms onto the functional space of the proteome. The latter can be defined as a set of all functions committed by proteins of an organism. Alternative splicing (AS) allows an organism to generate diverse mature RNA transcripts from a single mRNA strand and thus it could be one of the key mechanisms of increasing of functional complexity of the organism's proteome and a driving force of biological evolution. Thus, the projection of transcription units (TU) and alternative splice-variant (SV) forms onto proteome functional space could generate new types of relational networks (e.g. SV-protein function networks, SFN) and lead to discoveries of novel evolutionarily conservative functional modules. Such types of networks might provide new reliable characteristics of organism complexity and a better understanding of the evolutionary integration and plasticity of interconnection of genome-transcriptome-proteome functions. Results We use the InterPro and UniProt databases to attribute descriptive features (keywords) to protein sequences. UniProt database includes a controlled and curated vocabulary of specific descriptors or keywords. The keywords have been assigned to a protein sequence via conserved domains or via similarity with annotated sequences. Then we consider the unique combinations of keywords as the protein functional labels (FL), which characterize the biological functions of the given protein and construct the contingency tables and graphs providing the projections of transcription units (TU) and alternative splice-variants (SV) onto all FL of the proteome of a given organism. We constructed SFNs for organisms with different evolutionary history and levels of complexity, and performed detailed statistical parameterization of the networks. Conclusions The application of the algorithm to organisms with different evolutionary history and level of biological complexity (nematode, fruit fly, vertebrata) reveals that the parameters describing SFN correlate with the complexity of a given organism. Using statistical analysis of the links of the functional networks, we propose new features of evolution of protein function acquisition. We reveal a group of genes and corresponding functions, which could be attributed to an early conservative part of the cellular machinery essential for cell viability and survival. We identify and provide characteristics of functional switches in the polyform group of TUs in different organisms. Based on comparison of mouse and human SFNs, a role of alternative splicing as a necessary source of evolution towards more complex organisms is demonstrated. The entire set of FL across many organisms could be used as a draft of the catalogue of the functional space of the proteome world.

The entire protein domain set of the proteome of an organism we call the domainome. We define the list of domains in domainome, together with the numbers of their occurrences (links to proteins) found in the proteome to be the domain-to-protein linkage profile (DPLP). We estimated the DPLP of the proteomes of the 156 complete genomes represented in the InterPro database. This work presents several quantitative measures of the complexity of a proteome based on the DPLP. For each of the 156 studied genomes, we found two large sets of domains: D1, the domains that are not replicated within any protein of the proteome and D2, the domains that occur two or more times in at least one protein of the proteome. Statistics of the observed domain-to-protein links (DPLs) for set D1 and set D2 do not exhibit simple ‘scale-free network’ properties: for D1, the distribution of DPLs in proteome follows the Generalized Discrete Pareto function and for D2, the distribution of DPLs in proteome follows the inversed gamma probability function. Dynamical range of DPLs for D1 domains is larger than for D2 domains, and this range correlates with biological complexity of organism. D1 and D2 sets exhibit significant differences of molecular functions of the corresponding proteins, biological processes, and cellular components. The statistical distributions of the number of DPLs in the proteome and the estimates of the differences between the DPLPs for pairs of organisms are used as measures of relative biological complexity of the organisms. In particular, we show quantitatively the greater domain composition complexity of the human proteins relative to that of a mouse or a rat.

Integr8 is a new web portal for exploring the biology of organisms with completely deciphered genomes. For over 190 species, Integr8 provides access to general information, recent publications, and a detailed statistical overview of the genome and proteome of the organism. The preparation of this analysis is supported through Genome Reviews, a new database of bacterial and archaeal DNA sequences in which annotation has been upgraded (compared to the original submission) through the integration of data from many sources, including the EMBL Nucleotide Sequence Database, the UniProt Knowledgebase, InterPro, CluSTr, GOA and HOGENOM. Integr8 also allows the users to customize their own interactive analysis, and to download both customized and prepared datasets for their own use. Integr8 is available at http://www.ebi.ac.uk/integr8.

International audience; InterPro, an integrated documentation resource of protein families, domains and functional sites, was created to integrate the major protein signature databases. Currently, it includes PROSITE, Pfam, PRINTS, ProDom, SMART, TIGRFAMs, PIRSF and SUPERFAMILY. Signatures are manually integrated into InterPro entries that are curated to provide biological and functional information. Annotation is provided in an abstract, Gene Ontology mapping and links to specialized databases. New features of InterPro include extended protein match views, taxonomic range information and protein 3D structure data. One of the new match views is the InterPro Domain Architecture view, which shows the domain composition of protein matches. Two new entry types were introduced to better describe InterPro entries: these are active site and binding site. PIRSF and the structure-based SUPERFAMILY are the latest member databases to join InterPro, and CATH and PANTHER are soon to be integrated. InterPro release 8.0 contains 11 007 entries, representing 2573 domains, 8166 families, 201 repeats, 26 active sites, 21 binding sites and 20 post-translational modification sites. InterPro covers over 78% of all proteins in the Swiss-Prot and TrEMBL components of UniProt. The database is available for text- and sequence-based searches via a webserver (http://www.ebi.ac.uk/interpro), and for download by anonymous FTP (ftp://ftp.ebi.ac.uk/pub/databases/interpro).

Summary In the June 5th 2012 issue of Current Biology , Agoni et al. [1] reported finding 14 endogenous retrovirus (ERV) loci in the genome sequences of Neanderthal and/or Denisovan fossils (both ∼40,000 years old) that are not found in the human reference genome sequence. The authors [1] concluded that these retroviruses were infecting the germline of these archaic hominins at or subsequent to their divergence from modern humans (∼400,000 years ago). However, in our search for unfixed ERVs in the modern human population, we have found most of these loci. We explain this apparent contradiction using population genetic theory and suggest that it illustrates an important phenomenon for the study of transposable elements such as ERVs.

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology., An international team has systematically validated and annotated just over 21,000 human genes using full-length cDNA, thereby providing a valuable new resource for the human genetics community

A general overview of the protein sequence set for the mouse transcriptome produced during the FANTOM2 sequencing project is presented here. We applied different algorithms to characterize protein sequences derived from a nonredundant representative protein set (RPS) and a variant protein set (VPS) of the mouse transcriptome. The functional characterization and assignment of Gene Ontology terms was done by analysis of the proteome using InterPro. The Superfamily database analyses gave a detailed structural classification according to SCOP and provide additional evidence for the functional characterization of the proteome data. The MDS database analysis revealed new domains which are not presented in existing protein domain databases. Thus the transcriptome gives us a unique source of data for the detection of new functional groups. The data obtained for the RPS and VPS sets facilitated the comparison of different patterns of protein expression. A comparison of other existing mouse and human protein sequence sets (e.g., the International Protein Index) demonstrates the common patterns in mammalian proteomes. The analysis of the membrane organization within the transcriptome of multiple eukaryotes provides valuable statistics about the distribution of secretory and transmembrane proteins

The Proteome Analysis database (http://www.ebi.ac.uk/proteome/) has been developed by the Sequence Database Group at EBI utilizing existing resources and providing comparative analysis of the predicted protein coding sequences of the complete genomes of bacteria, archeae and eukaryotes. Three main projects are used, InterPro, CluSTr and GO Slim, to give an overview on families, domains, sites, and functions of the proteins from each of the complete genomes. Complete proteome analysis is available for a total of 89 proteome sets. A specifically designed application enables InterPro proteome comparisons for any one proteome against any other one or more of the proteomes in the database.

The applications of InterPro span a range of biologically important areas that includes automatic annotation of protein sequences and genome analysis. In automatic annotation of protein sequences InterPro has been utilised to provide reliable characterisation of sequences, identifying them as candidates for functional annotation. Rules based on the InterPro characterisation are stored and operated through a database called RuleBase. RuleBase is used as the main tool in the sequence database group at the EBI to apply automatic annotation to unknown sequences. The annotated sequences are stored and distributed in the TrEMBL protein sequence database. InterPro also provides a means to carry out statistical and comparative analyses of whole genomes. In the Proteome Analysis Database, InterPro analyses have been combined with other analyses based on CluSTr, the Gene Ontology (GO) and structural information on the proteins.

A principal challenge currently facing biologists is how to connect the complete DNA sequence of an organism to its development and behaviour. Large-scale targeted-deletions have been successful in defining gene functions in the single-celled yeast Saccharomyces cerevisiae, but comparable analyses have yet to be performed in an animal. Here we describe the use of RNA interference to inhibit the function of ∼86% of the 19,427 predicted genes of C. elegans. We identified mutant phenotypes for 1,722 genes, about two-thirds of which were not previously associated with a phenotype. We find that genes of similar functions are clustered in distinct, multi-megabase regions of individual chromosomes; genes in these regions tend to share transcriptional profiles. Our resulting data set and reusable RNAi library of 16,757 bacterial clones will facilitate systematic analyses of the connections among gene sequence, chromosomal location and gene function in C. elegans., R.S.K. was supported by a Howard Hughes Medical Institute Predoctoral Fellowship; A.G.F. by a US Army Breast Cancer Research Fellowship; Y.D., R.D., M.G., D.P.W. and P.Z. by the Wellcome Trust; G.P. by the Canadian Institute of Health Research and the Wellcome Trust; A.K. by the European Molecular Biology Laboratory; N.L.B. by the European Molecular Biology Organization; S.M. by the Centro de Investigacion del Cancer; M.S. by a Swiss National Science Foundation fellowship and J.A. by a Wellcome Trust Senior Research Fellowship.

InterPro (http://www.ebi.ac.uk/interpro/) is an integrated documentation resource for protein families, domains and sites, developed initially as a means of rationalizing the complementary efforts of the PROSITE, PRINTS, Pfam and ProDom database projects. It is a useful resource that aids the functional classification of proteins. Almost 90% of the actinopterygii protein sequences from SWISS-PROT and TrEMBL can be classified using InterPro. Over 30% of the actinopterygii protein sequences currently in SWISS-PROT and TrEMBL are of mitochondrial origin, the majority of which belong to the cytochrome b/b6 family. InterPro also gives insights into the domain composition of the classified proteins and has applications in the functional classification of newly determined sequences lacking biochemical characterization, and in comparative genome analysis. A comparison of the actinopterygii protein sequences against the sequences of other eukaryotes confirms the high representation of eukaryotic protein kinase in the organisms studied. The comparisons also show that, based on InterPro families, the trans-species evolution of MHC class I and II molecules in mammals and teleost fish can be recognized.

Signature databases are vital tools for identifying distant relationships in novel sequences and hence for inferring protein function. InterPro is an integrated documentation resource for protein families, domains and functional sites, which amalgamates the efforts of the PROSITE, PRINTS, Pfam and ProDom database projects. Each InterPro entry includes a functional description, annotation, literature references and links back to the relevant member database(s). Release 2.0 of InterPro (October 2000) contains over 3000 entries, representing families, domains, repeats and sites of post-translational modification encoded by a total of 6804 different regular expressions, profiles, fingerprints and Hidden Markov Models. Each InterPro entry lists all the matches against SWISS-PROT and TrEMBL (more than 1 000 000 hits from 462 500 proteins in SWISS-PROT and TrEMBL). The database is accessible for text- and sequence-based searches at http://www.ebi.ac.uk/interpro/. Questions can be emailed to interhelp@ebi.ac.uk.

The SWISS-PROT group at EBI has developed the Proteome Analysis Database utilising existing resources and providing comparative analysis of the predicted protein coding sequences of the complete genomes of bacteria, archaea and eukaryotes (http://www.ebi.ac.uk/proteome/). The two main projects used, InterPro and CluSTr, give a new perspective on families, domains and sites and cover 31–67% (InterPro statistics) of the proteins from each of the complete genomes. CluSTr covers the three complete eukaryotic genomes and the incomplete human genome data. The Proteome Analysis Database is accompanied by a program that has been designed to carry out InterPro proteome comparisons for any one proteome against any other one or more of the proteomes in the database.

The recent discovery of a variety of somatic splicesomal mutations in the myelodysplastic syndromes (MDS) has revealed a new leukaemogenic pathway involving spliceosomal dysfunction. Pre-mRNA splicing proceeds by way of two phosphoester transfer reactions and is catalyzed by the spliceosome, which consists of the U1, U2, U4/U6, and U5 small nuclear ribonucleoproteins (snRNPs) and numerous non-snRNP proteins. The snRNPs are involved in recognising short conserved sequences of the pre-mRNA, including the 5′ and 3′ splice sites and the branch site, and in positioning the reactive nucleotides for catalysis. The spliceosome is a dynamic molecular machine, undergoing several major structural rearrangements during its functional cycle. Mutation of the Splicing Factor 3B, subunit 1 (SF3B1) gene is common in MDS, occurring in over 70% of patients whose disease is characterised by ring sideroblasts (RARS). The close association between SF3B1 mutation and ring sideroblasts is consistent with a causal relationship, and makes this the first gene to be strongly associated with a specific feature of MDS. Sf3b1 heterozygous knockout mice show the presence of ringed sideroblasts. In order to investigate the role of SF3B1 haploinsufficiency in MDS we have silenced SF3B1 using siRNA in the myeloid cell lines K562, TF-1, SKM1, HeL and OCIM2. Cell growth was impaired in all the cell lines with SF3B1 knockdown. Using Flow Cytometry, cell cycle analysis showed a significant increase in cells in the sub-G0 phase as well as G2/M arrest in the cell lines. We also observed impaired erythroid differentiation in hemin treated K562 and TF-1 cell lines with SF3B1 knockdown. Gene expression profiling (GEP) was performed in two cell lines with SF3B1 knockdown (K562 and TF1). Deregulated pathways and gene ontology categories included cell cycle regulation and alternative splicing using Ingenuity Pathway Analysis. We next performed Gene Set Enrichment Analysis (GSEA). The GSEA showed a significant enrichment of nonsense-mediated mRNA decay (NMD) genes that were up-regulated in cells with SF3B1 knockdown, suggesting NMD activation following SF3B1 silencing. We used Human Exon-Junction arrays (Affymetrix) to evaluate global transcript exon usage in the K562 and TF1 cell lines with SF3B1 knockdown. We observed significant differential exon usage in genes involved in RNA degradation, spliceosome, cell cycle and apoptosis. We further observed aberrant splicing of the candidate gene ABCB7 showing exon skipping and TP53 gene showing exon skipping as well as intron retention. We have investigated the changes in the transcriptome in CD34+ cells from MDS patients with SF3B1 mutation by RNA sequencing and found many genes showing significant differential exon usage including CCND1, EIF3B, FKBP1A, BCL2 and RB1. Using Ingenuity Pathway Analysis we identified alternative splicing pattern of genes involved in cell cycle, RNA processing, mTOR signalling and P53 signalling pathways. We have studied CD34+ cells from MDS patients with SF3B1 mutation in vitro and observed impairment in cell growth compared to CD34+ cells from healthy controls or from MDS patients without splicing mutations. In colony forming assays we observed a decrease in the number of erythroid or myeloid colonies derived from CD34+ cells of patients with SF3B1 mutation compared to patient CD34+ cells without splicing factor mutation. The identification of SF3B1 downstream targets in SF3B1 mutant and wild-type erythroid and myeloid colonies from MDS patients is in progress using RNA sequencing. Our results show that knockdown of SF3B1 in haematopoietic cell lines results in impaired cell growth, deregulated global gene expression and aberrant splicing. Studies of the haematopoietic progenitor CD34+ cells of patients with SF3B1 mutation show impaired cell growth and erythroid differentiation as well as deregulation of many pathways including the cell cycle and RNA processing. The identification of the key target genes affected by the common splicing mutations in MDS is critical to our understanding of how the mutations contribute to the pathogenesis of this disorder. Disclosures: Maciejewski: NIH: Research Funding; Aplastic anemia&MDS International Foundation: Research Funding.

The Tn5 insertion into the genome of Rhizobium leguminosarum bv viciae VF39, resulting in non-mucoid growth and formation of non-N2-fixing nodule-like structures on Vicia faba plants, was mapped within a 1.4-kb EcoRV-SacI fragment. Nucleotide sequence analysis revealed an ORF (pss4) of 263 amino acids (aa). Three transcription start points (tsp) were determined. Two of them were localized upstream from the first GTG codon; the third tsp was mapped in front of the second putative start codon (GTG) corresponding to Val64 of the Pss4 aa sequence. The expression of pss4 in a T7 RNA polymerase/promoter system produced a single approx. 29-kDa protein. Pss4 reveals similarity to several proteins involved in polysaccharide biosynthesis in various Rhizobium species. A nearly complete homology was found with PssA from Rl biovar phaseoli 8002 [Borthakur et al., Mol. Gen. Genet. 213 (1988) 155–162], except that Pss4 has an additional 63 aa on its N terminus.

Motivation: The SWISS-PROT group at the EBI has developed the Proteome Analysis Database utilizing existing resources and providing comprehensive and integrated comparative analysis of the predicted protein coding sequencesof the complete genomes of bacteria, archaea and eukaryotes. The Proteome Analysis Database is accompanied by a program that has been designed to carry out interactive InterPro proteome comparisons for any one proteome against any other one or more of the proteomes in the database. Availability: http://www.ebi.ac.uk/proteome/comparisons.html Contact: alex@ebi.ac.uk; proteome@ebi.ac.uk * To whom all correspondence should be addressed.