Sensitive high-performance liquid chromatographic determination of ionic drugs in biological fluids with short-wavelength ultraviolet detection using column switching combined with ion-pair chromatography: application to basic compounds

Meningeal lymphatic dysfunction exacerbates traumatic brain injury pathogenesis

Traumatic brain injury (TBI) is a leading global cause of death and disability. Here we demonstrate in an experimental mouse model of TBI that mild forms of brain trauma cause severe deficits in meningeal lymphatic drainage that begin within hours and last out to at least one month post-injury. To investigate a mechanism underlying impaired lymphatic function in TBI, we examined how increased intracranial pressure (ICP) influences the meningeal lymphatics. We demonstrate that increased ICP can contribute to meningeal lymphatic dysfunction. Moreover, we show that pre-existing lymphatic dysfunction before TBI leads to increased neuroinflammation and negative cognitive outcomes. Finally, we report that rejuvenation of meningeal lymphatic drainage function in aged mice can ameliorate TBI-induced gliosis. These findings provide insights into both the causes and consequences of meningeal lymphatic dysfunction in TBI and suggest that therapeutics targeting the meningeal lymphatic system may offer strategies to treat TBI.

Microbial interactions in the mosquito gut determine Serratia colonization and blood-feeding propensity

How microbe–microbe interactions dictate microbial complexity in the mosquito gut is unclear. Previously we found that, Serratia, a gut symbiont that alters vector competence and is being considered for vector control, poorly colonized Aedes aegypti yet was abundant in Culex quinquefasciatus reared under identical conditions. To investigate the incompatibility between Serratia and Ae. aegypti, we characterized two distinct strains of Serratia marcescens from Cx. quinquefasciatus and examined their ability to infect Ae. aegypti. Both Serratia strains poorly infected Ae. aegypti, but when microbiome homeostasis was disrupted, the prevalence and titers of Serratia were similar to the infection in its native host. Examination of multiple genetically diverse Ae. aegypti lines found microbial interference to S. marcescens was commonplace, however, one line of Ae. aegypti was susceptible to infection. Microbiome analysis of resistant and susceptible lines indicated an inverse correlation between Enterobacteriaceae bacteria and Serratia, and experimental co-infections in a gnotobiotic system recapitulated the interference phenotype. Furthermore, we observed an effect on host behavior; Serratia exposure to Ae. aegypti disrupted their feeding behavior, and this phenotype was also reliant on interactions with their native microbiota. Our work highlights the complexity of host–microbe interactions and provides evidence that microbial interactions influence mosquito behavior.

Genome diversity of Chinese indigenous chicken and the selective signatures in Chinese gamecock chicken

Gamecock chickens are one of the earliest recorded birds in China, and have accumulated some unique morphological and behavioral signatures such as large body size, muscularity and aggressive behavior, whereby being excellent breeding materials and a good model for studying bird muscular development and behavior. In this study, we sequenced 126 chicken genomes from 19 populations, including four commercial chicken breeds that are commonly farmed in China, 13 nationwide Chinese typical indigenous chicken breeds (including two Chinese gamecock breeds), one red jungle fowl from Guangxi Province of China and three gamecock chickens from Laos. Combined with 31 published chicken genomes from three populations, a comparative genomics analysis was performed across 157 chickens. We found a severe confounding effect on potential cold adaptation exerted by introgression from commercial chickens into Chinese indigenous chickens, and argued that the genetic introgression from commercial chickens into indigenous chickens should be seriously considered for identifying selection footprint in indigenous chickens. LX gamecock chickens might have played a core role in recent breeding and conservation of other Chinese gamecock chickens. Importantly, AGMO (Alkylglycerol monooxygenase) and CPZ (Carboxypeptidase Z) might be crucial for determining the behavioral pattern of gamecock chickens, while ISPD (Isoprenoid synthase domain containing) might be essential for the muscularity of gamecock chickens. Our results can further the understanding of the evolution of Chinese gamecock chickens, especially the genetic basis of gamecock chickens revealed here was valuable for us to better understand the mechanisms underlying the behavioral pattern and the muscular development in chicken.

Suprachiasmatic VIP neurons are required for normal circadian rhythmicity and comprised of molecularly distinct subpopulations

The hypothalamic suprachiasmatic (SCN) clock contains several neurochemically defined cell groups that contribute to the genesis of circadian rhythms. Using cell-specific and genetically targeted approaches we have confirmed an indispensable role for vasoactive intestinal polypeptide-expressing SCN (SCN^VIP) neurons, including their molecular clock, in generating the mammalian locomotor activity (LMA) circadian rhythm. Optogenetic-assisted circuit mapping revealed functional, di-synaptic connectivity between SCN^VIP neurons and dorsomedial hypothalamic neurons, providing a circuit substrate by which SCN^VIP neurons may regulate LMA rhythms. In vivo photometry revealed that while SCN^VIP neurons are acutely responsive to light, their activity is otherwise behavioral state invariant. Single-nuclei RNA-sequencing revealed that SCN^VIP neurons comprise two transcriptionally distinct subtypes, including putative pacemaker and non-pacemaker populations. Altogether, our work establishes necessity of SCN^VIP neurons for the LMA circadian rhythm, elucidates organization of circadian outflow from and modulatory input to SCN^VIP cells, and demonstrates a subpopulation-level molecular heterogeneity that suggests distinct functions for specific SCN^VIP subtypes.

Cell groups in the hypothalamic suprachiasmatic clock contribute to the genesis of circadian rhythms. The authors identified two populations of vasoactive intestinal polypeptide-expressing neurons in the suprachiasmatic nucleus which regulate locomotor circadian rhythm in mice.

AmpliconReconstructor integrates NGS and optical mapping to resolve the complex structures of focal amplifications

Oncogene amplification, a major driver of cancer pathogenicity, is often mediated through focal amplification of genomic segments. Recent results implicate extrachromosomal DNA (ecDNA) as the primary driver of focal copy number amplification (fCNA) - enabling gene amplification, rapid tumor evolution, and the rewiring of regulatory circuitry. Resolving an fCNA's structure is a first step in deciphering the mechanisms of its genesis and the fCNA's subsequent biological consequences. We introduce a computational method, AmpliconReconstructor (AR), for integrating optical mapping (OM) of long DNA fragments (>150 kb) with next-generation sequencing (NGS) to resolve fCNAs at single-nucleotide resolution. AR uses an NGS-derived breakpoint graph alongside OM scaffolds to produce high-fidelity reconstructions. After validating its performance through multiple simulation strategies, AR reconstructed fCNAs in seven cancer cell lines to reveal the complex architecture of ecDNA, a breakage-fusion-bridge and other complex rearrangements. By reconstructing the rearrangement signatures associated with an fCNA's generative mechanism, AR enables a more thorough understanding of the origins of fCNAs.

Analysis of muntjac deer genome and chromatin architecture reveals rapid karyotype evolution

Closely related muntjac deer show striking karyotype differences. Here we describe chromosome-scale genome assemblies for Chinese and Indian muntjacs, Muntiacus reevesi (2n = 46) and Muntiacus muntjak vaginalis (2n = 6/7), and analyze their evolution and architecture. The genomes show extensive collinearity with each other and with other deer and cattle. We identified numerous fusion events unique to and shared by muntjacs relative to the cervid ancestor, confirming many cytogenetic observations with genome sequence. One of these M. muntjak fusions reversed an earlier fission in the cervid lineage. Comparative Hi-C analysis showed that the chromosome fusions on the M. muntjak lineage altered long-range, three-dimensional chromosome organization relative to M. reevesi in interphase nuclei including A/B compartment structure. This reshaping of multi-megabase contacts occurred without notable change in local chromatin compaction, even near fusion sites. A few genes involved in chromosome maintenance show evidence for rapid evolution, possibly associated with the dramatic changes in karyotype.

