A 20-nucleotide element in the intestinal fatty acid binding protein gene modulates its cell lineage-specific, differentiation-dependent, and cephalocaudal patterns of expression in transgenic mice.

A sequence of epithelial cell proliferation, allocation to four principal lineages, migration-associated differentiation, and cell loss occurs along the crypt-villus axis of the mouse intestine. The sequence is completed in a few days and is recapitulated throughout the life-span of the animal. We have used an intestine-specific fatty acid binding protein gene, Fabpi, as a model for studying regulation of gene expression in this unique developmental system. Promoter mapping studies in transgenic mice identified a 20-bp cis-acting element (5'-AGGTGGAAGCCATCACACTT-3') that binds small intestinal nuclear proteins and participates in the control of Fabpi's cephalocaudal, differentiation-dependent, and cell lineage-specific patterns of expression. Immunocytochemical studies using confocal and electron microscopy indicate that it does so by acting as a suppressor of gene expression in the distal small intestine/colon, as a suppressor of gene activation in proliferating and nonproliferating cells located in the crypts of Lieberkühn, and as a suppressor of expression in the growth factor and defensin-producing Paneth cell lineage. The 20-bp domain has no obvious sequence similarities to known transcription factor binding sites. The three functions modulated by this compact element represent the types of functions required to establish and maintain the intestine's remarkably complex spatial patterns of gene expression. The transgenes described in this report also appear to be useful in characterizing the crypt's stem cell hierarchy.

Phylogenetic plant community structure along elevation is lineage specific

The trend of closely related taxa to retain similar environmental preferences mediated by inherited traits suggests that several patterns observed at the community scale originate from longer evolutionary processes. While the effects of phylogenetic relatedness have been previously studied within a single genus or family, lineage-specific effects on the ecological processes governing community assembly have rarely been studied for entire communities or flora. Here, we measured how community phylogenetic structure varies across a wide elevation gradient for plant lineages represented by thirty-five families, using a co-occurrence index and net relatedness index (NRI). We propose a framework that analyses each lineage separately and reveals the trend of ecological assembly at tree nodes. We found prevailing phylogenetic clustering for more ancient nodes and overdispersion in more recent tree nodes. Closely related species may thus rapidly evolve new environmental tolerances to radiate into distinct communities, while older lineages likely retain inherent environmental tolerances to occupy communities in similar environments, either through efficient dispersal mechanisms or the exclusion of older lineages with more divergent environmental tolerances. Our study illustrates the importance of disentangling the patterns of community assembly among lineages to better interpret the ecological role of traits. It also sheds light on studies reporting absence of phylogenetic signal, and opens new perspectives on the analysis of niche and trait conservatism across lineages.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

K-r/K-c but not d(N)/d(S) correlates positively with body mass in birds, raising implications for inferring lineage-specific selection

Background: The ratio of the rates of non-synonymous and synonymous substitution (d(N)/d(S)) is commonly used to estimate selection in coding sequences. It is often suggested that, all else being equal, d(N)/d(S) should be lower in populations with large effective size (Ne) due to increased efficacy of purifying selection. As N-e is difficult to measure directly, life history traits such as body mass, which is typically negatively associated with population size, have commonly been used as proxies in empirical tests of this hypothesis. However, evidence of whether the expected positive correlation between body mass and d(N)/d(S) is consistently observed is conflicting. Results: Employing whole genome sequence data from 48 avian species, we assess the relationship between rates of molecular evolution and life history in birds. We find a negative correlation between dN/dS and body mass, contrary to nearly neutral expectation. This raises the question whether the correlation might be a method artefact. We therefore in turn consider non-stationary base composition, divergence time and saturation as possible explanations, but find no clear patterns. However, in striking contrast to d(N)/d(S), the ratio of radical to conservative amino acid substitutions (K-r/K-c) correlates positively with body mass. Conclusions: Our results in principle accord with the notion that non-synonymous substitutions causing radical amino acid changes are more efficiently removed by selection in large populations, consistent with nearly neutral theory. These findings have implications for the use of d(N)/d(S) and suggest that caution is warranted when drawing conclusions about lineage-specific modes of protein evolution using this metric.

Consequences of lineage-specific gene loss on functional evolution of surviving paralogs: ALDH1A and retinoic acid signaling in vertebrate genomes

