Circadian clock genes and sleep homeostasis.

Circadian and sleep-homeostatic processes both contribute to sleep timing and sleep structure. Elimination of circadian rhythms through lesions of the suprachiasmatic nuclei (SCN), the master circadian pacemaker, leads to fragmentation of wakefulness and sleep but does not eliminate the homeostatic response to sleep loss as indexed by the increase in EEG delta power. In humans, EEG delta power declines during sleep episodes nearly independently of circadian phase. Such observations have contributed to the prevailing notion that circadian and homeostatic processes are separate but recent data imply that this segregation may not extend to the molecular level. Here we summarize the criteria and evidence for a role for clock genes in sleep homeostasis. Studies in mice with targeted disruption for core circadian clock genes have revealed alterations in circadian rhythmicity as well as changes in sleep duration, sleep structure and EEG delta power. Clock-gene expression in brain areas outside the SCN, in particular the cerebral cortex, depends to a large extent on prior sleep-wake history. Evidence for effects of clock genes on sleep homeostasis has also been obtained in Drosophila and humans, pointing to a phylogenetically preserved pathway. These findings suggest that, while within the SCN clock genes are utilized to set internal time-of-day, in the forebrain the same feedback circuitry may be utilized to track time spent awake and asleep. The mechanisms by which clock-gene expression is coupled to the sleep-wake distribution could be through cellular energy charge whereby clock genes act as energy sensors. The data underscore the interrelationships between energy metabolism, circadian rhythmicity, and sleep regulation.

Genetics for clinicians: from candidate genes to whole genome scans (technological advances).

Human genetics has progressed at an unprecedented pace during the past 10 years. DNA microarrays currently allow screening of the entire human genome with high level of coverage and we are now entering the era of high-throughput sequencing. These remarkable technical advances are influencing the way medical research is conducted and have boosted our understanding of the structure of the human genome as well as of disease biology. In this context, it is crucial for clinicians to understand the main concepts and limitations of modern genetics. This review will describe key concepts in genetics, including the different types of genetic markers in the human genome, review current methods to detect DNA variation, describe major online public databases in genetics, explain key concepts in statistical genetics and finally present commonly used study designs in clinical and epidemiological research. This review will therefore concentrate on human genetic variation analysis.

Global GacA-steered control of cyanide and exoprotease production in Pseudomonas fluorescens involves specific ribosome binding sites.

The conserved two-component regulatory system GacS/GacA determines the expression of extracellular products and virulence factors in a variety of Gram-negative bacteria. In the biocontrol strain CHA0 of Pseudomonas fluorescens, the response regulator GacA is essential for the synthesis of extracellular protease (AprA) and secondary metabolites including hydrogen cyanide. GacA was found to exert its control on the hydrogen cyanide biosynthetic genes (hcnABC) and on the aprA gene indirectly via a posttranscriptional mechanism. Expression of a translational hcnA'-'lacZ fusion was GacA-dependent whereas a transcriptional hcnA-lacZ fusion was not. A distinct recognition site overlapping with the ribosome binding site appears to be primordial for GacA-steered regulation. GacA-dependence could be conferred to the Escherichia coli lacZ mRNA by a 3-bp substitution in the ribosome binding site. The gene coding for the global translational repressor RsmA of P. fluorescens was cloned. RsmA overexpression mimicked partial loss of GacA function and involved the same recognition site, suggesting that RsmA is a downstream regulatory element of the GacA control cascade. Mutational inactivation of the chromosomal rsmA gene partially suppressed a gacS defect. Thus, a central, GacA-dependent switch from primary to secondary metabolism may operate at the level of translation.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.

Sox genes in C. elegans: sox-2 role in postembryonic development

Report for the scientific sojourn carried out at the Columbia University, United States, from 2010 to 2012. Expression of SoxB genes correlates with the commitment of cells to a neural fate; however, the relevance of SoxB proteins in early vertebrate neurogenesis has been difficult to prove genetically due to embryonic lethality and presumed redundant functions. The nematode C. Elegants has only 5 sox genes: sox-2 and sox-3 form the SoxB group while sem-2, sox-4 and egl-13 belong to other Sox groups. Our results show that sox-2 and sem-2 are the sox genes expressed earliest and in a broader manner during embryogenesis, being expressed in several neuronal progenitors. sox-3, sox-4 and egl-13 are expressed in few cells during late embryogenesis, when most neurons are already born. Both sox-2 and sem-2 null mutants are early larval lethal but do not show neuronal specification defects during embryonic development as indicated by quantification of a panneuronal reporter. Potential redundancy or compensatory mechanisms between different sox genes have been ruled out, strongly suggesting that sox genes are not required for specification of embryonically-derived neurons. However, at the first larval stage there are still several blast cells that will give rise to different postembryonic lineages, which generate several neurons amongst other cell types. nterestingly, sox-2 is expressed in many of these progenitor cells. Using mosaic analysis we have so far identified neurons derived from two different postembryonic lineages which fail to be generated in C. elegans sox-2 mutants. These results support the idea that postembryonic progenitor competence is compromised in the absence of sox-2.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights.

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

