No abstract

While advances in genome sequencing technology make population-scale genomics a possibility, current approaches for analysis of these data rely upon parallelization strategies that have limited scalability, complex implementation and lack reproducibility. Churchill, a balanced regional parallelization strategy, overcomes these challenges, fully automating the multiple steps required to go from raw sequencing reads to variant discovery. Through implementation of novel deterministic parallelization techniques, Churchill allows computationally efficient analysis of a high-depth whole genome sample in less than two hours. The method is highly scalable, enabling full analysis of the 1000 Genomes raw sequence dataset in a week using cloud resources. http://churchill.nchri.org/.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-014-0577-x) contains supplementary material, which is available to authorized users.

Mucopolysaccharidosis (MPS) IIIB is a lysosomal storage disease with severe neurological manifestations due to alpha-N-acetylglucosaminidase (NaGlu) deficiency. The mechanism of neuropathology in MPS IIIB is unclear. This study investigates the role of immune responses in neurological disease of MPS IIIB in mice. By means of gene expression microarrays and real-time quantitative reverse transcriptase-polymerase chain reaction, we demonstrated significant up-regulation of numerous immune-related genes in MPS IIIB mouse brain involving a broad range of immune cells and molecules, including T cells, B cells, microglia/macrophages, complement, major histocompatibility complex class I, immunoglobulin, Toll-like receptors, and molecules essential for antigen presentation. The significantly enlarged spleen and lymph nodes in MPS IIIB mice were due to an increase in splenocytes/lymphocytes, and functional assays indicated that the T cells were activated. An autoimmune component to the disease was further suggested by the presence of putative autoantigen or autoantigens in brain extracts that reacted specifically with serum IgG from MPS IIIB mice. We also demonstrated for the first time that immunosuppression with prednisolone alone can significantly slow the central nervous system disease progression. Our data indicate that immune responses contribute greatly to the neuropathology of MPS IIIB and should be considered as an adjunct treatment in future therapeutic developments for optimal therapeutic effect.

Salmonella enterica serovar Typhimurium (Salmonella) is one of the most significant food-borne pathogens affecting both humans and agriculture. We have determined that Salmonella encodes an uptake and utilization pathway specific for a novel nutrient, fructose-asparagine (F-Asn), which is essential for Salmonella fitness in the inflamed intestine (modeled using germ-free, streptomycin-treated, ex-germ-free with human microbiota, and IL10−/− mice). The locus encoding F-Asn utilization, fra, provides an advantage only if Salmonella can initiate inflammation and use tetrathionate as a terminal electron acceptor for anaerobic respiration (the fra phenotype is lost in Salmonella SPI1− SPI2− or ttrA mutants, respectively). The severe fitness defect of a Salmonella fra mutant suggests that F-Asn is the primary nutrient utilized by Salmonella in the inflamed intestine and that this system provides a valuable target for novel therapies.

BackgroundBicuspid aortic valve (BAV) is the most common type of congenital heart disease with a population prevalence of 1-2%. While BAV is known to be highly heritable, mutations in single genes (such as GATA5 and NOTCH1) have been reported in few human BAV cases. Traditional gene sequencing methods are time and labor intensive, while next-generation high throughput sequencing remains costly for large patient cohorts and requires extensive bioinformatics processing. Here we describe an approach to targeted multi-gene sequencing with combinatorial pooling of samples from BAV patients.MethodsWe studied a previously described cohort of 78 unrelated subjects with echocardiogram-identified BAV. Subjects were identified as having isolated BAV or BAV associated with coarctation of aorta (BAV-CoA). BAV cusp fusion morphology was defined as right-left cusp fusion, right non-coronary cusp fusion, or left non-coronary cusp fusion. Samples were combined into 19 pools using a uniquely overlapping combinatorial design; a given mutation could be attributed to a single individual on the basis of which pools contained the mutation. A custom gene capture of 97 candidate genes was sequenced on the Illumina HiSeq 2000. Multistep bioinformatics processing was performed for base calling, variant identification, and in-silico analysis of putative disease-causing variants.ResultsTargeted capture identified 42 rare, non-synonymous, exonic variants involving 35 of the 97 candidate genes. Among these variants, in-silico analysis classified 33 of these variants as putative disease-causing changes. Sanger sequencing confirmed thirty-one of these variants, found among 16 individuals. There were no significant differences in variant burden among BAV fusion phenotypes or isolated BAV versus BAV-CoA. Pathway analysis suggests a role for the WNT signaling pathway in human BAV.ConclusionWe successfully developed a pooling and targeted capture strategy that enabled rapid and cost effective next generation sequencing of target genes in a large patient cohort. This approach identified a large number of putative disease-causing variants in a cohort of patients with BAV, including variants in 26 genes not previously associated with human BAV. The data suggest that BAV heritability is complex and polygenic. Our pooling approach saved over $39,350 compared to an unpooled, targeted capture sequencing strategy.

