Comparison of 12 reference genes for normalization of gene expression levels in Epstein-Barr virus-transformed lymphoblastoid cell lines and fibroblasts

Functional pathogenicity of ESRRB variant of uncertain significance contributes to hearing loss (DFNB35)

Advances in next-generation sequencing technologies have led to elucidation of sensorineural hearing loss genetics and associated clinical impacts. However, studies on the functional pathogenicity of variants of uncertain significance (VUS), despite their close association with clinical phenotypes, are lacking. Here we identified compound heterozygous variants in ESRRB transcription factor gene linked to DFNB35, specifically a novel splicing variant (NM_004452.4(ESRRB): c.397 + 2T>G) in trans with a missense variant (NM_004452.4(ESRRB): c.1144C>T p.(Arg382Cys)) whose pathogenicity remains unclear. The splicing variant (c.397 + 2T>G) caused exon 4 skipping, leading to premature stop codon formation and nonsense-mediated decay. The p.(Arg382Cys) variant was classified as a VUS due to its particularly higher allele frequency among East Asian population despite disease-causing in-silico predictions. However, functional assays showed that p.(Arg382Cys) variant disrupted key intramolecular interactions, leading to protein instability. This variant also reduced transcriptional activity and altered expression of downstream target genes essential for inner ear function, suggesting genetic contribution to disease phenotype. This study expanded the phenotypic and genotypic spectrum of ESRRB in DFNB35 and revealed molecular mechanisms underlying ESRRB-associated DFNB35. These findings suggest that variants with high allele frequencies can also possess functional pathogenicity, providing a breakthrough for cases where VUS, previously unexplored, could be reinterpreted by elucidating their functional roles and disease-causing characteristics.

Biallelic variants in TMEM222 cause a new autosomal recessive neurodevelopmental disorder

PURPOSE

To elucidate the novel molecular cause in families with a new autosomal recessive neurodevelopmental disorder.

METHODS

A combination of exome sequencing and gene matching tools was used to identify pathogenic variants in 17 individuals. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) and subcellular localization studies were used to characterize gene expression profile and localization.

RESULTS

Biallelic variants in the TMEM222 gene were identified in 17 individuals from nine unrelated families, presenting with intellectual disability and variable other features, such as aggressive behavior, shy character, body tremors, decreased muscle mass in the lower extremities, and mild hypotonia. We found relatively high TMEM222 expression levels in the human brain, especially in the parietal and occipital cortex. Additionally, subcellular localization analysis in human neurons derived from induced pluripotent stem cells (iPSCs) revealed that TMEM222 localizes to early endosomes in the synapses of mature iPSC-derived neurons.

CONCLUSION

Our findings support a role for TMEM222 in brain development and function and adds variants in the gene TMEM222 as a novel underlying cause of an autosomal recessive neurodevelopmental disorder.

De Novo Variants in SPOP Cause Two Clinically Distinct Neurodevelopmental Disorders

Recurrent somatic variants in SPOP are cancer specific; endometrial and prostate cancers result from gain-of-function and dominant-negative effects toward BET proteins, respectively. By using clinical exome sequencing, we identified six de novo pathogenic missense variants in SPOP in seven individuals with developmental delay and/or intellectual disability, facial dysmorphisms, and congenital anomalies. Two individuals shared craniofacial dysmorphisms, including congenital microcephaly, that were strikingly different from those of the other five individuals, who had (relative) macrocephaly and hypertelorism. We measured the effect of SPOP variants on BET protein amounts in human Ishikawa endometrial cancer cells and patient-derived cell lines because we hypothesized that variants would lead to functional divergent effects on BET proteins. The de novo variants c.362G>A (p.Arg121Gln) and c. 430G>A (p.Asp144Asn), identified in the first two individuals, resulted in a gain of function, and conversely, the c.73A>G (p.Thr25Ala), c.248A>G (p.Tyr83Cys), c.395G>T (p.Gly132Val), and c.412C>T (p.Arg138Cys) variants resulted in a dominant-negative effect. Our findings suggest that these opposite functional effects caused by the variants in SPOP result in two distinct and clinically recognizable syndromic forms of intellectual disability with contrasting craniofacial dysmorphisms.

De novo variants in CDK13 associated with syndromic ID/DD: Molecular and clinical delineation of 15 individuals and a further review

De novo variants in the gene encoding cyclin-dependent kinase 13 (CDK13) have been associated with congenital heart defects and intellectual disability (ID). Here, we present the clinical assessment of 15 individuals and report novel de novo missense variants within the kinase domain of CDK13. Furthermore, we describe 2 nonsense variants and a recurrent frame-shift variant. We demonstrate the synthesis of 2 aberrant CDK13 transcripts in lymphoblastoid cells from an individual with a splice-site variant. Clinical characteristics of the individuals include mild to severe ID, developmental delay, behavioral problems, (neonatal) hypotonia and a variety of facial dysmorphism. Congenital heart defects were present in 2 individuals of the current cohort, but in at least 42% of all known individuals. An overview of all published cases is provided and does not demonstrate an obvious genotype-phenotype correlation, although 2 individuals harboring a stop codons at the end of the kinase domain might have a milder phenotype. Overall, there seems not to be a clinically recognizable facial appearance. The variability in the phenotypes impedes an à vue diagnosis of this syndrome and therefore genome-wide or gene-panel driven genetic testing is needed. Based on this overview, we provide suggestions for clinical work-up and management of this recently described ID syndrome.

RNA sequencing of bipolar disorder lymphoblastoid cell lines implicates the neurotrophic factor HRP-3 in lithium's clinical efficacy

Objectives: Lithium remains the oldest and most effective treatment for mood stabilisation in bipolar disorder (BD), even though at least half of patients are only partially responsive or do not respond. This study aimed to identify biomarkers associated with lithium response in BD, based on comparing RNA sequencing information derived from lymphoblastoid cell lines (LCLs) of lithium-responsive (LR) versus lithium non-responsive (LNR) BD patients, to assess gene expression variations that might bear on treatment outcome. Methods: RNA sequencing was carried out on 24 LCLs from female BD patients (12 LR and 12 LNR) followed by qPCR validation in two additional independent cohorts (41 and 17 BD patients, respectively). Results: Fifty-six genes showed nominal differential expression comparing LR and LNR (FC ≥ |1.3|, P ≤ 0.01). The differential expression of HDGFRP3 and ID2 was validated by qPCR in the independent cohorts. Conclusions: We observed higher expression levels of HDGFRP3 and ID2 in BD patients who favourably respond to lithium. Both of these genes are involved in neurogenesis, and HDGFRP3 has been suggested to be a neurotrophic factor. Additional studies in larger BD cohorts are needed to confirm the potential of HDGFRP3 and ID2 expression levels in blood cells as tentative favourable lithium response biomarkers.

SIRT1, miR-132 and miR-212 link human longevity to Alzheimer's Disease

Alzheimer's Disease (AD) is the most common cause of dementia in the elderly. Centenarians - reaching the age of >100 years while maintaining good cognitive skills - seemingly have unique biological features allowing healthy aging and protection from dementia. Here, we studied the expression of SIRT1 along with miR-132 and miR-212, two microRNAs known to regulate SIRT1, in lymphoblastoid cell lines (LCLs) from 45 healthy donors aged 21 to 105 years and 24 AD patients, and in postmortem olfactory bulb and hippocampus tissues from 14 AD patients and 20 age-matched non-demented individuals. We observed 4.0-fold (P = 0.001) lower expression of SIRT1, and correspondingly higher expression of miR-132 (1.7-fold; P = 0.014) and miR-212 (2.1-fold; P = 0.036), in LCLs from AD patients compared with age-matched healthy controls. Additionally, SIRT1 expression was 2.2-fold (P = 0.001) higher in centenarian LCLs compared with LCLs from individuals aged 56-82 years; while centenarian LCLs miR-132 and miR-212 indicated 7.6-fold and 4.1-fold lower expression, respectively. Correlations of SIRT1, miR-132 and miR-212 expression with cognitive scores were observed for AD patient-derived LCLs and postmortem AD olfactory bulb and hippocampus tissues, suggesting that higher SIRT1 expression, possibly mediated by lower miR-132 and miR-212, may protect aged individuals from dementia and is reflected in their peripheral tissues.

