Chromosomal rearrangements and gene flow over time in an inter-specific hybrid zone of the Sorex araneus group.

Most hybrid zones have existed for hundreds or thousands of years but have generally been observed for only a short time period. Studies extending over periods long enough to track evolutionary changes in the zones or assess the ultimate outcome of hybridization are scarce. Here, we describe the evolution over time of the level of genetic isolation between two karyotypically different species of shrews (Sorex araneus and Sorex antinorii) at a hybrid zone located in the Swiss Alps. We first evaluated hybrid zone movement by contrasting patterns of gene flow and changes in cline parameters (centre and width) using 24 microsatellite loci, between two periods separated by 10 years apart. Additionally, we tested the role of chromosomal rearrangements on gene flow by analysing microsatellite loci located on both rearranged and common chromosomes to both species. We did not detect any movement of the hybrid zone during the period analysed, suggesting that the zone is a typical tension zone. However, the gene flow was significantly lower among the rearranged than the common chromosomes for the second period, whereas the difference was only marginally significant for the first period. This further supports the role of chromosomal rearrangements on gene flow between these taxa.

Real-time recording of circadian liver gene expression in freely moving mice reveals the phase-setting behavior of hepatocyte clocks.

The mammalian circadian timing system consists of a master pacemaker in the suprachiasmatic nucleus (SCN) in the hypothalamus, which is thought to set the phase of slave oscillators in virtually all body cells. However, due to the lack of appropriate in vivo recording technologies, it has been difficult to study how the SCN synchronizes oscillators in peripheral tissues. Here we describe the real-time recording of bioluminescence emitted by hepatocytes expressing circadian luciferase reporter genes in freely moving mice. The technology employs a device dubbed RT-Biolumicorder, which consists of a cylindrical cage with reflecting conical walls that channel photons toward a photomultiplier tube. The monitoring of circadian liver gene expression revealed that hepatocyte oscillators of SCN-lesioned mice synchronized more rapidly to feeding cycles than hepatocyte clocks of intact mice. Hence, the SCN uses signaling pathways that counteract those of feeding rhythms when their phase is in conflict with its own phase.

Regulation of the cardiac voltage-gated sodium channel Nav1.5 : regulation of Nav1.5 by ubiquitin-protein ligases of the Nedd4/Nedd4-like family and characterisation of two novel mutations in the gene encoding Nav1.5 (SCN5A) implicated in the Brugada syndrome triggered by fever

