Application of mixture models to detect differentially expressed genes

An important and common problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. As this problem concerns the selection of significant genes from a large pool of candidate genes, it needs to be carried out within the framework of multiple hypothesis testing. In this paper, we focus on the use of mixture models to handle the multiplicity issue. With this approach, a measure of the local FDR (false discovery rate) is provided for each gene. An attractive feature of the mixture model approach is that it provides a framework for the estimation of the prior probability that a gene is not differentially expressed, and this probability can subsequently be used in forming a decision rule. The rule can also be formed to take the false negative rate into account. We apply this approach to a well-known publicly available data set on breast cancer, and discuss our findings with reference to other approaches.

Using mixture models to detect differentially expressed genes

An important and common problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. As this problem concerns the selection of significant genes from a large pool of candidate genes, it needs to be carried out within the framework of multiple hypothesis testing. In this paper, we focus on the use of mixture models to handle the multiplicity issue. With this approach, a measure of the local false discovery rate is provided for each gene, and it can be implemented so that the implied global false discovery rate is bounded as with the Benjamini-Hochberg methodology based on tail areas. The latter procedure is too conservative, unless it is modified according to the prior probability that a gene is not differentially expressed. An attractive feature of the mixture model approach is that it provides a framework for the estimation of this probability and its subsequent use in forming a decision rule. The rule can also be formed to take the false negative rate into account.

Mixture models for detecting differentially expressed genes in microarrays

An important and common problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. As this problem concerns the selection of significant genes from a large pool of candidate genes, it needs to be carried out within the framework of multiple hypothesis testing. In this paper, we focus on the use of mixture models to handle the multiplicity issue. With this approach, a measure of the local FDR (false discovery rate) is provided for each gene. An attractive feature of the mixture model approach is that it provides a framework for the estimation of the prior probability that a gene is not differentially expressed, and this probability can subsequently be used in forming a decision rule. The rule can also be formed to take the false negative rate into account. We apply this approach to a well-known publicly available data set on breast cancer, and discuss our findings with reference to other approaches.

The dyslexia candidate locus on 2p12 is associated with general cognitive ability and white matter structure

Independent studies have shown that candidate genes for dyslexia and specific language impairment (SLI) impact upon reading/language-specific traits in the general population. To further explore the effect of disorder-associated genes on cognitive functions, we investigated whether they play a role in broader cognitive traits. We tested a panel of dyslexia and SLI genetic risk factors for association with two measures of general cognitive abilities, or IQ, (verbal and non-verbal) in the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (N>5,000). Only the MRPL19/C2ORF3 locus showed statistically significant association (minimum P = 0.00009) which was further supported by independent replications following analysis in four other cohorts. In addition, a fifth independent sample showed association between the MRPL19/C2ORF3 locus and white matter structure in the posterior part of the corpus callosum and cingulum, connecting large parts of the cortex in the parietal, occipital and temporal lobes. These findings suggest that this locus, originally identified as being associated with dyslexia, is likely to harbour genetic variants associated with general cognitive abilities by influencing white matter structure in localised neuronal regions.

Clonal composition of human adamantinomatous craniopharyngiomas and somatic mutation analyses of the patched (PTCH), Gsα and Gi2α genes

Craniopharyngioma is the most common childhood tumor and thought to arise from embryonic remnants of Rathke's pouch. The paucity of published data on the molecular basis of these tumors prompted us to examine 22 adamantinomatous craniopharyngiomas looking for genetic abnormalities. Using the X-linked polymorphic androgen receptor gene as a tool for X-chromosome inactivating analysis, we found that a subset of craniopharyngiomas are monoclonal and therefore are probably due to acquired somatic genetic defects. Thus, we investigated these tumours for mutations within three candidate genes, Gsα, Gi2α and patched (PTCH). Using single stranded conformational polymorphism (SSCP), denaturing gradient gel electrophoresis and direct sequencing, the presence of somatic mutations in these genes could not be demonstrated in any tumor. Our data indicate that a subset of craniopharyngiomas are monoclonal and the mutations in the PTCH, Gsα, and Gi2α contribute little if any to cranipharyngioma development.

Microarray analysis of expression of cell death-associated genes in spinal cord cells with cyclic tensile strain