The BRCA2 c.68-7T > A variant is not pathogenic: A model for clinical calibration of spliceogenicity

Although the spliceogenic nature of the BRCA2 c.68-7T > A variant has been demonstrated, its association with cancer risk remains controversial. In this study, we accurately quantified by real-time PCR and digital PCR (dPCR), the BRCA2 isoforms retaining or missing exon 3. In addition, the combined odds ratio for causality of the variant was estimated using genetic and clinical data, and its associated cancer risk was estimated by case-control analysis in 83,636 individuals. Co-occurrence in trans with pathogenic BRCA2 variants was assessed in 5,382 families. Exon 3 exclusion rate was 4.5-fold higher in variant carriers (13%) than controls (3%), indicating an exclusion rate for the c.68-7T > A allele of approximately 20%. The posterior probability of pathogenicity was 7.44 × 10-115 . There was neither evidence for increased risk of breast cancer (OR 1.03; 95% CI 0.86-1.24) nor for a deleterious effect of the variant when co-occurring with pathogenic variants. Our data provide for the first time robust evidence of the nonpathogenicity of the BRCA2 c.68-7T > A. Genetic and quantitative transcript analyses together inform the threshold for the ratio between functional and altered BRCA2 isoforms compatible with normal cell function. These findings might be exploited to assess the relevance for cancer risk of other BRCA2 spliceogenic variants.

CXorf56, a dendritic neuronal protein, identified as a new candidate gene for X-linked intellectual disability

Intellectual disability (ID) comprises a large group of heterogeneous disorders, often without a known molecular cause. X-linked ID accounts for 5-10% of male ID cases. We investigated a large, three-generation family with mild ID and behavior problems in five males and one female, with a segregation suggestive for X-linked inheritance. Linkage analysis mapped a disease locus to a 7.6 Mb candidate region on the X-chromosome (LOD score 3.3). Whole-genome sequencing identified a 2 bp insertion in exon 2 of the chromosome X open reading frame 56 gene (CXorf56), resulting in a premature stop codon. This insertion was present in all intellectually impaired individuals and carrier females. Additionally, X-inactivation status showed skewed methylation patterns favoring the inactivation of the mutated allele in the unaffected carrier females. We demonstrate that the insertion leads to nonsense-mediated decay and that CXorf56 mRNA expression is reduced in the impaired males and female. In murine brain slices and primary hippocampal neuronal cultures, CXorf56 protein was present and localized in the nucleus, cell soma, dendrites, and dendritic spines. Although no other families have been identified with pathogenic variants in CXorf56, these results suggest that CXorf56 is the causative gene in this family, and thus a novel candidate gene for X-linked ID with behavior problems.

RGS2 expression predicts amyloid-β sensitivity, MCI and Alzheimer's disease: genome-wide transcriptomic profiling and bioinformatics data mining

Alzheimer's disease (AD) is the most frequent cause of dementia. Misfolded protein pathological hallmarks of AD are brain deposits of amyloid-β (Aβ) plaques and phosphorylated tau neurofibrillary tangles. However, doubts about the role of Aβ in AD pathology have been raised as Aβ is a common component of extracellular brain deposits found, also by in vivo imaging, in non-demented aged individuals. It has been suggested that some individuals are more prone to Aβ neurotoxicity and hence more likely to develop AD when aging brains start accumulating Aβ plaques. Here, we applied genome-wide transcriptomic profiling of lymphoblastoid cells lines (LCLs) from healthy individuals and AD patients for identifying genes that predict sensitivity to Aβ. Real-time PCR validation identified 3.78-fold lower expression of RGS2 (regulator of G-protein signaling 2; P=0.0085) in LCLs from healthy individuals exhibiting high vs low Aβ sensitivity. Furthermore, RGS2 showed 3.3-fold lower expression (P=0.0008) in AD LCLs compared with controls. Notably, RGS2 expression in AD LCLs correlated with the patients' cognitive function. Lower RGS2 expression levels were also discovered in published expression data sets from postmortem AD brain tissues as well as in mild cognitive impairment and AD blood samples compared with controls. In conclusion, Aβ sensitivity phenotyping followed by transcriptomic profiling and published patient data mining identified reduced peripheral and brain expression levels of RGS2, a key regulator of G-protein-coupled receptor signaling and neuronal plasticity. RGS2 is suggested as a novel AD biomarker (alongside other genes) toward early AD detection and future disease modifying therapeutics.

De novo mutations in PLXND1 and REV3L cause Möbius syndrome

Möbius syndrome (MBS) is a neurological disorder that is characterized by paralysis of the facial nerves and variable other congenital anomalies. The aetiology of this syndrome has been enigmatic since the initial descriptions by von Graefe in 1880 and by Möbius in 1888, and it has been debated for decades whether MBS has a genetic or a non-genetic aetiology. Here, we report de novo mutations affecting two genes, PLXND1 and REV3L in MBS patients. PLXND1 and REV3L represent totally unrelated pathways involved in hindbrain development: neural migration and DNA translesion synthesis, essential for the replication of endogenously damaged DNA, respectively. Interestingly, analysis of Plxnd1 and Rev3l mutant mice shows that disruption of these separate pathways converge at the facial branchiomotor nucleus, affecting either motoneuron migration or proliferation. The finding that PLXND1 and REV3L mutations are responsible for a proportion of MBS patients suggests that de novo mutations in other genes might account for other MBS patients.

lt has been debated for decades if there is a genetic aetiology underlying Möbius syndrome, a neurological disorder characterized by facial paralysis. Here Tomas-Roca et al. use exome sequencing and identify de novo mutations in PLXND1 and REV3L, representing converging pathways in hindbrain development.

Latent Membrane Protein LMP2A Impairs Recognition of EBV-Infected Cells by CD8+ T Cells

The common pathogen Epstein-Barr virus (EBV) transforms normal human B cells and can cause cancer. Latent membrane protein 2A (LMP2A) of EBV supports activation and proliferation of infected B cells and is expressed in many types of EBV-associated cancer. It is not clear how latent EBV infection and cancer escape elimination by host immunity, and it is unknown whether LMP2A can influence the interaction of EBV-infected cells with the immune system. We infected primary B cells with EBV deleted for LMP2A, and established lymphoblastoid cell lines (LCLs). We found that CD8+ T cell clones showed higher reactivity against LMP2A-deficient LCLs compared to LCLs infected with complete EBV. We identified several potential mediators of this immunomodulatory effect. In the absence of LMP2A, expression of some EBV latent antigens was elevated, and cell surface expression of MHC class I was marginally increased. LMP2A-deficient LCLs produced lower amounts of IL-10, although this did not directly affect CD8+ T cell recognition. Deletion of LMP2A led to several changes in the cell surface immunophenotype of LCLs. Specifically, the agonistic NKG2D ligands MICA and ULBP4 were increased. Blocking experiments showed that NKG2D activation contributed to LCL recognition by CD8+ T cell clones. Our results demonstrate that LMP2A reduces the reactivity of CD8+ T cells against EBV-infected cells, and we identify several relevant mechanisms.

