Large scale assessment of regulatory evolution and transcriptome complexity in mammals

AbstractIn addition to genetic changes affecting the function of gene products, changes in gene expression have been suggested to underlie many or even most of the phenotypic differences among mammals. However, detailed gene expression comparisons were, until recently, restricted to closely related species, owing to technological limitations. Thus, we took advantage of the latest technologies (RNA-Seq) to generate extensive qualitative and quantitative transcriptome data for a unique collection of somatic and germline tissues from representatives of all major mammalian lineages (placental mammals, marsupials and monotremes) and birds, the evolutionary outgroup.In the first major project of my thesis, we performed global comparative analyses of gene expression levels based on these data. Our analyses provided fundamental insights into the dynamics of transcriptome change during mammalian evolution (e.g., the rate of expression change across species, tissues and chromosomes) and allowed the exploration of the functional relevance and phenotypic implications of transcription changes at a genome-wide scale (e.g., we identified numerous potentially selectively driven expression switches).In a second project of my thesis, which was also based on the unique transcriptome data generated in the context of the first project we focused on the evolution of alternative splicing in mammals. Alternative splicing contributes to transcriptome complexity by generating several transcript isoforms from a single gene, which can, thus, perform various functions. To complete the global comparative analysis of gene expression changes, we explored patterns of alternative splicing evolution. This work uncovered several general and unexpected patterns of alternative splicing evolution (e.g., we found that alternative splicing evolves extremely rapidly) as well as a large number of conserved alternative isoforms that may be crucial for the functioning of mammalian organs.Finally, the third and final project of my PhD consisted in analyzing in detail the unique functional and evolutionary properties of the testis by exploring the extent of its transcriptome complexity. This organ was previously shown to evolve rapidly both at the phenotypic and molecular level, apparently because of the specific pressures that act on this organ and are associated with its reproductive function. Moreover, my analyses of the amniote tissue transcriptome data described above, revealed strikingly widespread transcriptional activity of both functional and nonfunctional genomic elements in the testis compared to the other organs. To elucidate the cellular source and mechanisms underlying this promiscuous transcription in the testis, we generated deep coverage RNA-Seq data for all major testis cell types as well as epigenetic data (DNA and histone methylation) using the mouse as model system. The integration of these complete dataset revealed that meiotic and especially post-meiotic germ cells are the major contributors to the widespread functional and nonfunctional transcriptome complexity of the testis, and that this "promiscuous" spermatogenic transcription is resulting, at least partially, from an overall transcriptionally permissive chromatin state. We hypothesize that this particular open state of the chromatin results from the extensive chromatin remodeling that occurs during spermatogenesis which ultimately leads to the replacement of histones by protamines in the mature spermatozoa. Our results have important functional and evolutionary implications (e.g., regarding new gene birth and testicular gene expression evolution).Generally, these three large-scale projects of my thesis provide complete and massive datasets that constitute valuables resources for further functional and evolutionary analyses of mammalian genomes.

Premature replacement of mu with alpha immunoglobulin chains impairs lymphopoiesis and mucosal homing but promotes plasma cell maturation.

Sequentially along B cell differentiation, the different classes of membrane Ig heavy chains associate with the Ig alpha/Ig beta heterodimer within the B cell receptor (BCR). Whether each Ig class conveys specific signals adapted to the corresponding differentiation stage remains debated. We investigated the impact of the forced expression of an IgA-class receptor throughout murine B cell differentiation by knocking in the human C alpha Ig gene in place of the S mu region. Despite expression of a functional BCR, homozygous mutant mice showed a partial developmental blockade at the pro-B/pre-BI and large pre-BII cell stages, with decreased numbers of small pre-BII cells. Beyond this stage, peripheral B cell compartments of reduced size developed and allowed specific antibody responses, whereas mature cells showed constitutive activation and a strong commitment to plasma cell differentiation. Secreted IgA correctly assembled into polymers, associated with the murine J chain, and was transported into secretions. In heterozygous mutants, cells expressing the IgA allele competed poorly with those expressing IgM from the wild-type allele and were almost undetectable among peripheral B lymphocytes, notably in gut-associated lymphoid tissues. Our data indicate that the IgM BCR is more efficient in driving early B cell education and in mucosal site targeting, whereas the IgA BCR appears particularly suited to promoting activation and differentiation of effector plasma cells.

Transcatheter aortic valve replacement for degenerative bioprosthetic surgical valves: results from the global valve-in-valve registry.

