An integrated genetic and functional analysis of the role of type II transmembrane serine proteases (TMPRSSs) in hearing loss.

Building on our discovery that mutations in the transmembrane serine protease, TMPRSS3, cause nonsyndromic deafness, we have investigated the contribution of other TMPRSS family members to the auditory function. To identify which of the 16 known TMPRSS genes had a strong likelihood of involvement in hearing function, three types of biological evidence were examined: 1) expression in inner ear tissues; 2) location in a genomic interval that contains a yet unidentified gene for deafness; and 3) evaluation of hearing status of any available Tmprss knockout mouse strains. This analysis demonstrated that, besides TMPRSS3, another TMPRSS gene was essential for hearing and, indeed, mice deficient for Hepsin (Hpn) also known as Tmprss1 exhibited profound hearing loss. In addition, TMPRSS2, TMPRSS5, and CORIN, also named TMPRSS10, showed strong likelihood of involvement based on their inner ear expression and mapping position within deafness loci PKSR7, DFNB24, and DFNB25, respectively. These four TMPRSS genes were then screened for mutations in affected members of the DFNB24 and DFNB25 deafness families, and in a cohort of 362 sporadic deaf cases. This large mutation screen revealed numerous novel sequence variations including three potential pathogenic mutations in the TMPRSS5 gene. The mutant forms of TMPRSS5 showed reduced or absent proteolytic activity. Subsequently, TMPRSS genes with evidence of involvement in deafness were further characterized, and their sites of expression were determined. Tmprss1, 3, and 5 proteins were detected in spiral ganglion neurons. Tmprss3 was also present in the organ of Corti. TMPRSS1 and 3 proteins appeared stably anchored to the endoplasmic reticulum membranes, whereas TMPRSS5 was also detected at the plasma membrane. Collectively, these results provide evidence that TMPRSS1 and TMPRSS3 play and TMPRSS5 may play important and specific roles in hearing.

Construction and analysis of a new conditional ferritin H knockout mouse strain/