Abstract: The genesis of the cardiac action potential, which accounts for the cardiac contraction, is due to the sodium current INa mediated by the voltage-gated sodium channel Nav1.5. Several cardiac arrhythmias such as the Brugada syndrome are known te be caused by mutations in SCN5A, the gene encoding Nav1.5. Studies of these mutations allowed a better understanding of biophysical and functional properties of Nav1.5. However, only few investigations have been performed in order to understand the regulation of Nav1.5. During my thesis, I investigated different mechanisms of regulation of Nav1.5 using a heterologous expression system, HEK293 cells, coupled with a technique of sodium current recording: the patch clamp in whole cell configuration. In previous studies it has been shown that an enzyme of the Nedd4 family (Nedd4-2) regulates an epithelial sodium channel via the interaction with PY-motifs present in the latter. Interestingly, Nav1.5 contains a similar PY-motif, which motivated us to study the role of Nedd4-2 expressed in heart for the regulation of Nav1.5. In a second study, we investigated the implication of two Nav1.5 mutants, which were either less functional or net functional (Nav1.5 R535X and Nav1.5 L325R respectively) implied in the genesis of the Brugada syndrome by fever. Our results established two mechanisms implied in Nav1.5 regulation. The first one implies that following the interaction between the PY-motif of Nav1.5 and Nedd4- 2 Nav1.5 is ubiquitinated by Nedd4-2. This ubiquitination leads to the internalization of Nav1 .5. The second mechanism is a phenomenon called the "dominant negative" effect of Nav1.5 L325R on Nay1.5 where the decrease of 'Na is potentially due to the retention of Nav1.5 by Nav1.5 L325R in an undefined intracellular compartment. These studies defined two mechanisms of Nav1.5 regulation, which could play an important role for the genesis of cardiac arrhythmias where molecular processes are still poorly understood. Résumé La genèse du potentiel d'action cardiaque, permettant la contraction cardiaque, est due au courant sodique INa issu des canaux sodiques cardiaques dépendants du voltage Nav1.5. Nombreuses arythmies cardiaques telles que le syndrome de Brugada sont connues pour être liées à des mutations du gène SCN5A, codant pour Nav1.5. L'étude de ces mutations a permis une meilleure compréhension des propriétés structurelles et fonctionnelles de Nav1.5 et leurs implications dans la genèse de ces pathologies. Néanmoins peu d'études ont été menées afin de comprendre les mécanismes de régulation de Nav1.5. Mon travail de thèse a consisté à étudier des mécanismes de régulation de Nav1.5 en utilisant un système d'expression hétérologue, les cellules HEK293, couplé à une technique d'enregistrement des courants sodiques, le "patch clamp" en configuration cellule entière. La présence sur Nav1.5 d'un motif-PY similaire à ceux nécessaires pour la régulation d'un canal épithélial sodique par une enzyme de la famille de Nedd4, nous a amenée à étudier le rôle de ces ubiquitine-ligases, en particulier Nedd4-2, dans la régulation de Nav1.5. La seconde étude s'est intéressée aux conséquences de deux mutations de SCN5A codant pour deux mutants peu ou pas fonctionnels (Nav1.5 L325R et Nav1.5 R535X respectivement) retrouvées chez des patients présentant un syndrome de Brugada exacerbé par un état fébrile. Nos résultats ont permis d'établir deux mécanismes de régulation de Nav1.5 L'un par Nedd4-2 qui implique rubiquitination de Nav1.5 par cette ligase suite à l'interaction entre le motif-PY de Nav1.5 et Nedd4-2. Cette modification déclenche l'internalisation du canal impliquée dans la diminution d'INa. Le second mécanisme quant à lui est un effet "dominant négatif" de Nav1.5 L325R sur Nav1.5 aboutissant à une diminution d'INa suite à la séquestration intracellulaire potentielle de Nav1.5 par Nav1.5 L325R. Ces études ont mis en évidence deux mécanismes de régulation de Nav1.5 pouvant jouer un rôle majeur dans la genèse et/ou l'accentuation des arythmies cardiaques dont les processus moléculaires au sein des cardiomyocytes, impliquant des modifications du courant sodiques, sont encore mal compris. Résumé destiné à un large public La dépolarisation électrique de la membrane des cellules cardiaques permet la contraction du coeur. La génèse de cette activité électrique est due au courant sodique issu d'un type de canal à sodium situé dans la membrane des cellules cardiaques. De nombreuses pathologies provoquant des troubles du rythme cardiaque sont issues de mutations du gène qui code pour ce canal à sodium. Ces canaux mutants, entrainant diverses pathologies cardiaques telles que le syndrome de Brugada, ont été largement étudiées. Néanmoins, peu de travaux ont été réalisés sur les mécanismes de régulation de ce canal à sodium non muté. Mon travail de thèse a consisté à étudier certains des mécanismes de régulation de ce canal à sodium en utilisant une technique permettant l'enregistrement des courants sodiques issus de l'expression de ces canaux à sodium à la membrane de cellules mammifères. La présence sur ce canal à sodium d'une structure spécifique, similaire à celle nécessaire pour la régulation d'un canal épithélial à sodium par une enzyme appelée Nedd4-2, nous a amenée à étudier le rôle de cette enzyme dans la régulation de ce canal à sodium. La seconde étude s'est intéressée aux rôles de deux mutations du gène codant pour ce canal à sodium retrouvées chez des patients présentant un syndrome de Brugada exacerbé par la fièvre. Nos résultats nous ont permis d'établir deux mécanismes de régulation de ce canal à sodium diminuant le courant sodique l'un par l'action de l'enzyme Nedd4-2, suite à son interaction avec ce canal, qui modifie ce canal à sodium (ubiquitination) diminuant de ce fait la densité membranaire du canal. L'autre par un mécanisme suggérant un effet négatif de l'un des canaux mutants sur l'expression à la membrane du canal à sodium non muté. Ces études ont mis en évidence deux mécanismes de régulation de ce canal à sodium pouvant jouer un rôle majeur dans la genèse et/ou l'accentuation des troubles du rythme cardiaques dont les mécanismes cellulaires sont encore incompris.

Gene copy number polymorphisms in an arbuscular mycorrhizal fungal population.

Gene copy number polymorphism was studied in a population of the arbuscular mycorrhizal fungus Glomus intraradices by using a quantitative PCR approach on four different genomic regions. Variation in gene copy number was found for a pseudogene and for three ribosomal genes, providing conclusive evidence for a widespread occurrence of macromutational events in the population.

