An inherited duplication at the gene p21 Protein-Activated Kinase 7 (PAK7) is a risk factor for psychosis.

Identifying rare, highly penetrant risk mutations may be an important step in dissecting the molecular etiology of schizophrenia. We conducted a gene-based analysis of large (>100 kb), rare copy-number variants (CNVs) in the Wellcome Trust Case Control Consortium 2 (WTCCC2) schizophrenia sample of 1564 cases and 1748 controls all from Ireland, and further extended the analysis to include an additional 5196 UK controls. We found association with duplications at chr20p12.2 (P = 0.007) and evidence of replication in large independent European schizophrenia (P = 0.052) and UK bipolar disorder case-control cohorts (P = 0.047). A combined analysis of Irish/UK subjects including additional psychosis cases (schizophrenia and bipolar disorder) identified 22 carriers in 11 707 cases and 10 carriers in 21 204 controls [meta-analysis Cochran-Mantel-Haenszel P-value = 2 × 10(-4); odds ratio (OR) = 11.3, 95% CI = 3.7, ∞]. Nineteen of the 22 cases and 8 of the 10 controls carried duplications starting at 9.68 Mb with similar breakpoints across samples. By haplotype analysis and sequencing, we identified a tandem ~149 kb duplication overlapping the gene p21 Protein-Activated Kinase 7 (PAK7, also called PAK5) which was in linkage disequilibrium with local haplotypes (P = 2.5 × 10(-21)), indicative of a single ancestral duplication event. We confirmed the breakpoints in 8/8 carriers tested and found co-segregation of the duplication with illness in two additional family members of one of the affected probands. We demonstrate that PAK7 is developmentally co-expressed with another known psychosis risk gene (DISC1) suggesting a potential molecular mechanism involving aberrant synapse development and plasticity.

Identification of mechanosensitive genes during skeletal development: alteration of genes associated with cytoskeletal rearrangement and cell signalling pathways.

BACKGROUND: Mechanical stimulation is necessary for regulating correct formation of the skeleton. Here we test the hypothesis that mechanical stimulation of the embryonic skeletal system impacts expression levels of genes implicated in developmentally important signalling pathways in a genome wide approach. We use a mutant mouse model with altered mechanical stimulation due to the absence of limb skeletal muscle (Splotch-delayed) where muscle-less embryos show specific defects in skeletal elements including delayed ossification, changes in the size and shape of cartilage rudiments and joint fusion. We used Microarray and RNA sequencing analysis tools to identify differentially expressed genes between muscle-less and control embryonic (TS23) humerus tissue. RESULTS: We found that 680 independent genes were down-regulated and 452 genes up-regulated in humeri from muscle-less Spd embryos compared to littermate controls (at least 2-fold; corrected p-value ≤0.05). We analysed the resulting differentially expressed gene sets using Gene Ontology annotations to identify significant enrichment of genes associated with particular biological processes, showing that removal of mechanical stimuli from muscle contractions affected genes associated with development and differentiation, cytoskeletal architecture and cell signalling. Among cell signalling pathways, the most strongly disturbed was Wnt signalling, with 34 genes including 19 pathway target genes affected. Spatial gene expression analysis showed that both a Wnt ligand encoding gene (Wnt4) and a pathway antagonist (Sfrp2) are up-regulated specifically in the developing joint line, while the expression of a Wnt target gene, Cd44, is no longer detectable in muscle-less embryos. The identification of 84 genes associated with the cytoskeleton that are down-regulated in the absence of muscle indicates a number of candidate genes that are both mechanoresponsive and potentially involved in mechanotransduction, converting a mechanical stimulus into a transcriptional response. CONCLUSIONS: This work identifies key developmental regulatory genes impacted by altered mechanical stimulation, sheds light on the molecular mechanisms that interpret mechanical stimulation during skeletal development and provides valuable resources for further investigation of the mechanistic basis of mechanoregulation. In particular it highlights the Wnt signalling pathway as a potential point of integration of mechanical and molecular signalling and cytoskeletal components as mediators of the response.

Delineating the genetic heterogeneity of ALS using targeted high-throughput sequencing.

BACKGROUND: Over 100 genes have been implicated in the aetiology of amyotrophic lateral sclerosis (ALS). A detailed understanding of their independent and cumulative contributions to disease burden may help guide various clinical and research efforts. METHODS: Using targeted high-throughput sequencing, we characterised the variation of 10 Mendelian and 23 low penetrance/tentative ALS genes within a population-based cohort of 444 Irish ALS cases (50 fALS, 394 sALS) and 311 age-matched and geographically matched controls. RESULTS: Known or potential high-penetrance ALS variants were identified within 17.1% of patients (38% of fALS, 14.5% of sALS). 12.8% carried variants of Mendelian disease genes (C9orf72 8.78%; SETX 2.48%; ALS2 1.58%; FUS 0.45%; TARDBP 0.45%; OPTN 0.23%; VCP 0.23%. ANG, SOD1, VAPB 0%), 4.7% carried variants of low penetrance/tentative ALS genes and 9.7% (30% of fALS, 7.1% of sALS) carried previously described ALS variants (C9orf72 8.78%; FUS 0.45%; TARDBP 0.45%). 1.6% of patients carried multiple known/potential disease variants, including all identified carriers of an established ALS variant (p<0.01); TARDBP:c.859G>A(p.[G287S]) (n=2/2 sALS). Comparison of our results with those from studies of other European populations revealed significant differences in the spectrum of disease variation (p=1.7×10(-4)). CONCLUSIONS: Up to 17% of Irish ALS cases may carry high-penetrance variants within the investigated genes. However, the precise nature of genetic susceptibility differs significantly from that reported within other European populations. Certain variants may not cause disease in isolation and concomitant analysis of disease genes may prove highly important.

