Resolving the fibrotic niche of human liver cirrhosis at single-cell level.

Liver cirrhosis is a major cause of death worldwide and is characterized by extensive fibrosis. There are currently no effective antifibrotic therapies available. To obtain a better understanding of the cellular and molecular mechanisms involved in disease pathogenesis and enable the discovery of therapeutic targets, here we profile the transcriptomes of more than 100,000 single human cells, yielding molecular definitions for non-parenchymal cell types that are found in healthy and cirrhotic human liver. We identify a scar-associated TREM2+CD9+ subpopulation of macrophages, which expands in liver fibrosis, differentiates from circulating monocytes and is pro-fibrogenic. We also define ACKR1+ and PLVAP+ endothelial cells that expand in cirrhosis, are topographically restricted to the fibrotic niche and enhance the transmigration of leucocytes. Multi-lineage modelling of ligand and receptor interactions between the scar-associated macrophages, endothelial cells and PDGFRα+ collagen-producing mesenchymal cells reveals intra-scar activity of several pro-fibrogenic pathways including TNFRSF12A, PDGFR and NOTCH signalling. Our work dissects unanticipated aspects of the cellular and molecular basis of human organ fibrosis at a single-cell level, and provides a conceptual framework for the discovery of rational therapeutic targets in liver cirrhosis.

ELIXIR-UK role in bioinformatics training at the national level and across ELIXIR.

ELIXIR-UK is the UK node of ELIXIR, the European infrastructure for life science data. Since its foundation in 2014, ELIXIR-UK has played a leading role in training both within the UK and in the ELIXIR Training Platform, which coordinates and delivers training across all ELIXIR members. ELIXIR-UK contributes to the Training Platform's coordination and supports the development of training to address key skill gaps amongst UK scientists. As part of this work it acts as a conduit for nationally-important bioinformatics training resources to promote their activities to the ELIXIR community. ELIXIR-UK also leads ELIXIR's flagship Training Portal, TeSS, which collects information about a diverse range of training and makes it easily accessible to the community. ELIXIR-UK also works with others to provide key digital skills training, partnering with the Software Sustainability Institute to provide Software Carpentry training to the ELIXIR community and to establish the Data Carpentry initiative, and taking a lead role amongst national stakeholders to deliver the StaTS project - a coordinated effort to drive engagement with training in statistics.

Additive effects of oxytocin receptor gene polymorphisms on reward circuitry in youth with autism.

Several common alleles in the oxytocin receptor gene (OXTR) are associated with altered brain function in reward circuitry in neurotypical adults and may increase risk for autism spectrum disorders (ASD). Yet, it is currently unknown how variation in the OXTR relates to brain functioning in individuals with ASD, and, critically, whether neural endophenotypes vary as a function of aggregate genetic risk. Here, for we believe the first time, we use a multi-locus approach to examine how genetic variation across several OXTR single-nucleotide polymorphisms (SNPs) affect functional connectivity of the brain's reward network. Using data from 41 children with ASD and 41 neurotypical children, we examined functional connectivity of the nucleus accumbens (NAcc) - a hub of the reward network - focusing on how connectivity varies with OXTR risk-allele dosage. Youth with ASD showed reduced NAcc connectivity with other areas in the reward circuit as a function of increased OXTR risk-allele dosage, as well as a positive association between risk-allele dosage and symptom severity, whereas neurotypical youth showed increased NAcc connectivity with frontal brain regions involved in mentalizing. In addition, we found that increased NAcc-frontal cortex connectivity in typically developing youth was related to better scores on a standardized measure of social functioning. Our results indicate that cumulative genetic variation on the OXTR impacts reward system connectivity in both youth with ASD and neurotypical controls. By showing differential genetic effects on neuroendophenotypes, these pathways elucidate mechanisms of vulnerability versus resilience in carriers of disease-associated risk alleles.

A genome-wide association study implicates the APOE locus in nonpathological cognitive ageing.