Epstein-Barr virus (EBV) is carried by most humans. It can cause several types of cancer. In healthy infected people, EBV persists for life in a "latent" state in white blood cells called B cells. For infected persons to remain healthy, it is crucial that they harbor CD8-positive "killer" T cells that recognize and destroy precancerous EBV-infected cells. However, this protection is imperfect, because the virus is not eliminated from the body, and the danger of EBV-associated cancer remains. How does the virus counteract CD8+ T cell control? Here we study the effects of latent membrane protein 2A (LMP2A), which is an important viral molecule because it is present in several types of EBV-associated cancers, and in latently infected cells in healthy people. We show that LMP2A counteracts the recognition of EBV-infected B cells by antiviral killer cells. We found a number of mechanisms that are relevant to this effect. Notably, LMP2A disturbs expression of molecules on B cells that interact with NKG2D, a molecule on the surface of CD8+ T cells that aids their activation. In this way, LMP2A weakens important immune responses against EBV. Similar mechanisms may operate in different types of LMP2A-expressing cancers caused by EBV.

A novel single nucleotide splice site mutation in FHL1 confirms an Emery-Dreifuss plus phenotype with pulmonary artery hypoplasia and facial dysmorphology

We describe a Danish family with an, until recently, unknown X-linked disease with muscular dystrophy (MD), facial dysmorphology and pulmonary artery hypoplasia. One patient died suddenly before age 20 and another was resuscitated from cardiac arrest at the age of 28. Linkage analysis pointed to a region of 25 Mb from 123.6 Mb to 148.4 Mb on chromosome X containing over 100 genes. Exome sequencing identified a single nucleotide splice site mutation c.502-2A > T, which is located 5' to exon 6 in the gene encoding four and a half LIM domain 1 (FHL1) protein. FHL1 expresses three main splice variants, known as FHL1A, FHL1B and FHL1C. In healthy individuals, FHL1A is the predominant splice variant and is mainly found in skeletal and cardiac muscle. The FHL1 transcript profiles from two affected individuals were investigated in skin fibroblasts with quantitative real-time PCR. This demonstrated loss of isoform A and B, and an almost 200-fold overexpression of isoform C confirming that lack of FHL1A and overexpression of FHL1C results in an extended phenotype of EDMD as recently shown by Tiffin et al. [2013].

TDP2 protects transcription from abortive topoisomerase activity and is required for normal neural function

Topoisomerase II (TOP2) removes torsional stress from DNA and facilitates gene transcription by introducing transient DNA double-strand breaks (DSBs). Such DSBs are normally rejoined by TOP2 but on occasion can become abortive and remain unsealed. Here we identify homozygous mutations in the TDP2 gene encoding tyrosyl DNA phosphodiesterase-2, an enzyme that repairs 'abortive' TOP2-induced DSBs, in individuals with intellectual disability, seizures and ataxia. We show that cells from affected individuals are hypersensitive to TOP2-induced DSBs and that loss of TDP2 inhibits TOP2-dependent gene transcription in cultured human cells and in mouse post-mitotic neurons following abortive TOP2 activity. Notably, TDP2 is also required for normal levels of many gene transcripts in developing mouse brain, including numerous gene transcripts associated with neurological function and/or disease, and for normal interneuron density in mouse cerebellum. Collectively, these data implicate chromosome breakage by TOP2 as an endogenous threat to gene transcription and to normal neuronal development and maintenance.

Fibroblasts from patients affected by Pseudoxanthoma elasticum exhibit an altered PPi metabolism and are more responsive to pro-calcifying stimuli

BACKGROUND

Pseudoxanthoma elasticum (PXE) is a genetic disorder characterized by progressive calcification of soft connective tissues. The pathogenesis is still hard to pin down. In PXE dermal fibroblasts, in addition to impaired carboxylation of the vitamin K-dependent inhibitor matrix Gla protein (MGP), we have also demonstrated an up-regulation of alkaline phosphatase activity. In the light of these data we have suggested that both calcium and phosphate metabolism might be locally altered, both pathways acting in synergy on the occurrence of matrix calcification.

OBJECTIVE

This study aims to better explore if cultured PXE fibroblasts, compared to control cells, exhibit a modified inorganic pyrophosphate (PPi) metabolism and are more responsive to pro-calcifying stimuli.

METHODS

Primary human dermal fibroblasts isolated from healthy individuals and from PXE patients were cultured for different time points in standard and in pro-calcifying media. The expression of ANKH/ANKH, ENPP1/PC1, ALPL/TNAP, SPP1/OPN was evaluated by qRT-PCR and Western blot, respectively. TNAP activity was measured by spectrophotometric analyses, whereas calcification was investigated by light and electron microscopy as well as by micro-analytical techniques.

RESULTS

In the presence of pro-calcifying stimuli, dermal fibroblasts alter their phenotype favouring matrix mineralization. In particular, ENPP1/PC1 and SPP1/OPN expression, as well as TNAP activity, was differently expressed in control and in PXE fibroblasts. Moreover, in pathologic cells the ratio between factors favouring and reducing PPi availability exhibits a more pronounced shift towards a pro-calcifying balance.

CONCLUSION

PXE fibroblasts are more susceptible to pro-calcifying stimuli and in these cells an altered PPi metabolism contributes to matrix calcification.

Early presentation of cystic kidneys in a family with a homozygous INVS mutation

Nephronophthisis (NPHP) is an autosomal recessive cystic kidney disease that is the most frequent monogenic cause of end-stage renal disease in children. Infantile NPHP, often in combination with other features like situs inversus, are commonly caused by mutations in the INVS gene. INVS encodes the ciliary protein inversin, and mutations induce dysfunction of the primary cilia. In this article, we present a family with two severely affected fetuses that were aborted after discovery of grossly enlarged cystic kidneys by ultrasonography before 22 weeks gestation. Exome sequencing showed that the fetuses were homozygous for a previously unreported nonsense mutation, resulting in a truncation in the IQ1 domain of inversin. This mutation induces nonsense-mediated RNA decay, as suggested by a reduced RNA level in fibroblasts derived from the fetus. However, a significant amount of mutant INVS RNA was present in these fibroblasts, yielding mutant inversin protein that was mislocalized. In control fibroblasts, inversin was present in the ciliary axoneme as well as at the basal body, whereas in the fibroblasts from the fetus, inversin could only be detected at the basal body. The phenotype of both fetuses is partly characteristic of infantile NPHP and Potter sequence. We also identified that the fetuses had mild skeletal abnormalities, including shortening and bowing of long bones, which may expand the phenotypic spectrum associated with INVS mutations.

A 3-base pair deletion, c.9711_9713del, in DMD results in intellectual disability without muscular dystrophy

We have identified a deletion of 3 base pairs in the dystrophin gene (DMD), c.9711_9713del, in a family with nonspecific X-linked intellectual disability (ID) by sequencing of the exons of 86 known X-linked ID genes. This in-frame deletion results in the deletion of a single-amino-acid residue, Leu3238, in the brain-specific isoform Dp71 of dystrophin. Linkage analysis supported causality as the mutation was present in the 7.6 cM linkage interval on Xp22.11-Xp21.1 with a maximum positive LOD score of 2.41 (MRX85 locus). Molecular modeling predicts that the p.(Leu3238del) deletion results in the destabilization of the C-terminal domain of dystrophin and hence reduces the ability to interact with β-dystroglycan. Correspondingly, Dp71 protein levels in lymphoblastoid cells from the index patient are 6.7-fold lower than those in control cell lines (P=0.08). Subsequent determination of the creatine kinase levels in blood of the index patient showed a mild but significant elevation in serum creatine kinase, which is in line with impaired dystrophin function. In conclusion, we have identified the first DMD mutation in Dp71 that results in ID without muscular dystrophy.

Validation of microarray data in human lymphoblasts shows a role of the ubiquitin-proteasome system and NF-kB in the pathogenesis of Down syndrome

BACKGROUND

Down syndrome (DS) is a complex disorder caused by the trisomy of either the entire, or a critical region of chromosome 21 (21q22.1-22.3). Despite representing the most common cause of mental retardation, the molecular bases of the syndrome are still largely unknown.

METHODS