Objective Complement-mediated vasculopathy of muscle and skin are clinical features of juvenile dermatomyositis (JDM). We assess gene copy-number variations (CNVs) for complement C4 and its isotypes, C4A and C4B, in genetic risks and pathogenesis of JDM. Methods The study population included 105 JDM patients and 500 healthy European Americans. Gene copy-numbers (GCNs) for total C4, C4A, C4B and HLA-DRB1 genotypes were determined by Southern blots and PCRs. Processed activation product C4d bound to erythrocytes (E-C4d) was measured by flow cytometry. Global gene-expression microarrays were performed in 19 JDM and 7 controls using PAXgene-blood RNA. Differential expression levels for selected genes were validated by qPCR. Results Significantly lower GCNs and differences in distribution of GCN groups for total C4 and C4A were observed between JDM and controls. Lower GCN of C4A in JDM remained among HLA DR3-positive subjects (p=0.015). Homozygous or heterozygous C4A-deficiency was present in 40.0% of JDM compared to 18.2% of controls [odds ratio (OR)=3.00 (1.87–4.79), p=8.2x10−6]. JDM had higher levels of E-C4d than controls (p=0.004). In JDM, C4A-deficient subjects had higher levels of E-C4d (p=0.0003) and higher frequency of elevated levels of multiple serum muscle enzymes at diagnosis (p=0.004). Microarray profiling of blood RNA revealed upregulation of type I Interferon-stimulated genes and lower abundance of transcripts for T-cell and chemokine function genes in JDM, but this was less prominent among C4A-deficient or DR3-positive patients. Conclusions Complement C4A-deficiency appears to be an important factor for the genetic risk and pathogenesis of JDM, particularly in patients with a DR3-positive background.

Regular bouts of physical activity may cause changes in gene expression that accumulate over time and ultimately affect phenotypes, such as body weight, blood lipid profile and tumour development. Furthermore, acute activity may affect gene expression and phenotypes differently depending on whether the individual is regularly inactive or active. One-month-old male Sprague-Dawley rats (n = 72) were equally divided into SED (standard laboratory cage, n = 24), PA (large activity box, n = 24) and EX groups (exercise wheel inside standard cage, n = 24). At 3 months of age, half the animals from each group were killed at rest and the other half following 30 min of physical activity. The RNA was extracted from cardiac tissue, and microarray analysis was performed on 27,000 genes. Select gene results were validated using quantitative PCR. No gene expression differences occurred when comparing all 3-month-old groups at rest. A relatively small percentage of genes (1.9%) were differentially expressed (P < 0.05) following acute swimming activity in all groups, but only 37 unique and identifiable genes reached or exceeded twofold differences in expression. The genes Atf3, Fos, Apold1 and Pxdn were expressed differently among SED, PA and EX following acute activity, with a clear separation of the magnitude in gene expression with SED > PA > EX. Differences in gene expression levels in young physically inactive and active animals following acute activity have different regulatory roles in gene networks that affect health-related phenotypes.

The approximately 25,000 genes in mammalian genomes can be transcribed at different levels. Measurements of gene expression for ten thousands of genes in parallel give the most comprehensive picture of steady-state levels of transcripts and is used in basic and applied research. Microarrays are the most frequently used technology for genome-wide expression profiling; from the various available microarray platforms, Affymetrix GeneChips are most frequently used for expression profiling and over 3,000 scientific publications describe results of this technology. In medical research, expression profiling by microarrays holds great promises for better understanding of diseases, identification of new therapeutic targets, and subclassification of diseases to identify individualized treatment strategies.