The contribution of de novo coding mutations to autism spectrum disorder.

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

Genomic analysis reveals a potential role for cell cycle perturbation in HCV-mediated apoptosis of cultured hepatocytes.

The mechanisms of liver injury associated with chronic HCV infection, as well as the individual roles of both viral and host factors, are not clearly defined. However, it is becoming increasingly clear that direct cytopathic effects, in addition to immune-mediated processes, play an important role in liver injury. Gene expression profiling during multiple time-points of acute HCV infection of cultured Huh-7.5 cells was performed to gain insight into the cellular mechanism of HCV-associated cytopathic effect. Maximal induction of cell-death-related genes and appearance of activated caspase-3 in HCV-infected cells coincided with peak viral replication, suggesting a link between viral load and apoptosis. Gene ontology analysis revealed that many of the cell-death genes function to induce apoptosis in response to cell cycle arrest. Labeling of dividing cells in culture followed by flow cytometry also demonstrated the presence of significantly fewer cells in S-phase in HCV-infected relative to mock cultures, suggesting HCV infection is associated with delayed cell cycle progression. Regulation of numerous genes involved in anti-oxidative stress response and TGF-beta1 signaling suggest these as possible causes of delayed cell cycle progression. Significantly, a subset of cell-death genes regulated during in vitro HCV infection was similarly regulated specifically in liver tissue from a cohort of HCV-infected liver transplant patients with rapidly progressive fibrosis. Collectively, these data suggest that HCV mediates direct cytopathic effects through deregulation of the cell cycle and that this process may contribute to liver disease progression. This in vitro system could be utilized to further define the cellular mechanism of this perturbation.

Impaired cholesterol biosynthesis in a neuronal cell line persistently infected with measles virus.

Measles virus remains a substantial cause of morbidity and mortality, producing acute infection with a potential for development of viral persistence. To study the events underlying acute and persistent measles virus infection, we performed a global transcriptional analysis on murine neuroblastoma cells that were acutely or persistently infected with measles virus. In general, we found that acute infection induced significantly more gene expression changes than did persistent infection. A functional enrichment analysis to identify which host pathways were perturbed during each of these infections identified several pathways related to cholesterol biosynthesis, including cholesterol metabolic processes, hydroxymethylglutaryl-coenzyme A (CoA) reductase activity, and acetyl-CoA C-acetyltransferase activity. We also found that measles virus colocalized to lipid rafts in both acute and persistent infection models and that the majority of genes associated with cholesterol synthesis were downregulated in persistent infection relative to acute infection, suggesting a possible link with the defective viral budding in persistent infection. Further, we found that pharmacological inhibition of cholesterol synthesis resulted in the inhibition of viral budding during acute infection. In summary, persistent measles viral infection was associated with decreased cholesterol synthesis, a lower abundance of cholesterol and lipid rafts in the cell membrane, and inhibition of giant-cell formation and release of viral progeny.

MAID : an effect size based model for microarray data integration across laboratories and platforms.

BACKGROUND: Gene expression profiling has the potential to unravel molecular mechanisms behind gene regulation and identify gene targets for therapeutic interventions. As microarray technology matures, the number of microarray studies has increased, resulting in many different datasets available for any given disease. The increase in sensitivity and reliability of measurements of gene expression changes can be improved through a systematic integration of different microarray datasets that address the same or similar biological questions. RESULTS: Traditional effect size models can not be used to integrate array data that directly compare treatment to control samples expressed as log ratios of gene expressions. Here we extend the traditional effect size model to integrate as many array datasets as possible. The extended effect size model (MAID) can integrate any array datatype generated with either single or two channel arrays using either direct or indirect designs across different laboratories and platforms. The model uses two standardized indices, the standard effect size score for experiments with two groups of data, and a new standardized index that measures the difference in gene expression between treatment and control groups for one sample data with replicate arrays. The statistical significance of treatment effect across studies for each gene is determined by appropriate permutation methods depending on the type of data integrated. We apply our method to three different expression datasets from two different laboratories generated using three different array platforms and two different experimental designs. Our results indicate that the proposed integration model produces an increase in statistical power for identifying differentially expressed genes when integrating data across experiments and when compared to other integration models. We also show that genes found to be significant using our data integration method are of direct biological relevance to the three experiments integrated. CONCLUSION: High-throughput genomics data provide a rich and complex source of information that could play a key role in deciphering intricate molecular networks behind disease. Here we propose an extension of the traditional effect size model to allow the integration of as many array experiments as possible with the aim of increasing the statistical power for identifying differentially expressed genes.

High-density rhesus macaque oligonucleotide microarray design using early-stage rhesus genome sequence information and human genome annotations.