To better understand the pathogenesis of DS, we analyzed the genome-wide transcription profiles of lymphoblastoid cell lines (LCLs) from six DS and six euploid individuals and investigated differential gene expression and pathway deregulation associated with trisomy 21. Connectivity map and PASS-assisted exploration were used to identify compounds whose molecular signatures counteracted those of DS lymphoblasts and to predict their therapeutic potential. An experimental validation in DS LCLs and fetal fibroblasts was performed for the most deregulated GO categories, i.e. the ubiquitin mediated proteolysis and the NF-kB cascade.

RESULTS

We show, for the first time, that the level of protein ubiquitination is reduced in human DS cell lines and that proteasome activity is increased in both basal conditions and oxidative microenvironment. We also provide the first evidence that NF-kB transcription levels, a paradigm of gene expression control by ubiquitin-mediated degradation, is impaired in DS due to reduced IkB-alfa ubiquitination, increased NF-kB inhibitor (IkB-alfa) and reduced p65 nuclear fraction. Finally, the DSCR1/DYRK1A/NFAT genes were analysed. In human DS LCLs, we confirmed the presence of increased protein levels of DSCR1 and DYRK1A, and showed that the levels of the transcription factor NFATc2 were decreased in DS along with a reduction of its nuclear translocation upon induction of calcium fluxes.

CONCLUSIONS

The present work offers new perspectives to better understand the pathogenesis of DS and suggests a rationale for innovative approaches to treat some pathological conditions associated to DS.

Genome-wide transcriptomic variations of human lymphoblastoid cell lines: insights from pairwise gene-expression correlations

Aims: Human lymphoblastoid cell lines (LCLs) are a rich resource of information on human interindividual genomic, transcriptomic, proteomic and phenomic variations, and are therefore gaining popularity for pharmacogenomic studies. In the present study we demonstrate that genome-wide transcriptomic data from a small LCL panel from unrelated individuals is sufficient for detecting pairs of genes that exhibit highly correlated expression levels and may thus convey insights about coregulated genes. Materials & methods: RNA samples were prepared from LCLs representing 12 unrelated healthy adult female Caucasian donors. Transcript levels were determined with the Affymetrix Human Gene arrays. Expression-level correlations were searched using Partek® Genomics Suite™ and the R environment. Sequences of detected correlated gene pairs were compared for shared conserved 3´-UTR miRNA binding. Results: Most of the approximately 33,000 transcripts covered by the Affymetrix arrays showed closely similar expression levels in LCLs from unrelated donors. However, the expression levels of some transcripts showed large inter-individual variations. When comparing the expression levels of each of the top 1000 genes showing the largest interindividual expression variations against the others, two sets containing 156 and 4438 correlated gene pairs with false-discovery rates of 0.01 and 0.05 were detected, respectively. Similar analysis of another gene-expression data set from LCLs (GSE11582) indicated that 61 and 39% of identified pairs matched the pairs detected from our transcriptomic data, respectively. Shared conserved 3´-UTR miRNA binding sites were noted for 14–17% of the top 100 gene pairs, suggesting that regulation by miRNA may contribute to their coordinated expression. Conclusion: Probing genome-wide transcriptomic data sets of LCLs from unrelated individuals may detect coregulated genes, adding insights on cellular regulation by miRNAs.

Original submitted 11 July 2012; Revision submitted 4 September 2012

The EBV immunoevasins vIL-10 and BNLF2a protect newly infected B cells from immune recognition and elimination

Lifelong persistence of Epstein-Barr virus (EBV) in infected hosts is mainly owed to the virus' pronounced abilities to evade immune responses of its human host. Active immune evasion mechanisms reduce the immunogenicity of infected cells and are known to be of major importance during lytic infection. The EBV genes BCRF1 and BNLF2a encode the viral homologue of IL-10 (vIL-10) and an inhibitor of the transporter associated with antigen processing (TAP), respectively. Both are known immunoevasins in EBV's lytic phase. Here we describe that BCRF1 and BNLF2a are functionally expressed instantly upon infection of primary B cells. Using EBV mutants deficient in BCRF1 and BNLF2a, we show that both factors contribute to evading EBV-specific immune responses during the earliest phase of infection. vIL-10 impairs NK cell mediated killing of infected B cells, interferes with CD4+ T-cell activity, and modulates cytokine responses, while BNLF2a reduces antigen presentation and recognition of newly infected cells by EBV-specific CD8+ T cells. Together, both factors significantly diminish the immunogenicity of EBV-infected cells during the initial, pre-latent phase of infection and may improve the establishment of a latent EBV infection in vivo.

B4GALT1-congenital disorders of glycosylation presents as a non-neurologic glycosylation disorder with hepatointestinal involvement

The clinical phenotype of congenital disorders of glycosylation is heterogeneous, mostly including a severe neurological involvement and multisystem disease. We identified a novel patient with a galactosyltransferase deficiency with mild hepatopathy and coagulation anomalies, but normal psychomotor development. The tissue-specific expression of the defective B4GALT1 gene correlated with the clinical phenotype.

Genetic CD21 deficiency is associated with hypogammaglobulinemia

BACKGROUND

Complement receptor 2 (CR2/CD21) is part of the B-cell coreceptor and expressed by mature B cells and follicular dendritic cells. CD21 is a receptor for C3d-opsonized immune complexes and enhances antigen-specific B-cell responses.

OBJECTIVE

Genetic inactivation of the murine CR2 locus results in impaired humoral immune responses. Here we report the first case of a genetic CD21 deficiency in human subjects.

METHODS

CD21 protein expression was analyzed by means of flow cytometry and Western blotting. CD21 transcripts were quantified by using real-time PCR. The CD21 gene was sequenced. Wild-type and mutant CD21 cDNA expression was studied after transfection of 293T cells. Binding of EBV-gp350 or C3d-containing immune complexes and induction of calcium flux in CD21-deficient B cells were analyzed by means of flow cytometry. Antibody responses to protein and polysaccharide vaccines were measured.

RESULTS

A 28-year-old man presented with recurrent infections, reduced class-switched memory B cells, and hypogammaglobulinemia. CD21 receptor expression was undetectable. Binding of C3d-containing immune complexes and EBV-gp350 to B cells was severely reduced. Sequence analysis revealed a compound heterozygous deleterious mutation in the CD21 gene. Functional studies with anti-immunoglobulin- and C3d-containing immune complexes showed a complete loss of costimulatory activity of C3d in enhancing suboptimal B-cell receptor stimulation. Vaccination responses to protein antigens were normal, but the response to pneumococcal polysaccharide vaccination was moderately impaired.

CONCLUSIONS

Genetic CD21 deficiency adds to the molecular defects observed in human subjects with hypogammaglobulinemia.

The lytic phase of epstein-barr virus requires a viral genome with 5-methylcytosine residues in CpG sites

Epstein-Barr virus (EBV) is a human herpesvirus which has been studied intensively for its role in certain human tumors. It also serves as a model of herpesviral latency because it establishes an immediate, latent infection in human B cells. When EBV infects quiescent, primary B cells it induces their continuous proliferation to yield growth-transformed B-cell lines in vitro. The lytic or productive phase of EBV's life cycle is induced by the expression of the viral BZLF1 gene in latently infected cells. The BZLF1 protein is a transactivator, which selectively binds to two classes of distinct DNA sequence motifs. One class is similar to the motifs that are bound by members of the AP-1 transcription factor family to which BZLF1 belongs. The second class, which contains CpG motifs, is predominant in viral promoters of early lytic genes and is BZLF1's preferred or exclusive target sequence when methylated. The BZLF1 gene is transiently expressed in newly infected B cells but fails to induce EBV's lytic cycle, potentially because the virion DNA is unmethylated. Here we report that the lack of 5-methylcytosine residues in CpG sites of virion DNA prevents the expression of essential lytic genes indispensable for viral DNA amplification during productive infection. This finding indicates that BZLF1 transactivates these promoters in a methylation-dependent fashion and explains how progeny virus synthesis is abrogated in newly infected B cells. Our data also reveal that viral lytic DNA synthesis precludes CpG methylation of virion DNA during EBV's lytic, productive cycle, which can be overcome by the ectopic expression of a prokaryotic cytosine methyltransferase to yield CpG-methylated virion DNA. Upon infection of B cells, randomly CpG-methylated virion DNA induces high expression of essential lytic genes in contrast to virion DNA free of 5-methylcytosine residues. Our data suggest that unmethylated virion DNA is part of EBV's strategy to prevent the viral lytic phase in newly infected B cells, allowing it to establish its characteristic latent infection in them.