Currently favored sampling practices for tumor sequencing can produce optimal results in the clinical setting

Tumor heterogeneity is a consequence of clonal evolution, resulting in a fractal-like architecture with spatially separated main clones, sub-clones and single-cells. As sequencing an entire tumor is not feasible, we ask the question whether there is an optimal clinical sampling strategy that can handle heterogeneity and hypermutations? Here, we tested the effect of sample size, pooling strategy as well as sequencing depth using whole-exome sequencing of ovarian tumor specimens paired with normal blood samples. Our study has an emphasis on clinical application—hence we compared single biopsy, combined local biopsies and combined multi-regional biopsies. Our results show that sequencing from spatially neighboring regions show similar genetic compositions, with few private mutations. Pooling samples from multiple distinct regions of the primary tumor did not increase the overall number of identified mutations but may increase the robustness of detecting clonal mutations. Hypermutating tumors are a special case, since increasing sample size can easily dilute sub-clonal private mutations below detection thresholds. In summary, we compared the effects of sampling strategies (single biopsy, multiple local samples, pooled global sample) on mutation detection by next generation sequencing. In view of the limitations of present tools and technologies, only one sequencing run per sample combined with high coverage (100–300 ×) sequencing is affordable and practical, regardless of the number of samples taken from the same patient.

Parallel reductive genome evolution in Desulfovibrio ectosymbionts independently acquired by Trichonympha protists in the termite gut

Several Trichonympha protist species in the termite gut have independently acquired Desulfovibrio ectosymbionts in apparently different stages of symbiosis. Here, we obtained the near-complete genome sequence of Desulfovibrio phylotype ZnDsv-02, which attaches to the surface of Trichonympha collaris cells, and compared it with a previously obtained genome sequence of 'Candidatus Desulfovibrio trichonymphae' phylotype Rs-N31, which is almost completely embedded in the cytoplasm of Trichonympha agilis. Single-nucleotide polymorphism analysis indicated that although Rs-N31 is almost clonal, the ZnDsv-02 population on a single host cell is heterogeneous. Despite these differences, the genome of ZnDsv-02 has been reduced to 1.6 Mb, which is comparable to that of Rs-N31 (1.4 Mb), but unlike other known ectosymbionts of protists with a genome similar in size to their free-living relatives. Except for the presence of a lactate utilization pathway, cell-adhesion components and anti-phage defense systems in ZnDsv-02, the overall gene-loss pattern between the two genomes is very similar, including the loss of genes responsive to environmental changes. Our study suggests that genome reduction can occur in ectosymbionts, even when they can be transmitted horizontally and obtain genes via lateral transfer, and that the symbiont genome size depends heavily on their role in the symbiotic system.

Exploring the genetic basis of fatty liver development in geese

Although geese possess an adaptive physiological capacity for lipid storage, few candidate genes contributing to this ability are characterised. By comparing the genomes of individuals with extremely high and low fatty liver weights (FLW), candidate genes were identified, including ARAP2, GABRE, and IL6. Single-nucleotide polymorphisms in or near these genes were significantly (p < 0.05) associated with carcass traits (FLW) and biochemical indexes (very-low-density lipoprotein and N-terminal procollagen III), suggesting contribution to trait variation. A common variant at the 5'-end of LCORL explained ~ 18% and ~ 26% of the phenotypic variance in body weight with/without overfeeding and had significant effects on FLW (p < 0.01). ZFF36L1, ARHGEF1 and IQCJ, involved in bile acid metabolism, blood pressure, and lipid concentration modulation, were also identified. The presence of highly divergent haplotypes within these genes suggested involvement in protection against negative effects from excessive lipids in the liver or circulatory system. Based on this and transcriptomic data, we concluded that geese hepatosteatosis results from severe imbalance between lipid accumulation and secretion, comparable to human non-alcohol fatty liver disease but involving other genes. Our results provided valuable insights into the genesis of geese fatty liver and detected potential target genes for treatment of lipid-related diseases.

Meta-analysis of GWAS in canola blackleg (Leptosphaeria maculans) disease traits demonstrates increased power from imputed whole-genome sequence

Blackleg disease causes yield losses in canola (Brassica napus L.). To identify resistance genes and genomic regions, genome-wide association studies (GWAS) of 585 diverse winter and spring canola accessions were performed using imputed whole-genome sequence (WGS) and transcriptome genotype-by-sequencing (GBSt). Blackleg disease phenotypes were collected across three years in six trials. GWAS were performed in several ways and their respective power was judged by the number of significant single nucleotide polymorphisms (SNP), the false discovery rate (FDR), and the percentage of SNP that validated in additional field trials in two subsequent years. WGS GWAS with 1,234,708 million SNP detected a larger number of significant SNP, achieved a lower FDR and a higher validation rate than GBSt with 64,072 SNP. A meta-analysis combining survival and average internal infection resulted in lower FDR but also lower validation rates. The meta-analysis GWAS identified 79 genomic regions (674 SNP) conferring potential resistance to L. maculans. While several GWAS signals localised in regions of known Rlm genes, fifty-three new potential resistance regions were detected. Seventeen regions had underlying genes with putative functions related to disease defence or stress response in Arabidopsis thaliana. This study provides insight into the genetic architecture and potential molecular mechanisms underlying canola L. maculans resistance.

Improved reference genome of the arboviral vector Aedes albopictus

BACKGROUND

The Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. With limited antiviral drugs and vaccines available, vector control is the primary approach to prevent mosquito-borne diseases. A reliable and accurate DNA sequence of the Ae. albopictus genome is essential to develop new approaches that involve genetic manipulation of mosquitoes.

RESULTS

We use long-read sequencing methods and modern scaffolding techniques (PacBio, 10X, and Hi-C) to produce AalbF2, a dramatically improved assembly of the Ae. albopictus genome. AalbF2 reveals widespread viral insertions, novel microRNAs and piRNA clusters, the sex-determining locus, and new immunity genes, and enables genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we build the first physical map for this species with 75% of the assembled genome anchored to the chromosomes.

CONCLUSION

The AalbF2 genome assembly represents the most up-to-date collective knowledge of the Ae. albopictus genome. These resources represent a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

Ancient Jomon genome sequence analysis sheds light on migration patterns of early East Asian populations

Anatomically modern humans reached East Asia more than 40,000 years ago. However, key questions still remain unanswered with regard to the route(s) and the number of wave(s) in the dispersal into East Eurasia. Ancient genomes at the edge of the region may elucidate a more detailed picture of the peopling of East Eurasia. Here, we analyze the whole-genome sequence of a 2,500-year-old individual (IK002) from the main-island of Japan that is  characterized with a typical Jomon culture. The phylogenetic analyses support multiple waves of migration, with IK002 forming a basal lineage to the East and Northeast Asian genomes examined, likely representing some of the earliest-wave migrants who went north from Southeast Asia to East Asia. Furthermore, IK002 shows strong genetic affinity with the indigenous Taiwan aborigines, which may support a coastal route of the Jomon-ancestry migration. This study highlights the power of ancient genomics to provide new insights into the complex history of human migration into East Eurasia.

Takashi Gakuhari, Shigeki Nakagome et al. report the genomic analysis on a 2.5 kya individual from the ancient Jomon culture in present-day Japan. Phylogenetic analysis with comparison to other Eurasian sequences suggests early migration patterns in Asia and provides insight into the genetic affinities between peoples of the region.

Whole-genome sequencing and gene network modules predict gemcitabine/carboplatin-induced myelosuppression in non-small cell lung cancer patients

