Induction of kidney-related gene programs through co-option of SALL1 in mole ovotestes

Changes in gene expression represent an important source of phenotypic innovation. Yet how such changes emerge and impact the evolution of traits remains elusive. Here, we explore the molecular mechanisms associated with the development of masculinizing ovotestes in female moles. By performing integrative analyses of epigenetic and transcriptional data in mole and mouse, we identified the co-option of SALL1 expression in mole ovotestes formation. Chromosome conformation capture analyses highlight a striking conservation of the 3D organization at the SALL1 locus, but an evolutionary divergence of enhancer activity. Interspecies reporter assays support the capability of mole-specific enhancers to activate transcription in urogenital tissues. Through overexpression experiments in transgenic mice, we further demonstrate the capability of SALL1 to induce kidney-related gene programs, which are a signature of mole ovotestes. Our results highlight the co-option of gene expression, through changes in enhancer activity, as a plausible mechanism for the evolution of traits.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

The mole genome reveals regulatory rearrangements associated with adaptive intersexuality

Linking genomic variation to phenotypical traits remains a major challenge in evolutionary genetics. In this study, we use phylogenomic strategies to investigate a distinctive trait among mammals: the development of masculinizing ovotestes in female moles. By combining a chromosome-scale genome assembly of the Iberian mole, Talpa occidentalis, with transcriptomic, epigenetic, and chromatin interaction datasets, we identify rearrangements altering the regulatory landscape of genes with distinct gonadal expression patterns. These include a tandem triplication involving CYP17A1, a gene controlling androgen synthesis, and an intrachromosomal inversion involving the pro-testicular growth factor gene FGF9, which is heterochronically expressed in mole ovotestes. Transgenic mice with a knock-in mole CYP17A1 enhancer or overexpressing FGF9 showed phenotypes recapitulating mole sexual features. Our results highlight how integrative genomic approaches can reveal the phenotypic impact of noncoding sequence changes.

International survey of neurosurgical anesthesia (iSonata) : An international survey of current practices in neurosurgical anesthesia

BACKGROUND

No standardized recommendations have been currently defined for anesthesia management of patients undergoing elective intracranial surgery. It can therefore be assumed that international clinical institutions have diverging approaches or standard operating procedures (SOP) which determine the type of general anesthesia, hemodynamic management, neuromuscular blockade, implementation of hypothermia and postoperative patient care.

OBJECTIVE

This international survey aimed to assess perioperative patient management during elective intracranial procedures. This survey was performed from February to October 2018 and 311 neurosurgical, maximum care centers across 19 European countries were contacted. The aim was to evaluate the anesthesia management to provide relevant data of neuroanesthesia practices across European centers. The survey differentiated between vascular and non-vascular as well as supratentorial and infratentorial procedures.

RESULTS

A total of 109 (35.0%) completed questionnaires from 15 European countries were analyzed. The results illustrated that total intravenous anesthesia was most commonly implemented during elective intracranial procedures (83.8%). All centers performed endotracheal intubation prior to major intracranial surgery (100%). Central venous lines were placed in 63.3% of cases. Moderate intraoperative hypothermia was carried out in 12.8% of the procedures, especially during vascular supratentorial and infratentorial surgery. A neuromuscular blockade during surgery was implemented in 74.1% of patients. Assessment of the neuromuscular junction was performed in 59.2% of cases, 76.7% of patients were immediately extubated in the operating room. 84.7% of these patients were directly transferred to a monitoring ward or an intensive care unit (ICU) and 55.1% of ventilated patients were transferred directly to an ICU.

CONCLUSION

The data demonstrate that many aspects of anesthesia management during elective intracranial surgery vary between European institutions. The data also suggest that a broad consensus exists regarding the implementation of total intravenous anesthesia, airway management (endotracheal intubation), the implementation of urinary catheters, large bore peripheral venous lines and the broad availability of cross-matched red blood cell concentrates. Nevertheless, anesthesia management (e.g. central venous catheterization, moderate hypothermia, neuromuscular monitoring) is still handled differently across many European institutions. A lack of standardized guidelines defining anesthetic management in patients undergoing intracranial procedures could explain this variability. Further studies could help establish optimal anesthesia management for these patients. This in turn could help in the development of national and international guidelines and SOPs which could define optimal management strategies for intracranial procedures.

Pathogenic variants in E3 ubiquitin ligase RLIM/RNF12 lead to a syndromic X-linked intellectual disability and behavior disorder

RLIM, also known as RNF12, is an X-linked E3 ubiquitin ligase acting as a negative regulator of LIM-domain containing transcription factors and participates in X-chromosome inactivation (XCI) in mice. We report the genetic and clinical findings of 84 individuals from nine unrelated families, eight of whom who have pathogenic variants in RLIM (RING finger LIM domain-interacting protein). A total of 40 affected males have X-linked intellectual disability (XLID) and variable behavioral anomalies with or without congenital malformations. In contrast, 44 heterozygous female carriers have normal cognition and behavior, but eight showed mild physical features. All RLIM variants identified are missense changes co-segregating with the phenotype and predicted to affect protein function. Eight of the nine altered amino acids are conserved and lie either within a domain essential for binding interacting proteins or in the C-terminal RING finger catalytic domain. In vitro experiments revealed that these amino acid changes in the RLIM RING finger impaired RLIM ubiquitin ligase activity. In vivo experiments in rlim mutant zebrafish showed that wild type RLIM rescued the zebrafish rlim phenotype, whereas the patient-specific missense RLIM variants failed to rescue the phenotype and thus represent likely severe loss-of-function mutations. In summary, we identified a spectrum of RLIM missense variants causing syndromic XLID and affecting the ubiquitin ligase activity of RLIM, suggesting that enzymatic activity of RLIM is required for normal development, cognition and behavior.

Mutation p.R356Q in the Collybistin Phosphoinositide Binding Site Is Associated With Mild Intellectual Disability

The recruitment of inhibitory GABA_A receptors to neuronal synapses requires a complex interplay between receptors, neuroligins, the scaffolding protein gephyrin and the GDP-GTP exchange factor collybistin (CB). Collybistin is regulated by protein-protein interactions at the N-terminal SH3 domain, which can bind neuroligins 2/4 and the GABA_AR α2 subunit. Collybistin also harbors a RhoGEF domain which mediates interactions with gephyrin and catalyzes GDP-GTP exchange on Cdc42. Lastly, collybistin has a pleckstrin homology (PH) domain, which binds phosphoinositides, such as phosphatidylinositol 3-phosphate (PI3P/PtdIns3P) and phosphatidylinositol 4-monophosphate (PI4P/PtdIns4P). PI3P located in early/sorting endosomes has recently been shown to regulate the postsynaptic clustering of gephyrin and GABA_A receptors and consequently the strength of inhibitory synapses in cultured hippocampal neurons. This process is disrupted by mutations in the collybistin gene (ARHGEF9), which cause X-linked intellectual disability (XLID) by a variety of mechanisms converging on disrupted gephyrin and GABA_A receptor clustering at central synapses. Here we report a novel missense mutation (chrX:62875607C>T, p.R356Q) in ARHGEF9 that affects one of the two paired arginine residues in the PH domain that were predicted to be vital for binding phosphoinositides. Functional assays revealed that recombinant collybistin CB3_SH3- ^R356Q was deficient in PI3P binding and was not able to translocate EGFP-gephyrin to submembrane microaggregates in an in vitro clustering assay. Expression of the PI3P-binding mutants CB3_SH3- ^R356Q and CB3_SH3- ^R356N/R357N in cultured hippocampal neurones revealed that the mutant proteins did not accumulate at inhibitory synapses, but instead resulted in a clear decrease in the overall number of synaptic gephyrin clusters compared to controls. Molecular dynamics simulations suggest that the p.R356Q substitution influences PI3P binding by altering the range of structural conformations adopted by collybistin. Taken together, these results suggest that the p.R356Q mutation in ARHGEF9 is the underlying cause of XLID in the probands, disrupting gephyrin clustering at inhibitory GABAergic synapses via loss of collybistin PH domain phosphoinositide binding.

