The evolution of gene expression levels in mammalian organs.

Changes in gene expression are thought to underlie many of the phenotypic differences between species. However, large-scale analyses of gene expression evolution were until recently prevented by technological limitations. Here we report the sequencing of polyadenylated RNA from six organs across ten species that represent all major mammalian lineages (placentals, marsupials and monotremes) and birds (the evolutionary outgroup), with the goal of understanding the dynamics of mammalian transcriptome evolution. We show that the rate of gene expression evolution varies among organs, lineages and chromosomes, owing to differences in selective pressures: transcriptome change was slow in nervous tissues and rapid in testes, slower in rodents than in apes and monotremes, and rapid for the X chromosome right after its formation. Although gene expression evolution in mammals was strongly shaped by purifying selection, we identify numerous potentially selectively driven expression switches, which occurred at different rates across lineages and tissues and which probably contributed to the specific organ biology of various mammals.

Solving ambiguities in contig assembly of Idiomarina loihiensis L2TR chromosome by in silico analyses.

Nucleotide composition analyses of bacterial genomes such as cumulative GC skew highlight the atypical, strongly asymmetric architecture of the recently published chromosome of Idiomarina loihiensis L2TR, suggesting that an inversion of a 600-kb chromosomal segment occurred. The presence of 3.4-kb inverted repeated sequences at the borders of the putative rearrangement supports this hypothesis. Reverting in silico this segment restores (1) a symmetric chromosome architecture; (2) the co-orientation of transcription of all rRNA operons with DNA replication; and (3) a better conservation of gene order between this chromosome and other gamma-proteobacterial ones. Finally, long-range PCRs encompassing the ends of the 600-kb segment reveal the existence of the reverted configuration but not of the published one. This demonstrates how cumulative nucleotide-skew analyses can validate genome assemblies.

Population consequences of environmental sex reversal.

When sex determination in a species is predominantly genetic but environmentally reversible, exposure to (anthropogenic) changes in the environment can lead to shifts in a population's sex ratio. Such scenarios may be common in many fishes and amphibians, yet their ramifications remain largely unexplored. We used a simple model to study the (short-term) population consequences of environmental sex reversal (ESR). We examined the effects on sex ratios, sex chromosome frequencies, and population growth and persistence after exposure to environmental forces with feminizing or masculinizing tendencies. When environmental feminization was strong, X chromosomes were driven to extinction. Analogously, extinction of normally male-linked genetic factors (e.g., Y chromosomes) was caused by continuous environmental masculinization. Although moderate feminization was beneficial for population growth in the absence of large viability effects, our results suggest that the consequences of ESR are generally negative in terms of population size and the persistence of sex chromosomes. Extreme sex ratios resulting from high rates of ESR also reduced effective population sizes considerably. This may limit any evolutionary response to the deleterious effects of ESR. Our findings suggest that ESR changes population growth and sex ratios in some counter-intuitive ways and can change the predominant factor in sex determination from genetic to fully environmental, often within only a few tens of generations. Populations that lose genetic sex determination may quickly go extinct if the environmental forces that cause sex reversal cease.

A sex-specific marker reveals male heterogamety in European tree frogs.

Most amphibians examined so far show undifferentiated sex chromosomes. The heterogametic sex's identity, usually revealed through indirect means, often varies among closely related species or even populations (as do sex-linkage groups), suggesting great evolutionary instability of the sex-determining genes. Here we take advantage of a sex-specific marker that amplifies in several related species of European tree frogs (Hyla arborea group) to disclose a homogeneous pattern of male heterogamety. Besides relevance for evolutionary studies of sex determination in amphibians, our results have potential for addressing practical issues in conservation biology because sex reversal by anthropogenic endocrine disruptors is considered one possible cause of amphibian decline.

Evolutionary dynamics of coding and non-coding transcriptomes.

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

Paper 4: EUROCAT statistical monitoring: identification and investigation of ten year trends of congenital anomalies in Europe.