BACKGROUND: Until recently, few genomic reagents specific for non-human primate research have been available. To address this need, we have constructed a macaque-specific high-density oligonucleotide microarray by using highly fragmented low-pass sequence contigs from the rhesus genome project together with the detailed sequence and exon structure of the human genome. Using this method, we designed oligonucleotide probes to over 17,000 distinct rhesus/human gene orthologs and increased by four-fold the number of available genes relative to our first-generation expressed sequence tag (EST)-derived array. RESULTS: We constructed a database containing 248,000 exon sequences from 23,000 human RefSeq genes and compared each human exon with its best matching sequence in the January 2005 version of the rhesus genome project list of 486,000 DNA contigs. Best matching rhesus exon sequences for each of the 23,000 human genes were then concatenated in the proper order and orientation to produce a rhesus "virtual transcriptome." Microarray probes were designed, one per gene, to the region closest to the 3' untranslated region (UTR) of each rhesus virtual transcript. Each probe was compared to a composite rhesus/human transcript database to test for cross-hybridization potential yielding a final probe set representing 18,296 rhesus/human gene orthologs, including transcript variants, and over 17,000 distinct genes. We hybridized mRNA from rhesus brain and spleen to both the EST- and genome-derived microarrays. Besides four-fold greater gene coverage, the genome-derived array also showed greater mean signal intensities for genes present on both arrays. Genome-derived probes showed 99.4% identity when compared to 4,767 rhesus GenBank sequence tag site (STS) sequences indicating that early stage low-pass versions of complex genomes are of sufficient quality to yield valuable functional genomic information when combined with finished genome information from a closely related species. CONCLUSION: The number of different genes represented on microarrays for unfinished genomes can be greatly increased by matching known gene transcript annotations from a closely related species with sequence data from the unfinished genome. Signal intensity on both EST- and genome-derived arrays was highly correlated with probe distance from the 3' UTR, information often missing from ESTs yet present in early-stage genome projects.

Distinct cellular responses differentiating alcohol- and hepatitis C virus-induced liver cirrhosis.

BACKGROUND: Little is known at the molecular level concerning the differences and/or similarities between alcohol and hepatitis C virus induced liver disease. Global transcriptional profiling using oligonucleotide microarrays was therefore performed on liver biopsies from patients with cirrhosis caused by either chronic alcohol consumption or chronic hepatitis C virus (HCV). RESULTS: Global gene expression patterns varied significantly depending upon etiology of liver disease, with a greater number of differentially regulated genes seen in HCV-infected patients. Many of the gene expression changes specifically observed in HCV-infected cirrhotic livers were expectedly associated with activation of the innate antiviral immune response. We also compared severity (CTP class) of cirrhosis for each etiology and identified gene expression patterns that differentiated ethanol-induced cirrhosis by class. CTP class A ethanol-cirrhotic livers showed unique expression patterns for genes implicated in the inflammatory response, including those related to macrophage activation and migration, as well as lipid metabolism and oxidative stress genes. CONCLUSION: Stages of liver cirrhosis could be differentiated based on gene expression patterns in ethanol-induced, but not HCV-induced, disease. In addition to genes specifically regulating the innate antiviral immune response, mechanisms responsible for differentiating chronic liver damage due to HCV or ethanol may be closely related to regulation of lipid metabolism and to effects of macrophage activation on deposition of extracellular matrix components.

A 52-kb deletion in the SOST-MEOX1 intergenic region on 17q12-q21 is associated with van Buchem disease in the Dutch population.

Van Buchem disease is an autosomal recessive sclerosing bone dysplasia characterized by skeletal hyperostosis, overgrowth of the mandible, and a liability to entrapment of the seventh and eighth cranial nerves. The genetic determinant maps to chromosome 17q12-q21. We refined the critical interval to the < 1-Mb region between D17S2250 and D17S2253 in 15 affected individuals, all of whom shared a common disease haplotype. Furthermore, we report here the identification of a 52-kb deletion located within the interval and encompassing D17S1789 that is 100% concordant with the disorder. Although the deletion itself does not appear to disrupt the coding region of any known or novel gene(s), the closest flanking genes are MEOX1 on the proximal side, and SOST on the distal side of the deletion. MEOX1 is known to be important for the development of the axial skeleton, whereas the SOST gene is the determinant of sclerosteosis, a disorder that shares many features with van Buchem disease, thus raising the possibility that van Buchem disease results from dysregulation of the expression of one or both of these genes.

Proteomic profiling of human liver biopsies: hepatitis C virus-induced fibrosis and mitochondrial dysfunction.