Rare GABRA3 variants are associated with epileptic seizures, encephalopathy and dysmorphic features

Genetic epilepsies are caused by mutations in a range of different genes, many of them encoding ion channels, receptors or transporters. While the number of detected variants and genes increased dramatically in the recent years, pleiotropic effects have also been recognized, revealing that clinical syndromes with various degrees of severity arise from a single gene, a single mutation, or from different mutations showing similar functional defects. Accordingly, several genes coding for GABAA receptor subunits have been linked to a spectrum of benign to severe epileptic disorders and it was shown that a loss of function presents the major correlated pathomechanism. Here, we identified six variants in GABRA3 encoding the α3-subunit of the GABAA receptor. This gene is located on chromosome Xq28 and has not been previously associated with human disease. Five missense variants and one microduplication were detected in four families and two sporadic cases presenting with a range of epileptic seizure types, a varying degree of intellectual disability and developmental delay, sometimes with dysmorphic features or nystagmus. The variants co-segregated mostly but not completely with the phenotype in the families, indicating in some cases incomplete penetrance, involvement of other genes, or presence of phenocopies. Overall, males were more severely affected and there were three asymptomatic female mutation carriers compared to only one male without a clinical phenotype. X-chromosome inactivation studies could not explain the phenotypic variability in females. Three detected missense variants are localized in the extracellular GABA-binding NH2-terminus, one in the M2-M3 linker and one in the M4 transmembrane segment of the α3-subunit. Functional studies in Xenopus laevis oocytes revealed a variable but significant reduction of GABA-evoked anion currents for all mutants compared to wild-type receptors. The degree of current reduction correlated partially with the phenotype. The microduplication disrupted GABRA3 expression in fibroblasts of the affected patient. In summary, our results reveal that rare loss-of-function variants in GABRA3 increase the risk for a varying combination of epilepsy, intellectual disability/developmental delay and dysmorphic features, presenting in some pedigrees with an X-linked inheritance pattern.

EIF2S3 Mutations Associated with Severe X-Linked Intellectual Disability Syndrome MEHMO

Impairment of translation initiation and its regulation within the integrated stress response (ISR) and related unfolded-protein response has been identified as a cause of several multisystemic syndromes. Here, we link MEHMO syndrome, whose genetic etiology was unknown, to this group of disorders. MEHMO is a rare X-linked syndrome characterized by profound intellectual disability, epilepsy, hypogonadism and hypogenitalism, microcephaly, and obesity. We have identified a C-terminal frameshift mutation (Ile465Serfs) in the EIF2S3 gene in three families with MEHMO syndrome and a novel maternally inherited missense EIF2S3 variant (c.324T>A; p.Ser108Arg) in another male patient with less severe clinical symptoms. The EIF2S3 gene encodes the γ subunit of eukaryotic translation initiation factor 2 (eIF2), crucial for initiation of protein synthesis and regulation of the ISR. Studies in patient fibroblasts confirm increased ISR activation due to the Ile465Serfs mutation and functional assays in yeast demonstrate that the Ile465Serfs mutation impairs eIF2γ function to a greater extent than tested missense mutations, consistent with the more severe clinical phenotype of the Ile465Serfs male mutation carriers. Thus, we propose that more severe EIF2S3 mutations cause the full MEHMO phenotype, while less deleterious mutations cause a milder form of the syndrome with only a subset of the symptoms.

Is applanation tonometry a reliable method for monitoring blood pressure in morbidly obese patients undergoing bariatric surgery?

BACKGROUND

The aim of this study was to evaluate the validity of non-invasive continuous BP measurement by applanation tonometry in morbidly obese patients undergoing bariatric surgery.

METHODS

Arterial blood pressure (AP) was recorded intraoperatively both by applanation tonometry (AT) (T-Line 200pro, Tensys Medical(®), USA) and an arterial line (AL) after radial cannulation in obese patients undergoing bariatric surgery. Discrepancies between the two methods were assessed as bias, limits of agreement and percentage error. Mean, systolic, and diastolic arterial pressures were assessed (MAP, SAP, DAP respectively). Trending ability was assessed by concordance based on four-quadrant plotting.

RESULTS

Mean (sd) BMI of the 28 patients was 49.4 (9.7 kg m(-2)). A total of 201 907 time points were available for analysis. Bias for MAPAT compared with MAPAL was +3.97 mm Hg (SAPAT +3.45 mm Hg; DAPAT +3.66 mm Hg) with limits of agreement for MAPAT of -14.47 and +22.41 mm Hg (SAPAT -22.0 and +28.9 mm Hg; DAPAT -15.7 and +23.1 mm Hg). Percentage error for MAPAT was 23.5% (23.4% for SAPAT; 30.5% for DAPAT). Trending ability for MAP, SAP, and DAP revealed a concordance of 0.74, 0.72, and 0.71, respectively.

CONCLUSIONS

Continuous BP assessment by applanation tonometry is feasible in morbidly obese patients undergoing bariatric surgery. However, despite a low mean difference, 95% limits of agreement and trending ability indicate that the technology needs to be improved further, before being recommended for routine use in this group of patients.

Novel Missense Mutation A789V in IQSEC2 Underlies X-Linked Intellectual Disability in the MRX78 Family

Disease gene discovery in neurodevelopmental disorders, including X-linked intellectual disability (XLID) has recently been accelerated by next-generation DNA sequencing approaches. To date, more than 100 human X chromosome genes involved in neuronal signaling pathways and networks implicated in cognitive function have been identified. Despite these advances, the mutations underlying disease in a large number of XLID families remained unresolved. We report the resolution of MRX78, a large family with six affected males and seven affected females, showing X-linked inheritance. Although a previous linkage study had mapped the locus to the short arm of chromosome X (Xp11.4-p11.23), this region contained too many candidate genes to be analyzed using conventional approaches. However, our X-chromosome exome resequencing, bioinformatics analysis and inheritance testing revealed a missense mutation (c.C2366T, p.A789V) in IQSEC2, encoding a neuronal GDP-GTP exchange factor for Arf family GTPases (ArfGEF) previously implicated in XLID. Molecular modeling of IQSEC2 revealed that the A789V substitution results in the insertion of a larger side-chain into a hydrophobic pocket in the catalytic Sec7 domain of IQSEC2. The A789V change is predicted to result in numerous clashes with adjacent amino acids and disruption of local folding of the Sec7 domain. Consistent with this finding, functional assays revealed that recombinant IQSEC2^A789V was not able to catalyze GDP-GTP exchange on Arf6 as efficiently as wild-type IQSEC2. Taken together, these results strongly suggest that the A789V mutation in IQSEC2 is the underlying cause of XLID in the MRX78 family.