Résumé Le fer joue un rôle important dans la plupart des fonctions biologiques mais sa présence excessive provoque la production de molécules réactives d'oxygène (ROS) qui peuvent contribuer à diverses maladies. La protéine de stockage du fer, la ferritine H, capte l'excès en fer et le stocke sous forme non-toxique, ce qui empêche des dommages potentiels. La délétion de la ferritine H dans des souris knock-out a été essayée antérieurement, mais ces souris mouraient au stade précoce du développement embryonnaire. Pour étudier l'importance du fer, et en particulier son stockage dans la ferritine, et pour pouvoir mieux comprendre les fonctions de la ferritine H, nous avons créé un modèle de souris knock-out conditionnelles de la ferritine H, selon le système classique de Cre-LoxP. Le premier exon et la région du promoteur du gène de la ferritine H ont été entourés de sites loxP. La mortalité embryonnaire provoquée par la délétion constitutive du gène de la ferritine H a été confirmée en croisant nos souris avec des souris exprimant nestin-Cre1. En croisant nos souris avec des souris transgéniques Mx-Cre, nous avons observé que l'induction de Cre par injection de polyI-polyC provoque la délétion presque complète de la ferritine H dans le foie (> 99%) et la rate (> 88%). Ces tissus ont également perdu une grande partie de leur réserve de fer. Cette observation apporte pour la première fois la preuve in vivo que la ferritine H est indispensable pour le stockage du fer, que les fonctions de la ferritine H et de la ferritine L ne sont pas équivalentes, et que la ferritine L ne peut pas assumer seule la fonction de stockage du fer. Dans le foie des souris knock-out, l'expression de l'ARN messager de l'hepcidine a été induite après 10 jours. En même temps, l'expression de l'ARN messager des gènes codant pour des protéines de l'absorption de fer (DMT1, ferroportin, Dcytb1 et hephaestin) a été réprimée mais dans le duodénum seulement. L'expression d'hepcidine est inversément corrélée avec celle des gènes liés à l'absorption de fer. Cette observation corrobore des études antérieures. Mais, en plus, elle montre également que cette répression se produit seulement dans l'intestin. Nous pouvons ainsi tirer la conclusion suivante : ou bien l'hepcidine a un récepteur spécifique dans le duodénum ou bien les gènes liés à l'absorption de fer dans le duodénum ont un facteur spécifique de transcription sensible à l'hepcidine. Aucune répression de DMT1 et de ferroportin n'a été observée dans les macrophages de la rate après l'induction d'hepcidine. La délétion de ferritine H a entraîné une augmentation du taux de mortalité des cellules hépatiques, ainsi que des altérations dans l'architecture normale du tissu de la rate. Vu par l'immunohistologie, le nombre de lymphocytes B et T était réduit dans la rate, tendant à démontrer que la ferritine H et l'homéostase du fer jouent un rôle dans l'immunité. En conclusion, le modèle de souris knock-out conditionnelles de la ferritine H nous fournit un outil précieux pour l'étude in vivo du rôle joué par la ferritine dans l'homéostase du fer, dans les dommages créés par les ROS, ainsi que dans l'apoptose et l'immunité. Summary Iron plays an important role in most biological functions. However, excess of iron results in production of reactive oxygen species (ROS) which could substantially contribute to pathology of various diseases. Ferritin H scavenges excess of iron and stores it in non-toxic form and potentially prevents the damage. Fenitin H targeting in mice has been attempted before, however, straight knockout was lethal in early embryonic stage. To study the role of iron and its storage protein ferritin and to further elucidate ferritin H functions, we aimed at creating a conditional ferritin H knockout mouse model by classical Cre-LoxP system. First exon along with promoter region of the ferritin H gene was foxed. Embryonic lethality of the constitutive ferritin H deletion was confirmed by crossing the foxed mice with mice expressing nestin Cre-1 as transgene. Almost complete deletion was observed in liver (> 99%) and spleen (>88%) upon induction of Cre by injecting polyI-polyC in Fth Lox/Lox; MxCre mice. These tissues also lost substantial fraction of their iron stores. This provides first in vivo evidence that ferritin H is required for iron storage, ferritin H and L functions are not redundant and that ferritin L cannot perform iron storage function alone. Hepcidin mRNA expression was induced after 10 days in the livers of deleted mice and, simultaneously, mRNA expression of iron absorption related genes (DMT 1, ferroportin, Dcytb1 and hephaestin) was repressed in duodenum only. Hepcidin expression is inversely correlated with that of duodenal iron absorption related genes. This is in agreement with previous studies. However, we also show that this repression happens only in intestine. This leads to the conclusion that either hepcidin has a specific receptor in duodenum or the iron absorption related genes have duodenum specific transcription factor that is responsive to hepcidin. No repression of DMT1 and ferroportin was observed in spleen macrophages upon hepcidin induction. Ferritin H deletion showed increased cell death in liver and disruption of normal architecture of spleen. B lymphocytes were reduced in spleen on immunohistology which point towards a role of ferritin H and iron homeostasis in immunity. In conclusion, ferritin H conditional knockout mouse model provides us with an invaluable tool to study the in vivo role of ferritin H in iron homeostasis, ROS mediated damage, apoptosis and immunity.

Analysis of c-Myc function in mouse epidermis