Genome duplications increase genetic diversity and may facilitate the evolution of gene subfunctions. Little attention, however, has focused on the evolutionary impact of lineage-specific gene loss. Here, we show that identifying lineage-specific gene loss after genome duplication is important for understanding the evolution of gene subfunctions in surviving paralogs and for improving functional connectivity among human and model organism genomes. We examine the general principles of gene loss following duplication, coupled with expression analysis of the retinaldehyde dehydrogenase Aldh1a gene family during retinoic acid signaling in eye development as a case study. Humans have three ALDH1A genes, but teleosts have just one or two. We used comparative genomics and conserved syntenies to identify loss of ohnologs (paralogs derived from genome duplication) and to clarify uncertain phylogenies. Analysis showed that Aldh1a1 and Aldh1a2 form a clade that is sister to Aldh1a3-related genes. Genome comparisons showed secondarily loss of aldh1a1 in teleosts, revealing that Aldh1a1 is not a tetrapod innovation and that aldh1a3 was recently lost in medaka, making it the first known vertebrate with a single aldh1a gene. Interestingly, results revealed asymmetric distribution of surviving ohnologs between co-orthologous teleost chromosome segments, suggesting that local genome architecture can influence ohnolog survival. We propose a model that reconstructs the chromosomal history of the Aldh1a family in the ancestral vertebrate genome, coupled with the evolution of gene functions in surviving Aldh1a ohnologs after R1, R2, and R3 genome duplications. Results provide evidence for early subfunctionalization and late subfunction-partitioning and suggest a mechanistic model based on altered regulation leading to heterochronic gene expression to explain the acquisition or modification of subfunctions by surviving ohnologs that preserve unaltered ancestral developmental programs in the face of gene loss.

Lineage-specific transcription factor aberrations in AML

Transcription factors play a key role in the commitment of hematopoietic stem cells to differentiate into specific lineages [78]. This is particularly important in that a block in terminal differentiation is the key contributing factor in acute leukemias. This general theme of the role of transcription factors in differentiation may also extend to other tissues, both in terms of normal development and cancer. Consistent with the role of transcription factors in hematopoietic lineage commitment is the frequent finding of aberrations in transcription factors in AML patients. Here, we intend to review recent findings on aberrations in lineage-restricted transcription factors as observed in patients with acute myeloid leukemia (AML).

Arabinosylguanine: Metabolism, action, and lineage-specific cytotoxicity in human leukemia cells

9-β-D-arabinosylguanine (ara-G), an analogue of deoxyguanosine, has demonstrated T-lymphoblast selective anti-leukemia activity both in vitro and in vivo in cell lines and primary cells and in phase I investigations. The present work was initiated to identify factors that result in this selectivity. ^ The cytotoxicity of ara-G is manifest only after its phosphorylation. Experiments using cell lines transfected to overexpress specific nucleoside kinases demonstrated that the phosphorylation of ara-G to its monophosphate is by both cytoplasmic deoxycytidine kinase and mitochondria) deoxyguanosine kinase. Ara-G monophosphate is converted to its 5′-triphosphate (ara-GTP) in cells by these kinases and then incorporated into DNA. Mechanistic studies demonstrated that incorporation of ara-GTP into DNA was a necessary event for the induction of cell death. ^ Pharmacokinetic and pharmacodynamic studies utilizing three human acute leukemia cell lines, CEM (T-lymphoblastic), Raji (B-lymphoblastic), and ML-1 (myeloid) were performed. CEM cells were most sensitive to ara-G-induced inhibition of colony formation, accumulated ara-GTP at a faster rate and to a greater degree than either Raji or ML-1, but incorporated the lowest number of ara-G molecules into DNA. The position of incorporation was internal and similar in all cell lines. The terminal elimination phase of ara-GTP was >24 h and similar in these cells. Comparisons between inhibition of colony formation and ara-GTP incorporation into DNA demonstrated that while within a cell line there was correlation among these parameters, between cell lines there was no relationship between number of incorporated ara-G molecules and ara-G(TP)-mediated toxicity suggesting that there were additional factors. ^ The expression of membrane bound Fas and Fast was unchanged in all cell lines. In contrast, there was a 2-fold increase in soluble Fast, which was found exclusively in CEM cells. Ara-G-mediated apoptosis in CEM occurred from all phases of the cell cycle and was abrogated partially by Fas antagonist antibodies. These data suggest that Fas-mediated cell death due to the liberation of sFasL may be responsible for the hypersensitivity to ara-G manifested by immature T-cells such as CEM. The role of Fas in ara-G induced death of acute T-lymphoblastic leukemia cells during therapy needs to be tested. ^

Minimal DNA sequences that control the cell lineage-specific expression of the pro$\alpha$2(I) collagen promoter in transgenic mice