X-exome sequencing of 405 unresolved families identifies seven novel intellectual disability genes

X-linked intellectual disability (XLID) is a clinically and genetically heterogeneous disorder. During the past two decades in excess of 100 X-chromosome ID genes have been identified. Yet, a large number of families mapping to the X-chromosome remained unresolved suggesting that more XLID genes or loci are yet to be identified. Here, we have investigated 405 unresolved families with XLID. We employed massively parallel sequencing of all X-chromosome exons in the index males. The majority of these males were previously tested negative for copy number variations and for mutations in a subset of known XLID genes by Sanger sequencing. In total, 745 X-chromosomal genes were screened. After stringent filtering, a total of 1297 non-recurrent exonic variants remained for prioritization. Co-segregation analysis of potential clinically relevant changes revealed that 80 families (20%) carried pathogenic variants in established XLID genes. In 19 families, we detected likely causative protein truncating and missense variants in 7 novel and validated XLID genes (CLCN4, CNKSR2, FRMPD4, KLHL15, LAS1L, RLIM and USP27X) and potentially deleterious variants in 2 novel candidate XLID genes (CDK16 and TAF1). We show that the CLCN4 and CNKSR2 variants impair protein functions as indicated by electrophysiological studies and altered differentiation of cultured primary neurons from Clcn4(-/-) mice or after mRNA knock-down. The newly identified and candidate XLID proteins belong to pathways and networks with established roles in cognitive function and intellectual disability in particular. We suggest that systematic sequencing of all X-chromosomal genes in a cohort of patients with genetic evidence for X-chromosome locus involvement may resolve up to 58% of Fragile X-negative cases.

Increased STAG2 dosage defines a novel cohesinopathy with intellectual disability and behavioral problems

Next generation genomic technologies have made a significant contribution to the understanding of the genetic architecture of human neurodevelopmental disorders. Copy number variants (CNVs) play an important role in the genetics of intellectual disability (ID). For many CNVs, and copy number gains in particular, the responsible dosage-sensitive gene(s) have been hard to identify. We have collected 18 different interstitial microduplications and 1 microtriplication of Xq25. There were 15 affected individuals from 6 different families and 13 singleton cases, 28 affected males in total. The critical overlapping region involved the STAG2 gene, which codes for a subunit of the cohesin complex that regulates cohesion of sister chromatids and gene transcription. We demonstrate that STAG2 is the dosage-sensitive gene within these CNVs, as gains of STAG2 mRNA and protein dysregulate disease-relevant neuronal gene networks in cells derived from affected individuals. We also show that STAG2 gains result in increased expression of OPHN1, a known X-chromosome ID gene. Overall, we define a novel cohesinopathy due to copy number gain of Xq25 and STAG2 in particular.

A Novel Mutation in RPL10 (Ribosomal Protein L10) Causes X-Linked Intellectual Disability, Cerebellar Hypoplasia, and Spondylo-Epiphyseal Dysplasia

RPL10 encodes ribosomal protein L10 (uL16), a highly conserved multifunctional component of the large ribosomal subunit, involved in ribosome biogenesis and function. Using X-exome resequencing, we identified a novel missense mutation (c.191C>T; p.(A64V)) in the N-terminal domain of the protein, in a family with two affected cousins presenting with X-linked intellectual disability, cerebellar hypoplasia, and spondylo-epiphyseal dysplasia (SED). We assessed the impact of the mutation on the translational capacity of the cell using yeast as model system. The mutation generates a functional ribosomal protein, able to complement the translational defects of a conditional lethal mutation of yeast rpl10. However, unlike previously reported mutations, this novel RPL10 missense mutation results in an increase in the actively translating ribosome population. Our results expand the mutational and clinical spectrum of RPL10 identifying a new genetic cause of SED and highlight the emerging role of ribosomal proteins in the pathogenesis of neurodevelopmental disorders.

Variants in CUL4B are associated with cerebral malformations

Variants in cullin 4B (CUL4B) are a known cause of syndromic X-linked intellectual disability. Here, we describe an additional 25 patients from 11 families with variants in CUL4B. We identified nine different novel variants in these families and confirmed the pathogenicity of all nontruncating variants. Neuroimaging data, available for 15 patients, showed the presence of cerebral malformations in ten patients. The cerebral anomalies comprised malformations of cortical development (MCD), ventriculomegaly, and diminished white matter volume. The phenotypic heterogeneity of the cerebral malformations might result from the involvement of CUL-4B in various cellular pathways essential for normal brain development. Accordingly, we show that CUL-4B interacts with WDR62, a protein in which variants were previously identified in patients with microcephaly and a wide range of MCD. This interaction might contribute to the development of cerebral malformations in patients with variants in CUL4B.

THOC2 Mutations Implicate mRNA-Export Pathway in X-Linked Intellectual Disability

Export of mRNA from the cell nucleus to the cytoplasm is essential for protein synthesis, a process vital to all living eukaryotic cells. mRNA export is highly conserved and ubiquitous. Mutations affecting mRNA and mRNA processing or export factors, which cause aberrant retention of mRNAs in the nucleus, are thus emerging as contributors to an important class of human genetic disorders. Here, we report that variants in THOC2, which encodes a subunit of the highly conserved TREX mRNA-export complex, cause syndromic intellectual disability (ID). Affected individuals presented with variable degrees of ID and commonly observed features included speech delay, elevated BMI, short stature, seizure disorders, gait disturbance, and tremors. X chromosome exome sequencing revealed four missense variants in THOC2 in four families, including family MRX12, first ascertained in 1971. We show that two variants lead to decreased stability of THOC2 and its TREX-complex partners in cells derived from the affected individuals. Protein structural modeling showed that the altered amino acids are located in the RNA-binding domains of two complex THOC2 structures, potentially representing two different intermediate RNA-binding states of THOC2 during RNA transport. Our results show that disturbance of the canonical molecular pathway of mRNA export is compatible with life but results in altered neuronal development with other comorbidities.

Comprehensive genomic profiles of small cell lung cancer

We have sequenced the genomes of 110 small cell lung cancers (SCLC), one of the deadliest human cancers. In nearly all the tumours analysed we found bi-allelic inactivation of TP53 and RB1, sometimes by complex genomic rearrangements. Two tumours with wild-type RB1 had evidence of chromothripsis leading to overexpression of cyclin D1 (encoded by the CCND1 gene), revealing an alternative mechanism of Rb1 deregulation. Thus, loss of the tumour suppressors TP53 and RB1 is obligatory in SCLC. We discovered somatic genomic rearrangements of TP73 that create an oncogenic version of this gene, TP73Δex2/3. In rare cases, SCLC tumours exhibited kinase gene mutations, providing a possible therapeutic opportunity for individual patients. Finally, we observed inactivating mutations in NOTCH family genes in 25% of human SCLC. Accordingly, activation of Notch signalling in a pre-clinical SCLC mouse model strikingly reduced the number of tumours and extended the survival of the mutant mice. Furthermore, neuroendocrine gene expression was abrogated by Notch activity in SCLC cells. This first comprehensive study of somatic genome alterations in SCLC uncovers several key biological processes and identifies candidate therapeutic targets in this highly lethal form of cancer.