Cognitive decline is a feared aspect of growing old. It is a major contributor to lower quality of life and loss of independence in old age. We investigated the genetic contribution to individual differences in nonpathological cognitive ageing in five cohorts of older adults. We undertook a genome-wide association analysis using 549 692 single-nucleotide polymorphisms (SNPs) in 3511 unrelated adults in the Cognitive Ageing Genetics in England and Scotland (CAGES) project. These individuals have detailed longitudinal cognitive data from which phenotypes measuring each individual's cognitive changes were constructed. One SNP--rs2075650, located in TOMM40 (translocase of the outer mitochondrial membrane 40 homolog)--had a genome-wide significant association with cognitive ageing (P=2.5 × 10(-8)). This result was replicated in a meta-analysis of three independent Swedish cohorts (P=2.41 × 10(-6)). An Apolipoprotein E (APOE) haplotype (adjacent to TOMM40), previously associated with cognitive ageing, had a significant effect on cognitive ageing in the CAGES sample (P=2.18 × 10(-8); females, P=1.66 × 10(-11); males, P=0.01). Fine SNP mapping of the TOMM40/APOE region identified both APOE (rs429358; P=3.66 × 10(-11)) and TOMM40 (rs11556505; P=2.45 × 10(-8)) as loci that were associated with cognitive ageing. Imputation and conditional analyses in the discovery and replication cohorts strongly suggest that this effect is due to APOE (rs429358). Functional genomic analysis indicated that SNPs in the TOMM40/APOE region have a functional, regulatory non-protein-coding effect. The APOE region is significantly associated with nonpathological cognitive ageing. The identity and mechanism of one or multiple causal variants remain unclear.

Phospholipase Czeta, the trigger of egg activation in mammals, is present in a non-mammalian species.

The activation of the egg to begin development into an embryo is triggered by a sperm-induced increase in intracellular egg Ca2+. There has been much controversy about how the sperm induces this fundamental developmental event, but recent studies suggest that, in mammals, egg activation is triggered by a testis-specific phospholipase C: PLCzeta. Since the discovery of PLCzeta, it has been unclear whether its role in triggering egg activation is common to all vertebrates, or is confined to mammals. Here, we demonstrate for the first time that PLCzeta is present in a non-mammalian vertebrate. Using genomic and cDNA databases, we have identified the cDNA encoding a PLCzeta orthologue in the domestic chicken that, like the mammalian isoforms, is a testis-specific gene. The chicken PLCzeta cDNA is 2152 bp in size and encodes an open reading frame of 639 amino acids. When injected into mouse oocytes, chicken PLCzeta cRNA triggers Ca2+ oscillations, indicating that it has functional properties similar to those of mammalian PLCzeta. Our findings suggest that PLCzeta may have a universal role in triggering egg activation in vertebrates.

Stromal cell-derived receptor 2 and cytochrome b561 are functional ferric reductases.

Iron has a variety of functions in cellular organisms ranging from electron transport and DNA synthesis to adenosine triphosphate (ATP) and neurotransmitter synthesis. Failure to regulate the homeostasis of iron can lead to cognition and demyelination disorders when iron levels are deficient, and to neurodegenerative disorders when iron is in excess. In this study we show that three members of the b561 family of predicted ferric reductases, namely mouse cytochrome b561 and mouse and fly stromal cell-derived receptor 2 (SDR2), have ferric reductase activity. Given that a fourth member, duodenal cytochrome b (Dcytb), has previously been shown to be a ferric reductase, it is likely that all remaining members of this family also exhibit this activity. Furthermore, we show that the rat sdr2 message is predominantly expressed in the liver and kidney, with low expression in the duodenum. In hypotransferrinaemic (hpx) mice, sdr2 expression in the liver and kidney is reduced, suggesting that it may be regulated by iron. Moreover, we demonstrate the presence of mouse sdr2 in the choroid plexus and in the ependymal cells lining the four ventricles, through in situ hybridization analysis.