BACKGROUND: As part of EUROCAT's surveillance of congenital anomalies in Europe, a statistical monitoring system has been developed to detect recent clusters or long-term (10 year) time trends. The purpose of this article is to describe the system for the identification and investigation of 10-year time trends, conceived as a "screening" tool ultimately leading to the identification of trends which may be due to changing teratogenic factors.METHODS: The EUROCAT database consists of all cases of congenital anomalies including livebirths, fetal deaths from 20 weeks gestational age, and terminations of pregnancy for fetal anomaly. Monitoring of 10-year trends is performed for each registry for each of 96 non-independent EUROCAT congenital anomaly subgroups, while Pan-Europe analysis combines data from all registries. The monitoring results are reviewed, prioritized according to a prioritization strategy, and communicated to registries for investigation. Twenty-one registries covering over 4 million births, from 1999 to 2008, were included in monitoring in 2010.CONCLUSIONS: Significant increasing trends were detected for abdominal wall anomalies, gastroschisis, hypospadias, Trisomy 18 and renal dysplasia in the Pan-Europe analysis while 68 increasing trends were identified in individual registries. A decreasing trend was detected in over one-third of anomaly subgroups in the Pan-Europe analysis, and 16.9% of individual registry tests. Registry preliminary investigations indicated that many trends are due to changes in data quality, ascertainment, screening, or diagnostic methods. Some trends are inevitably chance phenomena related to multiple testing, while others seem to represent real and continuing change needing further investigation and response by regional/national public health authorities.

Prognostic factors in early-stage leiomyosarcoma of the uterus.

Uterine leiomyosarcomas (LMSs) are rare cancers representing less than 1% of all uterine malignancies. Clinical International Federation of Gynecology and Obstetrics (FIGO) stage is the most important prognostic factor. Other significant prognostic factors, especially for early stages, are difficult to establish because most of the published studies have included localized and extra-pelvian sarcomas. The aim of our study was to search for significant prognostic factors in clinical stage I and II uterine LMS. The pathologic features of 108 uterine LMS including 72 stage I and II lesions were reviewed using standardized criteria. The prognostic significance of different pathologic features was assessed. The median follow-up in the whole group was 64 months (range, 6-223 months). The 5-year overall survival (OS) and metastasis-free interval and local relapse-free interval rates in the whole group and early-stage group (FIGO stages I and II) were 40% and 57%, 42% and 50%, 56% and 62%, respectively. Clinical FIGO stage was the most important prognostic factor for OS in the whole group (P = 4 x 10). In the stage I and II group, macroscopic circumscription was the most significant factor predicting OS (P = 0.001). In the same group, mitotic score and vascular invasion were associated with metastasis-free interval (P = 0.03 and P = 0.04, respectively). Uterine LMSs diagnosed using standardized criteria have a poor prognosis, and clinical FIGO stage is an ominous prognostic factor. In early-stage LMS, pathologic features such as mitotic score, vascular invasion, and tumor circumscription significantly impact patient outcome.

Prediction of 18-month survival in patients with primary myelodysplastic syndrome. A regression model and scoring system based on the combination of chromosome findings and the Bournemouth score.

The predictive potential of six selected factors was assessed in 72 patients with primary myelodysplastic syndrome using univariate and multivariate logistic regression analysis of survival at 18 months. Factors were age (above median of 69 years), dysplastic features in the three myeloid bone marrow cell lineages, presence of chromosome defects, all metaphases abnormal, double or complex chromosome defects (C23), and a Bournemouth score of 2, 3, or 4 (B234). In the multivariate approach, B234 and C23 proved to be significantly associated with a reduction in the survival probability. The similarity of the regression coefficients associated with these two factors means that they have about the same weight. Consequently, the model was simplified by counting the number of factors (0, 1, or 2) present in each patient, thus generating a scoring system called the Lausanne-Bournemouth score (LB score). The LB score combines the well-recognized and easy-to-use Bournemouth score (B score) with the chromosome defect complexity, C23 constituting an additional indicator of patient outcome. The predicted risk of death within 18 months calculated from the model is as follows: 7.1% (confidence interval: 1.7-24.8) for patients with an LB score of 0, 60.1% (44.7-73.8) for an LB score of 1, and 96.8% (84.5-99.4) for an LB score of 2. The scoring system presented here has several interesting features. The LB score may improve the predictive value of the B score, as it is able to recognize two prognostic groups in the intermediate risk category of patients with B scores of 2 or 3. It has also the ability to identify two distinct prognostic subclasses among RAEB and possibly CMML patients. In addition to its above-described usefulness in the prognostic evaluation, the LB score may bring new insights into the understanding of evolution patterns in MDS. We used the combination of the B score and chromosome complexity to define four classes which may be considered four possible states of myelodysplasia and which describe two distinct evolutional pathways.

