Functional divergence in the Arabidopsis LOB-domain gene family

The Arabidopsis LOB-domain (LBD) gene family is composed by 43 members divided in two classes based on amino acid conservation within the LOB-domain. The LOB domain is known to be responsible for DNA binding and protein-protein interactions. There is very little functional information available for most genes in the LBD family and many lbd single mutants do not exhibit conspicuous phenotypes. One plausible explanation for the limited loss-of-function phenotypes observed in this family is that LBD genes exhibit significant functional redundancy. Here we discuss an example of one phylogenetic subgroup of the LBD family, in which genes that are closely related based on phylogeny exhibit distinctly different expression patterns and do not have overlapping functions. We discuss the challenges of using phylogenetic analyses to predict redundancy in gene families.     

Integrated analyses of copy number variations and gene expression in lung adenocarcinoma

Numerous efforts have been made to elucidate the etiology and improve the treatment of lung cancer, but the overall five-year survival rate is still only 15%. Identification of prognostic biomarkers for lung cancer using gene expression microarrays poses a major challenge in that very few overlapping genes have been reported among different studies. To address this issue, we have performed concurrent genome-wide analyses of copy number variation and gene expression to identify genes reproducibly associated with tumorigenesis and survival in non-smoking female lung adenocarcinoma. The genomic landscape of frequent copy number variable regions (CNVRs) in at least 30% of samples was revealed, and their aberration patterns were highly similar to several studies reported previously. Further statistical analysis for genes located in the CNVRs identified 475 genes differentially expressed between tumor and normal tissues (p<10⁻⁵). We demonstrated the reproducibility of these genes in another lung cancer study (p = 0.0034, Fisher''s exact test), and showed the concordance between copy number variations and gene expression changes by elevated Pearson correlation coefficients. Pathway analysis revealed two major dysregulated functions in lung tumorigenesis: survival regulation via AKT signaling and cytoskeleton reorganization. Further validation of these enriched pathways using three independent cohorts demonstrated effective prediction of survival. In conclusion, by integrating gene expression profiles and copy number variations, we identified genes/pathways that may serve as prognostic biomarkers for lung tumorigenesis.     

Mapping of multiple subunits of the neuronal nicotinic acetylcholine receptor to chromosome 15 in man and chromosome 9 in mouse

The alpha 3, alpha 5, and beta 4 genes (human gene symbols CHRNA3, CHRNA5, and CHRNB4 respectively; mouse gene symbols Acra-3, Acra-5, and Acrb-4, respectively) are members of the nicotinic acetylcholine receptor gene family and are clustered within a 68-kb segment of the rat genome (Boulter et al., 1990, J. Biol. Chem. 265:4472). By somatic cell hybrid analysis, three cDNAs corresponding to these genes were used to map the homologous loci to human chromosome 15 and to mouse chromosome 9. Linkage analysis using CEPH pedigrees showed that the CHRNA5 gene was closely linked to the following chromosome 15 loci: D15S46, D15S52, D15S28, D15S34, and D15S35. Using interspecies crosses in mice, the Acra-5 gene was found closely linked to the Mpi-1 locus. The mapping of these members of a neurotransmitter receptor gene family may facilitate the identification of relationships between the neurotransmitter receptors and murine or human phenotypes.     

Functional characterization of the Arabidopsis ubiquitin-specific protease gene family reveals specific role and redundancy of individual members in development

Ubiquitin-specific proteases (UBPs) are a highly conserved family of proteins in eukaryotes, and play critical roles in protein de-ubiquitination. Here we report a systematic genetic and expression profiling analysis of the UBP gene family in the Arabidopsis thaliana genome. Mutation analysis of 25 of the 27 member genes representing 13 of the 14 sub-families of the UBP gene family revealed that single-gene mutants of three genes in two sub-families exhibit visible phenotypes. Two of these three genes belonging to the UBP15 sub-family were selected for further characterization. The ubp15 mutants display narrower, serrated and flat rosette leaves, partially due to a defect in cell proliferation, as well as other phenotypes such as early flowering, weak apical dominance and reduced fertility, while the line over-expressing UBP15 shows opposite phenotypes. We demonstrated that UPB15 has UBP activity in vitro , and that this biochemical activity is essential for its in vivo function. A genetic interaction analysis among members of this sub-family revealed that UBP15 and UBP16, but not UBP17, have functional redundancy. Our data thus suggest that distinct UBPs, even within a closely related sub-family, can function in different developmental pathways. Although there are clearly functional redundancies among related sub-family members, those redundancies cannot be inferred simply based on the amino acid identity of the family members.     

