EHMT2 suppresses the variation of transcriptional switches in the mouse embryo

EHMT2 is the main euchromatic H3K9 methyltransferase. Embryos with zygotic, or maternal mutation in the Ehmt2 gene exhibit variable developmental delay. To understand how EHMT2 prevents variable developmental delay we performed RNA sequencing of mutant and somite stage-matched normal embryos at 8.5–9.5 days of gestation. Using four-way comparisons between delayed and normal embryos we clarified what it takes to be normal and what it takes to develop. We identified differentially expressed genes, for example Hox genes that simply reflected the difference in developmental progression of wild type and the delayed mutant uterus-mate embryos. By comparing wild type and zygotic mutant embryos along the same developmental window we detected a role of EHMT2 in suppressing variation in the transcriptional switches. We identified transcription changes where precise switching during development occurred only in the normal but not in the mutant embryo. At the 6-somite stage, gastrulation-specific genes were not precisely switched off in the Ehmt2−/− zygotic mutant embryos, while genes involved in organ growth, connective tissue development, striated muscle development, muscle differentiation, and cartilage development were not precisely switched on. The Ehmt2mat−/+ maternal mutant embryos displayed high transcriptional variation consistent with their variable survival. Variable derepression of transcripts occurred dominantly in the maternally inherited allele. Transcription was normal in the parental haploinsufficient wild type embryos despite their delay, consistent with their good prospects. Global profiling of transposable elements revealed EHMT2 targeted DNA methylation and suppression at LTR repeats, mostly ERVKs. In Ehmt2−/− embryos, transcription over very long distances initiated from such misregulated ‘driver’ ERVK repeats, encompassing a multitude of misexpressed ‘passenger’ repeats. In summary, EHMT2 reduced transcriptional variation of developmental switch genes and developmentally switching repeat elements at the six-somite stage embryos. These findings establish EHMT2 as a suppressor of transcriptional and developmental variation at the transition between gastrulation and organ specification.

Functionally Mature CD1c+ Dendritic Cells Preferentially Accumulate in the Inflammatory Arthritis Synovium

ObjectiveTo examine the role of synovial CD1c+DCs in patients with Inflammatory Arthritis (IA) with a specific focus on the transcriptional and maturation signatures that govern their function.MethodsRNA sequencing was performed on healthy control (HC) peripheral blood (PB), IA PB, and IA synovial fluid (SF) CD1c+DCs. Multiparametric flow-cytometry and SPICE analysis were used to examine site [SF and Synovial Tissue (ST) CD1c+DCs] and disease specific characteristics of CD1c+DCs, while functional assays such as antigen processing, activation, and MMP production were also performed.ResultsIncreased frequency of CD1c+DCs (p<0.01) with a concomitant increase in CD80, CCR7 (p<0.01), and CXCR3 (p<0.05) expression was identified in IA PB compared to HC PB. Enrichment of CD1c+DCs was identified in IA synovial tissue (ST) (p<0.01) and IA SF (p<0.0001) compared to IA PB, while RNAseq revealed distinct transcriptional variation between PB and SF CD1c+DCs. Flow cytometry revealed increased expression of CD83, CD80, PD-L1, and BTLA (all p<0.05) in IA SF CD1c+DCs compared to PB, while SPICE identified synovial cells with unique co-expression patterns, expressing multiple DC maturation markers simultaneously. Functionally, synovial CD1c+DCs are hyper-responsive to TLR7/8 ligation (p<0.05), have decreased antigen processing capacity (p=0.07), and display dysregulated production of MMPs. Finally, examination of both synovial CD1c+DCs and synovial CD141+DCs revealed distinct maturation and transcriptomic profiles.ConclusionSynovial CD1c+DCs accumulate in the inflamed IA synovium in a variety of distinct poly-maturational states, distinguishing them transcriptionally and functionally from CD1c+DCs in the periphery and synovial CD141+DCs.

Comparative transcriptomic analysis of maize ear heterosis during the inflorescence meristem differentiation stage

Heterosis is widely used in many crops; however, its genetic mechanisms are only partly understood. Here, we sampled inflorescence meristem (IM) ears from the single-segment substitution maize (Zea mays) line lx9801hlEW2b, containing a heterotic locus hlEW2b associated with ear width, the receptor parent lx9801, the test parent Zheng58, and their corresponding hybrids. After transcriptomic analysis, 1638 genes were identified in at least one hybrid with nonadditively expressed patterns and different expression levels between the two hybrids. In particular, 2263 (12.89%) and 2352 (14.65%) genes showed allele-specific expression (ASE) in Zheng58 × lx9801 and Zheng58 × lx9801hlEW2b, respectively. A functional analysis showed that these genes were enriched in development-related processes and biosynthesis and catabolism processes, which are potentially associated with heterosis. Additionally, nonadditive expression and ASE may fine-tune the expression levels of crucial genes (such as WUS and KNOX that control IM development) controlling auxin metabolism and ear development to optimal states, and transcriptional variation may play important roles in maize ear heterosis. The results provide new information that increases our understanding of the relationship between transcriptional variation and heterosis formation during maize ear development, which may be helpful in clarifying the genetic and molecular mechanisms of heterosis.