[Goal-directed hemodynamic therapy: Concepts, indications and risks]

Goal-directed hemodynamic therapy is becoming increasingly more interesting for anesthesiologists and intensive care physicians. Meta-analyses of studies evaluating perioperative therapy algorithms demonstrated a reduction of postoperative morbidity compared to the previous clinical practices. In this review article the basic concepts of goal-directed hemodynamic therapy and the principles of previously employed therapy algorithms are described and discussed. Furthermore, the questions of how these therapy strategies can be transferred into daily clinical practice and whether these therapeutic approaches might even bear risks for patients are elucidated.

Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions

Mammalian genomes are organized into megabase-scale topologically associated domains (TADs). We demonstrate that disruption of TADs can rewire long-range regulatory architecture and result in pathogenic phenotypes. We show that distinct human limb malformations are caused by deletions, inversions, or duplications altering the structure of the TAD-spanning WNT6/IHH/EPHA4/PAX3 locus. Using CRISPR/Cas genome editing, we generated mice with corresponding rearrangements. Both in mouse limb tissue and patient-derived fibroblasts, disease-relevant structural changes cause ectopic interactions between promoters and non-coding DNA, and a cluster of limb enhancers normally associated with Epha4 is misplaced relative to TAD boundaries and drives ectopic limb expression of another gene in the locus. This rewiring occurred only if the variant disrupted a CTCF-associated boundary domain. Our results demonstrate the functional importance of TADs for orchestrating gene expression via genome architecture and indicate criteria for predicting the pathogenicity of human structural variants, particularly in non-coding regions of the human genome.

New insights into Brunner syndrome and potential for targeted therapy

We report two families with Brunner syndrome living in one state of Australia. The first family had a predicted protein-truncating variant of monoamine oxidase A (MAOA) (p.S251KfsX2). Affected males had mild intellectual disability (ID), obsessive behaviour, limited friendships and were introverted and placid during clinical interview. The family disclosed episodic explosive aggression after a diagnosis was made. The second family had a missense variant in MAOA (p.R45W). Affected males had borderline-mild ID, attention deficit disorder and limited friendships. One had a history of explosive aggression in childhood and episodic symptoms of flushing, headaches and diarrhoea. Their carrier mother had normal intelligence but similar episodic symptoms. Characteristic biochemical abnormalities included high serum serotonin and urinary metanephrines and low urinary 5-hydroxyindoleacetic acid (5-HIAA) and vanillylmandelic acid (VMA). Symptomatic individuals in the second family had particularly high serotonin levels, and treatment with a serotonin reuptake inhibitor and dietary modification resulted in reversal of biochemical abnormalities, reduction of 'serotonergic' symptoms and behavioural improvement. Brunner syndrome should be considered as a cause of mild ID with paroxysmal behavioural symptoms. It can be screened for with serum/urine metanephrine and serotonin measurement. Cautious treatment with a serotonin reuptake inhibitor, dietary modifications and avoidance of medications contraindicated in patients on monoamine oxidase inhibitors can improve symptoms.

Identification of novel fusion genes in lung cancer using breakpoint assembly of transcriptome sequencing data