Genomic expression patterns distinguish long-term from short-term glioblastoma survivors: a preliminary feasibility study

We used microarray analysis to investigate associations between genotypic expression profiles and survival phenotypes in patients with primary glioblastoma (GBM). Tumor samples from 7 long-term glioblastoma survivors (>24 months) and 13 short-term survivors (<9 months) were analyzed to detect differential patterns of gene expression between these groups and to identify genotypic subclasses of glioblastomas that correlate with survival phenotypes. Five unsupervised and three supervised clustering algorithms consistently and accurately grouped the tumors into genotypic subgroups corresponding to the two clinical survival phenotypes. Three unique prospective mathematical classification algorithms were subsequently trained to use expression data to stratify unknown glioblastomas between survival groups and performed this task with 100% accuracy in validation studies. A set of 1478 genes with significant differential expression (p<0.01) between long-term and short-term survivors was identified, and additional mathematical filtering was used to isolate a 43-gene "fingerprint" that distinguished survival phenotypes. Differential regulation of a subset of these genes was confirmed using RT-PCR. Gene ontology analysis of the fingerprint demonstrated pathophysiologic functions for the gene products that are consistent with current models of tumor biology, suggesting that differential expression of these genes may contribute etiologically to the observed differences in survival. These results demonstrate that unique expression profiles characterize genotypic subsets of primary GBMs associated with differential survival phenotypes, and these profiles can be used in a prospective fashion to assign unknown tumors to survival groups. Future efforts will focus on building more robust classifiers and identifying additional subclasses of gliomas with phenotypic significance.     

Genome-wide identification and characterisation of F-box family in maize

F-box-containing proteins, as the key components of the protein degradation machinery, are widely distributed in higher plants and are considered as one of the largest known families of regulatory proteins. The F-box protein family plays a crucial role in plant growth and development and in response to biotic and abiotic stresses. However, systematic analysis of the F-box family in maize (Zea mays) has not been reported yet. In this paper, we identified and characterised the maize F-box genes in a genome-wide scale, including phylogenetic analysis, chromosome distribution, gene structure, promoter analysis and gene expression profiles. A total of 359 F-box genes were identified and divided into 15 subgroups by phylogenetic analysis. The F-box domain was relatively conserved, whereas additional motifs outside the F-box domain may indicate the functional diversification of maize F-box genes. These genes were unevenly distributed in ten maize chromosomes, suggesting that they expanded in the maize genome because of tandem and segmental duplication events. The expression profiles suggested that the maize F-box genes had temporal and spatial expression patterns. Putative cis-acting regulatory DNA elements involved in abiotic stresses were observed in maize F-box gene promoters. The gene expression profiles under abiotic stresses also suggested that some genes participated in stress responsive pathways. Furthermore, ten genes were chosen for quantitative real-time PCR analysis under drought stress and the results were consistent with the microarray data. This study has produced a comparative genomics analysis of the maize ZmFBX gene family that can be used in further studies to uncover their roles in maize growth and development.     

The Rice B-Box Zinc Finger Gene Family: Genomic Identification,Characterization, Expression Profiling and Diurnal Analysis

Background

The B-box (BBX) -containing proteins are a class of zinc finger proteins that contain one or two B-box domains and play important roles in plant growth and development. The Arabidopsis BBX gene family has recently been re-identified and renamed. However, there has not been a genome-wide survey of the rice BBX (OsBBX) gene family until now.

Methodology/Principal Findings