Cloning and developmental expression analysis of ltd-1, the Caenorhabditis elegans homologue of the mouse kyphoscoliosis (ky) gene.

We have characterized the developmental expression pattern of the Caenorhabditis elegans homologue of the mouse ky gene. The Ky protein has a putative key function in muscle development and has homologues in invertebrates, fungi and a cyanobacterium. The C. elegans Ky homologue gene has been named ltd-1 for LIM and transglutaminase domains gene. The LTD-1::GFP construct is expressed in developing hypodermal cells from the twofold stage embryo through adulthood. These data define the ltd-1 gene as a novel marker for C. elegans epithelial cell development.

Novel protein domains and repeats in Drosophila melanogaster: insights into structure, function, and evolution.

Sequence database searching methods such as BLAST, are invaluable for predicting molecular function on the basis of sequence similarities among single regions of proteins. Searches of whole databases however, are not optimized to detect multiple homologous regions within a single polypeptide. Here we have used the prospero algorithm to perform self-comparisons of all predicted Drosophila melanogaster gene products. Predicted repeats, and their homologs from all species, were analyzed further to detect hitherto unappreciated evolutionary relationships. Results included the identification of novel tandem repeats in the human X-linked retinitis pigmentosa type-2 gene product, repeated segments in cystinosin, associated with a defect in cystine transport, and 'nested' homologous domains in dysferlin, whose gene is mutated in limb girdle muscular dystrophy. Novel signaling domain families were found that may regulate the microtubule-based cytoskeleton and ubiquitin-mediated proteolysis, respectively. Two families of glycosyl hydrolases were shown to contain internal repetitions that hint at their evolution via a piecemeal, modular approach. In addition, three examples of fruit fly genes were detected with tandem exons that appear to have arisen via internal duplication. These findings demonstrate how completely sequenced genomes can be exploited to further understand the relationships between molecular structure, function, and evolution.

A new sequence motif linking lissencephaly, Treacher Collins and oral-facial-digital type 1 syndromes, microtubule dynamics and cell migration.

A previously unidentified sequence motif has been identified in the products of genes mutated in Miller-Dieker lissencephaly, Treacher Collins, oral-facial-digital type 1 and contiguous syndrome ocular albinism with late onset sensorineural deafness syndromes. An additional homologous motif was detected in a gene product fused to the fibroblast growth factor receptor type 1 in patients with an atypical stem cell myeloproliferative disorder. In total, over 100 eukaryotic intracellular proteins are shown to possess a LIS1 homology (LisH) motif, including several katanin p60 subunits, muskelin, tonneau, LEUNIG, Nopp140, aimless and numerous WD repeat-containing beta-propeller proteins. It is suggested that LisH motifs contribute to the regulation of microtubule dynamics, either by mediating dimerization, or else by binding cytoplasmic dynein heavy chain or microtubules directly. The predicted secondary structure of LisH motifs, and their occurrence in homologues of Gbeta beta-propeller subunits, suggests that they are analogues of Ggamma subunits, and might associate with the periphery of beta-propeller domains. The finding of LisH motifs in both treacle and Nopp140 reinforces previous observations of functional similarities between these nucleolar proteins. Uncharacterized LisH motif-containing proteins represent candidates for other diseases associated with aberrant microtubule dynamics and defects of cell migration, nucleokinesis or chromosome segregation.

Mutations in the fukutin-related protein gene (FKRP) cause a form of congenital muscular dystrophy with secondary laminin alpha2 deficiency and abnormal glycosylation of alpha-dystroglycan.