Equine herpesvirus-2 E10 gene product, but not its cellular homologue, activates NF-kappaB transcription factor and c-Jun N-terminal kinase.

We have previously reported on the death effector domain containing E8 gene product from equine herpesvirus-2, designated FLICE inhibitory protein (v-FLIP), and on its cellular homologue, c-FLIP, which inhibit the activation of caspase-8 by death receptors. Here we report on the structure and function of the E10 gene product of equine herpesvirus-2, designated v-CARMEN, and on its cellular homologue, c-CARMEN, which contain a caspase-recruiting domain (CARD) motif. c-CARMEN is highly homologous to the viral protein in its N-terminal CARD motif but differs in its C-terminal extension. v-CARMEN and c-CARMEN interact directly in a CARD-dependent manner yet reveal different binding specificities toward members of the tumor necrosis factor receptor-associated factor (TRAF) family. v-CARMEN binds to TRAF6 and weakly to TRAF3 and, upon overexpression, potently induces the c-Jun N-terminal kinase (JNK), p38, and nuclear factor (NF)-kappaB transcriptional pathways. c-CARMEN or truncated versions thereof do not appear to induce JNK and NF-kappaB activation by themselves, nor do they affect the JNK and NF-kappaB activating potential of v-CARMEN. Thus, in contrast to the cellular homologue, v-CARMEN may have additional properties in its unique C terminus that allow for an autonomous activator effect on NF-kappaB and JNK. Through activation of NF-kappaB, v-CARMEN may regulate the expression of the cellular and viral genes important for viral replication.

Phenotypic and molecular characterization of Bruck syndrome (Osteogenesis imperfecta with contractures of the large joints) caused by a recessive mutation in PLOD2

Le Syndrome de Bruck (Bruck Syndrome; BS) est une maladie autosomique récessive assemblant la combinaison inhabituelle de fragilité osseuse semblable à celle de l'Ostéogenèse Imparfaite (0I) avec des contractures congénitales tendineuses et cutanées des grandes articulations («ptérygia»). Les cas décrits jusqu'à ce jour mettent en évidence une grande hétérogénéité du tableau clinique, liée en partie au manque d'un diagnostic biochimique ou moléculaire. Nous savons que dans le BS les gènes codant pour le collagène 1 ne sont pas mutés, mais savons néanmoins, grâce à l'étude du collagène extrait de biopsies osseuses, qu'il y a un déficit d'hydroxylation des résidus de lysine dans les télopeptides du collagène 1 qui servent à la formation des liens intermoléculaires (crosslinks) et donc à la stabilisation des fibres de collagène. Un locus génétique du BS à été mappé sur 17q12, mais le gène responsable sur ce locus reste inconnu; plus récemment, deux mutations dans le gène de la lysyl hydroxylase 2 (PLOD2, position chromosomique 3q23-q24) ont été identifiées, démontrant l'hétérogénéité génétique du ES. La proportion de ES liée à 17p22 (BS type 1) et celle liée à une mutation dans PLOD2 (BS type 2) est encore incertaine et nous manquons de données sur la corrélation phenotype-génotype. Nous avons étudié le cas d'un garçon avec des contractures et des ptérygia dès la naissance, combinées à une ostéopénie sévère de type OI menant à des fractures multiples. Ses urines contenaient une quantité élevée d'hydroxyproline, indiquant un remaniement important du tissu osseux, mais peu de produits de dégradation des crosslinks du collagène, indiquant donc une réduction de la proportion de crosslinks dans le collagène in vivo. Nous avons pu démontrer chez lui la présence d'une nouvelle mutation homozygote dans le gène PLOD2 menant à une substitution Arg598His; les deux parents du sujet étaient hétérozygotes pour la mutation et celle-ci était absente dans notre population témoin. La mutation est adjacente aux deux mutations rapportées précédemment (Gly601Val et Thr608Ile), ce qui suggère la présence d'un ''hotspot'' mutationnel mais aussi d'une région de grande importance fonctionnelle sur PLOD2 : cette observation est importante pour la création d'inhibiteurs de PLOD2, recherchés en ce moment pour le traitement de la fibrose. La combinaison de ptérygia et de fragilité osseuse, comme illustrée par notre patient est apparemment contradictoire et donc difficilement explicable mais indique que l'hydroxylation des résidus lysyl des télopeptides est importante non seulement pour la stabilité osseuse mais aussi dans la morphogénèse et la formation des articulations dans la période prénatale. Finalement, la mesure des produits de dégradation du collagène dans l'urine et l'analyse de mutation de PLOD2 permet le diagnostic du syndrome de Bruck et permet de le différencier de l'Osteogénèse Imparfaite. -- Bruck syndrome (BS) is a recessively-inherited phenotypic disorder featuring the unusual combination of skeletal changes resembling osteogenesis imperfecta (0I) with congenital contractures of the large joints. Clinical heterogeneity is apparent in cases reported thus far. While the genes coding for collagen 1 chains are unaffected in BS, there is biochemical evidence for a defect in the hydroxylation of lysine residues in collagen 1 telopeptides. One BS locus has been mapped at 17p12, but more recently, two mutations in the lysyl hydroxylase 2 gene (PLOD2, 3q23-q24) have been identified in BS, showing genetic heterogeneity. The proportion of BS cases linked to 17p22 (BS type 1) or caused by mutations in PLOD2 (BS type 2) is still uncertain, and phenotypic correlations are lacking. We report on a boy who had congenital contractures with pterygia at birth and severe 0I-like osteopenia and multiple frac-tures. His urine contained high amounts of hydroxyproline but low amounts of collagen crosslinks degradation products; and he was shown to be homozygous for a novel mutation leading to an Arg598His substitution in PLOD2. The mutation is adjacent to the two mutations previously reported (Gly601Val and Thr608Ile), suggesting a functionally important hotspot in PLOD2. The combination of pterygia with bone fragility, as illustrated by this case, is difficult to explain; it suggests that telopeptide lysyl hydroxylation must be involved in prenatal joint formation and morphogenesis. Collagen degradation products in urine and mutation analysis ofPLOD2 maybe used to diagnose BS and differentiate it from M.