The origins of the Guinness stout yeast.

Beer is made via the fermentation of an aqueous extract predominantly composed of malted barley flavoured with hops. The transforming microorganism is typically a single strain of Saccharomyces cerevisiae, and for the majority of major beer brands the yeast strain is a unique component. The present yeast used to make Guinness stout brewed in Dublin, Ireland, can be traced back to 1903, but its origins are unknown. To that end, we used Illumina and Nanopore sequencing to generate whole-genome sequencing data for a total of 22 S. cerevisiae yeast strains: 16 from the Guinness collection and 6 other historical Irish brewing. The origins of the Guinness yeast were determined with a SNP-based analysis, demonstrating that the Guinness strains occupy a distinct group separate from other historical Irish brewing yeasts. Assessment of chromosome number, copy number variation and phenotypic evaluation of key brewing attributes established Guinness yeast-specific SNPs but no specific chromosomal amplifications. Our analysis also demonstrated the effects of yeast storage on phylogeny. Altogether, our results suggest that the Guinness yeast used today is related to the first deposited Guinness yeast; the 1903 Watling Laboratory Guinness yeast.

Global endometrial transcriptomic profiling: transient immune activation precedes tissue proliferation and repair in healthy beef cows.

BACKGROUND: All cows experience bacterial contamination and tissue injury in the uterus postpartum, instigating a local inflammatory immune response. However mechanisms that control inflammation and achieve a physiologically functioning endometrium, while avoiding disease in the postpartum cow are not succinctly defined. This study aimed to identify novel candidate genes indicative of inflammation resolution during involution in healthy beef cows. Previous histological analysis of the endometrium revealed elevated inflammation 15 days postpartum (DPP) which was significantly decreased by 30 DPP. The current study generated a genome-wide transcriptomic profile of endometrial biopsies from these cows at both time points using mRNA-Seq. The pathway analysis tool GoSeq identified KEGG pathways enriched by significantly differentially expressed genes at both time points. Novel candidate genes associated with inflammatory resolution were subsequently validated in additional postpartum animals using quantitative real-time PCR (qRT-PCR). RESULTS: mRNA-Seq revealed 1,107 significantly differentially expressed genes, 73 of which were increased 15 DPP and 1,034 were increased 30 DPP. Early postpartum, enriched immune pathways (adjusted P < 0.1) included the T cell receptor signalling pathway, graft-versus-host disease and cytokine-cytokine receptor interaction pathways. However 30 DPP, where the majority of genes were differentially expressed, the enrichment (adjusted P < 0.1) of tissue repair and proliferative activity pathways was observed. Nineteen candidate genes selected from mRNA-Seq results, were independently assessed by qRT-PCR in additional postpartum cows (5 animals) at both time points. SAA1/2, GATA2, IGF1, SHC2, and SERPINA14 genes were significantly elevated 30 DPP and are functionally associated with tissue repair and the restoration of uterine homeostasis postpartum. CONCLUSIONS: The results of this study reveal an early activation of the immune response which undergoes a temporal functional change toward tissue proliferation and regeneration during endometrial involution in healthy postpartum cows. These molecular changes mirror the activation and resolution of endometrial inflammation during involution previously classified by the degree of neutrophil infiltration. SAA1/2, GATA2, IGF1, SHC2, and SERPINA14 genes may become potential markers for resolution of endometrial inflammation in the postpartum cow.

Translation of 5' leaders is pervasive in genes resistant to eIF2 repression.

Eukaryotic cells rapidly reduce protein synthesis in response to various stress conditions. This can be achieved by the phosphorylation-mediated inactivation of a key translation initiation factor, eukaryotic initiation factor 2 (eIF2). However, the persistent translation of certain mRNAs is required for deployment of an adequate stress response. We carried out ribosome profiling of cultured human cells under conditions of severe stress induced with sodium arsenite. Although this led to a 5.4-fold general translational repression, the protein coding open reading frames (ORFs) of certain individual mRNAs exhibited resistance to the inhibition. Nearly all resistant transcripts possess at least one efficiently translated upstream open reading frame (uORF) that represses translation of the main coding ORF under normal conditions. Site-specific mutagenesis of two identified stress resistant mRNAs (PPP1R15B and IFRD1) demonstrated that a single uORF is sufficient for eIF2-mediated translation control in both cases. Phylogenetic analysis suggests that at least two regulatory uORFs (namely, in SLC35A4 and MIEF1) encode functional protein products.