BACKGROUND: Transcatheter aortic valve-in-valve implantation is an emerging therapeutic alternative for patients with a failed surgical bioprosthesis and may obviate the need for reoperation. We evaluated the clinical results of this technique using a large, worldwide registry. METHODS AND RESULTS: The Global Valve-in-Valve Registry included 202 patients with degenerated bioprosthetic valves (aged 77.7±10.4 years; 52.5% men) from 38 cardiac centers. Bioprosthesis mode of failure was stenosis (n=85; 42%), regurgitation (n=68; 34%), or combined stenosis and regurgitation (n=49; 24%). Implanted devices included CoreValve (n=124) and Edwards SAPIEN (n=78). Procedural success was achieved in 93.1% of cases. Adverse procedural outcomes included initial device malposition in 15.3% of cases and ostial coronary obstruction in 3.5%. After the procedure, valve maximum/mean gradients were 28.4±14.1/15.9±8.6 mm Hg, and 95% of patients had ≤+1 degree of aortic regurgitation. At 30-day follow-up, all-cause mortality was 8.4%, and 84.1% of patients were at New York Heart Association functional class I/II. One-year follow-up was obtained in 87 patients, with 85.8% survival of treated patients. CONCLUSIONS: The valve-in-valve procedure is clinically effective in the vast majority of patients with degenerated bioprosthetic valves. Safety and efficacy concerns include device malposition, ostial coronary obstruction, and high gradients after the procedure.

Evaluation of Bridge Replacement Alternatives for the County Bridge System, HR-365 August 1994

Recent reports have indicated that 23.5 percent of the nation's highway bridges are structurally deficient and 17.7 percent are functionally obsolete. A significant number of these bridges are on the Iowa county road system. The objective of the investigation described in this report was to identify, review and evaluate replacement bridges currently being used by various counties in Iowa and surrounding states. Iowa county engineers, county engineers in neighboring states as well as private manufacturers of bridge components, and regional precad prestressed concrete manufacturers were contacted to determine the most common replacement bridge types being used. Depending upon the findings of the review, possible improvements and/or new replacement bridge systems were to be proposed. A questionnaire was developed and sent to county engineers in Iowa and several counties in surrounding states. The results of the questionnaire showed that the most common replacement bridges in Iowa are the continuous concrete slab and prestressed concrete bridges. The primary reason these types are used is because of the availability of standard designs and because of their ease of maintenance. Counties seldom construct these types of bridges using their own labor forces, but instead contract the work. However, county forces are used to construct steel stringer, precast reinforced concrete and timber bridges. In general, 69 percent of the counties indicate an ability and willingness to use their own forces to design and construct relatively short span bridges (i.e., 40 A or less) provided the construction procedures are relatively simple. Several unique replacement bridge types used in Iowa that are constructed by county forces are documented and presented in this report. Sufficient details are provided to allow county engineers to determine if some of these bridges could be used to resolve some of their own replacement bridge problems. Where possible, cost information has also been provided. Each of these bridge types were evaluated for various criteria (e.g., cost effectiveness, conformance to AASI-ITO standards, range of sizes, etc.) by a panel of four Iowa county engineers; a summary of this critique is included. After evaluating the questionnaire responses from the counties and evaluating the various bridge replacement concepts currently in use, one new bridge replacement concept and one modification of a current Iowa county bridge replacement concept were developed. Both of these concepts would utilize county labor forces.

Mechanisms and evolutionary patterns of mammalian and avian dosage compensation.

As a result of sex chromosome differentiation from ancestral autosomes, male mammalian cells only contain one X chromosome. It has long been hypothesized that X-linked gene expression levels have become doubled in males to restore the original transcriptional output, and that the resulting X overexpression in females then drove the evolution of X inactivation (XCI). However, this model has never been directly tested and patterns and mechanisms of dosage compensation across different mammals and birds generally remain little understood. Here we trace the evolution of dosage compensation using extensive transcriptome data from males and females representing all major mammalian lineages and birds. Our analyses suggest that the X has become globally upregulated in marsupials, whereas we do not detect a global upregulation of this chromosome in placental mammals. However, we find that a subset of autosomal genes interacting with X-linked genes have become downregulated in placentals upon the emergence of sex chromosomes. Thus, different driving forces may underlie the evolution of XCI and the highly efficient equilibration of X expression levels between the sexes observed for both of these lineages. In the egg-laying monotremes and birds, which have partially homologous sex chromosome systems, partial upregulation of the X (Z in birds) evolved but is largely restricted to the heterogametic sex, which provides an explanation for the partially sex-biased X (Z) expression and lack of global inactivation mechanisms in these lineages. Our findings suggest that dosage reductions imposed by sex chromosome differentiation events in amniotes were resolved in strikingly different ways.