Repression of flagellar genes in exponential phase by CsgD and CpxR, two crucial modulators of Escherichia coli biofilm formation.

Escherichia coli adapts its lifestyle to the variations of environmental growth conditions, swapping between swimming motility or biofilm formation. The stationary-phase sigma factor RpoS is an important regulator of this switch, since it stimulates adhesion and represses flagellar biosynthesis. By measuring the dynamics of gene expression, we show that RpoS inhibits the transcription of the flagellar sigma factor, FliA, in exponential growth phase. RpoS also partially controls the expression of CsgD and CpxR, two transcription factors important for bacterial adhesion. We demonstrate that these two regulators repress the transcription of fliA, flgM, and tar and that this regulation is dependent on the growth medium. CsgD binds to the flgM and fliA promoters around their -10 promoter element, strongly suggesting direct repression. We show that CsgD and CpxR also affect the expression of other known modulators of cell motility. We propose an updated structure of the regulatory network controlling the choice between adhesion and motility.

The CNGCb and CNGCd genes from Physcomitrella patens moss encode for thermosensory calcium channels responding to fluidity changes in the plasma membrane.

Land plants need precise thermosensors to timely establish molecular defenses in anticipation of upcoming noxious heat waves. The plasma membrane-embedded cyclic nucleotide-gated Ca(2+) channels (CNGCs) can translate mild variations of membrane fluidity into an effective heat shock response, leading to the accumulation of heat shock proteins (HSP) that prevent heat damages in labile proteins and membranes. Here, we deleted by targeted mutagenesis the CNGCd gene in two Physcomitrella patens transgenic moss lines containing either the heat-inducible HSP-GUS reporter cassette or the constitutive UBI-Aequorin cassette. The stable CNGCd knockout mutation caused a hyper-thermosensitive moss phenotype, in which the heat-induced entry of apoplastic Ca(2+) and the cytosolic accumulation of GUS were triggered at lower temperatures than in wild type. The combined effects of an artificial membrane fluidizer and elevated temperatures suggested that the gene products of CNGCd and CNGCb are paralogous subunits of Ca(2+)channels acting as a sensitive proteolipid thermocouple. Depending on the rate of temperature increase, the duration and intensity of the heat priming preconditions, terrestrial plants may thus acquire an array of HSP-based thermotolerance mechanisms against upcoming, otherwise lethal, extreme heat waves.

Bioinformática: interfaz web para mostrar las relaciones de expresión entre grupos de genes

El projecte que es presenta en aquesta memòria està composat per un preproces i una aplicació web. L’objectiu del treball realizat és mostrar les relacions d’expresió entre grups de gens. A més de dissenyar l’aplicació web en : http://revolutionresearch.uab.es que permetrà l’estudi dels resultats obtinguts.

Targeted Gene Deletion and In Vivo Analysis of Putative Virulence Gene Function in the Pathogenic Dermatophyte Arthroderma benhamiae.

Dermatophytes cause the majority of superficial mycoses in humans and animals. However, little is known about the pathogenicity of this specialized group of filamentous fungi, for which molecular research has been limited thus far. During experimental infection of guinea pigs by the human pathogenic dermatophyte Arthroderma benhamiae, we recently detected the activation of the fungal gene encoding malate synthase AcuE, a key enzyme of the glyoxylate cycle. By the establishment of the first genetic system for A. benhamiae, specific ΔacuE mutants were constructed in a wild-type strain and, in addition, in a derivative in which we inactivated the nonhomologous end-joining pathway by deletion of the A. benhamiae KU70 gene. The absence of AbenKU70 resulted in an increased frequency of the targeted insertion of linear DNA by homologous recombination, without notably altering the monitored in vitro growth abilities of the fungus or its virulence in a guinea pig infection model. Phenotypic analyses of ΔacuE mutants and complemented strains depicted that malate synthase is required for the growth of A. benhamiae on lipids, major constituents of the skin. However, mutant analysis did not reveal a pathogenic role of the A. benhamiae enzyme in guinea pig dermatophytosis or during epidermal invasion of the fungus in an in vitro model of reconstituted human epidermis. The presented efficient system for targeted genetic manipulation in A. benhamiae, paired with the analyzed infection models, will advance the functional characterization of putative virulence determinants in medically important dermatophytes.

Network-guided analysis of genes with altered somatic copy number and gene expression reveals pathways commonly perturbed in metastatic melanoma.