Gemcitabine/carboplatin chemotherapy commonly induces myelosuppression, including neutropenia, leukopenia, and thrombocytopenia. Predicting patients at risk of these adverse drug reactions (ADRs) and adjusting treatments accordingly is a long-term goal of personalized medicine. This study used whole-genome sequencing (WGS) of blood samples from 96 gemcitabine/carboplatin-treated non-small cell lung cancer (NSCLC) patients and gene network modules for predicting myelosuppression. Association of genetic variants in PLINK found 4594, 5019, and 5066 autosomal SNVs/INDELs with p ≤ 1 × 10⁻³ for neutropenia, leukopenia, and thrombocytopenia, respectively. Based on the SNVs/INDELs we identified the toxicity module, consisting of 215 unique overlapping genes inferred from MCODE-generated gene network modules of 350, 345, and 313 genes, respectively. These module genes showed enrichment for differentially expressed genes in rat bone marrow, human bone marrow, and human cell lines exposed to carboplatin and gemcitabine (p < 0.05). Then using 80% of the patients as training data, random LASSO reduced the number of SNVs/INDELs in the toxicity module into a feasible prediction model consisting of 62 SNVs/INDELs that accurately predict both the training and the test (remaining 20%) data with high (CTCAE 3–4) and low (CTCAE 0–1) maximal myelosuppressive toxicity completely, with the receiver-operating characteristic (ROC) area under the curve (AUC) of 100%. The present study shows how WGS, gene network modules, and random LASSO can be used to develop a feasible and tested model for predicting myelosuppressive toxicity. Although the proposed model predicts myelosuppression in this study, further evaluation in other studies is required to determine its reproducibility, usability, and clinical effect.

Circadian regulation of chemotherapy-induced peripheral neuropathic pain and the underlying transcriptomic landscape

Growing evidence demonstrates circadian rhythms of pain hypersensitivity in various chronic disorders. In chemotherapy-induced peripheral neuropathy (CIPN), agents such as paclitaxel are known to elicit chronic neuropathic pain in cancer patients and seriously compromise their quality of life. Here, we report that the mechanical threshold for allodynia in paclitaxel-treated rats exhibited a robust circadian oscillation, reaching the nadir during the daytime (inactive phase). Using Per2::LucSV circadian reporter mice expressing a PER2::LUC fusion protein, we isolated dorsal root ganglia (DRG), the primary sensory cell body for peripheral nerve injury generated hypersensitivity, and monitored ex vivo reporter bioluminescence. We observed strong circadian reporter rhythms in DRG neurons which are highly entrainable by external cues. Paclitaxel treatment significantly lengthened DRG circadian periods, with little effects on the amplitude of oscillation. We further observed the core protein BMAL1 and PER2 in DRG neurons and satellite cells. Using DRG and dorsal horn (DH; another key structure for CIPN pain response) tissues from vehicle and paclitaxel treated rats, we performed RNA-sequencing and identified diurnal expression of core clock genes as well as clock-controlled genes in both sites. Interestingly, 20.1% and 30.4% of diurnal differentially expressed genes (DEGs) overlapped with paclitaxel-induced DEGs in the DRG and the DH respectively. In contrast, paclitaxel-induced DEGs displayed only a modest overlap between daytime and nighttime (Zeitgeber Time 8 and 20). Furthermore, paclitaxel treatment induced de novo diurnal DEGs, suggesting reciprocal interaction of circadian rhythms and chemotherapy. Our study therefore demonstrates a circadian oscillation of CIPN and its underlying transcriptomic landscape.

Native American ancestry significantly contributes to neuromyelitis optica susceptibility in the admixed Mexican population

Neuromyelitis Optica (NMO) is an autoimmune disease with a higher prevalence in non-European populations. Because the Mexican population resulted from the admixture between mainly Native American and European populations, we used genome-wide microarray, HLA high-resolution typing and AQP4 gene sequencing data to analyze genetic ancestry and to seek genetic variants conferring NMO susceptibility in admixed Mexican patients. A total of 164 Mexican NMO patients and 1,208 controls were included. On average, NMO patients had a higher proportion of Native American ancestry than controls (68.1% vs 58.6%; p = 5 × 10^-6). GWAS identified a HLA region associated with NMO, led by rs9272219 (OR = 2.48, P = 8 × 10^-10). Class II HLA alleles HLA-DQB1*03:01, -DRB1*08:02, -DRB1*16:02, -DRB1*14:06 and -DQB1*04:02 showed the most significant associations with NMO risk. Local ancestry estimates suggest that all the NMO-associated alleles within the HLA region are of Native American origin. No novel or missense variants in the AQP4 gene were found in Mexican patients with NMO or multiple sclerosis. To our knowledge, this is the first study supporting the notion that Native American ancestry significantly contributes to NMO susceptibility in an admixed population, and is consistent with differences in NMO epidemiology in Mexico and Latin America.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.

Patterns of DNA variation between the autosomes, the X chromosome and the Y chromosome in Bos taurus genome

The new ARS-UCD1.2 assembly of the bovine genome has considerable improvements over the previous assembly and thus more accurate identification of patterns of genetic variation can be achieved with it. We explored differences in genetic variation between autosomes, the X chromosome, and the Y chromosome. In particular, variant densities, annotations, lengths (only for InDels), nucleotide divergence, and Tajima’s D statistics between chromosomes were considered. Whole-genome DNA sequences of 217 individuals representing different cattle breeds were examined. The analysis included the alignment to the new reference genome and variant identification. 23,655,295 SNPs and 3,758,781 InDels were detected. In contrast to autosomes, both sex chromosomes had negative values of Tajima’s D and lower nucleotide divergence. That implies a correlation between nucleotide diversity and recombination rate, which is obviously reduced for sex chromosomes. Moreover, the accumulation of nonsynonymous mutations on the Y chromosome could be associated with loss of recombination. Also, the relatively lower effective population size for sex chromosomes leads to a lower expected density of variants.

Repeatable ecological dynamics govern the response of experimental communities to antibiotic pulse perturbation

In an era of pervasive anthropogenic ecological disturbances, there is a pressing need to understand the factors that constitute community response and resilience. A detailed understanding of disturbance response needs to go beyond associations and incorporate features of disturbances, species traits, rapid evolution and dispersal. Multispecies microbial communities that experience antibiotic perturbation represent a key system with important medical dimensions. However, previous microbiome studies on this theme have relied on high-throughput sequencing data from uncultured species without the ability to explicitly account for the role of species traits and immigration. Here, we serially passage a 34-species defined bacterial community through different levels of pulse antibiotic disturbance, manipulating the presence or absence of species immigration. To understand the ecological community response measured using amplicon sequencing, we combine initial trait data measured for each species separately and metagenome sequencing data revealing adaptive mutations during the experiment. We found that the ecological community response was highly repeatable within the experimental treatments, which could be attributed in part to key species traits (antibiotic susceptibility and growth rate). Increasing antibiotic levels were also coupled with an increasing probability of species extinction, making species immigration critical for community resilience. Moreover, we detected signals of antibiotic-resistance evolution occurring within species at the same time scale, leaving evolutionary changes in communities despite recovery at the species compositional level. Together, these observations reveal a disturbance response that presents as classic species sorting, but is nevertheless accompanied by rapid within-species evolution.

Wapl repression by Pax5 promotes V gene recombination by Igh loop extrusion

Nuclear processes, such as V(D)J recombination, are orchestrated by the three-dimensional organization of chromosomes at multiple levels, including compartments1 and topologically associated domains (TADs)2,3 consisting of chromatin loops4. TADs are formed by chromatin-loop extrusion5-7, which depends on the loop-extrusion function of the ring-shaped cohesin complex8-12. Conversely, the cohesin-release factor Wapl13,14 restricts loop extension10,15. The generation of a diverse antibody repertoire, providing humoral immunity to pathogens, requires the participation of all V genes in V(D)J recombination16, which depends on contraction of the 2.8-Mb-long immunoglobulin heavy chain (Igh) locus by Pax517,18. However, how Pax5 controls Igh contraction in pro-B cells remains unknown. Here we demonstrate that locus contraction is caused by loop extrusion across the entire Igh locus. Notably, the expression of Wapl is repressed by Pax5 specifically in pro-B and pre-B cells, facilitating extended loop extrusion by increasing the residence time of cohesin on chromatin. Pax5 mediates the transcriptional repression of Wapl through a single Pax5-binding site by recruiting the polycomb repressive complex 2 to induce bivalent chromatin at the Wapl promoter. Reduced Wapl expression causes global alterations in the chromosome architecture, indicating that the potential to recombine all V genes entails structural changes of the entire genome in pro-B cells.