Abstract The c-myc gene is one of the most frequently mutated oncogenes found in human tumors. c-Myc has been implicated in the regulation of various biological processes including cell cycle progression, cellular growth, differentiation, angiogenesis, immortalization and apoptosis. To assess the normal role of c-Myc in epithelial cell types in vitro and in vivo we have deleted the c-myc gene in keratinocytes and in the adult skin epidermis by conditional Cre/loxP mediated recombination. Similar to what we have previously shown in mouse embryonic fibroblasts acute elimination of c-Myc activity in cultured keratinocytes causes cells to cease proliferation and adapt a flat cell morphology. Mutant cells accumulate in a diploid Ki67neg stage, indicative of a quiescent Go stage. This demonstrates that c-Myc activity is essential to maintain keratinocytes in a productive cell cycle. In addition, mutant keratinocytes showed a defect in Ca2+ induced induction of the differentiation marker Keratin 1 suggesting a role for c-Myc during differentiation. To assess the in vivo role of c-Myc we used a tamoxifen inducible K5::CreERT transgene to delete the c-myc gene in the adult skin epidermis. Unexpectedly, despite strong c-Myc expression in the basal compartment it is not required for maintenance of the skin epidermis in the adult mouse. The epidermis appeared normal with respect to both proliferation and differentiation. In addition, no selection against c-Myc deficient epidermal cells occurred over many months, further confirming that c-Myc is dispensable for normal skin homeostasis. Even more surprising, TPA induced hyperproliferation also occurred in a c-Myc independent manner. Treatment of the skin with the mutagen DMBA prior to TPA is a classical way to induce papillomas by selecting for mutations that lead to dominant activation of the oncogene Ha-Ras. Most interestingly tumor formation was severely inhibited suggesting that tumor progression requires endogenous c-Myc. Further studies are required to address whether the role of c-Myc in the activation of telomerase or the Werner protein, or its role to induce angiogenesis is required for skin tumor progression, In conclusion, this work shows that while c-Myc is not required for maintenance or hyperplasia of mouse epidermis, it is essential for skin tumor progression in collaboration with Ras. Résumé Le gène c-myc est un des oncogènes les plus fréquemment mutés dans les tumeurs humaines. c-Myc est impliqué dans la régulation de processus biologiques variés, comme la progression du cycle cellulaire, la croissance cellulaire, la différenciation, l'angiogenèse, l'immortalisation et l'apoptose. Pour caractériser le rôle physiologique de c-Myc dans les cellules de type épithélial in vitro et in vivo, le gène c-myc a été délété dans des kératinocytes primaires et dans l'épiderme de peau de souris adultes par des recombinaisons conditionnelles (système Cre/loxP). De la même façon que dans les fibroblastes d'embryon de souris, l'élimination aiguë de l'activité de c-Myc dans les kératinocytes en culture primaire provoque l'arrêt de la prolifération des cellules et leur applatissement morphologique. Les cellules mutantes restent dans un stade diploïde Ki67neg, indiquant un stade quiescent Go. Cela démontre que l'activité de c-Myc est essentielle pour maintenir les kératinocytes dans le cycle cellulaire. De plus, les kératinocytes mutants montrent une déficience pour le marqueur de différenciation Kératine 1 au cours de la différenciation induite par le calcium, suggérant un rôle de c-Myc dans la différenciation cellulaire. Pour comprendre le rôle de c-Myc in vivo, le transgène K5::CreERT inductible par le tamoxifen a été utilisé pour déléter le gène c-inyc dans l'épiderme de souris adultes. Etonnemment, malgré une forte expression de c-Myc dans le compartiment basal de l'épiderme, ce gène n'est pas nécessaire pour la maintenance de l'épiderme de la peau chez la souris adulte. L'épiderme apparait normal avec une prolifération et une différenciation physiologique des cellules. De plus, il n'y a pas de sélection contre les cellules épidennales c-Myc déficientes après plusieurs mois, ce qui confirme que c-Myc n'est pas nécessaire pour l'homéostasie normale de la peau. Encore plus surprenant, une hyperprolifération est également induite par du TPA chez les souris mutantes, impliquant une voie de prolifération indépendante de c-Myc. Le traitement de la peau par le mutagène DMBA avant le traitement au TPA est une voie classique d'induction de papillomes, par sélection de mutations conduisant à l'activation de l'oncogène Ha-Ras. La formation des tumeurs est fortement inhibée chez les souris mutantes, suggérant que la progression des tumeurs nécessite la présence endogène de c-Myc. De nouvelles études sont nécessaires pour savoir si c-Myc a un rôle dans l'activation de la télomérase ou de la protéine de Werner, ou encore dans l'angiogénèse, qui sont nécessaires pour la progression tumorale. En conclusion, ce travail montre que même si c-Myc n'est pas nécessaire pour la maintenance ou l'hyperplasie de la peau de souris, il est essentiel pour la progression des tumeurs de la peau en collaboration avec Ras.

