Early antibiotic administration affects the gut barrier function and the immune response to oral antigen in suckling rats

Abstract: The increasingly high hygienic standards characterizing westernized societies correlate with an increasingly high prevalence of allergic disease. Initially based on these observations, the hygiene hypothesis postulates that reduced microbial stimulation during infancy impairs the immune system development and increases the risk of allergy. Moreover, there is increasing evidence that the crosstalk existing between the intestine and the resident microbiota is crucial for gut homeostasis. In particular, bacterial colonization of the gut affects the integrity of the gut barrier and stimulates the development of the gut associated immune tissue, both phenomena being essential for the immune system to mount a controlled response to food antigens. Therefore, alterations in the microbial colonization process, by compromising the barrier homeostasis, may increase the risk of food allergy. In this context, antibiotic treatment, frequently prescribed during infancy, affects gut colonization by bacteria. However, little is known about the impact of alterations in the colonization process on the maturation of the gut barrier and on the immunological response to oral antigens. The objective of this work was to determine the impact of a commercial antibiotic preparation employed in pediatric settings on the gut barrier status at the critical period of the suckling/weaning transition and to evaluate the physiological consequences of this treatment in terms of immune response to food antigens. We established an antibiotic-treated suckling rat model relevant to the pediatric population in terms of type, dose and route of administration of the antibiotic and of changes in the patterns of microbial colonization. Oral tolerance to a novel luminal antigen (ovalbumin) was impaired when the antigen was introduced during antibiotic treatment. These results paralleled to alterations in the intestinal permeability to macromolecules and reduced intestinal expression of genes coding for the major histocomptatibility complex II molecules, which suggest a reduced capacity of antigen handling and presentation in the intestine of the antibiotic-treated animals. In addition, low luminal IgA levels and reduced intestinal expression of genes coding for antimicrobial proteins suggest that protection against pathogens was reduced under antibiotic treatment. In conclusion, we observed in suckling rats that treatment with abroad-spectrum antibiotic commonly used in pediatric practices reduced the capacity of the immune system to develop tolerance. The impact of the antibiotic treatment on the immune response to the antigen-was likely mediated by the alterations of the gut microbiota, through impairment in the mechanisms of antigen handling and presentation. This work reinforces the body of data supporting a key role of the intestinal microbiota modulating the risk of allergy development and leads us to propose that the introduction of new food antigens should be avoided during antibiotic treatment in infants. Résumé: L'augmentation du niveau d'hygiène caractérisant les sociétés occidentales semble être fortement corrélée avec l'augmentation des cas d'allergie dans ces pays. De cette observation est née l'hypothèse qu'une diminution des stimuli microbiens pendant l'enfance modifie le développement du système immunitaire augmentant ainsi le risque d'allergie. En ce sens, un nombre croissant de données indiquent que les interactions existant entre l'intestin et les bactéries résidantes sont cruciales pour l'équilibre du système. En effet, la présence de bactéries dans l'intestin affecte l'intégrité de sa fonction de barrière et stimule le développement du système immunitaire intestinal. Ces deux paramètres étant essentiels à la mise en place d'une réponse contrôlée vis à vis d'un antigène reçu oralement, toute modification du processus naturel de colonisation compromettant l'équilibre intestinal pourrait augmenter le risque d'allergie. Les traitements aux antibiotiques, fréquemment prescrits en pédiatrie, modifient de façon conséquente le processus de colonisation bactérienne. Cependant peu de données existent concernant l'impact d'une altération du processus de colonisation sur la maturation de la barrière intestinale et de la réponse immunitaire dirigée contre un antigène. L'objectif de ce travail était de déterminer l'impact d'un antibiotique commercial et employé en pédiatrie sur l'état de la barrière intestinale au moment critique du sevrage et d'évaluer les conséquences physiologiques d'un tel traitement sur la réponse immune à un antigène alimentaire. Nous avons mis en place un modèle de rats allaités, traités à l'antibiotique, le plus proche possible des pratiques pédiatriques, en terme de nature, dose et voie d'administration de l'antibiotique. Nous avons constaté que l'établissement de la tolérance orale à un nouvel antigène (l'ovalbumine) est altéré quand celui-ci est donné pour la première fois au cours du traitement antibiotique. Ces résultats coïncident avec une diminution de la perméabilité intestinale aux macromolécules, ainsi qu'avec une diminution de l'expression des gènes codant pour les molécules du complexe majeur d'histocomptatibilité de classe II, suggérant une modification de l'apprêtement et de la présentation de l'antigène au niveau intestinal chez les rats traités à l'antibiotique. De plus, un faible taux d'IgA et une diminution de l'expression des gènes codant pour des protéines antimicrobiennes, observés après l'administration d'antibiotique, laissent à penser que la protection contre un pathogène est diminuée lors d'un traitement antibiotique. En conclusion, nous avons observé qu'un traitement antibiotique à large spectre d'activité, couramment utilisé en pédiatrie, réduit la capacité d'induction de la tolérance orale chez le rat allaité. L'impact du traitement antibiotique sur la réponse immune semble induite par l'altération de la flore intestinale via son effet sur les mécanismes d'apprêtement et de présentation de l'antigène. Ce travail renforce l'ensemble des données existantes qui accorde à la flore intestinale un rôle clef dans la modulation du risque de développement d'allergie et nous amène à recommander d'éviter l'introduction d'un nouvel aliment lorsqu'un enfant est traité aux antibiotiques.

Eptingiacea and Saturnaliacea (Radiolaria) from the middle Carnian of Turkey and some late Ladinian to early Norian samples from Oman and Alaska

This article is a taxonomic study of the radiolarian species of the superfamilies Eptingiacea and Saturnaliacea occurring in the middle Carnian fauna from the Koseyahya section, near the town of Elbistan, southeastern Turkey. This fauna is characteristic of the Tetraporobrachia haeckeli Radiolarian Zone as defined in Austria and later found also in Turkey and Oman. It comes from an 8 m thick succession of clayey/cherty limestones from the lower part of the section. In addition, a few species from the late Ladinian and Carnian from Oman and the early Norian from Alaska have also been included in this study, in order to improve some generic diagnoses and to show the diversity and evolutionary trends of some genera. 32 radiolarian species of which 22 are new are described and illustrated, and assigned to 16 genera of which three are new (Capnuchospyris, Veleptingium, and Triassolaguncula). The diagnoses of some species, genera, subfamilies and families have been revised, and the family Eptingiidae has been raised to the rank of superfamily.

Lipoxygenase 2 and the early steps of jasmonic acid biosynthesis in "Arabidopsis thaliana"

Résumé Les oxylipines, telles que l'acide jasmonique (AJ ou jasmonate), jouent un rôle central en réponse à la blessure et à la pathogenèse. De nombreuses études ont montré l'importance de la voie canonique du jasmonate lors de la défense des plantes. De plus, un précurseur cyclopentenone de l'AJ, l'acide oxo-phyto-dienoic (OPDA), a été impliqué comme jouant le rôle d'une molécule signal lors de la défense contre certains pathogènes. En utilisant des mutants bloqués dans la biosynthèse de l'acide jasmonique (aos) ou dans sa perception (coi1-1), nous avons cherché à définir dans quelle mesure l'OPDA joue un rôle de signal induisant l'expression génétique en réponse à la blessure chez Arabidopsis. A l'aide de puces à ADN (microarray), nous avons montré que les transcriptomes d'aos et de coi1-1 sont très semblables après blessure, ce qui suggère que les produits d'AOS sont tous perçus via COI1. Pourtant, lorsqu'on analyse les métabolites présents chez ces mutants, une différence est visible, puisque aos n'accumule pas d'AJ, alors que coi1-1 en accumule encore rapidement après blessure. Nous avons étudié la possibilité qu'un mécanisme de régulation post-traductionnelle sur la voie de biosynthèse du jasmonate explique l'accumulation d'AJ chez coi1-1 après blessure. La lipoxygenase 2 (LOX2) est la première enzyme impliquée dans la biosynthèse de l'AJ et est donc une cible potentielle d'un tel mécanisme. Un indice sur la manière dont l'activité LOX pourrait être régulée vient du mutant fou2 (pour fatty acid oxygenation upregcilated 2) dans lequel l'activité LOX ainsi que le niveau d'AJ sont constitutivement élevés. Cette mutation implique un flux de cation dans la régulation de la production de l'AJ. De plus, il a été montré que plusieurs LOXs, dans des organismes autres que des plantes, peuvent lier le calcium. Nous montrons que l'activité LOX requiert l'addition de cations divalents pour être maximale in vitro, et que non seulement le calcium mais aussi le magnésium joue ce rôle. De plus, nous caractérisons un mutant récessif de LOX2 chez Arabidopsis (lox2-1). Ces plantes sont fertiles, et une analyse quantitative montre qu'elles accumulent toujours un peu d'AJ après blessure. Ceci suggère que LOX2 n'est pas la seule LOX impliquée dans la synthèse d'AJ. Aussi les plantes lox2-1 ne sont pas plus sensibles que les plantes de type sauvage lorsqu'elles sont infectées par la moisissure Botrytis cinerea ou lorsqu'elles sont exposées à un détritivore, néanmoins elles sont plus sensibles lorsqu'elles sont offertes en nourriture à un insecte herbivore. Les insectes et les plantes ont co-évolué conjointement, ainsi une plante ne contenant qu'un niveau réduit d'AJ favorise l'insecte. La disponibilité d'un mutant avec un niveau intermédiaire d'AJ va permettre de mieux comprendre pourquoi les plantes produisent autant de jasmonate. Abstract Oxylipins such as jasmonic acid (JA) play central roles in the wound response and during pathogenesis and many studies have confirmed the important role of the canonical jasmonate pathway in plant defense. Moreover, the cyclopentenone precursor of JA, oxo-phytodienoic acid (OPDA), is also thought to function as a signaling molecule in defense towards some pathogens. Its action was reported to depend on a different signal pathway to JA. By using mutants blocked in the biosynthesis (aos) or perception (coil-1) of JA, we investigated to which extend OPDA works as signaling molecule to trigger gene expression in the wound response of Arabidopsis. Using microarrays, we showed that aos and coil-1 transcriptome are similar in response to wounding, suggesting that products of AOS are all perceived by COI1. However, we found a difference between the two mutants at the metobolomic level, since aos is devoid of JA, but coil-1 can still rapidly accumulate JA upon wounding. We investigated the possibility that the post-translational activation of JA biosynthesis could explain the fast accumulation of JA in coil-1 plants upon wounding. Lipoxygenase (LOX) 2 is the first enzyme implicated in JA synthesis and was thus chosen as a potential target for posttranslational regulation. A clue as to how LOX activity might be regulated came from the fatty acid oxygenation upregulated 2 (foul) mutant in which LOX activity and JA levels are elevated. The foul mutant implicates cations flux in the regulation of JA production, and several LOXs in organisms other than plants have been shown to bind calcium. We showed that Arabidopsis LOX requires divalent cations for full activity in vitro, and that not only calcium but also magnesium can play this role. Moreover, a single recessive mutant of AtLOX2 was characterized. These plants are fully fertile. Quantitative oxylipin analysis showed that lox2-1 can still accumulate some JA after wounding, which suggests that LOX2 is not the only LOX involved in JA biosynthesis. lox2-1 plants do not show altered susceptibility to the fungus Botrytis cinerea or to a detritivore, however, they are more susceptible to an insect herbivore. The insect and plants are closely co-evolved and a reduced ability to synthesize JA favors the insect. The availability of a lox2-1 mutant with intermediate JA levels will further help understanding why plants produce elevated JA levels.

Mothers' and Fathers' Early Relationship With Their Infant: Similar Yet Temporally Discordant Themes [La relation initiale des parents avec leur nourrisson: Thèmes similaires quoique temporellement discordants]

Objective: To examine first-time mothers' and fathers' themes in their relationship with their infant, how these themes change during the first four months postpartum, and similarities and differences in mothers' and fathers' themes. Participants: Eighteen first-time mother-father couples were separately interviewed at one; six and 16 weeks postpartum. Data Analysis: Audio-taped, transcribed interviews were analysed using a Grounded Theory approach. Results: Our findings reveal a common set of themes for mothers and fathers in relation to the infant : 1: Discovery, 2: Physical Proximity, 3: Emotional Closeness, 4: Initiation of Complementary Interactions and 5: Commitment to Love and Care. However, there was a striking lack of concordance between mothers and fathers for these themes at each point in time. Conclusions: Mothers' and fathers' experience of the early relationship with their infant is unique. Focussing on maternal as well as paternal ways of experiencing the early relationship with their infant sets the way to understanding early developing relationships in the family context.

Paper money and national distress: William Huskisson and the early theories of credit, speculation and crises

This paper examines the explanation of commercial crises offered by William Huskisson in 1810 in the wake of the debate on the Bullion Report. Huskisson argued that the suspension of convertibility made it possible to extend issues of paper currency beyond its proper limits. Such an expansion, being in the interest of all parties concerned, would actually take place and stimulate excessive speculations, which would eventually prove unsustainable and bring generalized ruin and distress. Although some elements of this explanations were not new (having been anticipated by writers sucha as James Currie in 1793, William Roscoe in 1793, William Anderson in 1797 and an anonymous in 1796), Huskisson's explanation is more systematic and better organized, and his emphasis on the endogenous character of the crisis and on the instability of the dynamics of trade and credit makes it an interesting foreshadower of the theories of crises that were advanced half a century later.

CYR61 and alphaVbeta5 integrin cooperate to promote invasion and metastasis of tumors growing in preirradiated stroma.

Radiotherapy is widely used to treat human cancer. Patients locally recurring after radiotherapy, however, have increased risk of metastatic progression and poor prognosis. The clinical management of postradiation recurrences remains an unresolved issue. Tumors growing in preirradiated tissues have an increased fraction of hypoxic cells and are more metastatic, a condition known as tumor bed effect. The transcription factor hypoxia inducible factor (HIF)-1 promotes invasion and metastasis of hypoxic tumors, but its role in the tumor bed effect has not been reported. Here, we show that tumor cells derived from SCCVII and HCT116 tumors growing in a preirradiated bed, or selected in vitro through repeated cycles of severe hypoxia, retain invasive and metastatic capacities when returned to normoxia. HIF activity, although facilitating metastatic spreading of tumors growing in a preirradiated bed, is not essential. Through gene expression profiling and gain- and loss-of-function experiments, we identified the matricellular protein CYR61 and alphaVbeta5 integrin as proteins cooperating to mediate these effects. The anti-alphaV integrin monoclonal antibody 17E6 and the small molecular alphaVbeta3/alphaVbeta5 integrin inhibitor EMD121974 suppressed invasion and metastasis induced by CYR61 and attenuated metastasis of tumors growing within a preirradiated field. These results represent a conceptual advance to the understanding of the tumor bed effect and identify CYR61 and alphaVbeta5 integrin as proteins that cooperate to mediate metastasis. They also identify alphaV integrin inhibition as a potential therapeutic approach for preventing metastasis in patients at risk for postradiation recurrences.

Hyperhidrosis: a new and often early symptom in Fabry disease. International experience and data from the Fabry Outcome Survey.

Hypohidrosis is a classic feature of Fabry disease; in contrast, hyperhidrosis has only been rarely described. The aim of the study is to characterise the baseline descriptive data on hyperhidrosis (frequency, age at onset, sex ratio and outcome with and without enzyme replacement therapy) in hemizygous male and heterozygous female patients with Fabry disease. We describe case histories of five patients with Fabry disease and hyperhidrosis seen at three different centres. We have also analysed a cohort of 21 paediatric patients in the UK and a large European cohort of patients enrolled in the Fabry Outcome Survey (FOS). Five patients (three female, two male) with hyperhidrosis were originally identified, although each had additional symptoms related to Fabry disease. The age at onset of hyperhidrosis was less than 18 years in four cases. In the cohort of 21 paediatric patients (12 female, nine male), one female had hyperhidrosis; the age at onset of this symptom was 11 years. In the FOS cohort, 66 of 714 patients with Fabry disease had hyperhidrosis (44 of 369 females, 11.9%; 22 of 345 males, 6.4%). The female predominance was observed in seven of nine countries from which data were analysed. Hyperhidrosis is an increasingly recognised feature of the Fabry disease phenotype. It is more prevalent in females than in males and often appears in childhood or adolescence. The efficacy of enzyme replacement therapy on this recently recognised symptom should be assessed.

Transforming personal experience and emotions through education to cultural diversity: an interplay between unicity and genericity

When subjects studied at school are close to societal discourses and to the students' social identities, when they have high emotional resonance, is it possible to enable the students to distance themselves from their emotions and personal experience, and to conceptualise them? Examining the relation between emotion and learning through the lens of socio-cultural psychology, the aim of our study was to shed light on "secondarisation" processes, that is, processes that transform personal experience and emotions into conceptualised forms of thinking. We analysed 85 video-recorded lessons in education for cultural diversity involving 12 teachers (of primary and secondary schools). Having identified episodes in which emotions were put into words or personal experience was reported, we analysed the use of pronouns (taken as indicators of secondarisation processes) and found a recurrent pattern: "the unicity-genericity routine". We illustrate the functioning of this routine with various excerpts taken from lessons in education for diversity taught in the classes of two teachers in primary school. The results show that the interplay between unicity and genericity works as a discursive resource for the development of secondarisation processes.

Congenital heart defects in Europe: prevalence and perinatal mortality, 2000 to 2005.

BACKGROUND: This study determines the prevalence of Congenital Heart Defects (CHD), diagnosed prenatally or in infancy, and fetal and perinatal mortality associated with CHD in Europe. METHODS AND RESULTS: Data were extracted from the European Surveillance of Congenital Anomalies central database for 29 population-based congenital anomaly registries in 16 European countries covering 3.3 million births during the period 2000 to 2005. CHD cases (n=26 598) comprised live births, fetal deaths from 20 weeks gestation, and terminations of pregnancy for fetal anomaly (TOPFA). The average total prevalence of CHD was 8.0 per 1000 births, and live birth prevalence was 7.2 per 1000 births, varying between countries. The total prevalence of nonchromosomal CHD was 7.0 per 1000 births, of which 3.6% were perinatal deaths, 20% prenatally diagnosed, and 5.6% TOPFA. Severe nonchromosomal CHD (ie, excluding ventricular septal defects, atrial septal defects, and pulmonary valve stenosis) occurred in 2.0 per 1000 births, of which 8.1% were perinatal deaths, 40% were prenatally diagnosed, and 14% were TOPFA (TOPFA range between countries 0% to 32%). Live-born CHD associated with Down syndrome occurred in 0.5 per 1000 births, with > 4-fold variation between countries. CONCLUSION: Annually in the European Union, we estimate 36 000 children are live born with CHD and 3000 who are diagnosed with CHD die as a TOFPA, late fetal death, or early neonatal death. Investing in primary prevention and pathogenetic research is essential to reduce this burden, as well as continuing to improve cardiac services from in utero to adulthood.

A non-interventional phase IV Belgian survey to assess the antiviral effectiveness of pegylated interferon-alpha-2b and ribavirin treatment according to the stage of liver fibrosis in previously untreated patients with genotype 1/4/5/6 chronic hepatitis C (PRACTICE).

BACKGROUND AND STUDY AIMS: This was an observational, non-interventional, multicenter, phase IV study, in patients with genotype 1/4/5/6 chronic hepatitis C (CHC). The primary objectives were to evaluate SVR in patients with no or minimal fibrosis (METAVIR F0-F1) versus well established fibrosis (F2-F4), and to estimate response on Weeks 12, 24 and 48 on treatment in previously untreated patients with genotypes 1/4/5/6 CHC. PATIENTS AND METHODS: 538 patients treated with pegylated interferon alfa 2b 1.5 mcg/kg in combination with ribavirin 800-1200 mg/day were enrolled in 55 sites in Belgium and Luxembourg, 505 being considered for the analysis. 40% of the patients were female and 60% male, the average age was 47.5 years, 10.5% were 65 or older. RESULTS: SVR was observed in 35% of the patients, EVR in 68%, of which pEVR in 33% and cEVR in 35%. SVR was observed in 43% of the low fibrosis group (F0, F1) and 30% of the high fibrosis group (F2, F3, F4) (p = 0.005). SVR rates were 34% for genotype 1, 37% for genotype 4, and 47% for genotype 5 (NS). Multivariate analysis showed that EVR and baseline METAVIR score are independent prognostic factors for SVR. CONCLUSIONS: This trial confirms that fibrosis stage and early viral response are the most important key-factors to predict sustained response, suggesting that the earlier patients are treated, the better the outcome. Non-invasive techniques enable us to closely monitor progression of fibrosis, allowing a better selection of patients for antiviral treatment in the DAA-era.

A novel and divergent role of granzyme a and B in resistance to helminth infection.

Granzyme (gzm) A and B, proteases of NK cells and T killer cells, mediate cell death, but also cleave extracellular matrices, inactivate intracellular pathogens, and induce cytokines. Moreover, macrophages, Th2 cells, regulatory T cells, mast cells, and B cells can express gzms. We recently reported gzm induction in human filarial infection. In this study, we show that in rodent filarial infection with Litomosoides sigmodontis, worm loads were significantly reduced in gzmA×B and gzmB knockout mice during the whole course of infection, but enhanced only early in gzmA knockout compared with wild-type mice. GzmA/B deficiency was associated with a defense-promoting Th2 cytokine and Ab shift, enhanced early inflammatory gene expression, and a trend of reduced alternatively activated macrophage induction, whereas gzmA deficiency was linked with reduced inflammation and a trend toward increased alternatively activated macrophages. This suggests a novel and divergent role for gzms in helminth infection, with gzmA contributing to resistance and gzmB promoting susceptibility.

Global transcriptional programs in peripheral nerve endoneurium and DRG are resistant to the onset of type 1 diabetic neuropathy in Ins2 mice.

While the morphological and electrophysiological changes underlying diabetic peripheral neuropathy (DPN) are relatively well described, the involved molecular mechanisms remain poorly understood. In this study, we investigated whether phenotypic changes associated with early DPN are correlated with transcriptional alterations in the neuronal (dorsal root ganglia [DRG]) or the glial (endoneurium) compartments of the peripheral nerve. We used Ins2(Akita/+) mice to study transcriptional changes underlying the onset of DPN in type 1 diabetes mellitus (DM). Weight, blood glucose and motor nerve conduction velocity (MNCV) were measured in Ins2(Akita/+) and control mice during the first three months of life in order to determine the onset of DPN. Based on this phenotypic characterization, we performed gene expression profiling using sciatic nerve endoneurium and DRG isolated from pre-symptomatic and early symptomatic Ins2(Akita/+) mice and sex-matched littermate controls. Our phenotypic analysis of Ins2(Akita/+) mice revealed that DPN, as measured by reduced MNCV, is detectable in affected animals already one week after the onset of hyperglycemia. Surprisingly, the onset of DPN was not associated with any major persistent changes in gene expression profiles in either sciatic nerve endoneurium or DRG. Our data thus demonstrated that the transcriptional programs in both endoneurial and neuronal compartments of the peripheral nerve are relatively resistant to the onset of hyperglycemia and hypoinsulinemia suggesting that either minor transcriptional alterations or changes on the proteomic level are responsible for the functional deficits associated with the onset of DPN in type 1 DM.

Phylogenetic relatedness and proboscis length contribute to structuring bumblebee communities in the extremes of abiotic and biotic gradients

Aim To improve our understanding of how biological communities assemble, we investigated changes in bumblebee communities in space along an elevation gradient. We assessed how much deterministic abiotic and biotic factors shape community assembly. We focused on proboscis length (influencing the species' dietary regime) and phylogenetic relatedness to investigate if competition and environmental filtering occur in more and less productive climates, respectively. Location Western Swiss Alps. Methods We recorded bumblebee species in 149 plots along a 1800-m wide elevation gradient. We contrasted two major clades of bumblebees, a short-tongued and a long-tongued clade. We calculated the phylogenetic and proboscis-length diversity of the bumblebee communities and compared these observed data with a random distribution to detect clustering likely to be caused by environmental filtering or overdispersion likely to be caused by competition. We compared the prevalence of clustered and overdispersed communities along the gradients of plant species richness (biotic) and temperature (abiotic). Results Under colder conditions, where plant species richness is lower and floral resources are scarcer, the clade with shorter proboscides prevails over the clade with longer proboscides, and communities are functionally and phylogenetic clustered. Under warmer conditions, we found phylogenetic but not functional overdispersion in communities. Main conclusions We show for the first time a strong correlation between phylogenetic relatedness, proboscis length and species distribution along temperature and plant richness gradients shaping bumblebee communities. The low temperatures and low levels of plant species richness limit the dispersal of the species from the long-tongued clade, which have more specialized diets, into high-elevation areas. Competition under warmer conditions may produce communities composed of less closely related species that share distinct ecological preferences. Our empirical results corroborate theoretical expectation as well as experiments on the prevalence of deterministic processes in the most severe and most productive parts of environmental gradients.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.