Intraspecific transcriptome variation and sex-biased expression in Anopheles arabiensis

The magnitude and functional patterns of intraspecific transcriptional variation in the anophelines, including those of sex-biased genes underlying sex-specific traits relevant for malaria transmission, remain understudied. As a result, how changes in expression levels drive adaptation in these species is poorly understood. We sequenced the female, male, and larval transcriptomes of three populations of Anopheles arabiensis from Burkina Faso. One-third of the genes were differentially expressed between populations, often involving insecticide resistance-related genes in a sample type-specific manner, and with the females showing the largest number of differentially expressed genes. At the genomic level, the X chromosome appears depleted of differentially expressed genes compared to the autosomes, chromosomes harbouring inversions do not exhibit evidence for enrichment of such genes, and genes that are top contributors to functional enrichment patterns of population differentiation tend to be clustered in the genome. Further, the magnitude of variation for the sex expression ratio across populations did not substantially differ between male- and female-biased genes, except for some populations in which male-limited expressed genes showed more variation than their female counterparts. In fact, female-biased genes exhibited a larger level of interpopulation variation than male-biased genes, both when assayed in males and females. Beyond uncovering extensive adaptive potential of transcriptional variation in An. arabiensis, our findings suggest that the evolutionary rate of changes in expression levels on the X chromosome exceeds that on the autosomes, while pointing to female-biased genes as the most variable component of the An. arabiensis transcriptome.

EHMT2 Controls Transcriptional Noise and the Developmental Switch after Gastrulation in the Mouse Embryo

Embryos that carry zygotic or parental mutations in Ehmt2, the gene encoding the main euchromatic histone H3K9 methyltransferase, EHMT2, exhibit variable developmental delay. We asked the question whether the delayed embryo is different transcriptionally from the normally developing embryo when they reach the same developmental stage. We collected embryos carrying a series of genetic deficiencies in the Ehmt2 gene and performed total RNA sequencing of somite stage-matched individual embryos. We applied novel four-way comparisons to detect differences between normal versus deficient embryos, and between 12-somite and 6-somite embryos. Importantly, we also identified developmental changes in transcription that only occur during the development of the normal embryo. We found that at the 6-somite stage, gastrulation-specific genes were not precisely turned off in the Ehmt2-/- embryos, and genes involved in organ growth, connective tissue development, striated muscle development, muscle differentiation, and cartilage development were not precisely switched on in the Ehmt2-/- embryos. Zygotic EHMT2 reduced transcriptional variation of developmental switch genes and at some repeat elements at the six-somite stage embryos. Maternal EHMT2-mutant embryos also displayed great transcriptional variation consistent with their variable survival, but transcription was normal in developmentally delayed parental haploinsufficient embryos, consistent with their good prospects. Global profiling of transposable elements in the embryo revealed that specific repeat classes responded to EHMT2. DNA methylation was specifically targeted by EHMT2 to LTR repeats, mostly ERVKs. Long noncoding transcripts initiated from those misregulated driver repeats in Ehmt2-/- embryos, and extended to several hundred kilobases, encompassing a multitude of additional, similarly misexpressed passenger repeats. These findings establish EHMT2 as an important regulator of the transition between gastrulation programs and organ specification programs and of variability.

Transcriptional Variability Associated With CRISPR-Mediated Gene Replacements at the Phytophthora sojae Avr1b-1 Locus

Transcriptional plasticity enables oomycetes to rapidly adapt to environmental challenges including emerging host resistance. For example, the soybean pathogen Phytophthora sojae can overcome resistance conferred by the host resistance gene Rps1b through natural silencing of its corresponding effector gene, Avr1b-1. With the Phytophthora CRISPR/Cas9 genome editing system, it is possible to generate site-specific knock-out (KO) and knock-in (KI) mutants and to investigate the biological functions of target genes. In this study, the Avr1b-1 gene was deleted from the P. sojae genome using a homology-directed recombination strategy that replaced Avr1b-1 with a gene encoding the fluorescent protein mCherry. As expected, all selected KO transformants gained virulence on Rps1b plants, while infection of plants lacking Rps1b was not compromised. When a sgRNA-resistant version of Avr1b-1 was reintroduced into the Avr1b-1 locus of an Avr1b KO transformant, KI transformants with a well-transcribed Avr1b-1 gene were unable to infect Rps1b-containing soybeans. However, loss of expression of the incoming Avr1b-1 gene was frequently observed in KI transformants, which resulted in these transformants readily infecting Rps1b soybeans. A similar variability in the expression levels of the incoming gene was observed with AVI- or mCherry-tagged Avr1b-1 constructs. Our results suggest that Avr1b-1 may be unusually susceptible to transcriptional variation.