Previous studies have described alterations in gene expression following spinal cord injury, but this response to mechanical stimuli is difficult to investigate in vivo. Therefore, we have investigated the effect of cyclic tensile strain on cultured spinal cord cells from E15 Sprague-Dawley rats. Microarray analysis of gene expression and categorization of identified genes were performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) systems. The application of cyclic tensile strain reduced the viability of cultured spinal cord cells significantly in a dose- and time-dependent manner. GO analysis identified candidate genes related to apoptosis (44) and to response to stimulus (17). KEGG analysis identified changes in the expression levels of 12 genes of the mitogen-activated protein kinase (MAPK) signaling pathway, which were confirmed to be upregulated and validated by RT-PCR analysis. Spinal cord cells undergo cell death in response to cyclic tensile strain, which were dose- and time-dependent, with upregulation of various genes, in particular of the MAPK pathway.

Integration of global SNP-based mapping and expression arrays reveals key regions, mechanisms, and genes important in the pathogenesis of multiple myeloma.

Multiple myeloma is characterized by genomic alterations frequently involving gains and losses of chromosomes. Single nucleotide polymorphism (SNP)-based mapping arrays allow the identification of copy number changes at the sub-megabase level and the identification of loss of heterozygosity (LOH) due to monosomy and uniparental disomy (UPD). We have found that SNP-based mapping array data and fluorescence in situ hybridization (FISH) copy number data correlated well, making the technique robust as a tool to investigate myeloma genomics. The most frequently identified alterations are located at 1p, 1q, 6q, 8p, 13, and 16q. LOH is found in these large regions and also in smaller regions throughout the genome with a median size of 1 Mb. We have identified that UPD is prevalent in myeloma and occurs through a number of mechanisms including mitotic nondisjunction and mitotic recombination. For the first time in myeloma, integration of mapping and expression data has allowed us to reduce the complexity of standard gene expression data and identify candidate genes important in both the transition from normal to monoclonal gammopathy of unknown significance (MGUS) to myeloma and in different subgroups within myeloma. We have documented these genes, providing a focus for further studies to identify and characterize those that are key in the pathogenesis of myeloma.

Identification of superior reference genes for data normalisation of expression studies via quantitative PCR in hybrid roses (Rosa hybrida)

Background: Gene expression studies are a prerequisite for understanding the biological function of genes. Because of its high sensitivity and easy use, quantitative PCR (qPCR) has become the gold standard for gene expression quantification. To normalise qPCR measurements between samples, the most prominent technique is the use of stably expressed endogenous control genes, the so called reference genes. However, recent studies show there is no universal reference gene for all biological questions. Roses are important ornamental plants for which there has been no evaluation of useful reference genes for gene expression studies. Results: We used three different algorithms (BestKeeper, geNorm and NormFinder) to validate the expression stability of nine candidate reference genes in different rose tissues from three different genotypes of Rosa hybrida and in leaves treated with various stress factors. The candidate genes comprised the classical "housekeeping genes" (Actin, EF-1α, GAPDH, Tubulin and Ubiquitin), and genes showing stable expression in studies in Arabidopsis (PP2A, SAND, TIP and UBC). The programs identified no single gene that showed stable expression under all of the conditions tested, and the individual rankings of the genes differed between the algorithms. Nevertheless the new candidate genes, specifically, PP2A and UBC, were ranked higher as compared to the other traditional reference genes. In general, Tubulin showed the most variable expression and should be avoided as a reference gene. Conclusions: Reference genes evaluated as suitable in experiments with Arabidopsis thaliana were stably expressed in roses under various experimental conditions. In most cases, these genes outperformed conventional reference genes, such as EF1-α and Tubulin. We identified PP2A, SAND and UBC as suitable reference genes, which in different combinations may be used for normalisation in expression analyses via qPCR for different rose tissues and stress treatments. However, the vast genetic variation found within the genus Rosa, including differences in ploidy levels, might also influence expression stability of reference genes, so that future research should also consider different genotypes and ploidy levels.

Prospection for genes expressed during wheat - Magnaporthe oryzae interaction.

2013

Expression of rice genes homologous of Arabidopsis genes previously related to drought tolerance.