Primary mammary angiosarcomas harbor frequent mutations in KDR and PIK3CA and show evidence of distinct pathogenesis

Angiosarcoma (AS) is the most frequent primary sarcoma of the breast but nevertheless remains uncommon, accounting for <0.05% of breast malignancies. Secondary mammary AS arise following radiation therapy for breast cancer, in contrast to primary AS which occur sporadically. Essentially all show aggressive clinical behavior independent of histologic grade and most are treated by mastectomy. MYC amplification is frequently identified in radiation-induced AS but only rarely in primary mammary AS (PMAS). As a heterogeneous group, AS from various anatomic sites have been shown to harbor recurrent alterations in TP53, MAP kinase pathway genes, and genes involved in angiogenic signaling including KDR (VEGFR2) and PTPRB. In part due to its rarity, the pathogenesis of PMAS has not been fully characterized. In this study, we examined the clinical, pathologic, and genomic features of ten cases of PMAS, including one patient with bilateral disease. Recurrent genomic alterations were identified in KDR (70%), PIK3CA/PIK3R1 (70%), and PTPRB (30%), each at higher frequencies than reported in AS across all sites. Six tumors harbored a KDR p.T771R hotspot mutation, and all seven KDR-mutant cases showed evidence suggestive of biallelism (four with loss of heterozygosity and three with two aberrations). Of the seven tumors with PI3K alterations, six harbored pathogenic mutations other than in the canonical PIK3CA residues which are most frequent in breast cancer. Three AS were hypermutated (≥10 mutations/megabase (Mb)); hypermutation was seen concurrent with KDR or PIK3CA mutations. The patient with bilateral disease demonstrated shared alterations, indicative of contralateral metastasis. No MYC or TP53 aberrations were detected in this series. Immunohistochemistry for VEGFR2 was unable to discriminate between KDR-mutant tumors and benign vascular lesions of the breast. These findings highlight the underrecognized frequency of KDR and PIK3CA mutation in PMAS, and a significant subset with hypermutation, suggesting a pathogenesis distinct from other AS.

Genomic profiling of primary and recurrent adult granulosa cell tumors of the ovary

Adult-type granulosa cell tumor (aGCT) is a rare malignant ovarian sex cord-stromal tumor, harboring recurrent FOXL2 c.C402G/p.C134W hotspot mutations in 97% of cases. These tumors are considered to have a favorable prognosis, however aGCTs have a tendency for local spread and late recurrences, which are associated with poor survival rates. We sought to determine the genetic alterations associated with aGCT disease progression. We subjected primary non-recurrent aGCTs (n = 7), primary aGCTs that subsequently recurred (n = 9) and their matched recurrences (n = 9), and aGCT recurrences without matched primary tumors (n = 10) to targeted massively parallel sequencing of ≥410 cancer-related genes. In addition, three primary non-recurrent aGCTs and nine aGCT recurrences were subjected to FOXL2 and TERT promoter Sanger sequencing analysis. All aGCTs harbored the FOXL2 C134W hotspot mutation. TERT promoter mutations were found to be significantly more frequent in recurrent (18/28, 64%) than primary aGCTs (5/19, 26%, p = 0.017). In addition, mutations affecting TP53, MED12, and TET2 were restricted to aGCT recurrences. Pathway annotation of altered genes demonstrated that aGCT recurrences displayed an enrichment for genetic alterations affecting cell cycle pathway-related genes. Analysis of paired primary and recurrent aGCTs revealed that TERT promoter mutations were either present in both primary tumors and matched recurrences or were restricted to the recurrence and absent in the respective primary aGCT. Clonal composition analysis of these paired samples further revealed that aGCTs display intra-tumor genetic heterogeneity and harbor multiple clones at diagnosis and relapse. We observed that in a subset of cases, recurrences acquired additional genetic alterations not present in primary aGCTs, including TERT, MED12, and TP53 mutations and CDKN2A/B homozygous deletions. Albeit harboring relatively simple genomes, our data provide evidence to suggest that aGCTs are genetically heterogeneous tumors and that TERT promoter mutations and/or genetic alterations affecting other cell cycle-related genes may be associated with disease progression and recurrences.

Genomic epidemiology of Mycobacterium tuberculosis in Santa Catarina, Southern Brazil

Mycobacterium tuberculosis (M.tb), the pathogen responsible for tuberculosis (TB) poses as the major cause of death among infectious diseases. The knowledge about the molecular diversity of M.tb enables the implementation of more effective surveillance and control measures and, nowadays, Whole Genome Sequencing (WGS) holds the potential to produce high-resolution epidemiological data in a high-throughput manner. Florianópolis, the state capital of Santa Catarina (SC) in south Brazil, shows a high TB incidence (46.0/100,000). Here we carried out a WGS-based evaluation of the M.tb strain diversity, drug-resistance and ongoing transmission in the capital metropolitan region. Resistance to isoniazid, rifampicin, streptomycin was identified respectively in 4.0% (n = 6), 2.0% (n = 3) and 1.3% (n = 2) of the 151 studied strains by WGS. Besides, resistance to pyrazinamide and ethambutol was detected in 0.7% (n = 1) and reistance to ethionamide and fluoroquinolone (FQ) in 1.3% (n = 2), while a single (0.7%) multidrug-resistant (MDR) strain was identified. SNP-based typing classified all isolates into M.tb Lineage 4, with high proportion of sublineages LAM (60.3%), T (16.4%) and Haarlem (7.9%). The average core-genome distance between isolates was 420.3 SNPs, with 43.7% of all isolates grouped across 22 genomic clusters thereby showing the presence of important ongoing TB transmission events. Most clusters were geographically distributed across the study setting which highlights the need for an urgent interruption of these large transmission chains. The data conveyed by this study shows the presence of important and uncontrolled TB transmission in the metropolitan area and provides precise data to support TB control measures in this region.

Occupancy maps of 208 chromatin-associated proteins in one human cell type

Transcription factors are DNA-binding proteins that have key roles in gene regulation^1,2. Genome-wide occupancy maps of transcriptional regulators are important for understanding gene regulation and its effects on diverse biological processes^3–6. However, only a minority of the more than 1,600 transcription factors encoded in the human genome has been assayed. Here we present, as part of the ENCODE (Encyclopedia of DNA Elements) project, data and analyses from chromatin immunoprecipitation followed by high-throughput sequencing (ChIP–seq) experiments using the human HepG2 cell line for 208 chromatin-associated proteins (CAPs). These comprise 171 transcription factors and 37 transcriptional cofactors and chromatin regulator proteins, and represent nearly one-quarter of CAPs expressed in HepG2 cells. The binding profiles of these CAPs form major groups associated predominantly with promoters or enhancers, or with both. We confirm and expand the current catalogue of DNA sequence motifs for transcription factors, and describe motifs that correspond to other transcription factors that are co-enriched with the primary ChIP target. For example, FOX family motifs are enriched in ChIP–seq peaks of 37 other CAPs. We show that motif content and occupancy patterns can distinguish between promoters and enhancers. This catalogue reveals high-occupancy target regions at which many CAPs associate, although each contains motifs for only a minority of the numerous associated transcription factors. These analyses provide a more complete overview of the gene regulatory networks that define this cell type, and demonstrate the usefulness of the large-scale production efforts of the ENCODE Consortium.