Using scRNA-seq to Identify Transcriptional Variation in the Malaria Parasite Ookinete Stage

The crossing of the mosquito midgut epithelium by the malaria parasite motile ookinete form represents the most extreme population bottleneck in the parasite life cycle and is a prime target for transmission blocking strategies. However, we have little understanding of the clonal variation that exists in a population of ookinetes in the vector, partially because the parasites are difficult to access and are found in low numbers. Within a vector, variation may result as a response to specific environmental cues or may exist independent of those cues as a potential bet-hedging strategy. Here we use single-cell RNA-seq to profile transcriptional variation in Plasmodium berghei ookinetes across different vector species, and between and within individual midguts. We then compare our results to low-input transcriptomes from individual Anopheles coluzzii midguts infected with the human malaria parasite Plasmodium falciparum. Although the vast majority of transcriptional changes in ookinetes are driven by development, we have identified candidate genes that may be responding to environmental cues or are clonally variant within a population. Our results illustrate the value of single-cell and low-input technologies in understanding clonal variation of parasite populations.

Cutting across structural and transcriptomic scales translates time across the lifespan in humans and chimpanzees

How the unique capacities of human cognition arose in evolution is a question of enduring interest. It is still unclear which developmental programmes are responsible for the emergence of the human brain. The inability to determine corresponding ages between humans and apes has hampered progress in detecting developmental programmes leading to the emergence of the human brain. I harness temporal variation in anatomical, behavioural and transcriptional variation to determine corresponding ages from fetal to postnatal development and ageing, between humans and chimpanzees. This multi-dimensional approach results in 137 corresponding time points across the lifespan, from embryonic day 44 to approximately 55 years of age, in humans and their equivalent ages in chimpanzees. I used these data to test whether developmental programmes, such as the timeline of prefrontal cortex (PFC) maturation, previously claimed to differ between humans and chimpanzees, do so once variation in developmental schedules is controlled for. I compared the maturation of frontal cortex projections from structural magnetic resonance (MR) scans and from temporal variation in the expression of genes used to track long-range projecting neurons (i.e. supragranular-enriched genes) in chimpanzees and humans. Contrary to what has been suggested, the timetable of PFC maturation is not unusually extended in humans. This dataset, which is the largest with which to determine corresponding ages across humans and chimpanzees, provides a rigorous approach to control for variation in developmental schedules and to identify developmental programmes responsible for unique features of the human brain.

Transcriptomic and proteomic intra-tumor heterogeneity of colorectal cancer varies depending on tumor location within the colorectum

Background Intra-tumor heterogeneity (ITH) of colorectal cancer (CRC) complicates molecular tumor classification, such as transcriptional subtyping. Differences in cellular states, biopsy cell composition, and tumor microenvironment may all lead to ITH. Here we analyze ITH at the transcriptomic and proteomic levels to ascertain whether subtype discordance between multiregional biopsies reflects relevant biological ITH or lack of classifier robustness. Further, we study the impact of tumor location on ITH. Methods Multiregional biopsies from stage II and III CRC tumors were analyzed by RNA sequencing (41 biopsies, 14 tumors) and multiplex immune protein analysis (89 biopsies, 29 tumors). CRC subtyping was performed using consensus molecular subtypes (CMS), CRC intrinsic subtypes (CRIS), and TUMOR types. ITH-scores and network maps were defined to determine the origin of heterogeneity. A validation cohort was used with one biopsy per tumor (162 tumors). Results Overall, inter-tumor transcriptional variation exceeded ITH, and subtyping calls were frequently concordant between multiregional biopsies. Still, some tumors had high transcriptional ITH and were classified discordantly. Subtyping of proximal MSS tumors were discordant for 50% of the tumors, this ITH was related to differences in the microenvironment. Subtyping of distal MSS tumors were less discordant, here the ITH was more cancer-cell related. The subtype discordancy reflected actual molecular ITH within the tumors. The relevance of the subtypes was reflected at protein level where several inflammation markers were significantly increased in immune related transcriptional subtypes, which was verified in an independent cohort (Wilcoxon rank sum test; p<0.05). Unsupervised hierarchical clustering of the protein data identified large ITH at protein level; as the multiregional biopsies clustered together for only 9 out of 29 tumors. Conclusion Our transcriptomic and proteomic analyses show that the tumor location along the colorectum influence the ITH of CRC, which again influence the concordance of subtyping.