Transcriptome analysis of the mouse E14.5 (TS23) developing humerus and differential expression in muscle-less mutant embryos lacking mechanical stimulation.

Mechanical stimulation is important for the correct formation of the skeleton. Splotch-delayed mutant embryos (Pax3 (Spd/Spd) ) that develop with no limb muscle and therefore no limb movement experience an altered mechanical environment resulting in specific defects in ossification and joint formation, particularly in the forelimb. To test the hypothesis that mechanical stimuli influence the regulation of genes important in skeletal development we generated a transcriptome profile of the developing humerus at Theiler stage 23 (TS23), and then identified differentially expressed genes in muscle-less mutant embryos compared to control littermates. Here we describe the experimental methods and analysis of the resulting data, publically available in the ArrayExpress database under E-MTAB-1745 (Transcriptome of control humerus), E-MTAB-1744 (Microarray; differential expression) and E-MTAB-1746 (RNA-sequencing; differential expression). Our data provide a resource for exploring the transcriptome that underlies skeletal development at TS23 in the mouse humerus. The interpretation and description of this data can be found in a recent publication in BMC Genomics [1]. This is a resource for exploring the molecular mechanisms that are involved in skeletal development and mechanotransduction.

Next-generation sequencing of the mitochondrial genome and association with IgA nephropathy in a renal transplant population.

Kidneys are highly aerobic organs that are critically dependent on the normal functioning of mitochondria. Genetic variations disrupting mitochondrial function are associated with multifactorial disorders including kidney disease. This study sequenced the entire mitochondrial genome in a renal transplant cohort of 64 individuals, using next-generation sequencing, to evaluate the association of genetic variants with IgA nephropathy and end-stage renal disease (ESRD, n = 100).

Development of strategies for SNP detection in RNA-seq data: application to lymphoblastoid cell lines and evaluation using 1000 Genomes data.

Next-generation RNA sequencing (RNA-seq) maps and analyzes transcriptomes and generates data on sequence variation in expressed genes. There are few reported studies on analysis strategies to maximize the yield of quality RNA-seq SNP data. We evaluated the performance of different SNP-calling methods following alignment to both genome and transcriptome by applying them to RNA-seq data from a HapMap lymphoblastoid cell line sample and comparing results with sequence variation data from 1000 Genomes. We determined that the best method to achieve high specificity and sensitivity, and greatest number of SNP calls, is to remove duplicate sequence reads after alignment to the genome and to call SNPs using SAMtools. The accuracy of SNP calls is dependent on sequence coverage available. In terms of specificity, 89% of RNA-seq SNPs calls were true variants where coverage is >10X. In terms of sensitivity, at >10X coverage 92% of all expected SNPs in expressed exons could be detected. Overall, the results indicate that RNA-seq SNP data are a very useful by-product of sequence-based transcriptome analysis. If RNA-seq is applied to disease tissue samples and assuming that genes carrying mutations relevant to disease biology are being expressed, a very high proportion of these mutations can be detected.

Multiplex target enrichment using DNA indexing for ultra-high throughput SNP detection.

Screening large numbers of target regions in multiple DNA samples for sequence variation is an important application of next-generation sequencing but an efficient method to enrich the samples in parallel has yet to be reported. We describe an advanced method that combines DNA samples using indexes or barcodes prior to target enrichment to facilitate this type of experiment. Sequencing libraries for multiple individual DNA samples, each incorporating a unique 6-bp index, are combined in equal quantities, enriched using a single in-solution target enrichment assay and sequenced in a single reaction. Sequence reads are parsed based on the index, allowing sequence analysis of individual samples. We show that the use of indexed samples does not impact on the efficiency of the enrichment reaction. For three- and nine-indexed HapMap DNA samples, the method was found to be highly accurate for SNP identification. Even with sequence coverage as low as 8x, 99% of sequence SNP calls were concordant with known genotypes. Within a single experiment, this method can sequence the exonic regions of hundreds of genes in tens of samples for sequence and structural variation using as little as 1 µg of input DNA per sample.

Gene-ontology enrichment analysis in two independent family-based samples highlights biologically plausible processes for autism spectrum disorders.

Recent genome-wide association studies (GWAS) have implicated a range of genes from discrete biological pathways in the aetiology of autism. However, despite the strong influence of genetic factors, association studies have yet to identify statistically robust, replicated major effect genes or SNPs. We apply the principle of the SNP ratio test methodology described by O'Dushlaine et al to over 2100 families from the Autism Genome Project (AGP). Using a two-stage design we examine association enrichment in 5955 unique gene-ontology classifications across four groupings based on two phenotypic and two ancestral classifications. Based on estimates from simulation we identify excess of association enrichment across all analyses. We observe enrichment in association for sets of genes involved in diverse biological processes, including pyruvate metabolism, transcription factor activation, cell-signalling and cell-cycle regulation. Both genes and processes that show enrichment have previously been examined in autistic disorders and offer biologically plausibility to these findings.