The identification and validation of candidate genes related to traits of interest is a time consuming and expensive process and the homology among genes from different species can facilitate the identification of genes of the target species from the genomic information of a model species. This study aimed to quantify the expression of homologous rice genes previously related to drought tolerance in Arabidopsis. Five genes (CPK6, PLDa, GluR2, CesA8, and EIN2) were identified in rice by the homology of the amino acid sequence between rice and Arabidopsis. The genotypes Douradão (drought tolerant) and Primavera (drought susceptible) were subjected to a water deficit experiment, and subsequently evaluated for gene expression by qPCR for the five homologous and Lsi1 genes. The qPCR analysis clearly showed that the five homologous genes were expressed in rice, which is an indication that these genes could preserve their function in rice as a response to drought. In Douradão, of the five homologous genes, all but OsGluR2 displayed an increase in the average expression in drought treatment when compared to the control, while in Primavera, the average expression of the five genes did not differ between the control and drought treatment. In Douradão, the OsPLDa1, which showed the higher expression level in drought in relation to the control (10.82), significantly increased the gene expression in the leaf and root tissues as a response to drought, in both vegetative and reproductive stages, whereas in Primavera, this gene was suppressed in both tissues and stages under drought. Therefore, the OsPLDa1 gene was the most important in relation to drought response and is an interesting candidate for further studies in developing rice cultivars that are more tolerant to this stress.

The genetic basis of human height : the role of estrogen

Height is a complex physical trait that displays strong heritability. Adult height is related to length of the long bones, which is determined by growth at the epiphyseal growth plate. Longitudinal bone growth occurs via the process of endochondral ossification, where bone forms over the differentiating cartilage template at the growth plate. Estrogen plays a major role in regulating longitudinal bone growth and is responsible for inducing the pubertal growth spurt and fusion of the epiphyseal growth plate. However, the mechanism by which estrogen promotes epiphyseal fusion is poorly understood. It has been hypothesised that estrogen functions to regulate growth plate fusion by stimulating chondrocyte apoptosis, angiogenesis and bone cell invasion in the growth plate. Another theory has suggested that estrogen exposure exhausts the proliferative capacity of growth plate chondrocytes, which accelerates the process of chondrocyte senescence, leading to growth plate fusion. The overall objective of this study was to gain a greater understanding of the molecular mechanisms behind estrogen-mediated growth and height attainment by examining gene regulation in chondrocytes and the role of some of these genes in normal height inheritance. With the heritability of height so well established, the initial hypothesis was that genetic variation in candidate genes associated with longitudinal bone growth would be involved in normal adult height variation. The height-related genes FGFR3, CBFA1, ER and CBFA1 were screened for novel polymorphisms using denaturing HPLC and RFLP analysis. In total, 24 polymorphisms were identified. Two SNPs in ER (rs3757323 C>T and rs1801132 G>C) were strongly associated with adult male height and displayed an 8 cm and 9 cm height difference between homozygous genotypes, respectively. The TC haplotype of these SNPs was associated with a 6 cm decrease in height and remarkably, no homozygous carriers of the TC haplotype were identified in tall subjects. No significant associations with height were found for polymorphisms in the FGFR3, CBFA1 or VDR genes. In the epiphyseal growth plate, chondrocyte proliferation, matrix synthesis and chondrocyte hypertrophy are all major contributors to long bone growth. As estrogen plays such a significant role in both growth and final height attainment, another hypothesis of this study was that estrogen exerted its effects in the growth plate by influencing chondrocyte proliferation and mediating the expression of chondrocyte marker genes. The examination of genes regulated by estrogen in chondrocyte-like cells aimed to identify potential regulators of growth plate fusion, which may further elucidate mechanisms involved in the cessation of linear growth. While estrogen did not dramatically alter the proliferation of the SW1353 cell line, gene expression experiments identified several estrogen regulated genes. Sixteen chondrocyte marker genes were examined in response to estrogen concentrations ranging from 10-12 M to 10-8 M over varying time points. Of the genes analysed, IHH, FGFR3, collagen II and collagen X were not readily detectable and PTHrP, GHR, ER, BMP6, SOX9 and TGF1 mRNAs showed no significant response to estrogen treatments. However, the expression of MMP13, CBFA1, BCL-2 and BAX genes were significantly decreased. Interestingly, the majority of estrogen regulated genes in SW1353 cells are expressed in the hypertrophic zone of the growth plate. Estrogen is also known to regulate systemic GH secretion and local GH action. At the molecular level, estrogen functions to inhibit GH action by negatively regulating GH signalling. GH treated SW1353 cells displayed increases in MMP9 mRNA expression (4.4-fold) and MMP13 mRNA expression (64-fold) in SW1353 cells. Increases were also detected in their respective proteins. Treatment with AG490, an established JAK2 inhibitor, blocked the GH mediated stimulation of both MMP9 and MMP13 mRNA expression. The application of estrogen and GH to SW1353 cells attenuated GH-stimulated MMP13 levels, but did not affect MMP9 levels. Investigation of GH signalling revealed that SW1353 cells have high levels of activated JAK2 and exposure to GH, estrogen, AG490 and other signalling inhibitors did not affect JAK2 phosphorylation. Interestingly, AG490 treatment dramatically decreased ERK2 signalling, although GH did stimulate ERK2 phosphorylation above control levels. AG490 also decreased CBFA1 expression, a transcription factor known to activate MMP9 and MMP13. Finally, GH and estrogen treatment increased expression of SOCS3 mRNA, suggesting that SOCS3 may regulate JAK/STAT signalling in SW1353 cells. The modulation of GH-mediated MMP expression by estrogen in SW1353 cells represents a potentially novel mechanism by which estrogen may regulate longitudinal bone growth. However, further investigation is required in order to elucidate the precise mechanisms behind estrogen and GH regulation of MMP13 expression in SW1353 cells. This study has provided additional evidence that estrogen and the ER gene are major factors in the regulation of growth and the determination of adult height. Newly identified polymorphisms in the ER gene not only contribute to our understanding of the genetic basis of human height, but may also be useful in association studies examining other complex traits. This study also identified several estrogen regulated genes and indicated that estrogen modifies the expression of genes which are primarily expressed in the hypertrophic region of the epiphyseal growth plate. Furthermore, synergistic studies incorporating GH and estrogen have revealed the ability of estrogen to attenuate the effects of GH on MMP13 expression, revealing potential pathways by which estrogen may modulate growth plate fusion, longitudinal bone growth and even arthritis.