Regulatory RNA as mediator in GacA/RsmA-dependent global control of exoproduct formation in Pseudomonas fluorescens CHA0.

In Pseudomonas fluorescens CHA0, an antagonist of root-pathogenic fungi, the GacS/GacA two-component system tightly controls the expression of antifungal secondary metabolites and exoenzymes at a posttranscriptional level, involving the RNA-binding protein and global regulator of secondary metabolism RsmA. This protein was purified from P. fluorescens, and RNA bound to it was converted to cDNA, which served as a probe to isolate the corresponding chromosomal locus, rsmZ. This gene encoded a regulatory RNA of 127 nucleotides and a truncated form lacking 35 nucleotides at the 3' end. Expression of rsmZ depended on GacA, increased with increasing population density, and was stimulated by the addition of a solvent-extractable extracellular signal produced by strain CHA0 at the end of exponential growth. This signal appeared to be unrelated to N-acyl-homoserine lactones. A conserved upstream element in the rsmZ promoter, but not the stress sigma factor RpoS, was involved in rsmZ expression. Overexpression of rsmZ effectively suppressed the negative effect of gacS and gacA mutations on target genes, i.e., hcnA (for hydrogen cyanide synthase) and aprA (for the major exoprotease). Mutational inactivation of rsmZ resulted in reduced expression of these target genes in the presence of added signal. Overexpression of rsmA had a similar, albeit stronger negative effect. These results support a model in which GacA upregulates the expression of regulatory RNAs, such as RsmZ of strain CHA0, in response to a bacterial signal. By a titration effect, RsmZ may then alleviate the repressing activity of RsmA on the expression of target mRNAs.

Crocidura cossyrensis Contoli, 1989 (Mammalia, Soridiae): karyotype, biochemical genetics and hybridization experiments

Crocidura cossyrensis Contoli, 1989 (Mammalia, Soricidae): karyotype, biochemical genetics and hybridization experiments. - The shrew Crocidura cossyrensis Contoli, 1989 from Pantelleria (I), a Mediterranean island 100 km south of Sicily and 70 km west from Tunisia, was investigated in order to understand its origin and its relationship with C. russula from Tunisia, Morocco and Switzerland. With the exception of a single heterozygote centric fusion, C. cossyrensis had a karyotype identical with that of C russula from Tunisia (2N = 42, NF = 70 to 72), but it differed from C russula from Morocco and Switzerland (2N = 42, NF = 60). The former have 5-6 pairs of chromosomes with small arms that are acrocentric in the latter. Genetic comparisons with allozyme data revealed small genetic distance (0.04) between C cossyrensis and C russula from Tunisia. In contrast, this eastern clade (Tunisia and Pantelleria) is separated from the western clade (Switzerland and Morocco) by a genetic distance of 0.14. A hybridization experiment between shrews from Pantelleria and Switzerland lead rapidly to an F1 generation. From 12 F1 hybrids that were backcrossed, females reproduced normally, but none of the males did so. Concluding from the results, C. cossyrensis from Pantelleria and C. russula cf. agilis from Tunisia belong to the same taxon that may have reached the differentiation of a biological species within the C. russula group. More geographic samples are needed to determine the definitive taxonomic positions of these shrews.