The pattern of expression of the pro$\alpha$2(I) collagen gene is highly tissue-specific in adult mice and shows its strongest expression in bones, tendons, and skin. Transgenic mice were generated harboring promoter fragments of the mouse pro$\alpha$2(I) collagen gene linked to the Escherichia coli $\beta$-galactosidase or firefly luciferase genes to examine the activity of these promoters during development. A region of the mouse pro$\alpha$2(I) collagen promoter between $-$2000 and +54 exhibited a pattern of $\beta$-galactosidase activity during embryonic development that corresponded to the expression pattern of the endogenous pro$\alpha$2(I) collagen gene as determined by in situ hybridization. A similar pattern of activity was also observed with much smaller promoter fragments containing either 500 or 350 bp of upstream sequence relative to the start of transcription. Embryonic regions expressing high levels of $\beta$-galactosidase activity included the valves of the developing heart, sclerotomes, meninges, limb buds, connective tissue fascia between muscle fibers, osteoblasts, tendon, periosteum, dermis, and peritoneal membranes. The pattern of $\beta$-galactosidase activity was similar to the extracellular immunohistochemical localization of transforming growth factor-$\beta$1 (TGF-$\beta$1). The $-$315 to $-$284 region of the pro$\alpha$2(I) collagen promoter was previously shown to mediate the stimulatory effects of TGF-$\beta$1 on the pro$\alpha$2(I) collagen promoter in DNA transfection experiments with cultured fibroblasts. A construct containing this sequence tandemly repeated 5$\sp\prime$ to both a very short $\alpha$2(I) collagen promoter ($-$40 to +54) and a heterologous minimal promoter showed preferential activity in tail and skin of 4-week old transgenic mice. The pattern of expression mimics that of the $-$350 to +54 pro$\alpha$2(I) collagen promoter linked to a luciferase reporter gene in transgenic mice. ^

Transcriptional regulation of cell lineage-specific expression of the murine type I collagen genes and of osteoblast differentiation

The aim of my project is to examine the mechanisms of cell lineage-specific transcriptional regulation of the two type I collagen genes by characterizing critical cis-acting elements and trans-acting factors. I hypothesize that the transcription factors that are involved in the cell lineage-specific expression of these genes may have a larger essential role in cell lineage commitment and differentiation. I first examined the proximal promoters of the proα1(I) and the proα2(I) collagen genes for cell type-specific DNA-protein interactions, using in vitro DNaseI and in vivo DMS footprinting. These experiments demonstrated that the cis-acting elements in these promoters are accessible to ubiquitous DNA-binding proteins in fibroblasts that express these genes, but not in other cells that do not express these genes. I speculate that in type I collagen-expressing cells, cell type-specific enhancer elements facilitate binding of ubiquitous proteins to the proximal promoters of these genes. Subsequently, examination of the upstream promoter of the proα(I) collagen gene by transgenic mice experiments delineated a 117 bp sequence (-1656 to -1540 bp) as the minimum element required for osteoblast-specific expression. This 117 bp element contained two segments that appeared to have different functions: (1) the A-segment, which was necessary to obtain osteoblast-specific expression and (2) the C-segment, which was dispensable for osteoblast-specific expression, but was necessary to obtain high-level expression. In experiments to identify trans-acting factors that bind to the 117 bp element, I have demonstrated that the cell lineage-restricted homeodomain proteins, Dlx2, Dlx5 and mHOX, bound to the A-segment and that the ubiquitous transcription factor, Sp1, bound to the C-segment of this element. These results suggested a model where the binding of cell lineage-restricted proteins to the A-segment and of ubiquitous proteins to the C-segment of the 117 bp element of the proα1 (I) collagen gene activated this gene in osteoblasts. These results, combined with additional evidence that Dlx2, Dlx5 and mHOX are probably involved in osteoblast differentiation, support my hypothesis that the transcription factors involved in osteoblast-specific expression of type I collagen genes may have essential role in osteoblast lineage commitment and differentiation. ^

A lineage-specific protein kinase crucial for myeloid maturation

To identify genes involved in macrophage development, we used the differential display technique and compared the gene expression profiles for human myeloid HL-60 leukemia cell lines susceptible and resistant to macrophage maturation. We identified a gene coding for a protein kinase, protein kinase X (PRKX), which was expressed in the maturation-susceptible, but not in the resistant, cell line. The expression of the PRKX gene was found to be induced during monocyte, macrophage, and granulocyte maturation of HL-60 cells. We also studied the expression of the PRKX gene in 12 different human tissues and transformed cell lines and found that, among these tissues and cell types, the PRKX gene is expressed only in blood. Among the blood cell lineages, the PRKX gene is specifically expressed in macrophages and granulocytes. Antisense inhibition of PRKX expression blocked terminal development in both the leukemic HL-60 cells and normal peripheral blood monocytes, implying that PRKX is a key mediator of macrophage and granulocyte maturation. Using the HL-60 cell variant deficient in protein kinase C-β (PKC-β) and several stable PKC-β transfectants, we found that PRKX gene expression is under control of PKC-β; hence PRKX is likely to act downstream of this PKC isozyme in the same signal transduction pathway leading to macrophage maturation.