The congenital muscular dystrophies (CMD) are a heterogeneous group of autosomal recessive disorders presenting in infancy with muscle weakness, contractures, and dystrophic changes on skeletal-muscle biopsy. Structural brain defects, with or without mental retardation, are additional features of several CMD syndromes. Approximately 40% of patients with CMD have a primary deficiency (MDC1A) of the laminin alpha2 chain of merosin (laminin-2) due to mutations in the LAMA2 gene. In addition, a secondary deficiency of laminin alpha2 is apparent in some CMD syndromes, including MDC1B, which is mapped to chromosome 1q42, and both muscle-eye-brain disease (MEB) and Fukuyama CMD (FCMD), two forms with severe brain involvement. The FCMD gene encodes a protein of unknown function, fukutin, though sequence analysis predicts it to be a phosphoryl-ligand transferase. Here we identify the gene for a new member of the fukutin protein family (fukutin related protein [FKRP]), mapping to human chromosome 19q13.3. We report the genomic organization of the FKRP gene and its pattern of tissue expression. Mutations in the FKRP gene have been identified in seven families with CMD characterized by disease onset in the first weeks of life and a severe phenotype with inability to walk, muscle hypertrophy, marked elevation of serum creatine kinase, and normal brain structure and function. Affected individuals had a secondary deficiency of laminin alpha2 expression. In addition, they had both a marked decrease in immunostaining of muscle alpha-dystroglycan and a reduction in its molecular weight on western blot analysis. We suggest these abnormalities of alpha-dystroglycan are caused by its defective glycosylation and are integral to the pathology seen in MDC1C.

CHROMA: consensus-based colouring of multiple alignments for publication.

CHROMA annotates multiple protein sequence alignments by consensus to produce formatted and coloured text suitable for incorporation into other documents for publication. The package is designed to be flexible and reliable, and has a simple-to-use graphical user interface running under Microsoft Windows. Both the executables and source code for CHROMA running under Windows and Linux (portable command-line only) are freely available at http://www.lg.ndirect.co.uk/chroma. Software enquiries should be directed to CHROMA@lg.ndirect.co.uk.

Sequence variation and disease in the wake of the draft human genome.

The sequencing phase of the human genome project will soon be over. In its wake, repertoires of sequence polymorphisms among the human population are being sampled and a battery of functional genomics projects, from gene and protein expression studies to whole proteome interaction experiments, are generating vast quantities of data. Now that the data, or the means to generate data, are available it is the application of this information in enhancing our understanding of biology that represents the next formidable challenge. Two prominent issues should be considered. First, existing data must be analysed using the best methods available. The prediction of enzymatic activity for bestrophin, whose gene is mutated in Best macular dystrophy, is described in this review. This is an example of the experimentally testable hypotheses that can result from such detailed and exhaustive analyses. Secondly, the torrents of data from high-throughput studies will need to be made more accessible to all using web-based resources that integrate and digest complementary data types. The internet sites that showcase the human genome sequence are blazing a new trail. Ultimately, the success of genome sequencing and functional genomics will be measured not by the quantity and accuracy of raw data generated, but how rapidly they can be harnessed to span the divide between genotype and phenotype.

Protein repeats: structures, functions, and evolution.

Internal repetition within proteins has been a successful strategem on multiple separate occasions throughout evolution. Such protein repeats possess regular secondary structures and form multirepeat assemblies in three dimensions of diverse sizes and functions. In general, however, internal repetition affords a protein enhanced evolutionary prospects due to an enlargement of its available binding surface area. Constraints on sequence conservation appear to be relatively lax, due to binding functions ensuing from multiple, rather than, single repeats. Considerable sequence divergence as well as the short lengths of sequence repeats mean that repeat detection can be a particularly arduous task. We also consider the conundrum of how multiple repeats, which show strong structural and functional interdependencies, ever evolved from a single repeat ancestor. In this review, we illustrate each of these points by referring to six prolific repeat types (repeats in beta-propellers and beta-trefoils and tetratricopeptide, ankyrin, armadillo/HEAT, and leucine-rich repeats) and in other less-prolific but nonetheless interesting repeats.