The ciliopathy-associated protein homologs RPGRIP1 and RPGRIP1L are linked to cilium integrity through interaction with Nek4 serine/threonine kinase

Recent studies have established ciliary dysfunction as the underlying cause of a broad range of multi-organ phenotypes, known as 'ciliopathies'. Ciliopathy-associated proteins have a common site of action in the cilium, however, their overall importance for ciliary function differs, as implied by the extreme variability in ciliopathy phenotypes. The aim of this study was to gain more insight in the function of two ciliopathy-associated protein homologs, RPGR interacting protein 1 (RPGRIP1) and RPGRIP1-like protein (RPGRIP1L). Mutations in RPGRIP1 lead to the eye-restricted disease Leber congenital amaurosis, while mutations in RPGRIP1L are causative for Joubert and Meckel syndrome, which affect multiple organs and are at the severe end of the ciliopathy spectrum. Using tandem affinity purification in combination with mass spectrometry, we identified Nek4 serine/threonine kinase as a prominent component of both the RPGRIP1- as well as the RPGRIP1L-associated protein complex. In ciliated cells, this kinase localized to basal bodies, while in ciliated organs, the kinase was predominantly detected at the ciliary rootlet. Down-regulation of NEK4 in ciliated cells led to a significant decrease in cilium assembly, pointing to a role for Nek4 in cilium dynamics. We now hypothesize that RPGRIP1 and RPGRIP1L function as cilium-specific scaffolds that recruit a Nek4 signaling network which regulates cilium stability. Our data are in line with previously established roles in the cilium of other members of the Nek protein family and define NEK4 as a ciliopathy candidate gene.

Genome-wide expression profiling of human lymphoblastoid cell lines identifies CHL1 as a putative SSRI antidepressant response biomarker

AIMS

Selective serotonin reuptake inhibitors (SSRIs) are the most commonly used class of antidepressants for treating major depression. However, approximately 30% of patients do not respond sufficiently to first-line antidepressant drug treatment and require alternative therapeutics. Genome-wide studies searching for SSRI response DNA biomarkers or studies of candidate serotonin-related genes so far have given inconclusive or contradictory results. Here, we present an alternative transcriptome-based genome-wide approach for searching antidepressant drug-response biomarkers by using drug-effect phenotypes in human lymphoblastoid cell lines (LCLs).

MATERIALS & METHODS

We screened 80 LCLs from healthy adult female individuals for growth inhibition by paroxetine. A total of 14 LCLs with reproducible high and low sensitivities to paroxetine (seven from each phenotypic group) were chosen for genome-wide expression profiling with commercial microarrays.

RESULTS

The most notable genome-wide transcriptome difference between LCLs displaying high versus low paroxetine sensitivities was a 6.3-fold lower (p = 0.0000256) basal expression of CHL1, a gene coding for a neuronal cell adhesion protein implicated in correct thalamocortical circuitry, schizophrenia and autism. The microarray findings were confirmed by real-time PCR (36-fold lower CHL1 expression levels in the high paroxetine sensitivity group). Several additional genes implicated in synaptogenesis or in psychiatric disorders, including ARRB1, CCL5, DDX60, DDX60L, ENDOD1, ENPP2, FLT1, GABRA4, GAP43, MCTP2 and SPRY2, also differed by more than 1.5-fold and a p-value of less than 0.005 between the two paroxetine sensitivity groups, as confirmed by real-time PCR experiments.

CONCLUSION

Genome-wide transcriptional profiling of in vitro phenotyped LCLs identified CHL1 and additional genes implicated in synaptogenesis and brain circuitry as putative SSRI response biomarkers. This method might be used as a preliminary tool for searching for potential depression treatment biomarkers.

A novel cerebello-ocular syndrome with abnormal glycosylation due to abnormalities in dolichol metabolism

Cerebellar hypoplasia and slowly progressive ophthalmological symptoms are common features in patients with congenital disorders of glycosylation type I. In a group of patients with congenital disorders of glycosylation type I with unknown aetiology, we have previously described a distinct phenotype with severe, early visual impairment and variable eye malformations, including optic nerve hypoplasia, retinal coloboma, congenital cataract and glaucoma. Some of the symptoms overlapped with the phenotype in other congenital disorders of glycosylation type I subtypes, such as vermis hypoplasia, anaemia, ichtyosiform dermatitis, liver dysfunction and coagulation abnormalities. We recently identified pathogenic mutations in the SRD5A3 gene, encoding steroid 5α-reductase type 3, in a group of patients who presented with this particular phenotype and a common metabolic pattern. Here, we report on the clinical, genetic and metabolic features of 12 patients from nine families with cerebellar ataxia and congenital eye malformations diagnosed with SRD5A3-congenital disorders of glycosylation due to steroid 5α-reductase type 3 defect. This enzyme is necessary for the reduction of polyprenol to dolichol, the lipid anchor for N-glycosylation in the endoplasmic reticulum. Dolichol synthesis is an essential metabolic step in protein glycosylation. The current defect leads to a severely abnormal glycosylation state already in the early phase of the N-glycan biosynthesis pathway in the endoplasmic reticulum. We detected high expression of SRD5A3 in foetal brain tissue, especially in the cerebellum, consistent with the finding of the congenital cerebellar malformations. Based on the overlapping clinical, biochemical and genetic data in this large group of patients with congenital disorders of glycosylation, we define a novel syndrome of cerebellar ataxia associated with congenital eye malformations due to a defect in dolichol metabolism.

CDK19 is disrupted in a female patient with bilateral congenital retinal folds, microcephaly and mild mental retardation

Microcephaly, mental retardation and congenital retinal folds along with other systemic features have previously been reported as a separate clinical entity. The sporadic nature of the syndrome and lack of clear inheritance patterns pointed to a genetic heterogeneity. Here, we report a genetic analysis of a female patient with microcephaly, congenital bilateral falciform retinal folds, nystagmus, and mental retardation. Karyotyping revealed a de novo pericentric inversion in chromosome 6 with breakpoints in 6p12.1 and 6q21. Fluorescence in situ hybridization analysis narrowed down the region around the breakpoints, and the breakpoint at 6q21 was found to disrupt the CDK19 gene. CDK19 was found to be expressed in a diverse range of tissues including fetal eye and fetal brain. Quantitative PCR of the CDK19 transcript from Epstein–Barr virus-transformed lymphoblastoid cell lines of the patient revealed ~50% reduction in the transcript (p = 0.02), suggesting haploinsufficiency of the gene. cdk8, the closest orthologue of human CDK19 in Drosophila has been shown to play a major role in eye development. Conditional knock-down of Drosophila cdk8 in multiple dendrite (md) neurons resulted in 35% reduced dendritic branching and altered morphology of the dendritic arbour, which appeared to be due in part to a loss of small higher order branches. In addition, Cdk8 mutant md neurons showed diminished dendritic fields revealing an important role of the CDK19 orthologue in the developing nervous system of Drosophila. This is the first time the CDK19 gene, a component of the mediator co-activator complex, has been linked to a human disease.

IPO8 and FBXL10: new reference genes for gene expression studies in human adipose tissue