Cancer genomes frequently contain somatic copy number alterations (SCNA) that can significantly perturb the expression level of affected genes and thus disrupt pathways controlling normal growth. In melanoma, many studies have focussed on the copy number and gene expression levels of the BRAF, PTEN and MITF genes, but little has been done to identify new genes using these parameters at the genome-wide scale. Using karyotyping, SNP and CGH arrays, and RNA-seq, we have identified SCNA affecting gene expression ('SCNA-genes') in seven human metastatic melanoma cell lines. We showed that the combination of these techniques is useful to identify candidate genes potentially involved in tumorigenesis. Since few of these alterations were recurrent across our samples, we used a protein network-guided approach to determine whether any pathways were enriched in SCNA-genes in one or more samples. From this unbiased genome-wide analysis, we identified 28 significantly enriched pathway modules. Comparison with two large, independent melanoma SCNA datasets showed less than 10% overlap at the individual gene level, but network-guided analysis revealed 66% shared pathways, including all but three of the pathways identified in our data. Frequently altered pathways included WNT, cadherin signalling, angiogenesis and melanogenesis. Additionally, our results emphasize the potential of the EPHA3 and FRS2 gene products, involved in angiogenesis and migration, as possible therapeutic targets in melanoma. Our study demonstrates the utility of network-guided approaches, for both large and small datasets, to identify pathways recurrently perturbed in cancer.

SGP-1: prediction and validation of homologous genes based on sequence alignments

Conventional methods of gene prediction rely on the recognition of DNA-sequence signals, the coding potential or the comparison of a genomic sequence with a cDNA, EST, or protein database. Reasons for limited accuracy in many circumstances are species-specific training and the incompleteness of reference databases. Lately, comparative genome analysis has attracted increasing attention. Several analysis tools that are based on human/mouse comparisons are already available. Here, we present a program for the prediction of protein-coding genes, termed SGP-1 (Syntenic Gene Prediction), which is based on the similarity of homologous genomic sequences. In contrast to most existing tools, the accuracy of SGP-1 depends little on species-specific properties such as codon usage or the nucleotide distribution. SGP-1 may therefore be applied to nonstandard model organisms in vertebrates as well as in plants, without the need for extensive parameter training. In addition to predicting genes in large-scale genomic sequences, the program may be useful to validate gene structure annotations from databases. To this end, SGP-1 output also contains comparisons between predicted and annotated gene structures in HTML format. The program can be accessed via a Web server at http://soft.ice.mpg.de/sgp-1. The source code, written in ANSI C, is available on request from the authors.

Assembling genes from predicted exons in linear time with dynamic programming

In a number of programs for gene structure prediction in higher eukaryotic genomic sequences, exon prediction is decoupled from gene assembly: a large pool of candidate exons is predicted and scored from features located in the query DNA sequence, and candidate genes are assembled from such a pool as sequences of nonoverlapping frame-compatible exons. Genes are scored as a function of the scores of the assembled exons, and the highest scoring candidate gene is assumed to be the most likely gene encoded by the query DNA sequence. Considering additive gene scoring functions, currently available algorithms to determine such a highest scoring candidate gene run in time proportional to the square of the number of predicted exons. Here, we present an algorithm whose running time grows only linearly with the size of the set of predicted exons. Polynomial algorithms rely on the fact that, while scanning the set of predicted exons, the highest scoring gene ending in a given exon can be obtained by appending the exon to the highest scoring among the highest scoring genes ending at each compatible preceding exon. The algorithm here relies on the simple fact that such highest scoring gene can be stored and updated. This requires scanning the set of predicted exons simultaneously by increasing acceptor and donor position. On the other hand, the algorithm described here does not assume an underlying gene structure model. Indeed, the definition of valid gene structures is externally defined in the so-called Gene Model. The Gene Model specifies simply which gene features are allowed immediately upstream which other gene features in valid gene structures. This allows for great flexibility in formulating the gene identification problem. In particular it allows for multiple-gene two-strand predictions and for considering gene features other than coding exons (such as promoter elements) in valid gene structures.

SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteinsusually misidentified and misannotated in sequencedatabases. The presence of an in-frame UGA (stop)codon in the coding sequence of selenoproteingenes precludes their identification and correctannotation. The in-frame UGA codons are recodedto cotranslationally incorporate selenocysteine,a rare selenium-containing amino acid. The developmentof ad hoc experimental and, more recently,computational approaches have allowed the efficientidentification and characterization of theselenoproteomes of a growing number of species.Today, dozens of selenoprotein families have beendescribed and more are being discovered in recentlysequenced species, but the correct genomic annotationis not available for the majority of thesegenes. SelenoDB is a long-term project that aims toprovide, through the collaborative effort of experimentaland computational researchers, automaticand manually curated annotations of selenoproteingenes, proteins and SECIS elements. Version 1.0 ofthe database includes an initial set of eukaryoticgenomic annotations, with special emphasis on thehuman selenoproteome, for immediate inspectionby selenium researchers or incorporation into moregeneral databases. SelenoDB is freely available athttp://www.selenodb.org.

Identifying protein-coding genes in genomic sequences

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.