On the evolution of protein folds: are similar motifs in different protein folds the result of convergence, insertion, or relics of an ancient peptide world?

This paper presents and discusses evidence suggesting how the diversity of domain folds in existence today might have evolved from peptide ancestors. We apply a structure similarity detection method to detect instances where localized regions of different protein folds contain highly similar sequences and structures. Results of performing an all-on-all comparison of known structures are described and compared with other recently published findings. The numerous instances of local sequence and structure similarities within different protein folds, together with evidence from proteins containing sequence and structure repeats, argues in favor of the evolution of modern single polypeptide domains from ancient short peptide ancestors (antecedent domain segments (ADSs)). In this model, ancient protein structures were formed by self-assembling aggregates of short polypeptides. Subsequently, and perhaps concomitantly with the evolution of higher fidelity DNA replication and repair systems, single polypeptide domains arose from the fusion of ADSs genes. Thus modern protein domains may have a polyphyletic origin.

Domain homologues of dopamine beta-hydroxylase and ferric reductase: roles for iron metabolism in neurodegenerative disorders?

One of the defining characteristics of neurodegenerative diseases, including Parkinson's, Alzheimer's and Huntington's diseases, is abnormal accumulations of iron, specifically in affected areas. Following injection of iron in rat brains, a relatively selective lesion of dopamine neurons, similar to parkinsonism, occurs. These observations indicate that Fe(II)-mediated generation of free radical species, by the Fenton reaction, might contribute to the pathoetiology of these diseases. Iron is known to possess multiple roles in the biosynthesis of catecholamines in dopaminergic neurons. These include, as Fe(II), facilitating the production of dopamine from phenylalanine by tyrosine hydroxylase, and as heme, assisting the recycling of ascorbate by cytochrome b-561 required for the generation of norepinephrine from dopamine by dopamine beta-hydroxylase. In this study, it is demonstrated that a human and mouse gene product, stromal cell-derived receptor 2, is a homologue of cytochrome b-561 and duodenal cytochrome b, and is thus predicted to be active as a ferric reductase. Moreover, this protein also contains a domain homologous to the N-terminal regulatory region of dopamine beta-hydroxylase. These findings from sequence analysis lead to a prediction that stromal cell-derived receptor 2 is a catecholamine-regulated ferric reductase active in the brain. Dysfunction of cytochrome b-561 or stromal cell-derived receptor 2, therefore, might predispose individuals to abnormal accumulation of Fe(III) and/or generation of cytotoxic free radicals as a consequence of a rapid cycling between Fe(III) and Fe(II). The hypothesis that aberrant ferric reductase activities are involved in the progression of neurodegenerative diseases should open up new avenues of research, and possibly therapy, for these devastating diseases.

Genome cartography through domain annotation.

The evolutionary history of eukaryotic proteins involves rapid sequence divergence, addition and deletion of domains, and fusion and fission of genes. Although the protein repertoires of distantly related species differ greatly, their domain repertoires do not. To account for the great diversity of domain contexts and an unexpected paucity of ortholog conservation, we must categorize the coding regions of completely sequenced genomes into domain families, as well as protein families.

Issues in predicting protein function from sequence.

Identifying homologues, defined as genes that arose from a common evolutionary ancestor, is often a relatively straightforward task, thanks to recent advances made in estimating the statistical significance of sequence similarities found from database searches. The extent by which homologues possess similarities in function, however, is less amenable to statistical analysis. Consequently, predicting function by homology is a qualitative, rather than quantitative, process and requires particular care to be taken. This review focuses on the various approaches that have been developed to predict function from the scale of the atom to that of the organism. Similarities in homologues' functions differ considerably at each of these different scales and also vary for different domain families. It is argued that due attention should be paid to all available clues to function, including orthologue identification, conservation of particular residue types, and the co-occurrence of domains in proteins. Pitfalls in database searching methods arising from amino acid compositional bias and database size effects are also discussed.