Housekeeping genes frequently used in gene expression studies are highly regulated in human adipose tissue. To ensure a correct interpretation of results, it is critical to select appropriate reference genes. Subcutaneous (SC) and omental (OM) adipose tissue expression was analyzed from lean and obese subjects using whole genome complementary DNA (cDNA) microarrays to identify stably expressed genes and commercial TaqMan low density arrays (LDAs), with 16 common control genes. The best candidate gene from microarrays analysis was F-box and leucine-rich repeat protein-10 (FBXL10) (fold-change 10(-3) P < 0.01), an ubiquitous nucleolar protein evolutionarily conserved. Hypoxanthine phosphoribosyltransferase 1 (HPRT1) and importin 8 (IPO8), were the best reference genes among the 16 genes in the LDAs with coefficient of variation (CV) of 4.51 and 4.55%, respectively. However, when the LDAs data were further analyzed by the geNorm and NormFinder softwares, IPO8, a nuclear protein mediating import of proteins, was the first and the third better reference gene, respectively. IPO8 and FBXL10 were further validated by real-time PCR in additional OM and SC fat samples and primary cultured preadipocytes. According to their CV, IPO8 resulted more suitable than FBXL10 in both adipose tissue depots and SC preadipocytes, whereas FBXL10 performed better than IPO8 in OM cultured preadipocytes. Both genes expression levels did not change throughout adipogenesis. Thus, we provide clear evidence that IPO8 and FBXL10 are good candidates to use as reference genes in gene expression studies in human OM and SC adipose tissues as well as differentiated primary preadipocytes.

Recurrent deletion of ZNF630 at Xp11.23 is not associated with mental retardation

ZNF630 is a member of the primate-specific Xp11 zinc finger gene cluster that consists of six closely related genes, of which ZNF41, ZNF81, and ZNF674 have been shown to be involved in mental retardation. This suggests that mutations of ZNF630 might influence cognitive function. Here, we detected 12 ZNF630 deletions in a total of 1,562 male patients with mental retardation from Brazil, USA, Australia, and Europe. The breakpoints were analyzed in 10 families, and in all cases they were located within two segmental duplications that share more than 99% sequence identity, indicating that the deletions resulted from non-allelic homologous recombination. In 2,121 healthy male controls, 10 ZNF630 deletions were identified. In total, there was a 1.6-fold higher frequency of this deletion in males with mental retardation as compared to controls, but this increase was not statistically significant (P-value = 0.174). Conversely, a 1.9-fold lower frequency of ZNF630 duplications was observed in patients, which was not significant either (P-value = 0.163). These data do not show that ZNF630 deletions or duplications are associated with mental retardation.

Gene expression analysis of calcineurin isoforms in T-lymphocytes--a method applied on kidney-transplant recipients

BACKGROUND

Tacrolimus exerts its immunosuppressive effect through inhibition of the intracellular enzyme calcineurin phosphatase (CaN). In this study, we set-up a validated real-time PCR method to measure the gene expression of the two major isoforms of the catalytic subunit of CaN in T-lymphocytes.

METHODS

20 stable kidney-transplant recipients, 10 early kidney-transplant recipients and 10 healthy non-medicated subjects had blood drawn and T-lymphocytes were isolated using E-rosette gradient centrifugation method. The cell counts were analyzed by DNA quantification using Hoeschst 33285. Gene expressions were analyzed using real-time PCR for CaN Aalpha, CaN Abeta and the reference genes CD3E and PPIB.

RESULTS

The real-time PCR method was found to be with high efficiencies and low intra- and inter-assay variabilities. No statistically significant differences were found in the gene expression levels of the two reference genes among the three groups. The two major isoforms of CaN A were expressed in equal amounts in the T-lymphocytes.

CONCLUSION

We found no significant difference in the reference genes between the three groups, but looking at the data there was a trend towards an up-regulation of CD3E. PPIB appears to be the more stable of the two reference genes tested in our study.

Toll-like receptor agonists synergistically increase proliferation and activation of B cells by epstein-barr virus

Epstein-Barr virus (EBV) efficiently drives proliferation of human primary B cells in vitro, a process relevant for human diseases such as infectious mononucleosis and posttransplant lymphoproliferative disease. Human B-cell proliferation is also driven by ligands of Toll-like receptors (TLRs), notably viral or bacterial DNA containing unmethylated CpG dinucleotides, which triggers TLR9. Here we quantitatively investigated how TLR stimuli influence EBV-driven B-cell proliferation and expression of effector molecules. CpG DNA synergistically increased EBV-driven proliferation and transformation, T-cell costimulatory molecules, and early production of interleukin-6. CpG DNA alone activated only memory B cells, but CpG DNA enhanced EBV-mediated transformation of both memory and naive B cells. Ligands for TLR2 or TLR7/8 or whole bacteria had a weaker but still superadditive effect on B-cell transformation. Additionally, CpG DNA facilitated the release of transforming virus by established EBV-infected lymphoblastoid cell lines. These results suggest that the proliferation of EBV-infected B cells and their capability to interact with immune effector cells may be directly influenced by components of bacteria or other microbes present at the site of infection.

OFD1 is mutated in X-linked Joubert syndrome and interacts with LCA5-encoded lebercilin

We ascertained a multi-generation Malaysian family with Joubert syndrome (JS). The presence of asymptomatic obligate carrier females suggested an X-linked recessive inheritance pattern. Affected males presented with mental retardation accompanied by postaxial polydactyly and retinitis pigmentosa. Brain MRIs showed the presence of a "molar tooth sign," which classifies this syndrome as classic JS with retinal involvement. Linkage analysis showed linkage to Xpter-Xp22.2 and a maximum LOD score of 2.06 for marker DXS8022. Mutation analysis revealed a frameshift mutation, p.K948NfsX8, in exon 21 of OFD1. In an isolated male with JS, a second frameshift mutation, p.E923KfsX3, in the same exon was identified. OFD1 has previously been associated with oral-facial-digital type 1 (OFD1) syndrome, a male-lethal X-linked dominant condition, and with X-linked recessive Simpson-Golabi-Behmel syndrome type 2 (SGBS2). In a yeast two-hybrid screen of a retinal cDNA library, we identified OFD1 as an interacting partner of the LCA5-encoded ciliary protein lebercilin. We show that X-linked recessive mutations in OFD1 reduce, but do not eliminate, the interaction with lebercilin, whereas X-linked dominant OFD1 mutations completely abolish binding to lebercilin. In addition, recessive mutations in OFD1 did not affect the pericentriolar localization of the recombinant protein in hTERT-RPE1 cells, whereas this localization was lost for dominant mutations. These findings offer a molecular explanation for the phenotypic spectrum observed for OFD1 mutations; this spectrum now includes OFD1 syndrome, SGBS2, and JS.

Regulation of MYCN expression in human neuroblastoma cells

Background

Amplification of the MYCN gene in neuroblastoma (NB) is associated with a poor prognosis. However, MYCN-amplification does not automatically result in higher expression of MYCN in children with NB. We hypothesized that the discrepancy between MYCN gene expression and prognosis in these children might be explained by the expression of either MYCN-opposite strand (MYCNOS) or the shortened MYCN-isoform (ΔMYCN) that was recently identified in fetal tissues. Both MYCNOS and ΔMYCN are potential inhibitors of MYCN either at the mRNA or at the protein level.

Methods

Expression of MYCN, MYCNOS and ΔMYCN was measured in human NB tissues of different stages. Transcript levels were quantified using a real-time reverse transcriptase polymerase chain reaction assay (QPCR). In addition, relative expression of these three transcripts was compared to the number of MYCN copies, which was determined by genomic real-time PCR (gQPCR).

Results