ChIP–seq and CETCh–seq data are used to analyse binding maps for 208 transcription factors and other chromatin-associated proteins in a single human cell type, providing a comprehensive catalogue of the transcription factor landscape and gene regulatory networks in these cells.

Mechanisms of stretch-mediated skin expansion at single-cell resolution

The ability of the skin to grow in response to stretching has been exploited in reconstructive surgery¹. Although the response of epidermal cells to stretching has been studied in vitro^2,3, it remains unclear how mechanical forces affect their behaviour in vivo. Here we develop a mouse model in which the consequences of stretching on skin epidermis can be studied at single-cell resolution. Using a multidisciplinary approach that combines clonal analysis with quantitative modelling and single-cell RNA sequencing, we show that stretching induces skin expansion by creating a transient bias in the renewal activity of epidermal stem cells, while a second subpopulation of basal progenitors remains committed to differentiation. Transcriptional and chromatin profiling identifies how cell states and gene-regulatory networks are modulated by stretching. Using pharmacological inhibitors and mouse mutants, we define the step-by-step mechanisms that control stretch-mediated tissue expansion at single-cell resolution in vivo.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

A chromosome-level genome assembly of the wild rice Oryza rufipogon facilitates tracing the origins of Asian cultivated rice

Oryza rufipogon Griff. is a wild progenitor of the Asian cultivated rice Oryza sativa. To better understand the genomic diversity of the wild rice, high-quality reference genomes of O. rufipogon populations are needed, which also facilitate utilization of the wild genetic resources in rice breeding. In this study, we generated a chromosome-level genome assembly of O. rufipogon using a combination of short-read sequencing, single-molecule sequencing, BioNano and Hi-C platforms. The genome sequence (399.8 Mb) was assembled into 46 scaffolds on the 12 chromosomes, with contig N50 and scaffold N50 of 13.2 Mb and 20.3 Mb, respectively. The genome contains 36,520 protein-coding genes, and 49.37% of the genome consists of repetitive elements. The genome has strong synteny with those of the O. sativa subspecies indica and japonica, but containing some large structural variations. Evolutionary analysis unveiled the polyphyletic origins of O. sativa, in which the japonica and indica genome formations involved different divergent O. rufipogon (including O. nivara) lineages, accompanied by introgression of genomic regions between japonica and indica. This high-quality reference genome provides insight on the genome evolution of the wild rice and the origins of the O. sativa subspecies, and valuable information for basic research and rice breeding.

Differences in local population history at the finest level: the case of the Estonian population

Several recent studies detected fine-scale genetic structure in human populations. Hence, groups conventionally treated as single populations harbour significant variation in terms of allele frequencies and patterns of haplotype sharing. It has been shown that these findings should be considered when performing studies of genetic associations and natural selection, especially when dealing with polygenic phenotypes. However, there is little understanding of the practical effects of such genetic structure on demography reconstructions and selection scans when focusing on recent population history. Here we tested the impact of population structure on such inferences using high-coverage (~30×) genome sequences of 2305 Estonians. We show that different regions of Estonia differ in both effective population size dynamics and signatures of natural selection. By analyzing identity-by-descent segments we also reveal that some Estonian regions exhibit evidence of a bottleneck 10-15 generations ago reflecting sequential episodes of wars, plague and famine, although this signal is virtually undetected when treating Estonia as a single population. Besides that, we provide a framework for relating effective population size estimated from genetic data to actual census size and validate it on the Estonian population. This approach may be widely used both to cross-check estimates based on historical sources as well as to get insight into times and/or regions with no other information available. Our results suggest that the history of human populations within the last few millennia can be highly region specific and cannot be properly studied without taking local genetic structure into account.

Murine models of IDH-wild-type glioblastoma exhibit spatial segregation of tumor initiation and manifestation during evolution

Recent characterization of spatiotemporal genomic architecture of IDH-wild-type multifocal glioblastomas (M-GBMs) suggests a clinically unobserved common-ancestor (CA) with a less aggressive phenotype, generating highly genetically divergent malignant gliomas/GBMs in distant brain regions. Using serial MRI/3D-reconstruction, whole-genome sequencing and spectral karyotyping-based single-cell phylogenetic tree building, we show two distinct types of tumor evolution in p53-mutant driven mouse models. Malignant gliomas/GBMs grow as a single mass (Type 1) and multifocal masses (Type 2), respectively, despite both exhibiting loss of Pten/chromosome 19 (chr19) and PI3K/Akt activation with sub-tetraploid/4N genomes. Analysis of early biopsied and multi-segment tumor tissues reveals no evidence of less proliferative diploid/2N lesions in Type 1 tumors. Strikingly, CA-derived relatively quiescent tumor precursors with ancestral diploid/2N genomes and normal Pten/chr19 are observed in the subventricular zone (SVZ), but are distantly segregated from multi focal Type 2 tumors. Importantly, PI3K/Akt inhibition by Rictor/mTORC2 deletion blocks distant dispersal, restricting glioma growth in the SVZ.

BRCA1-associated structural variations are a consequence of polymerase theta-mediated end-joining

Failure to preserve the integrity of the genome is a hallmark of cancer. Recent studies have revealed that loss of the capacity to repair DNA breaks via homologous recombination (HR) results in a mutational profile termed BRCAness. The enzymatic activity that repairs HR substrates in BRCA-deficient conditions to produce this profile is currently unknown. We here show that the mutational landscape of BRCA1 deficiency in C. elegans closely resembles that of BRCA1-deficient tumours. We identify polymerase theta-mediated end-joining (TMEJ) to be responsible: knocking out polq-1 suppresses the accumulation of deletions and tandem duplications in brc-1 and brd-1 animals. We find no additional back-up repair in HR and TMEJ compromised animals; non-homologous end-joining does not affect BRCAness. The notion that TMEJ acts as an alternative to HR, promoting the genome alteration of HR-deficient cells, supports the idea that polymerase theta is a promising therapeutic target for HR-deficient tumours.

A comparison of proteomic, genomic, and osteological methods of archaeological sex estimation

Sex estimation of skeletons is fundamental to many archaeological studies. Currently, three approaches are available to estimate sex-osteology, genomics, or proteomics, but little is known about the relative reliability of these methods in applied settings. We present matching osteological, shotgun-genomic, and proteomic data to estimate the sex of 55 individuals, each with an independent radiocarbon date between 2,440 and 100 cal BP, from two ancestral Ohlone sites in Central California. Sex estimation was possible in 100% of this burial sample using proteomics, in 91% using genomics, and in 51% using osteology. Agreement between the methods was high, however conflicts did occur. Genomic sex estimates were 100% consistent with proteomic and osteological estimates when DNA reads were above 100,000 total sequences. However, more than half the samples had DNA read numbers below this threshold, producing high rates of conflict with osteological and proteomic data where nine out of twenty conditional DNA sex estimates conflicted with proteomics. While the DNA signal decreased by an order of magnitude in the older burial samples, there was no decrease in proteomic signal. We conclude that proteomics provides an important complement to osteological and shotgun-genomic sex estimation.

Cross-species oncogenic signatures of breast cancer in canine mammary tumors

Genomic and precision medicine research has afforded notable advances in human cancer treatment, yet applicability to other species remains uncertain. Through whole-exome and transcriptome analyses of 191 spontaneous canine mammary tumors (CMTs) that exhibit the archetypal features of human breast cancers, we found a striking resemblance of genomic characteristics including frequent PIK3CA mutations (43.1%), aberrations of the PI3K-Akt pathway (61.7%), and key genes involved in cancer initiation and progression. We also identified three gene expression-based CMT subtypes, one of which segregated with basal-like human breast cancer subtypes with activated epithelial-to-mesenchymal transition, low claudin expression, and unfavorable disease prognosis. A relative lack of ERBB2 amplification and Her2-enrichment subtype in CMT denoted species-specific molecular mechanisms. Taken together, our results elucidate cross-species oncogenic signatures for a better understanding of universal and context-dependent mechanisms in breast cancer development and provide a basis for precision diagnostics and therapeutics for domestic dogs.