An experimental test for lineage-specific position effects on alcohol dehydrogenase (Adh) genes in Drosophila

Independent transgene insertions differ in expression based on their location in the genome; these position effects are of interest because they reflect the influence of genome organization on gene regulation. Position effects also represent potentially insurmountable obstacles to the rigorous functional comparison of homologous genes from different species because (i) quantitative variation in expression of each gene across genomic positions (generalized position effects, or GPEs) may overwhelm differences between the genes of interest, or (ii) divergent genes may be differentially sensitive to position effects, reflecting unique interactions between each gene and its genomic milieu (lineage-specific position effects, or LSPEs). We have investigated both types of position-effect variation by applying our method of transgene coplacement, which allows comparisons of transgenes in the same position in the genome of Drosophila melanogaster. Here we report an experimental test for LSPE in Drosophila. The alcohol dehydrogenase (Adh) genes of D. melanogaster and Drosophila affinidisjuncta differ in both tissue distribution and amounts of ADH activity. Despite this striking regulatory divergence, we found a very high correlation in overall ADH activity between the genes of the two species when placed in the same genomic position as assayed in otherwise Adh-null adults and larvae. These results argue against the influence of LSPE for these sequences, although the effects of GPE are significant. Our new findings validate the coplacement approach and show that it greatly magnifies the power to detect differences in expression between transgenes. Transgene coplacement thus dramatically extends the range of functional and evolutionary questions that can be addressed by transgenic technology.

Development of prospective mathematics teachers’ professional noticing in a specific domain: proportional reasoning

The aim of this research is to identify aspects that support the development of prospective mathematics teachers’ professional noticing in a b-learning context. The study presented here investigates the extent to which prospective secondary mathematics teachers attend and interpret secondary school students’ proportional reasoning and decide how to respond. Results show that interactions in an on-line discussion improve prospective mathematics teachers’ ability to identify and interpret important aspects of secondary school students’ mathematical thinking.

Emergence of novel domains in proteins

Proteins are composed of a combination of discrete, well-defined, sequence domains, associated with specific functions that have arisen at different times during evolutionary history. The emergence of novel domains is related to protein functional diversification and adaptation. But currently little is known about how novel domains arise and how they subsequently evolve. To gain insights into the impact of recently emerged domains in protein evolution we have identified all human young protein domains that have emerged in approximately the past 550 million years. We have classified them into vertebrate-specific and mammalian-specific groups, and compared them to older domains. We have found 426 different annotated young domains, totalling 995 domain occurrences, which represent about 12.3% of all human domains. We have observed that 61.3% of them arose in newly formed genes, while the remaining 38.7% are found combined with older domains, and have very likely emerged in the context of a previously existing protein. Young domains are preferentially located at the N-terminus of the protein, indicating that, at least in vertebrates, novel functional sequences often emerge there. Furthermore, young domains show significantly higher non-synonymous to synonymous substitution rates than older domains using human and mouse orthologous sequence comparisons. This is also true when we compare young and old domains located in the same protein, suggesting that recently arisen domains tend to evolve in a less constrained manner than older domains. We conclude that proteins tend to gain domains over time, becoming progressively longer. We show that many proteins are made of domains of different age, and that the fastest evolving parts correspond to the domains that have been acquired more recently.

Unique features of odorant-binding proteins of the parasitoid wasp Nasonia vitripennis revealed by genome annotation and comparative analyses

Insects are the most diverse group of animals on the planet, comprising over 90% of all metazoan life forms, and have adapted to a wide diversity of ecosystems in nearly all environments. They have evolved highly sensitive chemical senses that are central to their interaction with their environment and to communication between individuals. Understanding the molecular bases of insect olfaction is therefore of great importance from both a basic and applied perspective. Odorant binding proteins (OBPs) are some of most abundant proteins found in insect olfactory organs, where they are the first component of the olfactory transduction cascade, carrying odorant molecules to the olfactory receptors. We carried out a search for OBPs in the genome of the parasitoid wasp Nasonia vitripennis and identified 90 sequences encoding putative OBPs. This is the largest OBP family so far reported in insects. We report unique features of the N. vitripennis OBPs, including the presence and evolutionary origin of a new subfamily of double-domain OBPs (consisting of two concatenated OBP domains), the loss of conserved cysteine residues and the expression of pseudogenes. This study also demonstrates the extremely dynamic evolution of the insect OBP family: (i) the number of different OBPs can vary greatly between species; (ii) the sequences are highly diverse, sometimes as a result of positive selection pressure with even the canonical cysteines being lost; (iii) new lineage specific domain arrangements can arise, such as the double domain OBP subfamily of wasps and mosquitoes.