In this study, we identified 30 rice BBX genes through a comprehensive bioinformatics analysis. Each gene was assigned a uniform nomenclature. We described the chromosome localizations, gene structures, protein domains, phylogenetic relationship, whole life-cycle expression profile and diurnal expression patterns of the OsBBX family members. Based on the phylogeny and domain constitution, the OsBBX gene family was classified into five subfamilies. The gene duplication analysis revealed that only chromosomal segmental duplication contributed to the expansion of the OsBBX gene family. The expression profile of the OsBBX genes was analyzed by Affymetrix GeneChip microarrays throughout the entire life-cycle of rice cultivar Zhenshan 97 (ZS97). In addition, microarray analysis was performed to obtain the expression patterns of these genes under light/dark conditions and after three phytohormone treatments. This analysis revealed that the expression patterns of the OsBBX genes could be classified into eight groups. Eight genes were regulated under the light/dark treatments, and eleven genes showed differential expression under at least one phytohormone treatment. Moreover, we verified the diurnal expression of the OsBBX genes using the data obtained from the Diurnal Project and qPCR analysis, and the results indicated that many of these genes had a diurnal expression pattern.

Conclusions/Significance

The combination of the genome-wide identification and the expression and diurnal analysis of the OsBBX gene family should facilitate additional functional studies of the OsBBX genes.     

An Integrative Analysis to Identify Driver Genes in Esophageal Squamous Cell Carcinoma

BackgroundFew driver genes have been well established in esophageal squamous cell carcinoma (ESCC). Identification of the genomic aberrations that contribute to changes in gene expression profiles can be used to predict driver genes.MethodsWe searched for driver genes in ESCC by integrative analysis of gene expression microarray profiles and copy number data. To narrow down candidate genes, we performed survival analysis on expression data and tested the genetic vulnerability of each genes using public RNAi screening data. We confirmed the results by performing RNAi experiments and evaluating the clinical relevance of candidate genes in an independent ESCC cohort.ResultsWe found 10 significantly recurrent copy number alterations accompanying gene expression changes, including loci 11q13.2, 7p11.2, 3q26.33, and 17q12, which harbored CCND1, EGFR, SOX2, and ERBB2, respectively. Analysis of survival data and RNAi screening data suggested that GRB7, located on 17q12, was a driver gene in ESCC. In ESCC cell lines harboring 17q12 amplification, knockdown of GRB7 reduced the proliferation, migration, and invasion capacities of cells. Moreover, siRNA targeting GRB7 had a synergistic inhibitory effect when combined with trastuzumab, an anti-ERBB2 antibody. Survival analysis of the independent cohort also showed that high GRB7 expression was associated with poor prognosis in ESCC.ConclusionOur integrative analysis provided important insights into ESCC pathogenesis. We identified GRB7 as a novel ESCC driver gene and potential new therapeutic target.     

Identification of genes associated with shell color in the black-lipped pearl oyster,Pinctada margaritifera

Background

Color polymorphism in the nacre of pteriomorphian bivalves is of great interest for the pearl culture industry. The nacreous layer of the Polynesian black-lipped pearl oyster Pinctada margaritifera exhibits a large array of color variation among individuals including reflections of blue, green, yellow and pink in all possible gradients. Although the heritability of nacre color variation patterns has been demonstrated by experimental crossing, little is known about the genes involved in these patterns. In this study, we identify a set of genes differentially expressed among extreme color phenotypes of P. margaritifera using a suppressive and subtractive hybridization (SSH) method comparing black phenotypes with full and half albino individuals.

Results

Out of the 358 and 346 expressed sequence tags (ESTs) obtained by conducting two SSH libraries respectively, the expression patterns of 37 genes were tested with a real-time quantitative PCR (RT-qPCR) approach by pooling five individuals of each phenotype. The expression of 11 genes was subsequently estimated for each individual in order to detect inter-individual variation. Our results suggest that the color of the nacre is partially under the influence of genes involved in the biomineralization of the calcitic layer. A few genes involved in the formation of the aragonite tablets of the nacre layer and in the biosynthesis chain of melanin also showed differential expression patterns. Finally, high variability in gene expression levels were observed within the black phenotypes.

Conclusions

Our results revealed that three main genetic processes were involved in color polymorphisms: the biomineralization of the nacreous and calcitic layers and the synthesis of pigments such as melanin, suggesting that color polymorphism takes place at different levels in the shell structure. The high variability of gene expression found within black phenotypes suggests that the present work should serve as a basis for future studies exploring more thoroughly the expression patterns of candidate genes within black phenotypes with different dominant iridescent colors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1776-x) contains supplementary material, which is available to authorized users.     

Genetic Modifiers of Neurofibromatosis Type 1-Associated Café-au-Lait Macule Count Identified Using Multi-platform Analysis