In vivo enzymatic activity of acetylCoA synthetase in skeletal muscle revealed by 13C turnover from hyperpolarized [1-13C]acetate to [1-13C]acetylcarnitine

BACKGROUND: Acetate metabolism in skeletal muscle is regulated by acetylCoA synthetase (ACS). The main function of ACS is to provide cells with acetylCoA, a key molecule for numerous metabolic pathways including fatty acid and cholesterol synthesis and the Krebs cycle. METHODS: Hyperpolarized [1-(13)C]acetate prepared via dissolution dynamic nuclear polarization was injected intravenously at different concentrations into rats. The (13)C magnetic resonance signals of [1-(13)C]acetate and [1-(13)C]acetylcarnitine were recorded in vivo for 1min. The kinetic rate constants related to the transformation of acetate into acetylcarnitine were deduced from the 3s time resolution measurements using two approaches, either mathematical modeling or relative metabolite ratios. RESULTS: Although separated by two biochemical transformations, a kinetic analysis of the (13)C label flow from [1-(13)C]acetate to [1-(13)C]acetylcarnitine led to a unique determination of the activity of ACS. The in vivo Michaelis constants for ACS were KM=0.35±0.13mM and Vmax=0.199±0.031μmol/g/min. CONCLUSIONS: The conversion rates from hyperpolarized acetate into acetylcarnitine were quantified in vivo and, although separated by two enzymatic reactions, these rates uniquely defined the activity of ACS. The conversion rates associated with ACS were obtained using two analytical approaches, both methods yielding similar results. GENERAL SIGNIFICANCE: This study demonstrates the feasibility of directly measuring ACS activity in vivo and, since the activity of ACS can be affected by various pathological states such as cancer or diabetes, the proposed method could be used to non-invasively probe metabolic signatures of ACS in diseased tissue.

Identification of a novel splice-site mutation in the CYP1A2 gene.

AIMS: To identify the molecular basis for a low CYP1A2 metabolic status, as determined by a caffeine phenotyping test, in a 71-year-old, nonsmoking, Caucasian woman who presented with very high clozapine concentrations despite being administered a standard dose of the drug. METHODS: The nucleotide sequence of the 7 exons, exon-intron boundaries and 5'-flanking region of the CYP1A2 gene was analysed by direct sequencing. RESULTS: Only one heterozygous point mutation was identified in the donor splice site of intron 6 (3534G > A) of CYP1A2. This mutation could cause abnormal RNA splicing and therefore lead to a truncated nonfunctional enzyme. No other carrier of this mutation was identified in a population of 100 unrelated healthy Caucasians. CONCLUSIONS: This is the first report of a splice-site mutation affecting the CYP1A2 gene. This polymorphism is a likely explanation for the low CYP1A2 activity associated with high clozapine concentrations in this patient.

Recurrent mutations in the CDKL5 gene: Genotype-phenotype relationships.