UNLABELLED: Liver biopsies from hepatitis C virus (HCV)-infected patients offer the unique opportunity to study human liver biology and disease in vivo. However, the low protein yields associated with these small samples present a significant challenge for proteomic analysis. In this study we describe the application of an ultrasensitive proteomics platform for performing robust quantitative proteomic studies on microgram amounts of HCV-infected human liver tissue from 15 patients at different stages of fibrosis. A high-quality liver protein database containing 5,920 unique protein identifications supported high throughput quantitative studies using (16)O/(18)O stable isotope labeling in combination with the accurate mass and time (AMT) tag approach. A total of 1,641 liver biopsy proteins were quantified, and analysis of variance (ANOVA) identified 210 proteins exhibiting statistically significant differences associated with fibrosis stage. Hierarchical clustering showed that biopsies representative of later fibrosis stages (for example, Batts-Ludwig stages 3-4) exhibited a distinct protein expression profile, indicating an apparent down-regulation of many proteins when compared with samples from earlier fibrosis stages (for example, Batts-Ludwig stages 0-2). Functional analysis of these signature proteins suggests that impairment of key mitochondrial processes including fatty acid oxidation and oxidative phosphorylation, and response to oxidative stress and reactive oxygen species occurs during advanced stage 3 to 4 fibrosis. CONCLUSION: The results reported here represent a significant advancement in clinical proteomics providing to our knowledge, the first demonstration of global proteomic alterations accompanying liver disease progression in patients chronically infected with HCV. Our findings contribute to a generally emerging theme associating oxidative stress and hepatic mitochondrial dysfunction with HCV pathogenesis.

Functional gene analysis of individual response to challenge of SIVmac239 in M. mulatta PBMC culture.

It has previously been shown in macaques that individual animals exhibit varying responses to challenge with the same strain of SIV. We attempted to elucidate these differences using functional genomics and correlate them to biological response. Unfractionated PBMC from three rhesus macaques were isolated, activated, and infected with SIVmac239. Interestingly, one of the three animals used for these experiments exhibited a completely unique response to infection relative to the other two. After repeated attempts to infect the PBMC from this animal, little or no infectivity was seen across the time points considered, and corresponding to this apparent lack of infection, few genes were seen to be differentially expressed when compared to mock-infected cells. For the remaining two animals, gene expression analysis showed that while they exhibited responses for the same groups of pathways, these responses included differences specific to the individual animal at the gene level. In instances where the patterns of differential gene expression differed between these animals, the genes being differentially expressed were associated with the same categories of biological process, mainly immune response and cell signaling. At the pathway level, these animals again exhibited similar responses that could be predicted based on the experimental conditions. Even in these expected results, the degree of response and the specific genes being regulated differed greatly from animal to animal. The differences in gene expression on an individual level have the potential to be used as markers in identification of animals suitable for lentiviral infection experiments. Our results highlight the importance of individual variation in response to viral challenge.

Polymorphisms in the sclerosteosis/van Buchem disease gene (SOST) region are associated with bone-mineral density in elderly whites.

Osteoporosis has a strong genetic component, but the genes involved are poorly defined. We studied whether the sclerosteosis/van Buchem disease gene (SOST) is an osteoporosis-risk gene by examining its association with bone-mineral density (BMD). Mutations in SOST result in sclerosteosis, and alterations in the SOST gene expression may be causal in the closely related van Buchem disease. We used a set of eight polymorphisms from the SOST gene region to genotype 1,939 elderly men and women from a large population-based prospective-cohort study of Dutch whites. A 3-bp insertion (f=0.38) in the presumed SOST promoter region (SRP3) was associated with decreased BMD in women at the femoral neck (FN) (P=.05) and lumbar spine (LS) (P=.01), with evidence of an allele-dose effect in the oldest age group (P=.006). Similarly, a G variant (f=0.40) in the van Buchem deletion region (SRP9) was associated with increased BMD in men at the FN (P=.007) and LS (P=.02). In both cases, differences between extreme genotypes reached 0.2 SD. We observed no genotype effects on fracture risk, for the 234 osteoporotic fractures validated during 8.2 years of follow-up and for the 146 vertebral prevalent fractures analyzed. We did not find association between any of several frequent haplotypes across the SOST gene region and BMD. We did find evidence of additive effects of SRP3 with the COLIA1 Sp1 polymorphism but not with haplotypes of 3' polymorphisms in the vitamin-D receptor gene. The SOST-COLIA1 additive effect increased with age and reached 0.5 SD difference in BMD at LS in the oldest age group (P=.02). The molecular mechanism whereby these moderate SOST genotype effects are mediated remains to be elucidated, but it is likely to involve differences in regulation of SOST gene expression.

A novel mutation in CD83 results in the development of a unique population of CD4+ T cells.

Using a mouse mutagenesis screen, we have identified CD83 as being critical for the development of CD4(+) T cells and for their function postactivation. CD11c(+) dendritic cells develop and function normally in mice with a mutated CD83 gene but CD4(+) T cell development is substantially reduced. Additionally, we now show that those CD4(+) cells that develop in a CD83 mutant animal fail to respond normally following allogeneic stimulation. This is at least in part due to an altered cytokine expression pattern characterized by an increased production of IL-4 and IL-10 and diminished IL-2 production. Thus, in addition to its role in selection of CD4(+) T cells, absence of CD83 results in the generation of cells with an altered activation and cytokine profile.