Genomic translocation events frequently underlie cancer development through generation of gene fusions with oncogenic properties. Identification of such fusion transcripts by transcriptome sequencing might help to discover new potential therapeutic targets. We developed TRUP (Tumor-specimen suited RNA-seq Unified Pipeline) (https://github.com/ruping/TRUP), a computational approach that combines split-read and read-pair analysis with de novo assembly for the identification of chimeric transcripts in cancer specimens. We apply TRUP to RNA-seq data of different tumor types, and find it to be more sensitive than alternative tools in detecting chimeric transcripts, such as secondary rearrangements in EML4-ALK-positive lung tumors, or recurrent inactivating rearrangements affecting RASSF8.

Mutations in RAB39B cause X-linked intellectual disability and early-onset Parkinson disease with α-synuclein pathology

Advances in understanding the etiology of Parkinson disease have been driven by the identification of causative mutations in families. Genetic analysis of an Australian family with three males displaying clinical features of early-onset parkinsonism and intellectual disability identified a ∼45 kb deletion resulting in the complete loss of RAB39B. We subsequently identified a missense mutation (c.503C>A [p.Thr168Lys]) in RAB39B in an unrelated Wisconsin kindred affected by a similar clinical phenotype. In silico and in vitro studies demonstrated that the mutation destabilized the protein, consistent with loss of function. In vitro small-hairpin-RNA-mediated knockdown of Rab39b resulted in a reduction in the density of α-synuclein immunoreactive puncta in dendritic processes of cultured neurons. In addition, in multiple cell models, we demonstrated that knockdown of Rab39b was associated with reduced steady-state levels of α-synuclein. Post mortem studies demonstrated that loss of RAB39B resulted in pathologically confirmed Parkinson disease. There was extensive dopaminergic neuron loss in the substantia nigra and widespread classic Lewy body pathology. Additional pathological features included cortical Lewy bodies, brain iron accumulation, tau immunoreactivity, and axonal spheroids. Overall, we have shown that loss-of-function mutations in RAB39B cause intellectual disability and pathologically confirmed early-onset Parkinson disease. The loss of RAB39B results in dysregulation of α-synuclein homeostasis and a spectrum of neuropathological features that implicate RAB39B in the pathogenesis of Parkinson disease and potentially other neurodegenerative disorders.

Involvement of the kinesin family members KIF4A and KIF5C in intellectual disability and synaptic function

INTRODUCTION

Kinesin superfamily (KIF) genes encode motor proteins that have fundamental roles in brain functioning, development, survival and plasticity by regulating the transport of cargo along microtubules within axons, dendrites and synapses. Mouse knockout studies support these important functions in the nervous system. The role of KIF genes in intellectual disability (ID) has so far received limited attention, although previous studies have suggested that many ID genes impinge on synaptic function.

METHODS

By applying next-generation sequencing (NGS) in ID patients, we identified likely pathogenic mutations in KIF4A and KIF5C. To further confirm the pathogenicity of these mutations, we performed functional studies at the level of synaptic function in primary rat hippocampal neurons.

RESULTS AND CONCLUSIONS

Four males from a single family with a disruptive mutation in the X-linked KIF4A (c.1489-8_1490delins10; p.?- exon skipping) showed mild to moderate ID and epilepsy. A female patient with a de novo missense mutation in KIF5C (c.11465A>C; p.(Glu237Lys)) presented with severe ID, epilepsy, microcephaly and cortical malformation. Knock-down of Kif4a in rat primary hippocampal neurons altered the balance between excitatory and inhibitory synaptic transmission, whereas the mutation in Kif5c affected its protein function at excitatory synapses. Our results suggest that mutations in KIF4A and KIF5C cause ID by tipping the balance between excitatory and inhibitory synaptic excitability.

Frequent mutations in chromatin-remodelling genes in pulmonary carcinoids

Pulmonary carcinoids are rare neuroendocrine tumours of the lung. The molecular alterations underlying the pathogenesis of these tumours have not been systematically studied so far. Here we perform gene copy number analysis (n=54), genome/exome (n=44) and transcriptome (n=69) sequencing of pulmonary carcinoids and observe frequent mutations in chromatin-remodelling genes. Covalent histone modifiers and subunits of the SWI/SNF complex are mutated in 40 and 22.2% of the cases, respectively, with MEN1, PSIP1 and ARID1A being recurrently affected. In contrast to small-cell lung cancer and large-cell neuroendocrine lung tumours, TP53 and RB1 mutations are rare events, suggesting that pulmonary carcinoids are not early progenitor lesions of the highly aggressive lung neuroendocrine tumours but arise through independent cellular mechanisms. These data also suggest that inactivation of chromatin-remodelling genes is sufficient to drive transformation in pulmonary carcinoids.

CD74-NRG1 fusions in lung adenocarcinoma

We discovered a novel somatic gene fusion, CD74-NRG1, by transcriptome sequencing of 25 lung adenocarcinomas of never smokers. By screening 102 lung adenocarcinomas negative for known oncogenic alterations, we found four additional fusion-positive tumors, all of which were of the invasive mucinous subtype. Mechanistically, CD74-NRG1 leads to extracellular expression of the EGF-like domain of NRG1 III-β3, thereby providing the ligand for ERBB2-ERBB3 receptor complexes. Accordingly, ERBB2 and ERBB3 expression was high in the index case, and expression of phospho-ERBB3 was specifically found in tumors bearing the fusion (P < 0.0001). Ectopic expression of CD74-NRG1 in lung cancer cell lines expressing ERBB2 and ERBB3 activated ERBB3 and the PI3K-AKT pathway, and led to increased colony formation in soft agar. Thus, CD74-NRG1 gene fusions are activating genomic alterations in invasive mucinous adenocarcinomas and may offer a therapeutic opportunity for a lung tumor subtype with, so far, no effective treatment.

SIGNIFICANCE

CD74–NRG1 fusions may represent a therapeutic opportunity for invasive mucinous lung adenocarcinomas, a tumor with no effective treatment that frequently presents with multifocal unresectable disease.

Clinical and neurocognitive characterization of a family with a novel MED12 gene frameshift mutation

FG syndrome, Lujan syndrome, and Ohdo syndrome, the Maat-Kievit-Brunner type, have been described as distinct syndromes with overlapping non-specific features and different missense mutations of the MED12 gene have been reported in all of them. We report a family including 10 males and 1 female affected with profound non-specific intellectual disability (ID) which was linked to a 30-cM region extending from Xp11.21 (ALAS2) to Xq22.3 (COL4A5). Parallel sequencing of all X-chromosome exons identified a frameshift mutation (c.5898dupC) of MED12. Mutated mRNA was not affected by non-sense mediated RNA decay and induced an additional abnormal isoform due to activation of cryptic splice-sites in exon 41. Dysmorphic features common to most affected males were long narrow face, high forehead, flat malar area, high nasal bridge, and short philtrum. Language was absent or very limited. Most patients had a friendly personality. Cognitive impairment, varying from borderline to profound ID was similarly observed in seven heterozygous females. There was no correlation between cognitive function and X-chromosome inactivation profiles in blood cells. The severe degree of ID in male patients, as well as variable cognitive impairment in heterozygous females suggests that the duplication observed in the present family may have a more severe effect on MED12 function than missense mutations. In a cognitively impaired male from this family, who also presented with tall stature and dysmorphism and did not have the MED12 mutation, a 600-kb duplication at 17p13.3 including the YWHAE gene, was found in a mosaic state.

Comparison of an automated respiratory systolic variation test with dynamic preload indicators to predict fluid responsiveness after major surgery

BACKGROUND

Predicting the response of cardiac output to volume administration remains an ongoing clinical challenge. The objective of our study was to compare the ability to predict volume responsiveness of various functional measures of cardiac preload. These included pulse pressure variation (PPV), stroke volume variation (SVV), and the recently launched automated respiratory systolic variation test (RSVT) in patients after major surgery.

METHODS

In this prospective study, 24 mechanically ventilated patients after major surgery were enrolled. Three consecutive volume loading steps consisting of 300 ml 6% hydroxyethylstarch 130/0.4 were performed and cardiac index (CI) was assessed by transpulmonary thermodilution. Volume responsiveness was considered as positive if CI increased by >10%.

RESULTS

In total 72 volume loading steps were analysed, of which 41 showed a positive volume response. Receiver operating characteristic (ROC) curve analysis revealed an area under the curve (AUC) of 0.70 for PPV, 0.72 for SVV and 0.77 for RSVT. Areas under the curves of all variables did not differ significantly from each other (P>0.05). Suggested cut-off values were 9.9% for SVV, 10.1% for PPV, and 19.7° for RSVT as calculated by the Youden Index.

CONCLUSION

In predicting fluid responsiveness the new automated RSVT appears to be as accurate as established dynamic indicators of preload PPV and SVV in patients after major surgery. The automated RSVT is clinically easy to use and may be useful in guiding fluid therapy in ventilated patients.

Loss of function of KIAA2022 causes mild to severe intellectual disability with an autism spectrum disorder and impairs neurite outgrowth

Existence of a discrete new X-linked intellectual disability (XLID) syndrome due to KIAA2022 deficiency was questioned by disruption of KIAA2022 by an X-chromosome pericentric inversion in a XLID family we reported in 2004. Three additional families with likely pathogenic KIAA2022 mutations were discovered within the frame of systematic parallel sequencing of familial cases of XLID or in the context of routine array-CGH evaluation of sporadic intellectual deficiency (ID) cases. The c.186delC and c.3597dupA KIAA2022 truncating mutations were identified by X-chromosome exome sequencing, while array CGH discovered a 70 kb microduplication encompassing KIAA2022 exon 1 in the third family. This duplication decreased KIAA2022 mRNA level in patients' lymphocytes by 60%. Detailed clinical examination of all patients, including the two initially reported, indicated moderate-to-severe ID with autistic features, strabismus in all patients, with no specific dysmorphic features other than a round face in infancy and no structural brain abnormalities on magnetic resonance imaging (MRI). Interestingly, the patient with decreased KIAA2022 expression had only mild ID with severe language delay and repetitive behaviors falling in the range of an autism spectrum disorder (ASD). Since little is known about KIAA2022 function, we conducted morphometric studies in cultured rat hippocampal neurons. We found that siRNA-mediated KIAA2022 knockdown resulted in marked impairment in neurite outgrowth including both the dendrites and the axons, suggesting a major role for KIAA2022 in neuron development and brain function.

A noncoding, regulatory mutation implicates HCFC1 in nonsyndromic intellectual disability

The discovery of mutations causing human disease has so far been biased toward protein-coding regions. Having excluded all annotated coding regions, we performed targeted massively parallel resequencing of the nonrepetitive genomic linkage interval at Xq28 of family MRX3. We identified in the binding site of transcription factor YY1 a regulatory mutation that leads to overexpression of the chromatin-associated transcriptional regulator HCFC1. When tested on embryonic murine neural stem cells and embryonic hippocampal neurons, HCFC1 overexpression led to a significant increase of the production of astrocytes and a considerable reduction in neurite growth. Two other nonsynonymous, potentially deleterious changes have been identified by X-exome sequencing in individuals with intellectual disability, implicating HCFC1 in normal brain function.

Detecting genomic indel variants with exact breakpoints in single- and paired-end sequencing data using SplazerS

MOTIVATION

The reliable detection of genomic variation in resequencing data is still a major challenge, especially for variants larger than a few base pairs. Sequencing reads crossing boundaries of structural variation carry the potential for their identification, but are difficult to map.

RESULTS

Here we present a method for 'split' read mapping, where prefix and suffix match of a read may be interrupted by a longer gap in the read-to-reference alignment. We use this method to accurately detect medium-sized insertions and long deletions with precise breakpoints in genomic resequencing data. Compared with alternative split mapping methods, SplazerS significantly improves sensitivity for detecting large indel events, especially in variant-rich regions. Our method is robust in the presence of sequencing errors as well as alignment errors due to genomic mutations/divergence, and can be used on reads of variable lengths. Our analysis shows that SplazerS is a versatile tool applicable to unanchored or single-end as well as anchored paired-end reads. In addition, application of SplazerS to targeted resequencing data led to the interesting discovery of a complete, possibly functional gene retrocopy variant.

AVAILABILITY

SplazerS is available from http://www.seqan.de/projects/ splazers.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Next-generation sequencing identifies mutations of SMPX, which encodes the small muscle protein, X-linked, as a cause of progressive hearing impairment

In a Dutch family with an X-linked postlingual progressive hearing impairment, a critical linkage interval was determined to span a region of 12.9 Mb flanked by the markers DXS7108 and DXS7110. This interval overlaps with the previously described DFNX4 locus and contains 75 annotated genes. Subsequent next-generation sequencing (NGS) detected one variant within the linkage interval, a nonsense mutation in SMPX. SMPX encodes the small muscle protein, X-linked (SMPX). Further screening was performed on 26 index patients from small families for which X-linked inheritance of nonsyndromic hearing impairment (NSHI) was not excluded. We detected a frameshift mutation in SMPX in one of the patients. Segregation analysis of both mutations in the families in whom they were found revealed that the mutations cosegregated with hearing impairment. Although we show that SMPX is expressed in many different organs, including the human inner ear, no obvious symptoms other than hearing impairment were observed in the patients. SMPX had previously been demonstrated to be specifically expressed in striated muscle and, therefore, seemed an unlikely candidate gene for hearing impairment. We hypothesize that SMPX functions in inner ear development and/or maintenance in the IGF-1 pathway, the integrin pathway through Rac1, or both.

Functional analysis and identification of cis-regulatory elements of human chromosome 21 gene promoters

Given the inherent limitations of in silico studies relying solely on DNA sequence analysis, the functional characterization of mammalian promoters and associated cis-regulatory elements requires experimental support, which demands cloning and analysis of putative promoter regions. Focusing on human chromosome 21, we cloned 182 gene promoters of 2500 bp in length and conducted reporter gene assays on transfected-cell arrays. We found 56 promoters that were active in HEK293 cells, while another 49 promoters could be activated by treatment of cells with Trichostatin A or depletion of serum. We observed high correlations between promoter activities and endogenous transcript levels, RNA polymerase II occupancy, CpG islands and core promoter elements. Truncation of a subset of 62 promoters to ∼500 bp revealed that truncation rarely resulted in loss of activity, but rather in loss of responses to external stimuli, suggesting the presence of cis-regulatory response elements within distal promoter regions. In these regions, we found a strong enrichment of transcription factor binding sites that could potentially activate gene expression in the presence of stimuli. This study illustrates the modular functional architecture of chromosome 21 promoters and helps to reveal the complex mechanisms governing transcriptional regulation.

Erratum to: Mutation screening in 86 known X-linked mental retardation genes by droplet-based multiplex PCR and massive parallel sequencing

[This corrects the article DOI: 10.1007/s11568-010-9137-y.].

Mutation screening in 86 known X-linked mental retardation genes by droplet-based multiplex PCR and massive parallel sequencing

Massive parallel sequencing has revolutionized the search for pathogenic variants in the human genome, but for routine diagnosis, re-sequencing of the complete human genome in a large cohort of patients is still far too expensive. Recently, novel genome partitioning methods have been developed that allow to target re-sequencing to specific genomic compartments, but practical experience with these methods is still limited. In this study, we have combined a novel droplet-based multiplex PCR method and next generation sequencing to screen patients with X-linked mental retardation (XLMR) for mutations in 86 previously identified XLMR genes. In total, affected males from 24 large XLMR families were analyzed, including three in whom the mutations were already known. Amplicons corresponding to functionally relevant regions of these genes were sequenced on an Illumina/Solexa Genome Analyzer II platform. Highly specific and uniform enrichment was achieved: on average, 67.9% unambiguously mapped reads were derived from amplicons, and for 88.5% of the targeted bases, the sequencing depth was sufficient to reliably detect variations. Potentially disease-causing sequence variants were identified in 10 out of 24 patients, including the three mutations that were already known, and all of these could be confirmed by Sanger sequencing. The robust performance of this approach demonstrates the general utility of droplet-based multiplex PCR for parallel mutation screening in hundreds of genes, which is a prerequisite for the diagnosis of mental retardation and other disorders that may be due to defects of a wide variety of genes.

Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments

Alternative splicing, polyadenylation of pre-messenger RNA molecules and differential promoter usage can produce a variety of transcript isoforms whose respective expression levels are regulated in time and space, thus contributing specific biological functions. However, the repertoire of mammalian alternative transcripts and their regulation are still poorly understood. Second-generation sequencing is now opening unprecedented routes to address the analysis of entire transcriptomes. Here, we developed methods that allow the prediction and quantification of alternative isoforms derived solely from exon expression levels in RNA-Seq data. These are based on an explicit statistical model and enable the prediction of alternative isoforms within or between conditions using any known gene annotation, as well as the relative quantification of known transcript structures. Applying these methods to a human RNA-Seq dataset, we validated a significant fraction of the predictions by RT-PCR. Data further showed that these predictions correlated well with information originating from junction reads. A direct comparison with exon arrays indicated improved performances of RNA-Seq over microarrays in the prediction of skipped exons. Altogether, the set of methods presented here comprehensively addresses multiple aspects of alternative isoform analysis. The software is available as an open-source R-package called Solas at http://cmb.molgen.mpg.de/2ndGenerationSequencing/Solas/.

CpG-depleted promoters harbor tissue-specific transcription factor binding signals--implications for motif overrepresentation analyses

Motif overrepresentation analysis of proximal promoters is a common approach to characterize the regulatory properties of co-expressed sets of genes. Here we show that these approaches perform well on mammalian CpG-depleted promoter sets that regulate expression in terminally differentiated tissues such as liver and heart. In contrast, CpG-rich promoters show very little overrepresentation signal, even when associated with genes that display highly constrained spatiotemporal expression. For instance, while ∼50% of heart specific genes possess CpG-rich promoters we find that the frequently observed enrichment of MEF2-binding sites upstream of heart-specific genes is solely due to contributions from CpG-depleted promoters. Similar results are obtained for all sets of tissue-specific genes indicating that CpG-rich and CpG-depleted promoters differ fundamentally in their distribution of regulatory inputs around the transcription start site. In order not to dilute the respective transcription factor binding signals, the two promoter types should thus be treated as separate sets in any motif overrepresentation analysis.

PASTAA: identifying transcription factors associated with sets of co-regulated genes

MOTIVATION

A major challenge in regulatory genomics is the identification of associations between functional categories of genes (e.g. tissues, metabolic pathways) and their regulating transcription factors (TFs). While, for a limited number of categories, the regulating TFs are already known, still for many functional categories the responsible factors remain to be elucidated.

RESULTS

We put forward a novel method (PASTAA) for detecting transcriptions factors associated with functional categories, which utilizes the prediction of binding affinities of a TF to promoters. This binding strength information is compared to the likelihood of membership of the corresponding genes in the functional category under study. Coherence between the two ranked datasets is seen as an indicator of association between a TF and the category. PASTAA is applied primarily to the determination of TFs driving tissue-specific expression. We show that PASTAA is capable of recovering many TFs acting tissue specifically and, in addition, provides novel associations so far not detected by alternative methods. The application of PASTAA to detect TFs involved in the regulation of tissue-specific gene expression revealed a remarkable number of experimentally supported associations. The validated success for various datasets implies that PASTAA can directly be applied for the detection of TFs associated with newly derived gene sets.

AVAILABILITY

The PASTAA source code as well as a corresponding web interface is freely available at http://trap.molgen.mpg.de.

Gene identification and analysis of transcripts differentially regulated in fracture healing by EST sequencing in the domestic sheep

BACKGROUND

The sheep is an important model animal for testing novel fracture treatments and other medical applications. Despite these medical uses and the well known economic and cultural importance of the sheep, relatively little research has been performed into sheep genetics, and DNA sequences are available for only a small number of sheep genes.

RESULTS

In this work we have sequenced over 47 thousand expressed sequence tags (ESTs) from libraries developed from healing bone in a sheep model of fracture healing. These ESTs were clustered with the previously available 10 thousand sheep ESTs to a total of 19087 contigs with an average length of 603 nucleotides. We used the newly identified sequences to develop RT-PCR assays for 78 sheep genes and measured differential expression during the course of fracture healing between days 7 and 42 postfracture. All genes showed significant shifts at one or more time points. 23 of the genes were differentially expressed between postfracture days 7 and 10, which could reflect an important role for these genes for the initiation of osteogenesis.

CONCLUSION

The sequences we have identified in this work are a valuable resource for future studies on musculoskeletal healing and regeneration using sheep and represent an important head-start for genomic sequencing projects for Ovis aries, with partial or complete sequences being made available for over 5,800 previously unsequenced sheep genes.

T-STAG: resource and web-interface for tissue-specific transcripts and genes

T-STAG (tissue-specific transcripts and genes) is a resource and web-interface, designated to analyze tissue/tumor-specific expression patterns in human and mouse transcriptomes. It integrates our refined prediction of specific expression patterns both in genes as well as in individual isoforms with man–mouse orthology data. In combination with the features for combining/contrasting the genes expressed in different tissues, T-STAG implicates important biological applications, such as the detection of differentially expressed genes in tumors, the retrieval of orthologs with significant expression in the same tissue etc. Additionally, our refined categorization of expressed sequence tags (ESTs) according to the normalization of cDNA libraries allows searching for putative low-abundant transcripts. The results are tightly linked to our visualization tools, GeneNest (expression patterns of genes) and SpliceNest (gene structure and alternative splicing). The user-friendly interface of T-STAG offers a platform for comprehensive analysis of tissue and/or tumor-specific expression patterns revealed by the EST data. T-STAG is freely accessible at .

Intronic CA-repeat and CA-rich elements: a new class of regulators of mammalian alternative splicing

We have recently identified an intronic polymorphic CA-repeat region in the human endothelial nitric oxide synthase (eNOS) gene as an important determinant of the splicing efficiency, requiring specific binding of hnRNP L. Here, we analyzed the position requirements of this CA-repeat element, which revealed its potential role in alternative splicing. In addition, we defined the RNA binding specificity of hnRNP L by SELEX: not only regular CA repeats are recognized with high affinity but also certain CA-rich clusters. Therefore, we have systematically searched the human genome databases for CA-repeat and CA-rich elements associated with alternative 5' splice sites (5'ss), followed by minigene transfection assays. Surprisingly, in several specific human genes that we tested, intronic CA RNA elements could function either as splicing enhancers or silencers, depending on their proximity to the alternative 5'ss. HnRNP L was detected specifically bound to these diverse CA elements. These data demonstrated that intronic CA sequences constitute novel and widespread regulatory elements of alternative splicing.