Mutations in the cyclin-dependent kinase-like 5 gene (CDKL5) have been described in epileptic encephalopathies in females with infantile spasms with features that overlap with Rett syndrome. With more than 80 reported patients, the phenotype of CDKL5-related encephalopathy is well-defined. The main features consist of seizures starting before 6 months of age, severe intellectual disability with absent speech and hand stereotypies and deceleration of head growth, which resembles Rett syndrome. However, some clinical discrepancies suggested the influence of genetics and/or environmental factors. No genotype-phenotype correlation has been defined and thus there is a need to examine individual mutations. In this study, we analyzed eight recurrent CDKL5 mutations to test whether the clinical phenotype of patients with the same mutation is similar and whether patients with specific CDKL5 mutations have a milder phenotype than those with other CDKL5 mutations. Patients bearing missense mutations in the ATP binding site such as the p.Ala40Val mutation typically walked unaided, had normocephaly, better hand use ability, and less frequent refractory epilepsy when compared to girls with other CDKL5 mutations. In contrast, patients with mutations in the kinase domain (such as p.Arg59X, p.Arg134X, p.Arg178Trp/Pro/Gln, or c.145 + 2T > C) and frameshift mutations in the C-terminal region (such as c.2635_2636delCT) had a more severe phenotype with infantile spasms, refractory epileptic encephalopathy, absolute microcephaly, and inability to walk. It is important for clinicians to have this information when such patients are diagnosed. © 2012 Wiley Periodicals, Inc.

Effects of ploidy and sex-locus genotype on gene expression patterns in the fire ant Solenopsis invicta.

Males in many animal species differ greatly from females in morphology, physiology and behaviour. Ants, bees and wasps have a haplodiploid mechanism of sex determination whereby unfertilized eggs become males while fertilized eggs become females. However, many species also have a low frequency of diploid males, which are thought to develop from diploid eggs when individuals are homozygous at one or more sex determination loci. Diploid males are morphologically similar to haploids, though often larger and typically sterile. To determine how ploidy level and sex-locus genotype affect gene expression during development, we compared expression patterns between diploid males, haploid males and females (queens) at three developmental timepoints in Solenopsis invicta. In pupae, gene expression profiles of diploid males were very different from those of haploid males but nearly identical to those of queens. An unexpected shift in expression patterns emerged soon after adult eclosion, with diploid male patterns diverging from those of queens to resemble those of haploid males, a pattern retained in older adults. The finding that ploidy level effects on early gene expression override sex effects (including genes implicated in sperm production and pheromone production/perception) may explain diploid male sterility and lack of worker discrimination against them during development.

Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis.

BACKGROUND: Histologic grade in breast cancer provides clinically important prognostic information. However, 30%-60% of tumors are classified as histologic grade 2. This grade is associated with an intermediate risk of recurrence and is thus not informative for clinical decision making. We examined whether histologic grade was associated with gene expression profiles of breast cancers and whether such profiles could be used to improve histologic grading. METHODS: We analyzed microarray data from 189 invasive breast carcinomas and from three published gene expression datasets from breast carcinomas. We identified differentially expressed genes in a training set of 64 estrogen receptor (ER)-positive tumor samples by comparing expression profiles between histologic grade 3 tumors and histologic grade 1 tumors and used the expression of these genes to define the gene expression grade index. Data from 597 independent tumors were used to evaluate the association between relapse-free survival and the gene expression grade index in a Kaplan-Meier analysis. All statistical tests were two-sided. RESULTS: We identified 97 genes in our training set that were associated with histologic grade; most of these genes were involved in cell cycle regulation and proliferation. In validation datasets, the gene expression grade index was strongly associated with histologic grade 1 and 3 status; however, among histologic grade 2 tumors, the index spanned the values for histologic grade 1-3 tumors. Among patients with histologic grade 2 tumors, a high gene expression grade index was associated with a higher risk of recurrence than a low gene expression grade index (hazard ratio = 3.61, 95% confidence interval = 2.25 to 5.78; P < .001, log-rank test). CONCLUSIONS: Gene expression grade index appeared to reclassify patients with histologic grade 2 tumors into two groups with high versus low risks of recurrence. This approach may improve the accuracy of tumor grading and thus its prognostic value.

Genetic variability in Arbuscular Mycorrhizal Fungi : effect on gene transcription of "Oryza Sativa"