Family history of cancer: pooled analysis in the International Head and Neck Cancer Epidemiology Consortium.

Alcohol and tobacco consumption are well-recognized risk factors for head and neck cancer (HNC). Evidence suggests that genetic predisposition may also play a role. Only a few epidemiologic studies, however, have considered the relation between HNC risk and family history of HNC and other cancers. We pooled individual-level data across 12 case-control studies including 8,967 HNC cases and 13,627 controls. We obtained pooled odds ratios (OR) using fixed and random effect models and adjusting for potential confounding factors. All statistical tests were two-sided. A family history of HNC in first-degree relatives increased the risk of HNC (OR=1.7, 95% confidence interval, CI, 1.2-2.3). The risk was higher when the affected relative was a sibling (OR=2.2, 95% CI 1.6-3.1) rather than a parent (OR=1.5, 95% CI 1.1-1.8) and for more distal HNC anatomic sites (hypopharynx and larynx). The risk was also higher, or limited to, in subjects exposed to tobacco. The OR rose to 7.2 (95% CI 5.5-9.5) among subjects with family history, who were alcohol and tobacco users. A weak but significant association (OR=1.1, 95% CI 1.0-1.2) emerged for family history of other tobacco-related neoplasms, particularly with laryngeal cancer (OR=1.3, 95% CI 1.1-1.5). No association was observed for family history of nontobacco-related neoplasms and the risk of HNC (OR=1.0, 95% CI 0.9-1.1). Familial factors play a role in the etiology of HNC. In both subjects with and without family history of HNC, avoidance of tobacco and alcohol exposure may be the best way to avoid HNC.

Quantitative analysis of postnatal neurogenesis and neuron number in the macaque monkey dentate gyrus

The dentate gyrus is one of only two regions of the mammalian brain where substantial neurogenesis occurs postnatally. However, detailed quantitative information about the postnatal structural maturation of the primate dentate gyrus is meager. We performed design-based, stereological studies of neuron number and size, and volume of the dentate gyrus layers in rhesus macaque monkeys (Macaca mulatta) of different postnatal ages. We found that about 40% of the total number of granule cells observed in mature 5-10-year-old macaque monkeys are added to the granule cell layer postnatally; 25% of these neurons are added within the first three postnatal months. Accordingly, cell proliferation and neurogenesis within the dentate gyrus peak within the first 3 months after birth and remain at an intermediate level between 3 months and at least 1 year of age. Although granule cell bodies undergo their largest increase in size during the first year of life, cell size and the volume of the three layers of the dentate gyrus (i.e. the molecular, granule cell and polymorphic layers) continue to increase beyond 1 year of age. Moreover, the different layers of the dentate gyrus exhibit distinct volumetric changes during postnatal development. Finally, we observe significant levels of cell proliferation, neurogenesis and cell death in the context of an overall stable number of granule cells in mature 5-10-year-old monkeys. These data identify an extended developmental period during which neurogenesis might be modulated to significantly impact the structure and function of the dentate gyrus in adulthood.

An ecomorphological analysis of the determinants of mating success

Investigating the factors affecting the strength of sexual selection is important for understanding the evolution of sex-specific morphological and behavioural traits. Traditionally, sexual selection studies focus on male ornaments, although recent evidence indicates that sexual selection mechanisms also target organismal performance. In the present study, we investigated the role of sexually dimorphic morphological and performance traits of the common (viviparous) lizard (Zootoca vivipara, Jacquin 1787) with respect to determining mating behaviour. Using an experimental set-up controlling for size differences, we found that males with longer tails had a higher probability of mating a female. Unexpectedly, males with lower bite forces had an advantage over males with higher bite forces, whereas males with bigger heads copulated for a longer time with the female. This shows that predicting mating success is not straightforward and is sometimes counterintuitive because a longer tail appears to be beneficial, whereas biting harder is not, for male Z. vivipara in a male-female interaction context

Sequential analysis of trans-SNARE formation in intracellular membrane fusion.