Neurofibromatosis type 1 (NF1) is an autosomal dominant, monogenic disorder of dysregulated neurocutaneous tissue growth. Pleiotropy, variable expressivity and few NF1 genotype-phenotype correlates limit clinical prognostication in NF1. Phenotype complexity in NF1 is hypothesized to derive in part from genetic modifiers unlinked to the NF1 locus. In this study, we hypothesized that normal variation in germline gene expression confers risk for certain phenotypes in NF1. In a set of 79 individuals with NF1, we examined the association between gene expression in lymphoblastoid cell lines with NF1-associated phenotypes and sequenced select genes with significant phenotype/expression correlations. In a discovery cohort of 89 self-reported European-Americans with NF1 we examined the association between germline sequence variants of these genes with café-au-lait macule (CALM) count, a tractable, tumor-like phenotype in NF1. Two correlated, common SNPs (rs4660761 and rs7161) between DPH2 and ATP6V0B were significantly associated with the CALM count. Analysis with tiled regression also identified SNP rs4660761 as significantly associated with CALM count. SNP rs1800934 and 12 rare variants in the mismatch repair gene MSH6 were also associated with CALM count. Both SNPs rs7161 and rs4660761 (DPH2 and ATP6V0B) were highly significant in a mega-analysis in a combined cohort of 180 self-reported European-Americans; SNP rs1800934 (MSH6) was near-significant in a meta-analysis assuming dominant effect of the minor allele. SNP rs4660761 is predicted to regulate ATP6V0B, a gene associated with melanosome biology. Individuals with homozygous mutations in MSH6 can develop an NF1-like phenotype, including multiple CALMs. Through a multi-platform approach, we identified variants that influence NF1 CALM count.     

TTC4,a Novel Human Gene Containing the Tetratricopeptide Repeat and Mapping to the Region of Chromosome 1p31 That Is Frequently Deleted in Sporadic Breast Cancer

The 1p31 region shows loss of heterozygosity in up to 50% of human breast cancers, indicating the presence of a tumor suppressor gene in this location. We have mapped six novel ESTs to a 15-Mb contig of yeast artificial chromosomes spanning the critical region of 1p31. One of these ESTs was localized within the contig to the region most commonly undergoing loss of heterozygosity in breast cancer. The corresponding gene sequence for this EST was established by cDNA cloning and RACE procedures. This gene is 2 kb long and contains a tetratricopeptide repeat motif and a coiled-coil domain. This family of genes has been implicated in a wide variety of functions, including tumorigenesis. This is the fourth member of the human gene family, and so we have named this geneTTC4.Northern blot analysis demonstrates a ubiquitous pattern of gene expression that includes breast tissue. A preliminary screen of human breast cancer cell lines shows thatTTC4is expressed in all cases, but SSCP analysis of the coding region of this gene following RT-PCR failed to reveal any mutations. Clearly, because of its map location, a more extensive analysis is warranted to establish whether subtle mutations are present in breast cancers.     

Pathway deregulation and expression QTLs in response to Actinobacillus pleuropneumoniae infection in swine

Actinobacillus (A.) pleuropneumoniae is among the most important pathogens in pig. The agent causes severe economic losses due to decreased performance, the occurrence of acute or chronic pleuropneumonia, and an increase in death incidence. Since therapeutics cannot be used in a sustainable manner, and vaccination is not always available, new prophylactic measures are urgently needed. Recent research has provided evidence for a genetic predisposition in susceptibility to A. pleuropneumoniae in a Hampshire × German Landrace F2 family with 170 animals. The aim of the present study is to characterize the expression response in this family in order to unravel resistance and susceptibility mechanisms and to prioritize candidate genes for future fine mapping approaches. F2 pigs differed distinctly in clinical, pathological, and microbiological parameters after challenge with A. pleuropneumoniae. We monitored genome-wide gene expression from the 50 most and 50 least susceptible F2 pigs and identified 171 genes differentially expressed between these extreme phenotypes. We combined expression QTL analyses with network analyses and functional characterization using gene set enrichment analysis and identified a functional hotspot on SSC13, including 55 eQTL. The integration of the different results provides a resource for candidate prioritization for fine mapping strategies, such as TF, TFRC, RUNX1, TCN1, HP, CD14, among others.