New group of transmembrane proteins associated with desiccation tolerance in the anhydrobiotic midge Polypedilum vanderplanki

Larvae of the sleeping chironomid Polypedilum vanderplanki are known for their extraordinary ability to survive complete desiccation in an ametabolic state called "anhydrobiosis". The unique feature of P. vanderplanki genome is the presence of expanded gene clusters associated with anhydrobiosis. While several such clusters represent orthologues of known genes, there is a distinct set of genes unique for P. vanderplanki. These include Lea-Island-Located (LIL) genes with no known orthologues except two of LEA genes of P. vanderplanki, PvLea1 and PvLea3. However, PvLIL proteins lack typical features of LEA such as the state of intrinsic disorder, hydrophilicity and characteristic LEA_4 motif. They possess four to five transmembrane domains each and we confirmed membrane targeting for three PvLILs. Conserved amino acids in PvLIL are located in transmembrane domains or nearby. PvLEA1 and PvLEA3 proteins are chimeras combining LEA-like parts and transmembrane domains, shared with PvLIL proteins. We have found that PvLil genes are highly upregulated during anhydrobiosis induction both in larvae of P. vanderplanki and P. vanderplanki-derived cultured cell line, Pv11. Thus, PvLil are a new intriguing group of genes that are likely to be associated with anhydrobiosis due to their common origin with some LEA genes and their induction during anhydrobiosis.

A first genetic linage map construction and QTL mapping for growth traits in Larimichthys polyactis

Larimichthys polyactis is a commercially important marine fish species in Eastern Asia, yet very few genetic resources exist. In particular, genetic linkage maps are critical tools for genetic breeding. In this study, we generated a high resolution linkage map from a family of 110 individuals and their parents by resequencing the individuals. 3,802 effective SNPs were mapped to 24 linkage groups (LGs). The map spanned 2,567.39 cm, with an average marker interval of 0.66 cm. We used the map to conduct QTL analysis for growth traits, and found 31 markers were significantly associated with growth-related traits. Specifically, three SNPs were identified for total length, nineteen SNPs for body length, and nine SNPs for body weight. The identified SNPs could explain 15.2-22.6% of the phenotypic variation. SNPs associated with growth traits were distributed on LG6 and LG11, and candidate genes included, kif26b, bat1, gna1, gbgt1, and amfr, which may regulate growth. The linkage map and mapped QTLs would be useful for improving the quality of L. polyactis via marker-assisted selection.

Lactate production by Staphylococcus aureus biofilm inhibits HDAC11 to reprogramme the host immune response during persistent infection

Staphylococcus aureus (S. aureus) is a leading cause of biofilm-associated prosthetic joint infection (PJI), resulting in significant disability and prolonged treatment. It is known that host leukocyte IL-10 production is required for S. aureus biofilm persistence in PJI. A S. aureus bursa aurealis Tn library consisting of 1,952 non-essential genes was screened for mutants that failed to induce IL-10 in myeloid-derived suppressor cells (MDSCs), which identified a critical role for bacterial lactic acid biosynthesis. We generated a S. aureus ddh/ldh1/ldh2 triple Tn mutant that cannot produce D- or L-lactate. Co-culture of MDSCs or macrophages with ddh/ldh1/ldh2 mutant biofilm produced substantially less IL-10 compared with wild type S. aureus, which was also observed in a mouse model of PJI and led to reduced biofilm burden. Using MDSCs recovered from the mouse PJI model and in vitro leukocyte-biofilm co-cultures we show that bacterial-derived lactate inhibits histone deacetylase 11 (HDAC11), causing unchecked HDAC6 activity and increased histone 3 acetylation at the Il-10 promoter, resulting in enhanced Il-10 transcription in MDSCs and macrophages. Finally, we show that synovial fluid of patients with PJI contains elevated amounts of D-lactate and IL-10 compared with control subjects, and bacterial lactate increases IL-10 production by human monocyte-derived macrophages.

Biofilms are bacterial communities that are difficult to treat because of their tolerance to antibiotics and ability to evade immune-mediated clearance. Prosthetic joint infection (PJI), a devastating complication of arthroplasty, is characterized by biofilm formation. The current study has discovered a central role for lactic acid biosynthesis in S. aureus biofilm formation during PJI. Mechanistically, bacterial-derived lactate inhibits histone deacetylase 11 (HDAC11) activity, which causes extensive epigenetic changes at the promoters of numerous host genes, including the key anti-inflammatory cytokine Il-10. Indeed, IL-10 production by myeloid-derived suppressor cells (MDSCs) and macrophages is critical for biofilm persistence during PJI. HDAC11 inhibition by S. aureus lactate results in unchecked HDAC6 activity, a positive regulator of IL-10, thereby increasing IL-10 production by MDSCs and macrophages in vitro and in vivo. Similarly, S. aureus lactate promotes IL-10 production in human monocyte-derived macrophages following biofilm exposure. This study highlights how bacterial metabolism can influence the host immune response to promote infection persistence.

Biological rhythms in the deep-sea hydrothermal mussel Bathymodiolus azoricus

Biological rhythms are a fundamental property of life. The deep ocean covers 66% of our planet surface and is one of the largest biomes. The deep sea has long been considered as an arrhythmic environment because sunlight is totally absent below 1,000 m depth. In the present study, we have sequenced the temporal transcriptomes of a deep-sea species, the ecosystem-structuring vent mussel Bathymodiolus azoricus. We reveal that tidal cycles predominate in the transcriptome and physiology of mussels fixed directly at hydrothermal vents at 1,688 m depth at the Mid-Atlantic Ridge, whereas daily cycles prevail in mussels sampled after laboratory acclimation. We identify B. azoricus canonical circadian clock genes, and show that oscillations observed in deep-sea mussels could be either a direct response to environmental stimulus, or be driven endogenously by one or more biological clocks. This work generates in situ insights into temporal organisation in a deep-sea organism.

Little is known about gene expression of organisms in the deep sea, partially owing to constraints on sampling these organisms in situ. Here the authors circumvent this problem, fixing tissue of a deep-sea mussel at 1,688 m in depth, and later analyzing transcriptomes to reveal gene expression patterns showing tidal oscillations.

The chromatin remodeling enzyme Chd4 regulates genome architecture in the mouse brain

The development and function of the brain require tight control of gene expression. Genome architecture is thought to play a critical regulatory role in gene expression, but the mechanisms governing genome architecture in the brain in vivo remain poorly understood. Here, we report that conditional knockout of the chromatin remodeling enzyme Chd4 in granule neurons of the mouse cerebellum increases accessibility of gene regulatory sites genome-wide in vivo. Conditional knockout of Chd4 promotes recruitment of the architectural protein complex cohesin preferentially to gene enhancers in granule neurons in vivo. Importantly, in vivo profiling of genome architecture reveals that conditional knockout of Chd4 strengthens interactions among developmentally repressed contact domains as well as genomic loops in a manner that tightly correlates with increased accessibility, enhancer activity, and cohesin occupancy at these sites. Collectively, our findings define a role for chromatin remodeling in the control of genome architecture organization in the mammalian brain.

Massive haplotypes underlie ecotypic differentiation in sunflowers

Species often include multiple ecotypes that are adapted to different environments¹. However, it is unclear how ecotypes arise and how their distinctive combinations of adaptive alleles are maintained despite hybridization with non-adapted populations^2-4. Here, by resequencing 1,506 wild sunflowers from 3 species (Helianthus annuus, Helianthus petiolaris and Helianthus argophyllus), we identify 37 large (1-100 Mbp in size), non-recombining haplotype blocks that are associated with numerous ecologically relevant traits, as well as soil and climate characteristics. Limited recombination in these haplotype blocks keeps adaptive alleles together, and these regions differentiate sunflower ecotypes. For example, haplotype blocks control a 77-day difference in flowering between ecotypes of the silverleaf sunflower H. argophyllus (probably through deletion of a homologue of FLOWERING LOCUS T (FT)), and are associated with seed size, flowering time and soil fertility in dune-adapted sunflowers. These haplotypes are highly divergent, frequently associated with structural variants and often appear to represent introgressions from other-possibly now-extinct-congeners. These results highlight a pervasive role of structural variation in ecotypic adaptation.