Chromosomal rearrangement and genetic structure in the shrews of the Sorex araneus group

ABSTRACT The role of chromosomal rearrangements in the speciation process is much debated and many theoretical models have been developed. The shrews of the Sorex araneus group offer extraordinary opportunities to study the relationship between chromosomal variation and speciation. Indeed, this group of morphologically very similar species received a great deal of attention due to its karyotypic variability, which is mainly attributed to Robertsonian fusions. To explore the impact of karyotypic changes on genetic differentiation, we first studied the relationship between genetic and karyotypic structure among Alpine species and among chromosome races of the S. araneus group using Bayesian admixture analyses. The results of these analyses confirmed the taxonomic status of the studied species even though introgression can still be detected between species. Moreover, the strong spatial sub-structure highlighted the role of historical factors (e.g. geographical isolation) on genetic structure. Next, we studied gene flow at the chromosome level to address the question of the impact of chromosomal rearrangements on genetic differentiation. We used flow sorted chromosomes from three different karyotypic taxa of the S. araneus group to map microsatellite markers at the chromosóme arm level. We have been able to map 24 markers and to show that the karyotypic organisation of these taxa is well conserved, which suggests that these markers can be used for further inter-taxa studies. A general prediction of chromosomal speciation models is that genetic differentiation between two taxa should be larger across rearranged chromosomes than across chromosomes common to both taxa. We combined two approaches using mapped microsatellites to test this prediction. First, we studied the genetic differentiation among five shrew taxa placed at different evolutionary levels (i.e. within and among species). In this large scale study, we detected an overall significant difference in genetic structure between rearranged vs. common chromosomes. Moreover, this effect varied among pairwise comparisons, which allowed us to differentiate the role of the karyotypic complexity of hybrids and of the evolutionary divergence between taxa. Secondly, we compared the levels of gene flow measured across common vs. rearranged chromosomes in two karyotypically different hybrid zones (strong vs. low complexity of hybrids), which show similar levels of genetic structure. We detected a significantly stronger genetic structure across rearranged chromosomes in the hybrid zone showing the highest level of hybrid complexity. The large variance observed among loci suggested that other factors, such as the position of markers within the chromosome, also certainly affects genetic structure. In conclusion, our results strongly support the role of chromosomal rearrangements in the reproductive barrier and suggest their importance in speciation process of the S. araneus group. RESUME Le rôle des réarrangements chromosomiques dans les processus de spéciation est fortement débattu et de nombreux modèles théoriques ont été développés sur le sujet. Les musaraignes du groupe Sorex araneus présentent de nombreuses opportunités pour étudier les relations entre les variations chromosomiques et la spéciation. En effet, ce groupe d'espèces morphologiquement très proches a attiré l'attention des chercheurs en raison de sa variabilité caryotypique principalement attribuée à des fusions Robertsoniennes. Pour explorer l'impact des changements caryotypiques sur la différenciation génétique, nous avons tout d'abord étudié les relations entre la structure génétique et caryotypique de races chromosomiques et d'espèces alpine du groupe S. araneus en utilisant des analyses Bayesiennes d' « admixture ». Les résultats de ces analyses ont confirmé le statut taxonomique des espèces étudiées bien que nous ayons détecté de l'introgression entre espèces. L'observation d'une sous structure spatiale relativement forte souligne l'importance des facteurs historiques (telle que l'isolation géographique) sur la structure génétique de ce groupe. Ensuite, nous avons étudié le flux de gène au niveau des chromosomes pour aborder de manière directe la question de l'impact des réarrangements chromosomiques sur la différenciation génétique. En conséquence, nous avons utilisé des tris de chromosomes de trois taxons du groupe S. araneus pour localiser des marqueurs microsatellites au niveau du bras chromosomique. Au cours de cette étude, nous avons pu localiser 24 marqueurs et montrer une forte conservation dans l'organisation du caryotype de ces taxa. Ce résultat suggère que leur utilisation est appropriée pour des études entre taxa. Une prédiction générale à tous les modèles de spéciation chromosomique correspond à la plus grande différenciation génétique des chromosomes réarrangés que des chromosomes communs. Nous avons combiné deux approches utilisant des microsatellites localisés au niveau du bras chromosomique pour tester cette prédiction. Premièrement, nous avons étudié la différenciation génétique entre cinq taxa du groupe S. araneus se trouvant à des niveaux évolutifs différents (i.e. à l'intérieur et entre espèce). Au cours de cette étude, nous avons détecté une différenciation globale significativement plus élevée sur les chromosomes réarrangés. Cet effet varie entre les comparaisons, ce qui nous a permis de souligner le rôle de la complexité caryotypique des hybrides et du niveau de divergence évolutive entre taxa. Deuxièmement, nous avons comparé le flux de gènes des chromosomes communs et réarrangés dans deux zones d'hybridation caryotypiquement différentes (forte vs. Faible complexité des hybrides) mais présentant un niveau de différenciation génétique similaire. Ceci nous a permis de détecter une structure génétique significativement plus élevée sur les chromosomes réarrangés au centre de la zone d'hybridation présentant la plus grande complexité caryotypic. La forte variance observée entre loci souligne en outre le fait que d'autres facteurs, tel que la position du marqueur sur le chromosome, affectent probablement aussi la structure génétique mesurée. En conclusion, nos résultats supportent fortement le rôle des réarrangements chromosomiques dans la barrière reproductive entre espèces ainsi que leur importance dans les processus de spéciation des musaraignes du groupe S. araneus.

