Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

Exome sequencing identifies INPPL1 mutations as a cause of opsismodysplasia.

Opsismodysplasia (OPS) is a severe autosomal-recessive chondrodysplasia characterized by pre- and postnatal micromelia with extremely short hands and feet. The main radiological features are severe platyspondyly, squared metacarpals, delayed skeletal ossification, and metaphyseal cupping. In order to identify mutations causing OPS, a total of 16 cases (7 terminated pregnancies and 9 postnatal cases) from 10 unrelated families were included in this study. We performed exome sequencing in three cases from three unrelated families and only one gene was found to harbor mutations in all three cases: inositol polyphosphate phosphatase-like 1 (INPPL1). Screening INPPL1 in the remaining cases identified a total of 12 distinct INPPL1 mutations in the 10 families, present at the homozygote state in 7 consanguinous families and at the compound heterozygote state in the 3 remaining families. Most mutations (6/12) resulted in premature stop codons, 2/12 were splice site, and 4/12 were missense mutations located in the catalytic domain, 5-phosphatase. INPPL1 belongs to the inositol-1,4,5-trisphosphate 5-phosphatase family, a family of signal-modulating enzymes that govern a plethora of cellular functions by regulating the levels of specific phosphoinositides. Our finding of INPPL1 mutations in OPS, a severe spondylodysplastic dysplasia with major growth plate disorganization, supports a key and specific role of this enzyme in endochondral ossification.

A genomewide linkage study on suicidality in major depressive disorder confirms evidence for linkage to 2p12.

Suicidal behavior is commonly associated with depression. Twin studies indicate that both suicidality and major depressive disorder (MDD) are heritable. However, epidemiological evidence suggests that the inheritance of suicidality is likely to be independent of the underlying psychiatric disorder, implying a distinct genetic contribution to suicidality. We conducted a genomewide linkage search aiming to detect genomic loci that may harbor susceptibility genes contributing to risk for suicidality in recurrent MDD. Affected sibling pair (ASP) variance components analysis was performed using the Depression Network cohort of 971 ASPs. The quantitative trait measuring suicidality as a broad phenotype, encompassing ideation and suicide attempts, was established from Schedules for Clinical Assessment in Neuropsychiatry interview items. We examined 1,060 genotyped microsatellite markers with an average spacing of 3.3 cM. Empirical thresholds for linkage evidence were set by whole-genome simulations (LOD = 2.71 for genomewide significance, 1.71 for suggestive linkage). No genomewide significant findings were found. Marker D3S1234 on 3p14 achieved suggestive linkage and yielded a maximum LOD of 1.853 (P = 0.0017), loci 9p24.3 and 18q22-q23 achieved LOD scores >1.5. We found some support for linkage to 2p12 (LOD = 1.2, P = 0.0087) which was previously implicated in linkage studies of suicidality. Our follow-up meta-analysis of five studies showed strong linkage to this region (P = 2 × 10(-6) ). In conclusion, this study analyzed suicidality as a continuous trait in MDD. We found modest evidence for linkage on 3p14. Our meta-analysis supports previous evidence of linkage to suicidality on 2p12. Some candidate genes in these regions may plausibly be implicated in suicidality.

Low predictive value of elevated blood pressure during childhood and adolescence