Strengths and weaknesses of EST-based prediction of tissue-specific alternative splicing

Background

Alternative splicing contributes significantly to the complexity of the human transcriptome and proteome. Computational prediction of alternative splice isoforms are usually based on EST sequences that also allow to approximate the expression pattern of the related transcripts. However, the limited number of tissues represented in the EST data as well as the different cDNA construction protocols may influence the predictive capacity of ESTs to unravel tissue-specifically expressed transcripts.

Methods

We predict tissue and tumor specific splice isoforms based on the genomic mapping (SpliceNest) of the EST consensus sequences and library annotation provided in the GeneNest database. We further ascertain the potentially rare tissue specific transcripts as the ones represented only by ESTs derived from normalized libraries. A subset of the predicted tissue and tumor specific isoforms are then validated via RT-PCR experiments over a spectrum of 40 tissue types.

Results

Our strategy revealed 427 genes with at least one tissue specific transcript as well as 1120 genes showing tumor specific isoforms. While our experimental evaluation of computationally predicted tissue-specific isoforms revealed a high success rate in confirming the expression of these isoforms in the respective tissue, the strategy frequently failed to detect the expected restricted expression pattern. The analysis of putative lowly expressed transcripts using normalized cDNA libraries suggests that our ability to detect tissue-specific isoforms strongly depends on the expression level of the respective transcript as well as on the sensitivity of the experimental methods. Especially splice isoforms predicted to be disease-specific tend to represent transcripts that are expressed in a set of healthy tissues rather than novel isoforms.