Pervasive lesion segregation shapes cancer genome evolution

Cancers arise through the acquisition of oncogenic mutations and grow by clonal expansion1,2. Here we reveal that most mutagenic DNA lesions are not resolved into a mutated DNA base pair within a single cell cycle. Instead, DNA lesions segregate, unrepaired, into daughter cells for multiple cell generations, resulting in the chromosome-scale phasing of subsequent mutations. We characterize this process in mutagen-induced mouse liver tumours and show that DNA replication across persisting lesions can produce multiple alternative alleles in successive cell divisions, thereby generating both multiallelic and combinatorial genetic diversity. The phasing of lesions enables accurate measurement of strand-biased repair processes, quantification of oncogenic selection and fine mapping of sister-chromatid-exchange events. Finally, we demonstrate that lesion segregation is a unifying property of exogenous mutagens, including UV light and chemotherapy agents in human cells and tumours, which has profound implications for the evolution and adaptation of cancer genomes.

Histone H3.3 phosphorylation amplifies stimulation-induced transcription

Complex organisms can rapidly induce select genes in response to diverse environmental cues. This regulation occurs in the context of large genomes condensed by histone proteins into chromatin. The sensing of pathogens by macrophages engages conserved signalling pathways and transcription factors to coordinate the induction of inflammatory genes1-3. Enriched integration of histone H3.3, the ancestral histone H3 variant, is a general feature of dynamically regulated chromatin and transcription4-7. However, how chromatin is regulated at induced genes, and what features of H3.3 might enable rapid and high-level transcription, are unknown. The amino terminus of H3.3 contains a unique serine residue (Ser31) that is absent in 'canonical' H3.1 and H3.2. Here we show that this residue, H3.3S31, is phosphorylated (H3.3S31ph) in a stimulation-dependent manner along rapidly induced genes in mouse macrophages. This selective mark of stimulation-responsive genes directly engages the histone methyltransferase SETD2, a component of the active transcription machinery, and 'ejects' the elongation corepressor ZMYND118,9. We propose that features of H3.3 at stimulation-induced genes, including H3.3S31ph, provide preferential access to the transcription apparatus. Our results indicate dedicated mechanisms that enable rapid transcription involving the histone variant H3.3, its phosphorylation, and both the recruitment and the ejection of chromatin regulators.

Genomic signatures of domestication in Old World camels

Domestication begins with the selection of animals showing less fear of humans. In most domesticates, selection signals for tameness have been superimposed by intensive breeding for economical or other desirable traits. Old World camels, conversely, have maintained high genetic variation and lack secondary bottlenecks associated with breed development. By re-sequencing multiple genomes from dromedaries, Bactrian camels, and their endangered wild relatives, here we show that positive selection for candidate genes underlying traits collectively referred to as 'domestication syndrome' is consistent with neural crest deficiencies and altered thyroid hormone-based signaling. Comparing our results with other domestic species, we postulate that the core set of domestication genes is considerably smaller than the pan-domestication set - and overlapping genes are likely a result of chance and redundancy. These results, along with the extensive genomic resources provided, are an important contribution to understanding the evolutionary history of camels and the genomic features of their domestication.

A clinically and genomically annotated nerve sheath tumor biospecimen repository

Nerve sheath tumors occur as a heterogeneous group of neoplasms in patients with neurofibromatosis type 1 (NF1). The malignant form represents the most common cause of death in people with NF1, and even when benign, these tumors can result in significant disfigurement, neurologic dysfunction, and a range of profound symptoms. Lack of human tissue across the peripheral nerve tumors common in NF1 has been a major limitation in the development of new therapies. To address this unmet need, we have created an annotated collection of patient tumor samples, patient-derived cell lines, and patient-derived xenografts, and carried out high-throughput genomic and transcriptomic characterization to serve as a resource for further biologic and preclinical therapeutic studies. In this work, we release genomic and transcriptomic datasets comprised of 55 tumor samples derived from 23 individuals, complete with clinical annotation. All data are publicly available through the NF Data Portal and at http://synapse.org/jhubiobank.

Measurement(s)	gene expression • gene_variant
Technology Type(s)	RNA sequencing • exome sequencing • DNA sequencing
Factor Type(s)	tumor type
Sample Characteristic - Organism	Homo sapiens • Homo sapiens/Mus musculus xenograft

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12037599

Dissecting intratumour heterogeneity of nodal B-cell lymphomas at the transcriptional, genetic and drug-response levels

Tumour heterogeneity encompasses both the malignant cells and their microenvironment. While heterogeneity between individual patients is known to affect the efficacy of cancer therapy, most personalized treatment approaches do not account for intratumour heterogeneity. We addressed this issue by studying the heterogeneity of nodal B-cell lymphomas by single-cell RNA-sequencing and transcriptome-informed flow cytometry. We identified transcriptionally distinct malignant subpopulations and compared their drug-response and genomic profiles. Malignant subpopulations from the same patient responded strikingly differently to anti-cancer drugs ex vivo, which recapitulated subpopulation-specific drug sensitivity during in vivo treatment. Infiltrating T cells represented the majority of non-malignant cells, whose gene-expression signatures were similar across all donors, whereas the frequencies of T-cell subsets varied significantly between the donors. Our data provide insights into the heterogeneity of nodal B-cell lymphomas and highlight the relevance of intratumour heterogeneity for personalized cancer therapy.

Transcriptional regulation of strigolactone signalling in Arabidopsis

Plant hormones known as strigolactones control plant development and interactions between host plants and symbiotic fungi or parasitic weeds^1-4. In Arabidopsis thaliana and rice, the proteins DWARF14 (D14), MORE AXILLARY GROWTH 2 (MAX2), SUPPRESSOR OF MAX2-LIKE 6, 7 and 8 (SMXL6, SMXL7 and SMXL8) and their orthologues form a complex upon strigolactone perception and play a central part in strigolactone signalling^5-10. However, whether and how strigolactones activate downstream transcription remains largely unknown. Here we use a synthetic strigolactone to identify 401 strigolactone-responsive genes in Arabidopsis, and show that these plant hormones regulate shoot branching, leaf shape and anthocyanin accumulation mainly through transcriptional activation of the BRANCHED 1, TCP DOMAIN PROTEIN 1 and PRODUCTION OF ANTHOCYANIN PIGMENT 1 genes. We find that SMXL6 targets 729 genes in the Arabidopsis genome and represses the transcription of SMXL6, SMXL7 and SMXL8 by binding directly to their promoters, showing that SMXL6 serves as an autoregulated transcription factor to maintain the homeostasis of strigolactone signalling. These findings reveal an unanticipated mechanism through which a transcriptional repressor of hormone signalling can directly recognize DNA and regulate transcription in higher plants.

The genetic basis of sex determination in grapes

It remains a major challenge to identify the genes and mutations that lead to plant sexual differentiation. Here, we study the structure and evolution of the sex-determining region (SDR) in Vitis species. We report an improved, chromosome-scale Cabernet Sauvignon genome sequence and the phased assembly of nine wild and cultivated grape genomes. By resolving twenty Vitis SDR haplotypes, we compare male, female, and hermaphrodite haplotype structures and identify sex-linked regions. Coupled with gene expression data, we identify a candidate male-sterility mutation in the VviINP1 gene and potential female-sterility function associated with the transcription factor VviYABBY3. Our data suggest that dioecy has been lost during domestication through a rare recombination event between male and female haplotypes. This work significantly advances the understanding of the genetic basis of sex determination in Vitis and provides the information necessary to rapidly identify sex types in grape breeding programs.