SNARE complexes are required for membrane fusion in the endomembrane system. They contain coiled-coil bundles of four helices, three (Q(a), Q(b), and Q(c)) from target (t)-SNAREs and one (R) from the vesicular (v)-SNARE. NSF/Sec18 disrupts these cis-SNARE complexes, allowing reassembly of their subunits into trans-SNARE complexes and subsequent fusion. Studying these reactions in native yeast vacuoles, we found that NSF/Sec18 activates the vacuolar cis-SNARE complex by selectively displacing the vacuolar Q(a) SNARE, leaving behind a Q(bc)R subcomplex. This subcomplex serves as an acceptor for a Q(a) SNARE from the opposite membrane, leading to Q(a)-Q(bc)R trans-complexes. Activity tests of vacuoles with diagnostic distributions of inactivating mutations over the two fusion partners confirm that this distribution accounts for a major share of the fusion activity. The persistence of the Q(bc)R cis-complex and the formation of the Q(a)-Q(bc)R trans-complex are both sensitive to the Rab-GTPase inhibitor, GDI, and to mutations in the vacuolar tether complex, HOPS (HOmotypic fusion and vacuolar Protein Sorting complex). This suggests that the vacuolar Rab-GTPase, Ypt7, and HOPS restrict cis-SNARE disassembly and thereby bias trans-SNARE assembly into a preferred topology.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

A novel method for automated classification of epileptiform activity in the human electroencephalogram-based on independent component analysis.

Diagnosis of several neurological disorders is based on the detection of typical pathological patterns in the electroencephalogram (EEG). This is a time-consuming task requiring significant training and experience. Automatic detection of these EEG patterns would greatly assist in quantitative analysis and interpretation. We present a method, which allows automatic detection of epileptiform events and discrimination of them from eye blinks, and is based on features derived using a novel application of independent component analysis. The algorithm was trained and cross validated using seven EEGs with epileptiform activity. For epileptiform events with compensation for eyeblinks, the sensitivity was 65 +/- 22% at a specificity of 86 +/- 7% (mean +/- SD). With feature extraction by PCA or classification of raw data, specificity reduced to 76 and 74%, respectively, for the same sensitivity. On exactly the same data, the commercially available software Reveal had a maximum sensitivity of 30% and concurrent specificity of 77%. Our algorithm performed well at detecting epileptiform events in this preliminary test and offers a flexible tool that is intended to be generalized to the simultaneous classification of many waveforms in the EEG.

Epithelioid sarcoma: a clinicopathologic and immunohistochemical analysis of 106 cases from the French sarcoma group.

Epithelioid sarcoma (ES) is rare with a poor prognosis and for which a loss of INI1 expression has been recently reported. We report a study of 106 cases with clinical, histologic, and immunohistochemical data, including INI1 expression, and follow-up data. Of the 106 cases, 70 were the conventional subtype and 36 the large cell subtype. INI1 was negative in 86 cases (81.1%): 57 (81%) of 70 conventional and 29 (81%) of 36 large cell subtypes. Treatment modalities were available for 76 and follow-up for 80 patients. Of the 80 patients, 43 (54%) experienced metastasis and 25 (31%) died of the disease. Univariate analysis indicated that tumor size and mitotic index were significant for metastasis-free survival, whereas proximal location, tumor size, tumor multifocality, and mitotic index were significant for overall survival. Loss of expression of INI1 is frequent in the conventional and large cell subtypes of ES and can be used as a diagnostic marker, but it has no prognostic impact.

Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.

A 2D/3D image analysis system to track fluorescently labeled structures in rod-shaped cells: application to measure spindle pole asymmetry during mitosis.