Both ΔMYCN and MYCNOS are expressed in all NBs examined. In NBs with MYCN-amplification, these transcripts are significantly higher expressed. The ratio of MYCN:ΔMYCN expression was identical in all tested NBs. This indicates that ΔMYCN and MYCN are co-regulated, which suggests that ΔMYCN is not a regulator of MYCN in NB. However, the ratio of MYCNOS:MYCN expression is directly correlated with NB disease stage (p = 0.007). In the more advanced NB stages and NBs with MYCN-amplification, relatively more MYCNOS is present as compared to MYCN. Expression of the antisense gene MYCNOS might be relevant to the progression of NB, potentially by directly inhibiting MYCN transcription by transcriptional interference at the DNA level.

Conclusion

The MYCNOS:MYCN-ratio in NBs is significantly correlated with both MYCN-amplification and NB-stage. Our data indicate that in NB, MYCN expression levels might be influenced by MYCNOS but not by ΔMYCN.

Interleukin 17 acts in synergy with B cell-activating factor to influence B cell biology and the pathophysiology of systemic lupus erythematosus

Studies have suggested involvement of interleukin 17 (IL-17) in autoimmune diseases, although its effect on B cell biology has not been clearly established. Here we demonstrate that IL-17 alone or in combination with B cell-activating factor controlled the survival and proliferation of human B cells and their differentiation into immunoglobulin-secreting cells. This effect was mediated mainly through the nuclear factor-kappaB-regulated transcription factor Twist-1. In support of the relevance of our observations and the potential involvement of IL-17 in B cell biology, we found that the serum of patients with systemic lupus erythematosus had higher concentrations of IL-17 than did the serum of healthy people and that IL-17 abundance correlated with the disease severity of systemic lupus erythematosus.

Characterization of two ectrodactyly-associated translocation breakpoints separated by 2.5 Mb on chromosome 2q14.1-q14.2

Split hand-split foot malformation or ectrodactyly is a heterogeneous congenital defect of digit formation. The aim of this study is the mapping of the breakpoints and a detailed molecular characterization of the candidate genes for an isolated and syndromic form of ectrodactyly, both associated with de novo apparently balanced chromosome translocations involving the same chromosome 2 band, [t(2;11)(q14.2;q14.2)] and [t(2;4)(q14.1;q35)], respectively. Breakpoints were mapped by fluorescence in situ hybridization using bacterial artificial chromosome clones. Where possible, these breakpoints were further delimited. Candidate genes were screened for pathogenic mutations and the expression levels of two of them analysed. The isolated bilateral split foot malformation-associated chromosome 2 breakpoint was localized at 120.9 Mb, between the two main candidate genes, encoding GLI-Kruppel family member GLI2 and inhibin-betaB. The second breakpoint associated with holoprosencephaly, hypertelorism and ectrodactyly syndrome was mapped 2.5 Mb proximal at 118.4 Mb and the candidate genes identified from this region were the insulin-induced protein 2 and the homeobox protein engrailed-1. No clear pathogenic mutations were identified in any of these genes. The breakpoint between INHBB and GLI2 coincides with a previously identified translocation breakpoint associated with ectrodactyly. We propose a mechanism by which translocations in the 2q14.1-q14.2 region disrupt the specific arrangement of long-range regulatory elements that control the tight quantitative spatiotemporal expression of one or more genes from the breakpoint region.

Genomic selection of reference genes for real-time PCR in human myocardium

Background

Reliability of real-time PCR (RT-qPCR) data is dependent on the use of appropriate reference gene(s) for normalization. To date, no validated reference genes have been reported for normalizing gene expression in human myocardium. This study aimed to identify validated reference genes for use in gene expression studies of failed and non-failed human myocardium.

Methods

Bioinformatic analysis of published human heart gene expression arrays (195 failed hearts, 16 donor hearts) was used to identify 10 stable and abundant genes for further testing. The expression stability of these genes was investigated in 28 failed and 28 non-failed human myocardium samples by RT-qPCR using geNorm software.

Results

Signal recognition particle 14 kDa (SRP14), tumor protein, translationally-controlled 1 (TPT1) and eukaryotic elongation factor 1A1 (EEF1A1) were ranked the most stable genes. The commonly used reference gene, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was ranked the least stable of the genes tested. The normalization strategy was tested by comparing RT-qPCR data of both normalized and raw expression levels of brain natriuretic peptide precursor (NPPB), a gene known to be up-regulated in heart failure. Non-normalized levels of NPPB exhibited a marginally significant difference between failed and non-failed samples (p = 0.058). In contrast, normalized NPPB expression levels were significantly higher in heart-failed patients compared with controls (p = 0.023).

Conclusion

This study used publicly available gene array data to identify a strategy for normalization involving two reference genes in combination that may have broad application for accurate and reliable normalization of RT-qPCR data in failed and non-failed human myocardium.

Reference genes for quantitative real time PCR in UVB irradiated keratinocytes

Real time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is a sensitive and highly reproducible method often used for determining mRNA levels. To enable proper comparison of gene expression genes expressed at stabile levels within the cells in the studied experimental system need to be identified and used as reference. Ultraviolet B (UVB) radiation is an exogenous carcinogenic stimulus in keratinocytes, and UVB elicited changes have extensively been studied by qRT-PCR, yet a comparison of commonly used reference genes in UVB treatment is lacking. To find the best genes for compensating slight inter-sample variations in keratinocytes in UVB experiments and to understand the potential effects of improper reference gene (RG) selection we have analyzed the mRNA expression of 10 housekeeping genes in neonatal human epidermal keratinocytes (NHEK) after UVB treatment. The biological effect of the used UVB light source was validated by trypane blue exclusion, MTT and comet assays. 20-40mJ/cm(2) dose was chosen for the experiments. The stability of the 10 RGs was assessed by the GeNorm and Normfinder software tools. Regardless of their slightly different algorithm the programs found succinate dehydrogenase complex subunit A (SDHA) to be the best individual RG and SDHA and phosphoglycerate kinase-1 (PGK1) as the most suitable combination. Analysis of the expression of tumor necrosis factor alpha (TNFalpha) and vascular endothelial growth factor (VEGF) found that while the perception of changes in TNF-alpha, a gene undergoing marked upregulation after UVB irradiation is independent of the used RG, changes seen in the more modestly upregulated VEGF are greatly effected by reference gene selection. These findings highlight the importance of reference gene selection in UVB irradiation experiments, and provide evidence that using SDHA or the combination of SDHA and PGK1 as standards could be a reliable method for normalizing qRT-PCR results in keratinocytes after UVB treatment.

Integrated genomic and expression profiling in mantle cell lymphoma: identification of gene-dosage regulated candidate genes

Mantle cell lymphoma (MCL) is characterized by the t(11;14)(q13;q32) translocation and several other cytogenetic aberrations, including heterozygous loss of chromosomal arms 1p, 6q, 11q and 13q and/or gains of 3q and 8q. The common intervals of chromosomal imbalance have been narrowed down using array-comparative genomic hybridization (CGH). However, the chromosomal intervals still contain many genes potentially involved in MCL pathogeny. Combined analysis of tiling-resolution array-CGH with gene expression profiling on 11 MCL tumours enabled the identification of genomic alterations and their corresponding gene expression profiles. Only subsets of genes located within given cytogenetic anomaly-intervals showed a concomitant change in mRNA expression level. The genes that showed consistent correlation between DNA copy number and RNA expression levels are likely to be important in MCL pathology. Besides several 'anonymous genes', we also identified various fully annotated genes, whose gene products are involved in cyclic adenosine monophosphate-regulated pathways (PRKACB), DNA damage repair, maintenance of chromosome stability and prevention of rereplication (ATM, ERCC5, FBXO5), energy metabolism (such as genes that are involved in the synthesis of proteins encoded by the mitochondrial genome) and signal transduction (ARHGAP29). Deregulation of these gene products may interfere with the signalling pathways that are involved in MCL tumour development and maintenance.

Reference gene selection for real-time rtPCR in human epidermal keratinocytes

BACKGROUND