Background: Screening of elevated blood pressure (BP) in children has been advocated to early identify hypertension. However, identification of children with sustained elevated BP is challenging due to the high BP variability. The value of an elevated BP measure during childhood and adolescence for the prediction of future elevated BP is not well described. Objectives: We assessed the positive (PPV) and negative (NPV) predictive value of high BP for sustained elevated BP in cohorts of children of the Seychelles, a rapidly developing island state in the African region. Methods: Serial school-based surveys of weight, height, and BP were conducted yearly between 1998-2006 among all students of the country in four school grades (kindergarten [G0, mean age (SD): 5.5 (0.4) yr], G4 [9.2 (0.4) yr], G7 [12.5 (0.4) yr] and G10 (15.6 (0.5) yr]. We constituted three cohorts of children examined twice at 3-4 years interval: 4,557 children examined at G0 and G4, 6,198 at G4 and G7, and 6,094 at G7 and G10. The same automated BP measurement devices were used throughout the study. BP was measured twice at each exam and averaged. Obesity and elevated BP were defined using the CDC (BMI_95th sex-, and age-specific percentile) and the NHBPEP criteria (BP_95th sex-, age-, and height specific percentile), respectively. Results: Prevalence of obesity was 6.1% at G0, 7.1% at G4, 7.5% at G7, and 6.5% at G10. Prevalence of elevated BP was 10.2% at G0, 9.9% at G4, 7.1% at G7, and 8.7% at G10. Among children with elevated BP at initial exam, the PPV of keeping elevated BP was low but increased with age: 13% between G0 and G4, 19% between G4 and G7, and 27% between G7 and G10. Among obese children with elevated BP, the PPV was higher: 33%, 35% and 39% respectively. Overall, the probability for children with normal BP to remain in that category 3-4 years later (NPV) was 92%, 95%, and 93%, respectively. By comparison, the PPV for children initially obese to remain obese was much higher at 71%, 71%, and 62% (G7-G10), respectively. The NPV (i.e. the probability of remaining at normal weight) was 94%, 96%, and 98%, respectively. Conclusion: During childhood and adolescence, having an elevated BP at one occasion is a weak predictor of sustained elevated BP 3-4 years later. In obese children, it is a better predictor.

Suggestive linkage to chromosome 1q for bone mineral apparent density in Brazilian sister adolescents.

OBJECTIVE: To investigate linkage to chromosome 1q and 11q region for lumbar spine, femoral neck and total body BMD and volumetric BMD in Brazilian sister adolescents aged 10-20-year-old and 57 mothers. METHODS: We evaluated 161 sister pairs (n=329) aged 10-20 years old and 57 of their mothers in this study. Physical traits and lifestyle factors were collected as covariates for lumbar spine (LS), femoral neck (FN) and total body (TB) BMD and bone mineral apparent density (BMAD). We selected nine microsatellite markers in chromosome 1q region (spanning nearly 33cM) and eight in chromosome 11q region (spanning nearly 34cM) to perform linkage analysis. RESULTS: The highest LOD score values obtained from our data were in sister pairs LS BMAD analysis. Their values were: 1.32 (P<0.006), 2.61 (P<0.0002) and 2.44 (P<0.0004) in D1S218, D1S2640 and D1S2623 markers, respectively. No significant LOD score was found with LS and FN BMD/BMAD in chromosome 11q region. Only TB BMD showed significant linkage higher than 1.0 for chromosome 11q region in the markers D11S4191 and D11S937. DISCUSSION/CONCLUSIONS: Our results provided suggestive linkage for LS BMAD at D1S2640 marker in adolescent sister pairs and suggest a possible candidate gene (LHX4) related to adolescent LS BMAD in this region. These results reinforce chromosome 1q21-23 as a candidate region to harbor one or more bone formation/maintenance gene. In the other hand, it did not repeat for chromosome 11q12-13 in our population.

Defective tumor necrosis factor-alpha-dependent control of astrocyte glutamate release in a transgenic mouse model of Alzheimer disease.

The cytokine tumor necrosis factor-alpha (TNFalpha) induces Ca2+-dependent glutamate release from astrocytes via the downstream action of prostaglandin (PG) E2. By this process, astrocytes may participate in intercellular communication and neuromodulation. Acute inflammation in vitro, induced by adding reactive microglia to astrocyte cultures, enhances TNFalpha production and amplifies glutamate release, switching the pathway into a neurodamaging cascade (Bezzi, P., Domercq, M., Brambilla, L., Galli, R., Schols, D., De Clercq, E., Vescovi, A., Bagetta, G., Kollias, G., Meldolesi, J., and Volterra, A. (2001) Nat. Neurosci. 4, 702-710). Because glial inflammation is a component of Alzheimer disease (AD) and TNFalpha is overexpressed in AD brains, we investigated possible alterations of the cytokine-dependent pathway in PDAPP mice, a transgenic model of AD. Glutamate release was measured in acute hippocampal and cerebellar slices from mice at early (4-month-old) and late (12-month-old) disease stages in comparison with age-matched controls. Surprisingly, TNFalpha-evoked glutamate release, normal in 4-month-old PDAPP mice, was dramatically reduced in the hippocampus of 12-month-old animals. This defect correlated with the presence of numerous beta-amyloid deposits and hypertrophic astrocytes. In contrast, release was normal in cerebellum, a region devoid of beta-amyloid deposition and astrocytosis. The Ca2+-dependent process by which TNFalpha evokes glutamate release in acute slices is distinct from synaptic release and displays properties identical to those observed in cultured astrocytes, notably PG dependence. However, prostaglandin E2 induced normal glutamate release responses in 12-month-old PDAPP mice, suggesting that the pathology-associated defect involves the TNFalpha-dependent control of secretion rather than the secretory process itself. Reduced expression of DENN/MADD, a mediator of TNFalpha-PG coupling, might account for the defect. Alteration of this neuromodulatory astrocytic pathway is described here for the first time in relation to Alzheimer disease.

Imaging exocytosis and recycling of synaptic-like microvesicles in astrocytes.

Optical imaging techniques are well suited for following the dynamics of physiological processes in living cells. Total internal reflection fluorescence (TIRF) microscopy based on evanescent wave illumination (EWi) allows spectacular, real-time visualization of individual vesicle movements, fusions, and retrievals at the cell surface (i.e., within 100 nm of the plasma membrane). TIRF microscopy is an ideal approach for studying the properties of exocytosis and recycling in cultured astrocytes, particularly because these cells have a rather flat surface and contain secretory vesicles with sparse distribution. Among all populations of secretory vesicles, we focus here on synaptic-like microvesicles (SLMVs). We illustrate how TIRF microscopy using EWi is useful to study exocytosis and recycling of SLMVs at the single-vesicle level and, when combined with epifluorescence illumination (EPIi), can provide detailed information on the kinetics of exocytosis, endocytosis, and re-acidification at the whole-cell level.

Changes in arbuscular mycorrhizal fungal phenotypes and genotypes in response to plant species identity and phosphorus concentration.

* Arbuscular mycorrhizal fungi (AMF) are plant symbionts that improve floristic diversity and ecosystem productivity. Many AMF species are generalists with wide host ranges. Arbuscular mycorrhizal fungi individuals are heterokaryotic, and AMF populations are genetically diverse. Populations of AMF harbor two levels of genetic diversity on which selection can act, namely among individuals and within individuals. Whether environmental factors alter genetic diversity within populations is still unknown. * Here, we measured genetic changes and changes in fitness-related traits of genetically distinct AMF individuals from one field, grown with different concentrations of available phosphate or different host species. * We found significant genotype-by-environment interactions for AMF fitness traits in response to these treatments. Host identity had a strong effect on the fitness of different AMF, unearthing a specificity of response within Glomus intraradices. Arbuscular mycorrhizal fungi individuals grown in novel environments consistently showed a reduced presence of polymorphic genetic markers, providing some evidence for host or phosphate-induced genetic change in AMF. * Given that AMF individuals can form extensive hyphal networks colonizing different hosts simultaneously, contrasting habitats or soil properties may lead to evolution in the population. Local selection may alter the structure of AMF populations and maintain genetic diversity, potentially even within the hyphal network of one fungus.

Pituitary blastoma: a pathognomonic feature of germ-line DICER1 mutations.

Individuals harboring germ-line DICER1 mutations are predisposed to a rare cancer syndrome, the DICER1 Syndrome or pleuropulmonary blastoma-familial tumor and dysplasia syndrome [online Mendelian inheritance in man (OMIM) #601200]. In addition, specific somatic mutations in the DICER1 RNase III catalytic domain have been identified in several DICER1-associated tumor types. Pituitary blastoma (PitB) was identified as a distinct entity in 2008, and is a very rare, potentially lethal early childhood tumor of the pituitary gland. Since the discovery by our team of an inherited mutation in DICER1 in a child with PitB in 2011, we have identified 12 additional PitB cases. We aimed to determine the contribution of germ-line and somatic DICER1 mutations to PitB. We hypothesized that PitB is a pathognomonic feature of a germ-line DICER1 mutation and that each PitB will harbor a second somatic mutation in DICER1. Lymphocyte or saliva DNA samples ascertained from ten infants with PitB were screened and nine were found to harbor a heterozygous germ-line DICER1 mutation. We identified additional DICER1 mutations in nine of ten tested PitB tumor samples, eight of which were confirmed to be somatic in origin. Seven of these mutations occurred within the RNase IIIb catalytic domain, a domain essential to the generation of 5p miRNAs from the 5' arm of miRNA-precursors. Germ-line DICER1 mutations are a major contributor to PitB. Second somatic DICER1 "hits" occurring within the RNase IIIb domain also appear to be critical in PitB pathogenesis.

Host Genes Important to HIV Replication and Evolution.

Recent years have seen a significant increase in understanding of the host genetic and genomic determinants of susceptibility to HIV-1 infection and disease progression, driven in large part by candidate gene studies, genome-wide association studies, genome-wide transcriptome analyses, and large-scale in vitro genome screens. These studies have identified common variants in some host loci that clearly influence disease progression, characterized the scale and dynamics of gene and protein expression changes in response to infection, and provided the first comprehensive catalogs of genes and pathways involved in viral replication. Experimental models of AIDS and studies in natural hosts of primate lentiviruses have complemented and in some cases extended these findings. As the relevant technology continues to progress, the expectation is that such studies will increase in depth (e.g., to include host whole exome and whole genome sequencing) and in breadth (in particular, by integrating multiple data types).

Comparative functional analysis of two fibroblast growth factor receptor 1 (FGFR1) mutations affecting the same residue (R254W and R254Q) in isolated hypogonadotropic hypogonadism (IHH).

FGFR1 mutations have been identified in both Kallmann syndrome and normosmic HH (nIHH). To date, few mutations in the FGFR1 gene have been structurally or functionally characterized in vitro to identify molecular mechanisms that contribute to the disease pathogenesis. We attempted to define the in vitro functionality of two FGFR1 mutants (R254W and R254Q), resulting from two different amino acid substitutions of the same residue, and to correlate the in vitro findings to the patient phenotypes. Two unrelated GnRH deficient probands were found to harbor mutations in FGFR1 (R254W and R254Q). Mutant signaling activity and expression levels were evaluated in vitro and compared to a wild type (WT) receptor. Signaling activity was determined by a FGF2/FGFR1 dependent transcription reporter assay. Receptor total expression levels were assessed by Western blot and cell surface expression was measured by a radiolabeled antibody binding assay. The R254W maximal receptor signaling capacity was reduced by 45% (p<0.01) while R254Q activity was not different from WT. However, both mutants displayed diminished total protein expression levels (40 and 30% reduction relative to WT, respectively), while protein maturation was unaffected. Accordingly, cell surface expression levels of the mutant receptors were also significantly reduced (35% p<0.01 and 15% p<0.05, respectively). The p.R254W and p.R254Q are both loss-of-function mutations as demonstrated by their reduced overall and cell surface expression levels suggesting a deleterious effect on receptor folding and stability. It appears that a tryptophan substitution at R254 is more disruptive to receptor structure than the more conserved glutamine substitution. No clear correlation between the severity of in vitro loss-of-function and phenotypic presentation could be assigned.

Characterizing variation in mycorrhiza effect among diverse plant varieties.

Exploitation of arbuscular mycorrhizal fungi may be an important approach for development of reduced-input agriculture. We discuss the use of linear models to analyze variation in mycorrhiza response among diverse plant varieties in order to assess the value of mycorrhizas. Our approach allows elimination of variation linked to differences in plant performance in the absence of mycorrhizas and the selection of plant lines that might harbor genetic variation of use to improve the mycorrhizal symbiosis in agriculture. We illustrate our approach by applying it to previously published and to novel data. We suggest that in dealing with a relative trait such as mycorrhiza effect, the choice of measure used to quantify the trait greatly affects interpretation. In the plant populations under consideration, we find evidence for a greater potential to increase mycorrhiza benefit than previously suggested.

The Th2 lymphoproliferation developing in LatY136F mutant mice triggers polyclonal B cell activation and systemic autoimmunity.

Lat(Y136F) knock-in mice harbor a point mutation in Tyr(136) of the linker for activation of T cells and show accumulation of Th2 effector cells and IgG1 and IgE hypergammaglobulinemia. B cell activation is not a direct effect of the mutation on B cells since in the absence of T cells, mutant B cells do not show an activated phenotype. After adoptive transfer of linker for activation of T cell mutant T cells into wild-type, T cell-deficient recipients, recipient B cells become activated. We show in vivo and in vitro that the Lat(Y136F) mutation promotes T cell-dependent B cell activation leading to germinal center, memory, and plasma cell formation even in an MHC class II-independent manner. All the plasma and memory B cell populations found in physiological T cell-dependent B cell responses are found. Characterization of the abundant plasmablasts found in secondary lymphoid organs of Lat(Y136F) mice revealed the presence of a previously uncharacterized CD93-expressing subpopulation, whose presence was confirmed in wild-type mice after immunization. In Lat(Y136F) mice, B cell activation was polyclonal and not Ag-driven because the increase in serum IgG1 and IgE concentrations involved Abs and autoantibodies with different specificities equally. Although the noncomplement-fixing IgG1 and IgE are the only isotypes significantly increased in Lat(Y136F) serum, we observed early-onset systemic autoimmunity with nephritis showing IgE autoantibody deposits and severe proteinuria. These results show that Th2 cells developing in Lat(Y136F) mice can trigger polyclonal B cell activation and thereby lead to systemic autoimmune disease.

Hematopoietic stem cells reversibly switch from dormancy to self-renewal during homeostasis and repair.

Bone marrow hematopoietic stem cells (HSCs) are crucial to maintain lifelong production of all blood cells. Although HSCs divide infrequently, it is thought that the entire HSC pool turns over every few weeks, suggesting that HSCs regularly enter and exit cell cycle. Here, we combine flow cytometry with label-retaining assays (BrdU and histone H2B-GFP) to identify a population of dormant mouse HSCs (d-HSCs) within the lin(-)Sca1+cKit+CD150+CD48(-)CD34(-) population. Computational modeling suggests that d-HSCs divide about every 145 days, or five times per lifetime. d-HSCs harbor the vast majority of multilineage long-term self-renewal activity. While they form a silent reservoir of the most potent HSCs during homeostasis, they are efficiently activated to self-renew in response to bone marrow injury or G-CSF stimulation. After re-establishment of homeostasis, activated HSCs return to dormancy, suggesting that HSCs are not stochastically entering the cell cycle but reversibly switch from dormancy to self-renewal under conditions of hematopoietic stress.

SCRIB and PUF60 Are Primary Drivers of the Multisystemic Phenotypes of the 8q24.3 Copy-Number Variant.

Copy-number variants (CNVs) represent a significant interpretative challenge, given that each CNV typically affects the dosage of multiple genes. Here we report on five individuals with coloboma, microcephaly, developmental delay, short stature, and craniofacial, cardiac, and renal defects who harbor overlapping microdeletions on 8q24.3. Fine mapping localized a commonly deleted 78 kb region that contains three genes: SCRIB, NRBP2, and PUF60. In vivo dissection of the CNV showed discrete contributions of the planar cell polarity effector SCRIB and the splicing factor PUF60 to the syndromic phenotype, and the combinatorial suppression of both genes exacerbated some, but not all, phenotypic components. Consistent with these findings, we identified an individual with microcephaly, short stature, intellectual disability, and heart defects with a de novo c.505C>T variant leading to a p.His169Tyr change in PUF60. Functional testing of this allele in vivo and in vitro showed that the mutation perturbs the relative dosage of two PUF60 isoforms and, subsequently, the splicing efficiency of downstream PUF60 targets. These data inform the functions of two genes not associated previously with human genetic disease and demonstrate how CNVs can exhibit complex genetic architecture, with the phenotype being the amalgam of both discrete dosage dysfunction of single transcripts and also of binary genetic interactions.