BACKGROUND: The yeast Schizosaccharomyces pombe is frequently used as a model for studying the cell cycle. The cells are rod-shaped and divide by medial fission. The process of cell division, or cytokinesis, is controlled by a network of signaling proteins called the Septation Initiation Network (SIN); SIN proteins associate with the SPBs during nuclear division (mitosis). Some SIN proteins associate with both SPBs early in mitosis, and then display strongly asymmetric signal intensity at the SPBs in late mitosis, just before cytokinesis. This asymmetry is thought to be important for correct regulation of SIN signaling, and coordination of cytokinesis and mitosis. In order to study the dynamics of organelles or large protein complexes such as the spindle pole body (SPB), which have been labeled with a fluorescent protein tag in living cells, a number of the image analysis problems must be solved; the cell outline must be detected automatically, and the position and signal intensity associated with the structures of interest within the cell must be determined. RESULTS: We present a new 2D and 3D image analysis system that permits versatile and robust analysis of motile, fluorescently labeled structures in rod-shaped cells. We have designed an image analysis system that we have implemented as a user-friendly software package allowing the fast and robust image-analysis of large numbers of rod-shaped cells. We have developed new robust algorithms, which we combined with existing methodologies to facilitate fast and accurate analysis. Our software permits the detection and segmentation of rod-shaped cells in either static or dynamic (i.e. time lapse) multi-channel images. It enables tracking of two structures (for example SPBs) in two different image channels. For 2D or 3D static images, the locations of the structures are identified, and then intensity values are extracted together with several quantitative parameters, such as length, width, cell orientation, background fluorescence and the distance between the structures of interest. Furthermore, two kinds of kymographs of the tracked structures can be established, one representing the migration with respect to their relative position, the other representing their individual trajectories inside the cell. This software package, called "RodCellJ", allowed us to analyze a large number of S. pombe cells to understand the rules that govern SIN protein asymmetry. CONCLUSIONS: "RodCell" is freely available to the community as a package of several ImageJ plugins to simultaneously analyze the behavior of a large number of rod-shaped cells in an extensive manner. The integration of different image-processing techniques in a single package, as well as the development of novel algorithms does not only allow to speed up the analysis with respect to the usage of existing tools, but also accounts for higher accuracy. Its utility was demonstrated on both 2D and 3D static and dynamic images to study the septation initiation network of the yeast Schizosaccharomyces pombe. More generally, it can be used in any kind of biological context where fluorescent-protein labeled structures need to be analyzed in rod-shaped cells. AVAILABILITY: RodCellJ is freely available under http://bigwww.epfl.ch/algorithms.html, (after acceptance of the publication).

FGF21 signalling pathway and metabolic traits - genetic association analysis.

Fibroblast growth factor 21 (FGF21) is a novel master regulator of metabolic profile. The biological actions of FGF21 are elicited upon its klotho beta (KLB)-facilitated binding to FGF receptor 1 (FGFR1), FGFR2 and FGFR3. We hypothesised that common polymorphisms in the FGF21 signalling pathway may be associated with metabolic risk. At the screening stage, we examined associations between 63 common single-nucleotide polymorphisms (SNPs) in five genes of this pathway (FGF21, KLB, FGFR1, FGFR2, FGFR3) and four metabolic phenotypes (LDL cholesterol - LDL-C, HDL-cholesterol - HDL-C, triglycerides and body mass index) in 629 individuals from Silesian Hypertension Study (SHS). Replication analyses were performed in 5478 unrelated individuals of the Swiss CoLaus cohort (imputed genotypes) and in 3030 directly genotyped individuals of the German Myocardial Infarction Family Study (GerMIFS). Of 54 SNPs that met quality control criteria after genotyping in SHS, 4 (rs4733946 and rs7012413 in FGFR1; rs2071616 in FGFR2 and rs7670903 in KLB) showed suggestive association with LDL-C (P=0.0006, P=0.0013, P=0.0055, P=0.011, respectively) and 1 (rs2608819 in KLB) was associated with body mass index (P=0.011); all with false discovery rate q<0.5. Of these, only one FGFR2 polymorphism (rs2071616) showed replicated association with LDL-C in both CoLaus (P=0.009) and men from GerMIFS (P=0.017). The direction of allelic effect of rs2071616 upon LDL-C was consistent in all examined populations. These data show that common genetic variations in FGFR2 may be associated with LDL-C in subjects of white European ancestry.

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.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.