Efficacy of enzyme replacement therapy in Fabry disease.

Enzyme replacement therapy has recently been introduced to treat Fabry disease, a rare X-linked lysosomal storage disorder. The disease occurs due to deficient activity of alpha-galactosidase A, leading to progressive accumulation of globotriaosylceramide in multiple organs and tissues. Renal, cardiac and cerebrovascular manifestations of the disease result in premature death in both hemizygous males and heterozygous females. This paper outlines the clinical signs, symptoms and diagnosis of Fabry disease, and the development of the two available enzyme replacement therapies -- agalsidase alfa and agalsidase beta. Agalsidase alfa and agalsidase beta are produced in a human cell line and in Chinese hamster ovary cells, respectively, resulting in products with the same amino acid sequence as the native human enzyme, but with different patterns of glycosylation. Correct post-translational glycosylation is important in terms of the pharmacokinetics, biodistribution, clinical efficacy and tolerability of genetically engineered protein therapeutics. Differences in glycosylation, which may affect immunogenicity and mannose-6-phosphate receptor-mediated cellular internalisation of administered enzyme, possibly account for the differences in dosing, clinical effects and safety profiles reported for agalsidase alfa and agalsidase beta.

Monophyly of beta-tubulin and H+-ATPase gene variants in Glomus intraradices: consequences for molecular evolutionary studies of AM fungal genes.

Arbuscular mycorrhizal fungi (AMF) are an ecologically important group of fungi. Previous studies showed the presence of divergent copies of beta-tubulin and V-type vacuolar H+-ATPase genes in AMF genomes and suggested horizontal gene transfer from host plants or mycoparasites to AMF. We sequenced these genes from DNA isolated from an in vitro cultured isolate of Glomus intraradices that was free of any obvious contaminants. We found two highly variable beta-tubulin sequences and variable H+-ATPase sequences. Despite this high variation, comparison of the sequences with those in gene banks supported a glomeromycotan origin of G. intraradices beta-tubulin and H+-ATPase sequences. Thus, our results are in sharp contrast with the previously reported polyphyletic origin of those genes. We present evidence that some highly divergent sequences of beta-tubulin and H+-ATPase deposited in the databases are likely to be contaminants. We therefore reject the prediction of horizontal transfer to AMF genomes. High differences in GC content between glomeromycotan sequences and sequences grouping in other lineages are shown and we suggest they can be used as an indicator to detect such contaminants. H+-ATPase phylogeny gave unexpected results and failed to resolve fungi as a natural group. beta-Tubulin phylogeny supported Glomeromeromycota as sister group of the Chytridiomycota. Contrasts between our results and trees previously generated using rDNA sequences are discussed.

Evolutionary dynamics of molecular markers during local adaptation: a case study in Drosophila subobscura

Background: Natural selection and genetic drift are major forces responsible for temporal genetic changes in populations. Furthermore, these evolutionary forces may interact with each other. Here we study the impact of an ongoing adaptive process at the molecular genetic level by analyzing the temporal genetic changes throughout 40 generations of adaptation to a common laboratory environment. Specifically, genetic variability, population differentiation and demographic structure were compared in two replicated groups of Drosophila subobscura populations recently sampled from different wild sources. Results: We found evidence for a decline in genetic variability through time, along with an increase in genetic differentiation between all populations studied. The observed decline in genetic variability was higher during the first 14 generations of laboratory adaptation. The two groups of replicated populations showed overall similarity in variability patterns. Our results also revealed changing demographic structure of the populations during laboratory evolution, with lower effective population sizes in the early phase of the adaptive process. One of the ten microsatellites analyzed showed a clearly distinct temporal pattern of allele frequency change, suggesting the occurrence of positive selection affecting the region around that particular locus. Conclusion: Genetic drift was responsible for most of the divergence and loss of variability between and within replicates, with most changes occurring during the first generations of laboratory adaptation. We also found evidence suggesting a selective sweep, despite the low number of molecular markers analyzed. Overall, there was a similarity of evolutionary dynamics at the molecular level in our laboratory populations, despite distinct genetic backgrounds and some differences in phenotypic evolution.

Shaping mechanisms of metal specificity in a family of metazoan Metallothioneins: evolutionary differentiation of Mollusc MTs.