Conclusions

We propose to combine the computational prediction of alternative splice isoforms with experimental validation for efficient delineation of an accurate set of tissue-specific transcripts.

Genome wide identification and classification of alternative splicing based on EST data

MOTIVATION

Alternative splicing is currently seen to explain the vast disparity between the number of predicted genes in the human genome and the highly diverse proteome. The mapping of expressed sequences tag (EST) consensus sequences derived from the GeneNest database onto the genome provides an efficient way of predicting exon-intron boundaries, gene structure and alternative splicing events. However, the alternative splicing events are obscured by a large number of putatively artificial exon boundaries arising due to genomic contamination or alignment errors. The current work describes a methodology to associate quality values to the predicted exon-intron boundaries. High quality exon-intron boundaries are used to predict constitutive and alternative splicing ranked by confidence values, aiming to facilitate large-scale analysis of alternative splicing and splicing in general.

RESULTS

Applying the current methodology, constitutive splicing is observed in 33,270 EST clusters, out of which 45% are alternatively spliced. The classification derived from the computed confidence values for 17 of these splice events frequently correlate (15/17) with RT-PCR experiments performed for 40 different tissue samples. As an application of the confidence measure, an evaluation of distribution of alternative splicing revealed that majority of variants correspond to the coding regions of the genes. However, still a significant fraction maps to non-coding regions, thereby indicating a functional relevance of alternative splicing in untranslated regions.