Grapevine is one of a few ancestrally dioecious crops that are reverted to hermaphroditism during domestication. Here, the authors identify candidate genes related to male- and female-sterility in grapes and describe the genetic process that led to hermaphroditism during domestication.

ASCL1-regulated DARPP-32 and t-DARPP stimulate small cell lung cancer growth and neuroendocrine tumour cell proliferation

Background

Small cell lung cancer (SCLC) is the most aggressive form of lung cancer, and new molecular insights are necessary for prognostic and therapeutic advances.

Methods

Dopamine and cAMP-regulated phosphoprotein, Mr 32000 (DARPP-32) and its N-terminally truncated splice variant, t-DARPP, were stably overexpressed or ablated in human DMS-53 and H1048 SCLC cells. Functional assays and immunoblotting were used to assess how DARPP-32 isoforms regulate SCLC cell growth, proliferation, and apoptosis. DARPP-32-modulated SCLC cells were orthotopically injected into the lungs of SCID mice to evaluate how DARPP-32 and t-DARPP regulate neuroendocrine tumour growth. Immunostaining for DARPP-32 proteins was performed in SCLC patient-derived specimens. Bioinformatics analysis and subsequent transcription assays were used to determine the mechanistic basis of DARPP-32-regulated SCLC growth.

Results

We demonstrate in mice that DARPP-32 and t-DARPP promote SCLC growth through increased Akt/Erk-mediated proliferation and anti-apoptotic signalling. DARPP-32 isoforms are overexpressed in SCLC patient-derived tumour tissue, but undetectable in physiologically normal lung. Achaete-scute homologue 1 (ASCL1) transcriptionally activates DARPP-32 isoforms in human SCLC cells.

Conclusions

We reveal new regulatory mechanisms of SCLC oncogenesis that suggest DARPP-32 isoforms may represent a negative prognostic indicator for SCLC and serve as a potential target for the development of new therapies.

The characterization of Mediator 12 and 13 as conditional positive gene regulators in Arabidopsis

Mediator 12 (MED12) and MED13 are components of the Mediator multi-protein complex, that facilitates the initial steps of gene transcription. Here, in an Arabidopsis mutant screen, we identify MED12 and MED13 as positive gene regulators, both of which contribute broadly to morc1 de-repressed gene expression. Both MED12 and MED13 are preferentially required for the expression of genes depleted in active chromatin marks, a chromatin signature shared with morc1 re-activated loci. We further discover that MED12 tends to interact with genes that are responsive to environmental stimuli, including light and radiation. We demonstrate that light-induced transient gene expression depends on MED12, and is accompanied by a concomitant increase in MED12 enrichment during induction. In contrast, the steady-state expression level of these genes show little dependence on MED12, suggesting that MED12 is primarily required to aid the expression of genes in transition from less-active to more active states.

Mediator is a multiprotein complex required to activate gene transcription by RNAPII. Here, the authors report that MED12 and MED13 are conditional positive regulators that facilitate the expression of genes depleted in active chromatin marks and the induction of gene expression in response to environmental stimuli in Arabidopsis.

Analysis of cell-associated DENV RNA by oligo(dT) primed 5' capture scRNAseq

Dengue is one of the most widespread vector-borne viral diseases in the world. However, the size, heterogeneity, and temporal dynamics of the cell-associated viral reservoir during acute dengue virus (DENV) infection remains unclear. In this study, we analyzed cells infected in vitro with DENV and PBMC from an individual experiencing a natural DENV infection utilizing 5’ capture single cell RNA sequencing (scRNAseq). Both positive- and negative-sense DENV RNA was detected in reactions containing either an oligo(dT) primer alone, or in reactions supplemented with a DENV-specific primer. The addition of a DENV-specific primer did not increase the total amount of DENV RNA captured or the fraction of cells identified as containing DENV RNA. However, inclusion of a DENV-specific cDNA primer did increase the viral genome coverage immediately 5’ to the primer binding site. Furthermore, while the majority of intracellular DENV sequence captured in this analysis mapped to the 5’ end of the viral genome, distinct patterns of enhanced coverage within the DENV polyprotein coding region were observed. The 5’ capture scRNAseq analysis of PBMC not only recapitulated previously published reports by detecting virally infected memory and naïve B cells, but also identified cell-associated genomic variants not observed in contemporaneous serum samples. These results demonstrate that oligo(dT) primed 5’ capture scRNAseq can detect DENV RNA and quantify virus-infected cells in physiologically relevant conditions, and provides insight into viral sequence variability within infected cells.

The transcriptional repressor Blimp1/PRDM1 regulates the maternal decidual response in mice

The transcriptional repressor Blimp1 controls cell fate decisions in the developing embryo and adult tissues. Here we describe Blimp1 expression and functional requirements within maternal uterine tissues during pregnancy. Expression is robustly up-regulated at early post-implantation stages in the primary decidual zone (PDZ) surrounding the embryo. Conditional inactivation results in defective formation of the PDZ barrier and abnormal trophectoderm invasion. RNA-Seq analysis demonstrates down-regulated expression of genes involved in cell adhesion and markers of decidualisation. In contrast, genes controlling immune responses including IFNγ are up-regulated. ChIP-Seq experiments identify candidate targets unique to the decidua as well as those shared across diverse cell types including a highly conserved peak at the Csf-1 gene promoter. Interestingly Blimp1 inactivation results in up-regulated Csf1 expression and macrophage recruitment into maternal decidual tissues. These results identify Blimp1 as a critical regulator of tissue remodelling and maternal tolerance during early stages of pregnancy.

Simultaneous quantification of protein-DNA interactions and transcriptomes in single cells with scDam&T-seq

Protein-DNA interactions are essential for establishing cell type-specific chromatin architecture and gene expression. We recently developed scDam&T-seq, a multi-omics method that can simultaneously quantify protein-DNA interactions and the transcriptome in single cells. The method effectively combines two existing methods: DNA adenine methyltransferase identification (DamID) and CEL-Seq2. DamID works through the tethering of a protein of interest (POI) to the Escherichia coli DNA adenine methyltransferase (Dam). Upon expression of this fusion protein, DNA in proximity to the POI is methylated by Dam and can be selectively digested and amplified. CEL-Seq2, in contrast, makes use of poly-dT primers to reverse transcribe mRNA, followed by linear amplification through in vitro transcription. scDam&T-seq is the first technique capable of providing a combined readout of protein-DNA contact and transcription from single-cell samples. Once suitable cell lines have been established, the protocol can be completed in 5 d, with a throughput of hundreds to thousands of cells. The processing of raw sequencing data takes an additional 1-2 d. Our method can be used to understand the transcriptional changes a cell undergoes upon the DNA binding of a POI. It can be performed in any laboratory with access to FACS, robotic and high-throughput-sequencing facilities.

Colibactin DNA-damage signature indicates mutational impact in colorectal cancer

The mucosal epithelium is a common target of damage by chronic bacterial infections and the accompanying toxins, and most cancers originate from this tissue. We investigated whether colibactin, a potent genotoxin¹ associated with certain strains of Escherichia coli², creates a specific DNA-damage signature in infected human colorectal cells. Notably, the genomic contexts of colibactin-induced DNA double-strand breaks were enriched for an AT-rich hexameric sequence motif, associated with distinct DNA-shape characteristics. A survey of somatic mutations at colibactin target sites of several thousand cancer genomes revealed notable enrichment of this motif in colorectal cancers. Moreover, the exact double-strand-break loci corresponded with mutational hot spots in cancer genomes, reminiscent of a trinucleotide signature previously identified in healthy colorectal epithelial cells³. The present study provides evidence for the etiological role of colibactin in human cancer.