Linkage to chromosome 2q32.2-q33.3 in familial serrated neoplasia (Jass syndrome)

Abstract Causative genetic variants have to date been identified for only a small proportion of familial colorectal cancer (CRC). While conditions such as Familial Adenomatous Polyposis and Lynch syndrome have well defined genetic causes, the search for variants underlying the remainder of familial CRC is plagued by genetic heterogeneity. The recent identification of families with a heritable predisposition to malignancies arising through the serrated pathway (familial serrated neoplasia or Jass syndrome) provides an opportunity to study a subset of familial CRC in which heterogeneity may be greatly reduced. A genome-wide linkage screen was performed on a large family displaying a dominantly-inherited predisposition to serrated neoplasia genotyped using the Affymetrix GeneChip Human Mapping 10 K SNP Array. Parametric and nonparametric analyses were performed and resulting regions of interest, as well as previously reported CRC susceptibility loci at 3q22, 7q31 and 9q22, were followed up by finemapping in 10 serrated neoplasia families. Genome-wide linkage analysis revealed regions of interest at 2p25.2-p25.1, 2q24.3-q37.1 and 8p21.2-q12.1. Finemapping linkage and haplotype analyses identified 2q32.2-q33.3 as the region most likely to harbour linkage, with heterogeneity logarithm of the odds (HLOD) 2.09 and nonparametric linkage (NPL) score 2.36 (P = 0.004). Five primary candidate genes (CFLAR, CASP10, CASP8, FZD7 and BMPR2) were sequenced and no segregating variants identified. There was no evidence of linkage to previously reported loci on chromosomes 3, 7 and 9.

Massively parallel sequencing and analysis of expressed sequence tags in a successful invasive plant