AVAILABILITY

The predicted alternative splice variants are visualized in the SpliceNest database at http://splicenest.molgen.mpg.de

A DNA microarray for fission yeast: minimal changes in global gene expression after temperature shift

Completion of the fission yeast genome sequence has opened up possibilities for post-genomic approaches. We have constructed a DNA microarray for genome-wide gene expression analysis in fission yeast. The microarray contains DNA fragments, PCR-amplified from a genomic DNA template, that represent > 99% of the 5000 or so annotated fission yeast genes, as well as a number of control sequences. The GenomePRIDE software used attempts to design similarly sized DNA fragments corresponding to gene regions within single exons, near the 3'-end of genes that lack homology to other fission yeast genes. To validate the design and utility of the array, we studied expression changes after a 2 h temperature shift from 25 degrees C to 36 degrees C, conditions widely used when studying temperature-sensitive mutants. Obligingly, the vast majority of genes do not change more than two-fold, supporting the widely held view that temperature-shift experiments specifically reveal phenotypes associated with temperature-sensitive mutants. However, we did identify a small group of genes that showed a reproducible change in expression. Importantly, most of these corresponded to previously characterized heat-shock genes, whose expression has been reported to change after more extreme temperature shifts than those used here. We conclude that the DNA microarray represents a useful resource for fission yeast researchers as well as the broader yeast community, since it will facilitate comparison with the distantly related budding yeast, Saccharomyces cerevisiae. To maximize the utility of this resource, the array and its component parts are fully described in On-line Supplementary Information and are also available commercially.

Genome-wide RNAi analysis of growth and viability in Drosophila cells

A crucial aim upon completion of whole genome sequences is the functional analysis of all predicted genes. We have applied a high-throughput RNA-interference (RNAi) screen of 19,470 double-stranded (ds) RNAs in cultured cells to characterize the function of nearly all (91%) predicted Drosophila genes in cell growth and viability. We found 438 dsRNAs that identified essential genes, among which 80% lacked mutant alleles. A quantitative assay of cell number was applied to identify genes of known and uncharacterized functions. In particular, we demonstrate a role for the homolog of a mammalian acute myeloid leukemia gene (AML1) in cell survival. Such a systematic screen for cell phenotypes, such as cell viability, can thus be effective in characterizing functionally related genes on a genome-wide scale.

An integrated gene annotation and transcriptional profiling approach towards the full gene content of the Drosophila genome

A novel Drosophila microarray constructed on the basis of an integrated in silico/wet biology approach provides evidence for the transcription of approximately 2,600 additional genes. Validation indicates a lower limit of 2,000 novel annotations, thus raising the number of genes that make a fly.

Background

While the genome sequences for a variety of organisms are now available, the precise number of the genes encoded is still a matter of debate. For the human genome several stringent annotation approaches have resulted in the same number of potential genes, but a careful comparison revealed only limited overlap. This indicates that only the combination of different computational prediction methods and experimental evaluation of such in silico data will provide more complete genome annotations. In order to get a more complete gene content of the Drosophila melanogaster genome, we based our new D. melanogaster whole-transcriptome microarray, the Heidelberg FlyArray, on the combination of the Berkeley Drosophila Genome Project (BDGP) annotation and a novel ab initio gene prediction of lower stringency using the Fgenesh software.

Results

Here we provide evidence for the transcription of approximately 2,600 additional genes predicted by Fgenesh. Validation of the developmental profiling data by RT-PCR and in situ hybridization indicates a lower limit of 2,000 novel annotations, thus substantially raising the number of genes that make a fly.

Conclusions

The successful design and application of this novel Drosophila microarray on the basis of our integrated in silico/wet biology approach confirms our expectation that in silico approaches alone will always tend to be incomplete. The identification of at least 2,000 novel genes highlights the importance of gathering experimental evidence to discover all genes within a genome. Moreover, as such an approach is independent of homology criteria, it will allow the discovery of novel genes unrelated to known protein families or those that have not been strictly conserved between species.

Genome-scale design of PCR primers and long oligomers for DNA microarrays

During the last years, the demand for custom-made cDNA chips/arrays as well as whole genome chips is increasing rapidly. The efficient selection of gene-specific primers/oligomers is of the utmost importance for the successful production of such chips. We developed GenomePRIDE, a highly flexible and scalable software for designing primers/oligomers for large-scale projects. The program is able to generate either long oligomers (40-70 bases), or PCR primers for the amplification of gene-specific DNA fragments of user-defined length. Additionally, primers can be designed in-frame in order to facilitate large-scale cloning into expression vectors. Furthermore, GenomePRIDE can be adapted to specific applications such as the generation of genomic amplicon arrays or the design of fragments specific for alternative splice isoforms. We tested the performance of GenomePRIDE on the entire genomes of Listeria monocytogenes (1584 gene-specific PCRs, 48 long oligomers) as well as of eukaryotes such as Schizosaccharomyces pombe (5006 gene-specific PCRs), and Drosophila melanogaster (21 306 gene-specific PCRs). With its computing speed of 1000 primer pairs per hour and a PCR amplification success of 99%, GenomePRIDE represents an extremely cost- and time-effective program.

SYSTERS, GeneNest, SpliceNest: exploring sequence space from genome to protein

We have integrated the protein families from SYSTERS and the expressed sequence tag (EST) clusters from our database GeneNest with SpliceNest, a new database mapping EST contigs into genomic DNA. The SYSTERS protein sequence cluster set provides an automatically generated classification of all sequences of the SWISS-PROT, TrEMBL and PIR databases into disjoint protein family and superfamily clusters. GeneNest is a database and software package for producing and visualizing gene indices from ESTs and mRNAs. Currently, the database comprises gene indices of human, mouse, Arabidopsis thaliana and zebrafish. SpliceNest is a web-based graphical tool to explore gene structure, including alternative splicing, based on a mapping of the EST consensus sequences from GeneNest to the complete human genome. The integration of SYSTERS, GeneNest and SpliceNest into one framework now permits an overall exploration of the whole sequence space covering protein, mRNA and EST sequences, as well as genomic DNA. The databases are available for querying and browsing at http://cmb.molgen.mpg.de.