RNA interference represents a powerful tool with which to achieve suppression of specific target mRNA. Real-time rtPCR is a useful technique for assessing levels of mRNA expression. Critically, for real-time rtPCR to yield meaningful results, it is necessary to normalise expression of the gene of interest to stably expressed endogenous control genes.

OBJECTIVES

The study involved establishing expression profiles for 11 housekeeping genes in human epidermal keratinocyte cell lines determining their relative stability. Furthermore, the effect of the presence of shRNA on these expression profiles has been established.

METHODS

Keratinocytes were transfected using lipid-based transfection or AMAXA nucleofection. Real-time rtPCR was used to establish RNA expression levels. Data analysis was carried out using geNORM.

RESULTS

When using HaCaT or adult NHEK cells any combination of 8 of the housekeeping genes would be appropriate for normalisation. In contrast, with juvenile NHEK cells only 4 of the housekeeping genes were found to be sufficiently stable to be deemed appropriate for normalisation of expression data. Furthermore data demonstrated that the expression of housekeeping genes may sometimes be affected by the induction of the RNAi pathway.

CONCLUSIONS

The data obtained highlight the importance of characterising housekeeping gene expression profiles in each specific cell type prior to choosing a set of housekeeping genes for expression studies. The results from this study are applicable to researchers working with human epidermal keratinocytes and the experimental approach is more broadly applicable to any researcher carrying out real-time rtPCR.

Arts syndrome is caused by loss-of-function mutations in PRPS1

Arts syndrome is an X-linked disorder characterized by mental retardation, early-onset hypotonia, ataxia, delayed motor development, hearing impairment, and optic atrophy. Linkage analysis in a Dutch family and an Australian family suggested that the candidate gene maps to Xq22.1-q24. Oligonucleotide microarray expression profiling of fibroblasts from two probands of the Dutch family revealed reduced expression levels of the phosphoribosyl pyrophosphate synthetase 1 gene (PRPS1). Subsequent sequencing of PRPS1 led to the identification of two different missense mutations, c.455T-->C (p.L152P) in the Dutch family and c.398A-->C (p.Q133P) in the Australian family. Both mutations result in a loss of phosphoribosyl pyrophosphate synthetase 1 activity, as was shown in silico by molecular modeling and was shown in vitro by phosphoribosyl pyrophosphate synthetase activity assays in erythrocytes and fibroblasts from patients. This is in contrast to the gain-of-function mutations in PRPS1 that were identified previously in PRPS-related gout. The loss-of-function mutations of PRPS1 likely result in impaired purine biosynthesis, which is supported by the undetectable hypoxanthine in urine and the reduced uric acid levels in serum from patients. To replenish low levels of purines, treatment with S-adenosylmethionine theoretically could have therapeutic efficacy, and a clinical trial involving the two affected Australian brothers is currently underway.

Cis sequence effects on gene expression

Background

Sequence and transcriptional variability within and between individuals are typically studied independently. The joint analysis of sequence and gene expression variation (genetical genomics) provides insight into the role of linked sequence variation in the regulation of gene expression. We investigated the role of sequence variation in cis on gene expression (cis sequence effects) in a group of genes commonly studied in cancer research in lymphoblastoid cell lines. We estimated the proportion of genes exhibiting cis sequence effects and the proportion of gene expression variation explained by cis sequence effects using three different analytical approaches, and compared our results to the literature.

Results

We generated gene expression profiling data at N = 697 candidate genes from N = 30 lymphoblastoid cell lines for this study and used available candidate gene resequencing data at N = 552 candidate genes to identify N = 30 candidate genes with sufficient variance in both datasets for the investigation of cis sequence effects. We used two additive models and the haplotype phylogeny scanning approach of Templeton (Tree Scanning) to evaluate association between individual SNPs, all SNPs at a gene, and diplotypes, with log-transformed gene expression. SNPs and diplotypes at eight candidate genes exhibited statistically significant (p < 0.05) association with gene expression. Using the literature as a "gold standard" to compare 14 genes with data from both this study and the literature, we observed 80% and 85% concordance for genes exhibiting and not exhibiting significant cis sequence effects in our study, respectively.

Conclusion

Based on analysis of our results and the extant literature, one in four genes exhibits significant cis sequence effects, and for these genes, about 30% of gene expression variation is accounted for by cis sequence variation. Despite diverse experimental approaches, the presence or absence of significant cis sequence effects is largely supported by previously published studies.

Identification of new reference genes for the normalisation of canine osteoarthritic joint tissue transcripts from microarray data

BACKGROUND

Real-time reverse transcriptase quantitative polymerase chain reaction (real-time RT-qPCR) is the most accurate measure of gene expression in biological systems. The comparison of different samples requires the transformation of data through a process called normalisation. Reference or housekeeping genes are candidate genes which are selected on the basis of constitutive expression across samples, and allow the quantification of changes in gene expression. At present, no reference gene has been identified for any organism which is universally optimal for use across different tissue types or disease situations. We used microarray data to identify new reference genes generated from total RNA isolated from normal and osteoarthritic canine articular tissues (bone, ligament, cartilage, synovium and fat). RT-qPCR assays were designed and applied to each different articular tissue. Reference gene expression stability and ranking was compared using three different mathematical algorithms.

RESULTS

Twelve new potential reference genes were identified from microarray data. One gene (mitochondrial ribosomal protein S7 [MRPS7]) was stably expressed in all five of the articular tissues evaluated. One gene HIRA interacting protein 5 isoform 2 [HIRP5]) was stably expressed in four of the tissues evaluated. A commonly used reference gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was not stably expressed in any of the tissues evaluated. Most consistent agreement between rank ordering of reference genes was observed between Bestkeeper(c) and geNorm, although each method tended to agree on the identity of the most stably expressed genes and the least stably expressed genes for each tissue. New reference genes identified using microarray data normalised in a conventional manner were more stable than those identified by microarray data normalised by using a real-time RT-qPCR methodology.

CONCLUSION

Microarray data normalised by a conventional manner can be filtered using a simple stepwise procedure to identify new reference genes, some of which will demonstrate good measures of stability. Mitochondrial ribosomal protein S7 is a new reference gene worthy of investigation in other canine tissues and diseases. Different methods of reference gene stability assessment will generally agree on the most and least stably expressed genes, when co-regulation is not present.

Expression stability of commonly used reference genes in canine articular connective tissues

Background

The quantification of gene expression in tissue samples requires the use of reference genes to normalise transcript numbers between different samples. Reference gene stability may vary between different tissues, and between the same tissue in different disease states. We evaluated the stability of 9 reference genes commonly used in human gene expression studies. Real-time reverse transcription PCR and a mathematical algorithm were used to establish which reference genes were most stably expressed in normal and diseased canine articular tissues and two canine cell lines stimulated with lipolysaccaride (LPS).

Results

The optimal reference genes for comparing gene expression data between normal and diseased infrapatella fat pad were RPL13A and YWHAZ (M = 0.56). The ideal reference genes for comparing normal and osteoarthritic (OA) cartilage were RPL13A and SDHA (M = 0.57). The best reference genes for comparing normal and ruptured canine cranial cruciate ligament were B2M and TBP (M = 0.59). The best reference genes for normalising gene expression data from normal and LPS stimulated cell lines were SDHA and YWHAZ (K6) or SDHA and HMBS (DH82), which had expression stability (M) values of 0.05 (K6) and 0.07 (DH82) respectively. The number of reference genes required to reduce pairwise variation (V) to <0.20 was 4 for cell lines, 5 for cartilage, 7 for cranial cruciate ligament and 8 for fat tissue. Reference gene stability was not related to the level of gene expression.

Conclusion

The reference genes demonstrating the most stable expression within each different canine articular tissue were identified, but no single reference gene was identified as having stable expression in all different tissue types. This study underlines the necessity to select reference genes on the basis of tissue and disease specific expression profile evaluation and highlights the requirement for the identification of new reference genes with greater expression stability for use in canine articular tissue gene expression studies.