Background: The degree of metal binding specificity in metalloproteins such as metallothioneins (MTs) can be crucial for their functional accuracy. Unlike most other animal species, pulmonate molluscs possess homometallic MT isoforms loaded with Cu+ or Cd2+. They have, so far, been obtained as native metal-MT complexes from snail tissues, where they are involved in the metabolism of the metal ion species bound to the respective isoform. However, it has not as yet been discerned if their specific metal occupation is the result of a rigid control of metal availability, or isoform expression programming in the hosting tissues or of structural differences of the respective peptides determining the coordinative options for the different metal ions. In this study, the Roman snail (Helix pomatia) Cu-loaded and Cd-loaded isoforms (HpCuMT and HpCdMT) were used as model molecules in order t o elucidate the biochemical and evolutionary mechanisms permitting pulmonate MTs to achieve specificity for their cognate metal ion. Results: HpCuMT and HpCdMT were recombinantly synthesized in the presence of Cd2+, Zn2+ or Cu2+ and corresponding metal complexes analysed by electrospray mass spectrometry and circular dichroism (CD) and ultra violet-visible (UV-Vis) spectrophotometry. Both MT isoforms were only able to form unique, homometallic and stable complexes (Cd6-HpCdMT and Cu12-HpCuMT) with their cognate metal ions. Yeast complementation assays demonstrated that the two isoforms assumed metal-specific functions, in agreement with their binding preferences, in heterologous eukaryotic environments. In the snail organism, the functional metal specificity of HpCdMT and HpCuMT was contributed by metal-specific transcription programming and cell-specific expression. Sequence elucidation and phylogenetic analysis of MT isoforms from a number of snail species revealed that they possess an unspecific and two metal-specific MT isoforms, whose metal specificity was achieved exclusively by evolutionary modulation of non-cysteine amino acid positions. Conclusion: The Roman snail HpCdMT and HpCuMT isoforms can thus be regarded as prototypes of isoform families that evolved genuine metal-specificity within pulmonate molluscs. Diversification into these isoforms may have been initiated by gene duplication, followed by speciation and selection towards opposite needs for protecting copper-dominated metabolic pathways from nonessential cadmium. The mechanisms enabling these proteins to be metal-specific could also be relevant for other metalloproteins.

Comparative genomics of the odorant-binding and chemosensory protein gene families across the arthropoda: origin and evolutionary history of the chemosensory system

Chemoreception is a biological process essential for the survival of animals, as it allows the recognition of important volatile cues for the detection of food, egg-laying substrates, mates or predators, among other purposes. Furthermore, its role in pheromone detection may contribute to evolutionary processes such as reproductive isolation and speciation. This key role in several vital biological processes makes chemoreception a particularly interesting system for studying the role of natural selection in molecular adaptation. Two major gene families are involved in the perireceptor events of the chemosensory system: the odorant-binding protein (OBP) and chemosensory protein (CSP) families. Here, we have conducted an exhaustive comparative genomic analysis of these gene families in twenty Arthropoda species. We show that the evolution of the OBP and CSP gene families is highly dynamic, with a high number of gains and losses of genes, pseudogenes and independent origins of subfamilies. Taken together, our data clearly support the birth-and-death model for the evolution of these gene families with an overall high gene-turnover rate. Moreover, we show that the genome organization of the two families is significantly more clustered than expected by chance and, more important, that this pattern appears to be actively maintained across the Drosophila phylogeny. Finally, we suggest the homologous nature of the OBP and CSP gene families, dating back their MRCA (most recent common ancestor) to 380¿420 Mya, and we propose a scenario for the origin and diversification of these families.

Isolation of a novel monkey adenovirus reveals a new phylogenetic clade in the evolutionary history of simian adenoviruses

Adenoviruses of primates include human (HAdV) and simian (SAdV) isolates classified into 8 species (Human Adenovirus A to G, and Simian Adenovirus A). In this study, a novel adenovirus was isolated from a colony of cynomolgus macaques (Macaca fascicularis) and subcultured in VERO cells. Its complete genome was purified and a region encompassing the hexon gene, the protease gene, the DNA binding protein (DBP) and the 100 kDa protein was amplified by PCR and sequenced by primer walking. Sequence analysis of these four genes showed that the new isolate had 80% identity to other primate adenoviruses and lacked recombination events. The study of the evolutionary relationships of this new monkey AdV based on the combined sequences of the four genes supported a close relationship to SAdV-3 and SAdV-6, lineages isolated from Rhesus monkeys. The clade formed by these three types is separated from the remaining clades and establishes a novel branch that is related to species HAdV-A, F and G. However, the genetic distance corresponding to the newly isolated monkey AdV considerably differs from these as to belong to a new, not yet established species. Results presented here widen our knowledge on SAdV and represents an important contribution to the understanding of the evolutionary history of primate adenoviruses.