Tree of Sex: A database of sexual systems

The vast majority of eukaryotic organisms reproduce sexually, yet the nature of the sexual system and the mechanism of sex determination often vary remarkably, even among closely related species. Some species of animals and plants change sex across their lifespan, some contain hermaphrodites as well as males and females, some determine sex with highly differentiated chromosomes, while others determine sex according to their environment. Testing evolutionary hypotheses regarding the causes and consequences of this diversity requires interspecific data placed in a phylogenetic context. Such comparative studies have been hampered by the lack of accessible data listing sexual systems and sex determination mechanisms across the eukaryotic tree of life. Here, we describe a database developed to facilitate access to sexual system and sex chromosome information, with data on sexual systems from 11,038 plant, 705 fish, 173 amphibian, 593 non-avian reptilian, 195 avian, 479 mammalian, and 11,556 invertebrate species.

RNA-based gene duplication: mechanistic and evolutionary insights.

Gene copies that stem from the mRNAs of parental source genes have long been viewed as evolutionary dead-ends with little biological relevance. Here we review a range of recent studies that have unveiled a significant number of functional retroposed gene copies in both mammalian and some non-mammalian genomes. These studies have not only revealed previously unknown mechanisms for the emergence of new genes and their functions but have also provided fascinating general insights into molecular and evolutionary processes that have shaped genomes. For example, analyses of chromosomal gene movement patterns via RNA-based gene duplication have shed fresh light on the evolutionary origin and biology of our sex chromosomes.

Mutation hot spots in 5q31-linked corneal dystrophies.

Mutations in the BIGH3 gene on chromosome 5q31 cause four distinct autosomal dominant diseases of the human cornea: granular (Groenouw type I), Reis-Bücklers, lattice type I, and Avellino corneal dystrophies. All four diseases are characterized by both progressive accumulation of corneal deposits and eventual loss of vision. We have identified a specific recurrent missense mutation for each type of dystrophy, in 10 independently ascertained families. Genotype analysis with microsatellite markers surrounding the BIGH3 locus was performed in these 10 families and in 5 families reported previously. The affected haplotype could be determined in 10 of the 15 families and was different in each family. These data indicate that R555W, R124C, and R124H mutations occurred independently in several ethnic groups and that these mutations do not reflect a putative founder effect. Furthermore, this study confirms the specific importance of the R124 and R555 amino acids in the pathogenesis of autosomal dominant corneal dystrophies linked to 5q.