AbstractArbuscular Mycorrhizal Fungi (AMF) form obligate symbioses with the majority of land plants. These fungi influence the diversity and productivity of plants. AMF are unusual organisms, harbouring genetically different nuclei in a common cytoplasm (known as heterokaryosis). Genetic variability has been shown between AMF individuals coming from the same population. Recent findings showed that genetic exchange between genetically different AMF individuals was possible. Additionnaly, segregation was shown to occur at spore formation in AMF. These two processes were shown to increase genetic variability between AMF individuals.Because of the difficulty to study these organisms, almost nothing is known about the effect of intra-specific genetic variability in AMF on the plant transcriptome. The aim of this thesis was to bring insights into the effect of intra-specific genetic variability in AMF on plant gene transcription. We demonstrated that genetic exchange could influence expression of some symbiosis specific plant genes and the timing of the colonization of the fungi in plant roots. We also showed that segregation could have a large impact on plant gene transcription. Taken together, these results demonstrated that AMF intra-specific variability could profoundly affect the life of plants by altering various molecular pathways. Moreover, results obtained on rice open a field of research on AMF genetics in impromvment of growth in agricultural plants and should be taken into account for future experiments.RésuméLes champignons endomycorhiziens arbusculaires (CEA) forment une symbiose obligatoire avec la majorité des plantes sur terre. Ces champignons peuvent influencer la diversité et la productivité des plantes avec lesquelles ils forment la symbiose. Les CEA sont des organismes particuliers de part le fait qu'ils possèdent des noyaux génétiquement différents (appelés hétérocaryosis) dans un cytoplasme commun. Il a été montré qu'il existait de la variabilité génétique intra-specific chez les CEA. De plus, des études recentes ont montré que l'échange génétique chez les CEA était possible entre des individus génétiquement différents tout comme la ségrégation qui a aussi été démontrée au moment de la formation des nouvelles spores chez les CEA. Ces deux processus ont été montrés comme pouvant créer aussi de la variabilité génétique intra-specific.Du fait de la difficulté de travailler avec les CEA et à cause de la nouveauté de ces recherches, très peu de choses sont connues sur l'effet de l'échange génétique et de la ségrégation chez les CEA sur les plantes, et particulièrement au niveau moléculaire. Le but de cette thèse a été d'apporter la lumière sur les effets de la viariabilité génétique intra-specific chez les CEA, sur la transcription des gènes chez la plante. Nous avons pu montrer que l'échange génétique pouvait avoir des effets sur l'expression de gènes spécifiques à cette symbiose mais aussi pouvait influencer le timing de colonisation des racines de plantes par les CEA. Nous avons aussi montré que la ségrégation pouvait grandement influencer le transcriptome complet de la plante, et pas seulement les voies métaboliques spécifiques à la symbiose comme cela avait été montré auparavant.L'ensemble de ces résultats démontre l'importance de la variation intra-specific chez les CEA sur les plantes et leur implication sur leur cycle de vie en changeant l'expression de voies métaboliques. De plus, ces résultats obtenus sur le riz ouvrent un champ de recherches sur les plantes destinées à l'agriculture et devraient être pris en compte pour des expériences futures.

Additive functions in boolean models of gene regulatory network modules.

Gene-on-gene regulations are key components of every living organism. Dynamical abstract models of genetic regulatory networks help explain the genome's evolvability and robustness. These properties can be attributed to the structural topology of the graph formed by genes, as vertices, and regulatory interactions, as edges. Moreover, the actual gene interaction of each gene is believed to play a key role in the stability of the structure. With advances in biology, some effort was deployed to develop update functions in Boolean models that include recent knowledge. We combine real-life gene interaction networks with novel update functions in a Boolean model. We use two sub-networks of biological organisms, the yeast cell-cycle and the mouse embryonic stem cell, as topological support for our system. On these structures, we substitute the original random update functions by a novel threshold-based dynamic function in which the promoting and repressing effect of each interaction is considered. We use a third real-life regulatory network, along with its inferred Boolean update functions to validate the proposed update function. Results of this validation hint to increased biological plausibility of the threshold-based function. To investigate the dynamical behavior of this new model, we visualized the phase transition between order and chaos into the critical regime using Derrida plots. We complement the qualitative nature of Derrida plots with an alternative measure, the criticality distance, that also allows to discriminate between regimes in a quantitative way. Simulation on both real-life genetic regulatory networks show that there exists a set of parameters that allows the systems to operate in the critical region. This new model includes experimentally derived biological information and recent discoveries, which makes it potentially useful to guide experimental research. The update function confers additional realism to the model, while reducing the complexity and solution space, thus making it easier to investigate.