Economics of dialysis dependence following renal replacement therapy for critically ill acute kidney injury patients.

BACKGROUND: The obective of this study was to perform a cost-effectiveness analysis comparing intermittent with continuous renal replacement therapy (IRRT versus CRRT) as initial therapy for acute kidney injury (AKI) in the intensive care unit (ICU). METHODS: Assuming some patients would potentially be eligible for either modality, we modeled life year gained, the quality-adjusted life years (QALYs) and healthcare costs for a cohort of 1000 IRRT patients and a cohort of 1000 CRRT patients. We used a 1-year, 5-year and a lifetime horizon. A Markov model with two health states for AKI survivors was designed: dialysis dependence and dialysis independence. We applied Weibull regression from published estimates to fit survival curves for CRRT and IRRT patients and to fit the proportion of dialysis dependence among CRRT and IRRT survivors. We then applied a risk ratio reported in a large retrospective cohort study to the fitted CRRT estimates in order to determine the proportion of dialysis dependence for IRRT survivors. We conducted sensitivity analyses based on a range of differences for daily implementation cost between CRRT and IRRT (base case: CRRT day $632 more expensive than IRRT day; range from $200 to $1000) and a range of risk ratios for dialysis dependence for CRRT as compared with IRRT (from 0.65 to 0.95; base case: 0.80). RESULTS: Continuous renal replacement therapy was associated with a marginally greater gain in QALY as compared with IRRT (1.093 versus 1.078). Despite higher upfront costs for CRRT in the ICU ($4046 for CRRT versus $1423 for IRRT in average), the 5-year total cost including the cost of dialysis dependence was lower for CRRT ($37 780 for CRRT versus $39 448 for IRRT on average). The base case incremental cost-effectiveness analysis showed that CRRT dominated IRRT. This dominance was confirmed by extensive sensitivity analysis. CONCLUSIONS: Initial CRRT is cost-effective compared with initial IRRT by reducing the rate of long-term dialysis dependence among critically ill AKI survivors.

Decision Support Model for Assessing Archaeological Survey Needs for Bridge Replacement Projects in Iowa, 2006

The Bridges Decision Support Model is a geographic information system (GIS) that assembles existing data on archaeological sites, surveys, and their geologic contexts to assess the risk of bridge replacement projects encountering 13,000- to 150-year-old Native American sites. This project identifies critical variables for assessing prehistoric sites potential, examines the quality of available data about the variables, and applies the data to creating a decision support framework for use by the Iowa Department of Transportation (Iowa DOT) and others. An analysis of previous archaeological surveys indicates that subsurface testing to discover buried sites became increasingly common after 1980, but did not become routine until after the adoption of guidelines recommending such testing, in 1993. Even then, the average depth of testing has been relatively shallow. Alluvial deposits of sufficient age, deposited in depositional environments conducive to human habitation, are considerably thicker than archaeologists have routinely tested.

Isolation of a novel monkey adenovirus reveals a new phylogenetic clade in the evolutionary history of simian adenoviruses

Adenoviruses of primates include human (HAdV) and simian (SAdV) isolates classified into 8 species (Human Adenovirus A to G, and Simian Adenovirus A). In this study, a novel adenovirus was isolated from a colony of cynomolgus macaques (Macaca fascicularis) and subcultured in VERO cells. Its complete genome was purified and a region encompassing the hexon gene, the protease gene, the DNA binding protein (DBP) and the 100 kDa protein was amplified by PCR and sequenced by primer walking. Sequence analysis of these four genes showed that the new isolate had 80% identity to other primate adenoviruses and lacked recombination events. The study of the evolutionary relationships of this new monkey AdV based on the combined sequences of the four genes supported a close relationship to SAdV-3 and SAdV-6, lineages isolated from Rhesus monkeys. The clade formed by these three types is separated from the remaining clades and establishes a novel branch that is related to species HAdV-A, F and G. However, the genetic distance corresponding to the newly isolated monkey AdV considerably differs from these as to belong to a new, not yet established species. Results presented here widen our knowledge on SAdV and represents an important contribution to the understanding of the evolutionary history of primate adenoviruses.