Background Invasive species pose a significant threat to global economies, agriculture and biodiversity. Despite progress towards understanding the ecological factors associated with plant invasions, limited genomic resources have made it difficult to elucidate the evolutionary and genetic factors responsible for invasiveness. This study presents the first expressed sequence tag (EST) collection for Senecio madagascariensis, a globally invasive plant species. Methods We used pyrosequencing of one normalized and two subtractive libraries, derived from one native and one invasive population, to generate an EST collection. ESTs were assembled into contigs, annotated by BLAST comparison with the NCBI non-redundant protein database and assigned gene ontology (GO) terms from the Plant GO Slim ontologies. Key Results Assembly of the 221 746 sequence reads resulted in 12 442 contigs. Over 50 % (6183) of 12 442 contigs showed significant homology to proteins in the NCBI database, representing approx. 4800 independent transcripts. The molecular transducer GO term was significantly over-represented in the native (South African) subtractive library compared with the invasive (Australian) library. Based on NCBI BLAST hits and literature searches, 40 % of the molecular transducer genes identified in the South African subtractive library are likely to be involved in response to biotic stimuli, such as fungal, bacterial and viral pathogens. Conclusions This EST collection is the first representation of the S. madagascariensis transcriptome and provides an important resource for the discovery of candidate genes associated with plant invasiveness. The over-representation of molecular transducer genes associated with defence responses in the native subtractive library provides preliminary support for aspects of the enemy release and evolution of increased competitive ability hypotheses in this successful invasive. This study highlights the contribution of next-generation sequencing to better understanding the molecular mechanisms underlying ecological hypotheses that are important in successful plant invasions.

Deletion mapping suggests that the 1p22 melanoma susceptibility gene is a tumor suppressor localized to a 9-Mb interval

Loss of the short arm of chromosome 1 is frequently observed in many tumor types, including melanoma. We recently localized a third melanoma susceptibility locus to chromosome band 1p22. Critical recombinants in linked families localized the gene to a 15-Mb region between D1S430 and D1S2664. To map the locus more finely we have performed studies to assess allelic loss across the region in a panel of melanomas from 1p22-linked families, sporadic melanomas, and melanoma cell lines. Eighty percent of familial melanomas exhibited loss of heterozygosity (LOH) within the region, with a smallest region of overlapping deletions (SRO) of 9 Mb between D1S207 and D1S435. This high frequency of LOH makes it very likely that the susceptibility locus is a tumor suppressor. In sporadic tumors, four SROs were defined. SRO1 and SRO2 map within the critical recombinant and familial tumor region, indicating that one or the other is likely to harbor the susceptibility gene. However, SRO3 may also be significant because it overlaps with the markers with the highest 2-point LOD score (D1S2776), part of the linkage recombinant region, and the critical region defined in mesothelioma. The candidate genes PRKCL2 and GTF2B, within SRO2, and TGFBR3, CDC7, and EVI5, in a broad region encompassing SRO3, were screened in 1p22-linked melanoma kindreds, but no coding mutations were detected. Allelic loss in melanoma cell lines was significantly less frequent than in fresh tumors, indicating that this gene may not be involved late in progression, such as in overriding cellular senescence, necessary for the propagation of melanoma cells in culture.

Possible genetic predisposition to lymphedema after breast cancer

Background: Known risk factors for secondary lymphedema only partially explain who develops lymphedema following cancer, suggesting that inherited genetic susceptibility may influence risk. Moreover, identification of molecular signatures could facilitate lymphedema risk prediction prior to surgery or lead to effective drug therapies for prevention or treatment. Recent advances in the molecular biology underlying development of the lymphatic system and related congenital disorders implicate a number of potential candidate genes to explore in relation to secondary lymphedema. Methods and Results: We undertook a nested case-control study, with participants who had developed lymphedema after surgical intervention within the first 18 months of their breast cancer diagnosis serving as cases (n=22) and those without lymphedema serving as controls (n=98), identified from a prospective, population-based, cohort study in Queensland, Australia. TagSNPs that covered all known genetic variation in the genes SOX18, VEGFC, VEGFD, VEGFR2, VEGFR3, RORC, FOXC2, LYVE1, ADM and PROX1 were selected for genotyping. Multiple SNPs within three receptor genes, VEGFR2, VEGFR3 and RORC, were associated with lymphedema defined by statistical significance (p<0.05) or extreme risk estimates (OR<0.5 or >2.0). Conclusions: These provocative, albeit preliminary, findings regarding possible genetic predisposition to secondary lymphedema following breast cancer treatment warrant further attention for potential replication using larger datasets.