SAMBLASTER: fast duplicate marking and structural variant read extraction

The Lipid Metabolism as Target and Modulator of BOLD-100 Anticancer Activity: Crosstalk with Histone Acetylation

The leading first-in-class ruthenium-complex BOLD-100 currently undergoes clinical phase-II anticancer evaluation. Recently, BOLD-100 is identified as anti-Warburg compound. The present study shows that also deregulated lipid metabolism parameters characterize acquired BOLD-100-resistant colon and pancreatic carcinoma cells. Acute BOLD-100 treatment reduces lipid droplet contents of BOLD-100-sensitive but not -resistant cells. Despite enhanced glycolysis fueling lipid accumulation, BOLD-100-resistant cells reveal diminished lactate secretion based on monocarboxylate transporter 1 (MCT1) loss mediated by a frame-shift mutation in the MCT1 chaperone basigin. Glycolysis and lipid catabolism converge in the production of protein/histone acetylation substrate acetyl-coenzymeA (CoA). Mass spectrometric and nuclear magnetic resonance analyses uncover spontaneous cell-free BOLD-100-CoA adduct formation suggesting acetyl-CoA depletion as mechanism bridging BOLD-100-induced lipid metabolism alterations and histone acetylation-mediated gene expression deregulation. Indeed, BOLD-100 treatment decreases histone acetylation selectively in sensitive cells. Pharmacological targeting confirms histone de-acetylation as central mode-of-action of BOLD-100 and metabolic programs stabilizing histone acetylation as relevant Achilles' heel of acquired BOLD-100-resistant cell and xenograft models. Accordingly, histone gene expression changes also predict intrinsic BOLD-100 responsiveness. Summarizing, BOLD-100 is identified as epigenetically active substance acting via targeting several onco-metabolic pathways. Identification of the lipid metabolism as driver of acquired BOLD-100 resistance opens novel strategies to tackle therapy failure.

Evolutionary consequences of nascent multicellular life cycles

A key step in the evolutionary transition to multicellularity is the origin of multicellular groups as biological individuals capable of adaptation. Comparative work, supported by theory, suggests clonal development should facilitate this transition, although this hypothesis has never been tested in a single model system. We evolved 20 replicate populations of otherwise isogenic clonally reproducing 'snowflake' yeast (Δace2/∆ace2) and aggregative 'floc' yeast (GAL1p::FLO1 /GAL1p::FLO1) with daily selection for rapid growth in liquid media, which favors faster cell division, followed by selection for rapid sedimentation, which favors larger multicellular groups. While both genotypes adapted to this regime, growing faster and having higher survival during the group-selection phase, there was a stark difference in evolutionary dynamics. Aggregative floc yeast obtained nearly all their increased fitness from faster growth, not improved group survival; indicating that selection acted primarily at the level of cells. In contrast, clonal snowflake yeast mainly benefited from higher group-dependent fitness, indicating a shift in the level of Darwinian individuality from cells to groups. Through genome sequencing and mathematical modeling, we show that the genetic bottlenecks in a clonal life cycle also drive much higher rates of genetic drift-a result with complex implications for this evolutionary transition. Our results highlight the central role that early multicellular life cycles play in the process of multicellular adaptation.

Microhomology-mediated circular DNA formation from oligonucleosomal fragments during spermatogenesis

The landscape of extrachromosomal circular DNA (eccDNA) during mammalian spermatogenesis, as well as the biogenesis mechanism, remains to be explored. Here, we revealed widespread eccDNA formation in human sperms and mouse spermatogenesis. We noted that germline eccDNAs are derived from oligonucleosomal DNA fragmentation in cells likely undergoing cell death, providing a potential new way for quality assessment of human sperms. Interestingly, small-sized eccDNAs are associated with euchromatin, while large-sized ones are preferentially generated from heterochromatin. By comparing sperm eccDNAs with meiotic recombination hotspots and structural variations, we found that they are barely associated with de novo germline deletions. We further developed a bioinformatics pipeline to achieve nucleotide-resolution eccDNA detection even with the presence of microhomologous sequences that interfere with precise breakpoint identification. Empowered by our method, we provided strong evidence to show that microhomology-mediated end joining is the major eccDNA biogenesis mechanism. Together, our results shed light on eccDNA biogenesis mechanism in mammalian germline cells.

Clinical trio genome sequencing facilitates the interpretation of variants in cancer predisposition genes in paediatric tumour patients

The prevalence of pathogenic and likely pathogenic (P/LP) variants in genes associated with cancer predisposition syndromes (CPS) is estimated to be 8-18% for paediatric cancer patients. In more than half of the carriers, the family history is unsuspicious for CPS. Therefore, broad genetic testing could identify germline predisposition in additional children with cancer resulting in important implications for themselves and their families. We thus evaluated clinical trio genome sequencing (TGS) in a cohort of 72 paediatric patients with solid cancers other than retinoblastoma or CNS-tumours. The most prevalent cancer types were sarcoma (n = 26), neuroblastoma (n = 15), and nephroblastoma (n = 10). Overall, P/LP variants in CPS genes were identified in 18.1% of patients (13/72) and P/LP variants in autosomal-dominant CPS genes in 9.7% (7/72). Genetic evaluation would have been recommended for the majority of patients with P/LP variants according to the Jongmans criteria. Four patients (5.6%, 4/72) carried P/LP variants in autosomal-dominant genes known to be associated with their tumour type. With the immediate information on variant inheritance, TGS facilitated the identification of a de novo P/LP in NF1, a gonadosomatic mosaic in WT1 and two pathogenic variants in one patient (DICER1 and PALB2). TGS allows a more detailed characterization of structural variants with base-pair resolution of breakpoints which can be relevant for the interpretation of copy number variants. Altogether, TGS allows comprehensive identification of children with a CPS and supports the individualised clinical management of index patients and high-risk relatives.

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

Myeloid-specific KDM6B inhibition sensitizes glioblastoma to PD1 blockade

Glioblastoma (GBM) tumors are enriched in immune-suppressive myeloid cells and are refractory to immune checkpoint therapy (ICT). Targeting epigenetic pathways to reprogram the functional phenotype of immune-suppressive myeloid cells to overcome resistance to ICT remains unexplored. Single-cell and spatial transcriptomic analyses of human GBM tumors demonstrated high expression of an epigenetic enzyme-histone 3 lysine 27 demethylase (KDM6B)-in intratumoral immune-suppressive myeloid cell subsets. Importantly, myeloid cell-specific Kdm6b deletion enhanced proinflammatory pathways and improved survival in GBM tumor-bearing mice. Mechanistic studies showed that the absence of Kdm6b enhances antigen presentation, interferon response and phagocytosis in myeloid cells by inhibition of mediators of immune suppression including Mafb, Socs3 and Sirpa. Further, pharmacological inhibition of KDM6B mirrored the functional phenotype of Kdm6b-deleted myeloid cells and enhanced anti-PD1 efficacy. This study thus identified KDM6B as an epigenetic regulator of the functional phenotype of myeloid cell subsets and a potential therapeutic target for enhanced response to ICT.

Concerning the eXclusion in human genomics: the choice of sex chromosome representation in the human genome drastically affects the number of identified variants

Over the past 30 years, a community of scientists has pieced together every base pair of the human reference genome from telomere to telomere. Interestingly, most human genomics studies omit more than 5% of the genome from their analyses. Under "normal" circumstances, omitting any chromosome(s) from an analysis of the human genome would be a cause for concern, with the exception being sex chromosomes. Sex chromosomes in eutherians share an evolutionary origin as an ancestral pair of autosomes. In humans, they share 3 regions of high-sequence identity (∼98-100%), which, along with the unique transmission patterns of the sex chromosomes, introduce technical artifacts in genomic analyses. However, the human X chromosome bears numerous important genes, including more "immune response" genes than any other chromosome, which makes its exclusion irresponsible when sex differences across human diseases are widespread. To better characterize the possible effect of the inclusion/exclusion of the X chromosome on variants called, we conducted a pilot study on the Terra cloud platform to replicate a subset of standard genomic practices using both the CHM13 reference genome and the sex chromosome complement-aware reference genome. We compared the quality of variant calling, expression quantification, and allele-specific expression using these 2 reference genome versions across 50 human samples from the Genotype-Tissue Expression consortium annotated as females. We found that after correction, the whole X chromosome (100%) can generate reliable variant calls, allowing for the inclusion of the whole genome in human genomics analyses as a departure from the status quo of omitting the sex chromosomes from empirical and clinical genomics studies.

Integrating chromatin conformation information in a self-supervised learning model improves metagenome binning

Metagenome binning is a key step, downstream of metagenome assembly, to group scaffolds by their genome of origin. Although accurate binning has been achieved on datasets containing multiple samples from the same community, the completeness of binning is often low in datasets with a small number of samples due to a lack of robust species co-abundance information. In this study, we exploited the chromatin conformation information obtained from Hi-C sequencing and developed a new reference-independent algorithm, Metagenome Binning with Abundance and Tetra-nucleotide frequencies-Long Range (metaBAT-LR), to improve the binning completeness of these datasets. This self-supervised algorithm builds a model from a set of high-quality genome bins to predict scaffold pairs that are likely to be derived from the same genome. Then, it applies these predictions to merge incomplete genome bins, as well as recruit unbinned scaffolds. We validated metaBAT-LR's ability to bin-merge and recruit scaffolds on both synthetic and real-world metagenome datasets of varying complexity. Benchmarking against similar software tools suggests that metaBAT-LR uncovers unique bins that were missed by all other methods. MetaBAT-LR is open-source and is available at https://bitbucket.org/project-metabat/metabat-lr.

A scaffolded and annotated reference genome of giant kelp (Macrocystis pyrifera)

Macrocystis pyrifera (giant kelp), is a brown macroalga of great ecological importance as a primary producer and structure-forming foundational species that provides habitat for hundreds of species. It has many commercial uses (e.g. source of alginate, fertilizer, cosmetics, feedstock). One of the limitations to exploiting giant kelp's economic potential and assisting in giant kelp conservation efforts is a lack of genomic tools like a high quality, contiguous reference genome with accurate gene annotations. Reference genomes attempt to capture the complete genomic sequence of an individual or species, and importantly provide a universal structure for comparison across a multitude of genetic experiments, both within and between species. We assembled the giant kelp genome of a haploid female gametophyte de novo using PacBio reads, then ordered contigs into chromosome level scaffolds using Hi-C. We found the giant kelp genome to be 537 MB, with a total of 35 scaffolds and 188 contigs. The assembly N50 is 13,669,674 with GC content of 50.37%. We assessed the genome completeness using BUSCO, and found giant kelp contained 94% of the BUSCO genes from the stramenopile clade. Annotation of the giant kelp genome revealed 25,919 genes. Additionally, we present genetic variation data based on 48 diploid giant kelp sporophytes from three different Southern California populations that confirms the population structure found in other studies of these populations. This work resulted in a high-quality giant kelp genome that greatly increases the genetic knowledge of this ecologically and economically vital species.

BiocMAP: a Bioconductor-friendly, GPU-accelerated pipeline for bisulfite-sequencing data

BACKGROUND

Bisulfite sequencing is a powerful tool for profiling genomic methylation, an epigenetic modification critical in the understanding of cancer, psychiatric disorders, and many other conditions. Raw data generated by whole genome bisulfite sequencing (WGBS) requires several computational steps before it is ready for statistical analysis, and particular care is required to process data in a timely and memory-efficient manner. Alignment to a reference genome is one of the most computationally demanding steps in a WGBS workflow, taking several hours or even days with commonly used WGBS-specific alignment software. This naturally motivates the creation of computational workflows that can utilize GPU-based alignment software to greatly speed up the bottleneck step. In addition, WGBS produces raw data that is large and often unwieldy; a lack of memory-efficient representation of data by existing pipelines renders WGBS impractical or impossible to many researchers.

RESULTS

We present BiocMAP, a Bioconductor-friendly methylation analysis pipeline consisting of two modules, to address the above concerns. The first module performs computationally-intensive read alignment using Arioc, a GPU-accelerated short-read aligner. Since GPUs are not always available on the same computing environments where traditional CPU-based analyses are convenient, the second module may be run in a GPU-free environment. This module extracts and merges DNA methylation proportions-the fractions of methylated cytosines across all cells in a sample at a given genomic site. Bioconductor-based output objects in R utilize an on-disk data representation to drastically reduce required main memory and make WGBS projects computationally feasible to more researchers.

CONCLUSIONS

BiocMAP is implemented using Nextflow and available at http://research.libd.org/BiocMAP/ . To enable reproducible analysis across a variety of typical computing environments, BiocMAP can be containerized with Docker or Singularity, and executed locally or with the SLURM or SGE scheduling engines. By providing Bioconductor objects, BiocMAP's output can be integrated with powerful analytical open source software for analyzing methylation data.

Methods behind neoantigen prediction for personalized anticancer vaccines

Next to conventional cancer therapies, immunotherapies such as immune checkpoint inhibitors have broadened the cancer treatment landscape over the past decades. Recent advances in next generation sequencing and bioinformatics technologies have made it possible to identify a patient's own immunogenic neoantigens. These cancer neoantigens serve as important targets for personalized immunotherapy which has the benefit of being more active and effective in targeting cancer cells. This paper is a step-by-step guide discussing the different analyses and challenges encountered during in-silico neoantigen prediction. The protocol describes all the tools and steps required for the identification of immunogenic neoantigens.

Ribosome biogenesis disruption mediated chromatin structure changes revealed by SRAtac, a customizable end to end analysis pipeline for ATAC-seq

The nucleolus is a large nuclear body that serves as the primary site for ribosome biogenesis. Recent studies have suggested that it also plays an important role in organizing chromatin architecture. However, to establish a causal relationship between nucleolar ribosome assembly and chromatin architecture, genetic tools are required to disrupt nucleolar ribosome biogenesis. In this study, we used ATAC-seq to investigate changes in chromatin accessibility upon specific depletion of two ribosome biogenesis components, RPOA-2 and GRWD-1, in the model organism Caenorhabditis elegans. To facilitate the analysis of ATAC-seq data, we introduced two tools: SRAlign, an extensible NGS data processing workflow, and SRAtac, a customizable end-to-end ATAC-seq analysis pipeline. Our results revealed highly comparable changes in chromatin accessibility following both RPOA-2 and GRWD-1 perturbations. However, we observed a weak correlation between changes in chromatin accessibility and gene expression. While our findings corroborate the idea of a feedback mechanism between ribosomal RNA synthesis, nucleolar ribosome large subunit biogenesis, and chromatin structure during the L1 stage of C. elegans development, they also prompt questions regarding the functional impact of these alterations on gene expression.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Unravelling human hematopoietic progenitor cell diversity through association with intrinsic regulatory factors

Hematopoietic stem and progenitor cell (HSPC) transplantation is an essential therapy for hematological conditions, but finer definitions of human HSPC subsets with associated function could enable better tuning of grafts and more routine, lower-risk application. To deeply phenotype HSPCs, following a screen of 328 antigens, we quantified 41 surface proteins and functional regulators on millions of CD34+ and CD34- cells, spanning four primary human hematopoietic tissues: bone marrow, mobilized peripheral blood, cord blood, and fetal liver. We propose more granular definitions of HSPC subsets and provide new, detailed differentiation trajectories of erythroid and myeloid lineages. These aspects of our revised human hematopoietic model were validated with corresponding epigenetic analysis and in vitro clonal differentiation assays. Overall, we demonstrate the utility of using molecular regulators as surrogates for cellular identity and functional potential, providing a framework for description, prospective isolation, and cross-tissue comparison of HSPCs in humans.

DIS3 Variants are Associated With Primary Ovarian Insufficiency: Importance of Transcription/Translation in Oogenesis

CONTEXT

A genetic etiology accounts for the majority of unexplained primary ovarian insufficiency (POI).

OBJECTIVE

We hypothesized a genetic cause of POI for a sister pair with primary amenorrhea.

DESIGN

The study was an observational study. Subjects were recruited at an academic institution.

SUBJECTS

Subjects were sisters with primary amenorrhea caused by POI and their parents. Additional subjects included women with POI analyzed previously (n = 291). Controls were recruited for health in old age or were from the 1000 Genomes Project (total n = 233).

INTERVENTION

We performed whole exome sequencing, and data were analyzed using the Pedigree Variant Annotation, Analysis and Search Tool, which identifies genes harboring pathogenic variants in families. We performed functional studies in a Drosophila melanogaster model.

MAIN OUTCOME

Genes with rare pathogenic variants were identified.

RESULTS

The sisters carried compound heterozygous variants in DIS3. The sisters did not carry additional rare variants that were absent in publicly available datasets. DIS3 knockdown in the ovary of D. melanogaster resulted in lack of oocyte production and severe infertility.

CONCLUSIONS

Compound heterozygous variants in highly conserved amino acids in DIS3 and failure of oocyte production in a functional model suggest that mutations in DIS3 cause POI. DIS3 is a 3' to 5' exoribonuclease that is the catalytic subunit of the exosome involved in RNA degradation and metabolism in the nucleus. The findings provide further evidence that mutations in genes important for transcription and translation are associated with POI.

Analysis of SARS-CoV-2 Population Genetics from Samples Associated with Huanan Market and Early Cases Identifies Substitutions Associated with Future Variants of Concern

SARS-CoV-2 began spreading through human-to-human transmission first within China and then worldwide, with increasing sequence diversity associated with time and the further spread of the virus. The spillover events in the Huanan market were associated with two lineages of SARS-CoV-2 (lineages A and B). Infecting virus populations and those in infected individuals consist of a dominant genomic sequence and minor genomic variants; these latter populations can indicate sites on the genome that may be subject to mutational changes-either neutral or advantageous sites and those that act as a reservoir for future dominant variants-when placed under selection pressure. The earliest deposited sequences with human infections associated with the Huanan market shared very close homology with each other and were all lineage B. However, there were minor genomic variants present in each sample that encompassed synonymous and non-synonymous changes. Fusion sequences characteristic of defective RNA were identified that could potentially link transmission chains between individuals. Although all the individuals appeared to have lineage B as the dominant sequence, nucleotides associated with lineage A could be found at very low frequencies. Several substitutions (but not deletions) associated with much later variants of concern (VoCs) were already present as minor genomic variants. This suggests that low-frequency substitutions at the start of a pandemic could be a reservoir of future dominant variants and/or provide information on potential sites within the genome associated with future plasticity.

Whole genome sequencing of Canadian Saccharomyces cerevisiae strains isolated from spontaneous wine fermentations reveals a new Pacific West Coast Wine clade

Vineyards in wine regions around the world are reservoirs of yeast with oenological potential. Saccharomyces cerevisiae ferments grape sugars to ethanol and generates flavor and aroma compounds in wine. Wineries place a high-value on identifying yeast native to their region to develop a region-specific wine program. Commercial wine strains are genetically very similar due to a population bottleneck and in-breeding compared to the diversity of S. cerevisiae from the wild and other industrial processes. We have isolated and microsatellite-typed hundreds of S. cerevisiae strains from spontaneous fermentations of grapes from the Okanagan Valley wine region in British Columbia, Canada. We chose 75 S. cerevisiae strains, based on our microsatellite clustering data, for whole genome sequencing using Illumina paired-end reads. Phylogenetic analysis shows that British Columbian S. cerevisiae strains cluster into 4 clades: Wine/European, Transpacific Oak, Beer 1/Mixed Origin, and a new clade that we have designated as Pacific West Coast Wine. The Pacific West Coast Wine clade has high nucleotide diversity and shares genomic characteristics with wild North American oak strains but also has gene flow from Wine/European and Ecuadorian clades. We analyzed gene copy number variations to find evidence of domestication and found that strains in the Wine/European and Pacific West Coast Wine clades have gene copy number variation reflective of adaptations to the wine-making environment. The "wine circle/Region B", a cluster of 5 genes acquired by horizontal gene transfer into the genome of commercial wine strains is also present in the majority of the British Columbian strains in the Wine/European clade but in a minority of the Pacific West Coast Wine clade strains. Previous studies have shown that S. cerevisiae strains isolated from Mediterranean Oak trees may be the living ancestors of European wine yeast strains. This study is the first to isolate S. cerevisiae strains with genetic similarity to nonvineyard North American Oak strains from spontaneous wine fermentations.

The master growth regulator DELLA binding to histone H2A is essential for DELLA-mediated global transcription regulation

The DELLA genes, also known as 'Green Revolution' genes, encode conserved master growth regulators that control plant development in response to internal and environmental cues. Functioning as nuclear-localized transcription regulators, DELLAs modulate expression of target genes via direct protein-protein interaction of their carboxy-terminal GRAS domain with hundreds of transcription factors (TFs) and epigenetic regulators. However, the molecular mechanism of DELLA-mediated transcription reprogramming remains unclear. Here by characterizing new missense alleles of an Arabidopsis DELLA, repressor of ga1-3 (RGA), and co-immunoprecipitation assays, we show that RGA binds histone H2A via the PFYRE subdomain within its GRAS domain to form a TF-RGA-H2A complex at the target chromatin. Chromatin immunoprecipitation followed by sequencing analysis further shows that this activity is essential for RGA association with its target chromatin globally. Our results indicate that, although DELLAs are recruited to target promoters by binding to TFs via the LHR1 subdomain, DELLA-H2A interaction via the PFYRE subdomain is necessary to stabilize the TF-DELLA-H2A complex at the target chromatin. This study provides insights into the two distinct key modular functions in DELLA for its genome-wide transcription regulation in plants.

Phages are unrecognized players in the ecology of the oral pathogen Porphyromonas gingivalis

BACKGROUND

Porphyromonas gingivalis (hereafter "Pg") is an oral pathogen that has been hypothesized to act as a keystone driver of inflammation and periodontal disease. Although Pg is most readily recovered from individuals with actively progressing periodontal disease, healthy individuals and those with stable non-progressing disease are also colonized by Pg. Insights into the factors shaping the striking strain-level variation in Pg, and its variable associations with disease, are needed to achieve a more mechanistic understanding of periodontal disease and its progression. One of the key forces often shaping strain-level diversity in microbial communities is infection of bacteria by their viral (phage) predators and symbionts. Surprisingly, although Pg has been the subject of study for over 40 years, essentially nothing is known of its phages, and the prevailing paradigm is that phages are not important in the ecology of Pg.

RESULTS

Here we systematically addressed the question of whether Pg are infected by phages-and we found that they are. We found that prophages are common in Pg, they are genomically diverse, and they encode genes that have the potential to alter Pg physiology and interactions. We found that phages represent unrecognized targets of the prevalent CRISPR-Cas defense systems in Pg, and that Pg strains encode numerous additional mechanistically diverse candidate anti-phage defense systems. We also found that phages and candidate anti-phage defense system elements together are major contributors to strain-level diversity and the species pangenome of this oral pathogen. Finally, we demonstrate that prophages harbored by a model Pg strain are active in culture, producing extracellular viral particles in broth cultures.

CONCLUSION

This work definitively establishes that phages are a major unrecognized force shaping the ecology and intra-species strain-level diversity of the well-studied oral pathogen Pg. The foundational phage sequence datasets and model systems that we establish here add to the rich context of all that is already known about Pg, and point to numerous avenues of future inquiry that promise to shed new light on fundamental features of phage impacts on human health and disease broadly. Video Abstract.

OpenPBTA: The Open Pediatric Brain Tumor Atlas

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community.

Integrative profiling of extrachromosomal circular DNA in placenta and maternal plasma provides insights into the biology of fetal growth restriction and reveals potential biomarkers

Introduction: Fetal growth restriction (FGR) is a placenta-mediated pregnancy complication that predisposes fetuses to perinatal complications. Maternal plasma cell-free DNA harbors DNA originating from placental trophoblasts, which is promising for the prenatal diagnosis and prediction of pregnancy complications. Extrachromosomal circular DNA (eccDNA) is emerging as an ideal biomarker and target for several diseases. Methods: We utilized eccDNA sequencing and bioinformatic pipeline to investigate the characteristics and associations of eccDNA in placenta and maternal plasma, the role of placental eccDNA in the pathogenesis of FGR, and potential plasma eccDNA biomarkers of FGR. Results: Using our bioinformatics pipelines, we identified multi-chromosomal-fragment and single-fragment eccDNA in placenta, but almost exclusively single-fragment eccDNA in maternal plasma. Relative to that in plasma, eccDNA in placenta was larger and substantially more abundant in exons, untranslated regions, promoters, repetitive elements [short interspersed nuclear elements (SINEs)/Alu, SINEs/mammalian-wide interspersed repeats, long terminal repeats/endogenous retrovirus-like elements, and single recognition particle RNA], and transcription factor binding motifs. Placental multi-chromosomal-fragment eccDNA was enriched in confident enhancer regions predicted to pertain to genes in apoptosis, energy, cell growth, and autophagy pathways. Placental eccDNA-associated genes whose abundance differed between the FGR and control groups were associated with immunity-related gene ontology (GO) terms. The combined analysis of plasma and placental eccDNA-associated genes in the FGR and control groups led to the identification of potential biomarkers that were assigned to the GO terms of the epigenetic regulation of gene expression and nutrient-related processes, respectively. Conclusion: Together, our results highlight links between placenta functions and multi-chromosomal-fragment and single-fragment eccDNA. The integrative analysis of placental and plasma eccDNA confirmed the potential of these molecules as disease-specific biomarkers of FGR.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

Impacts of Sex Ratio Meiotic Drive on Genome Structure and Function in a Stalk-Eyed Fly

Stalk-eyed flies in the genus Teleopsis carry selfish genetic elements that induce sex ratio (SR) meiotic drive and impact the fitness of male and female carriers. Here, we assemble and describe a chromosome-level genome assembly of the stalk-eyed fly, Teleopsis dalmanni, to elucidate patterns of divergence associated with SR. The genome contains tens of thousands of transposable element (TE) insertions and hundreds of transcriptionally and insertionally active TE families. By resequencing pools of SR and ST males using short and long reads, we find widespread differentiation and divergence between XSR and XST associated with multiple nested inversions involving most of the SR haplotype. Examination of genomic coverage and gene expression data revealed seven X-linked genes with elevated expression and coverage in SR males. The most extreme and likely drive candidate involves an XSR-specific expansion of an array of partial copies of JASPer, a gene necessary for maintenance of euchromatin and associated with regulation of TE expression. In addition, we find evidence for rapid protein evolution between XSR and XST for testis expressed and novel genes, that is, either recent duplicates or lacking a Dipteran ortholog, including an X-linked duplicate of maelstrom, which is also involved in TE silencing. Overall, the evidence suggests that this ancient XSR polymorphism has had a variety of impacts on repetitive DNA and its regulation in this species.

An integrated cellular and molecular model of gastric neuroendocrine cancer evolution highlights therapeutic targets

Gastric neuroendocrine carcinomas (G-NEC) are aggressive malignancies with poorly understood biology and a lack of disease models. Here, we use genome sequencing to characterize the genomic landscapes of human G-NEC and its histologic variants. We identify global and subtype-specific alterations and expose hitherto unappreciated gains of MYC family members in a large part of cases. Genetic engineering and lineage tracing in mice delineate a model of G-NEC evolution, which defines MYC as a critical driver and positions the cancer cell of origin to the neuroendocrine compartment. MYC-driven tumors have pronounced metastatic competence and display defined signaling addictions, as revealed by large-scale genetic and pharmacologic screening of cell lines and organoid resources. We create global maps of G-NEC dependencies, highlight critical vulnerabilities, and validate therapeutic targets, including candidates for clinical drug repurposing. Our study gives comprehensive insights into G-NEC biology.

Extreme genome diversity and cryptic speciation in a harmful algal-bloom-forming eukaryote

Harmful algal blooms of the toxic haptophyte Prymnesium parvum are a recurrent problem in many inland and estuarine waters around the world. Strains of P. parvum vary in the toxins they produce and in other physiological traits associated with harmful algal blooms, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for 15 phylogenetically and geographically diverse strains of P. parvum, including Hi-C guided, near-chromosome-level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 to 845 Mbp. Strains included haploids, diploids, and polyploids, but not all differences in DNA content were due to variation in genome copy number. Haploid genome size between strains of different chemotypes differed by as much as 243 Mbp. Syntenic and phylogenetic analyses indicate that UTEX 2797, a common laboratory strain from Texas, is a hybrid that retains two phylogenetically distinct haplotypes. Investigation of gene families variably present across the strains identified several functional categories associated with metabolic and genome size variation in P. parvum, including genes for the biosynthesis of toxic metabolites and proliferation of transposable elements. Together, our results indicate that P. parvum comprises multiple cryptic species. These genomes provide a robust phylogenetic and genomic framework for investigations into the eco-physiological consequences of the intra- and inter-specific genetic variation present in P. parvum and demonstrate the need for similar resources for other harmful algal-bloom-forming morphospecies.

Neurite outgrowth deficits caused by rare PLXNB1 mutation in pediatric bipolar disorder

Pediatric bipolar disorder (PBD) is a severe mood dysregulation condition that affects 0.5-1% of children and teens in the United States. It is associated with recurrent episodes of mania and depression and an increased risk of suicidality. However, the genetics and neuropathology of PBD are largely unknown. Here, we used a combinatorial family-based approach to characterize cellular, molecular, genetic, and network-level deficits associated with PBD. We recruited a PBD patient and three unaffected family members from a family with a history of psychiatric illnesses. Using resting-state functional magnetic resonance imaging (rs-fMRI), we detected altered resting-state functional connectivity in the patient as compared to an unaffected sibling. Using transcriptomic profiling of patient and control induced pluripotent stem cell (iPSC)-derived telencephalic organoids, we found aberrant signaling in the molecular pathways related to neurite outgrowth. We corroborated the presence of neurite outgrowth deficits in patient iPSC-derived cortical neurons and identified a rare homozygous loss-of-function PLXNB1 variant (c.1360C>C; p.Ser454Arg) responsible for the deficits in the patient. Expression of wild-type PLXNB1, but not the variant, rescued neurite outgrowth in patient neurons, and expression of the variant caused the neurite outgrowth deficits in cortical neurons from PlxnB1 knockout mice. These results indicate that dysregulated PLXNB1 signaling may contribute to an increased risk of PBD and other mood dysregulation-related disorders by disrupting neurite outgrowth and functional brain connectivity. Overall, this study established and validated a novel family-based combinatorial approach for studying cellular and molecular deficits in psychiatric disorders and identified dysfunctional PLXNB1 signaling and neurite outgrowth as potential risk factors for PBD.

Comparative genomics of rainbow trout (Oncorhynchus mykiss): Is the genetic architecture of migratory behavior conserved among populations?

Rainbow trout (Oncorhynchus mykiss) are a partially migratory species wherein some individuals undergo long-distance anadromous migrations, and others stay as residents in their native freshwater streams. The decision to migrate is known to be highly heritable, and yet, the underlying genes and alleles associated with migration are not fully characterized. Here we used a pooled approach of whole-genome sequence data from migratory and resident trout of two native populations-Sashin Creek, Alaska and Little Sheep Creek, Oregon-to obtain a genome-wide perspective of the genetic architecture of resident and migratory life history. We calculated estimates of genetic differentiation, genetic diversity, and selection between the two phenotypes to locate regions of interest and then compared these associations between populations. We identified numerous genes and alleles associated with life history development in the Sashin Creek population with a notable area on chromosome 8 that may play a critical role in the development of the migratory phenotype. However, very few alleles appeared to be associated with life history development in the Little Sheep Creek system, suggesting population-specific genetic effects are likely important in the development of anadromy. Our results indicate that a migratory life history is not controlled by a singular gene or region but supports the idea that there are many independent ways for a migratory phenotype to emerge in a population. Therefore, conserving and promoting genetic diversity in migratory individuals is paramount to conserving these populations. Ultimately, our data add to a growing body of literature that suggests that population-specific genetic effects, likely mediated through environmental variation, contribute to life history development in rainbow trout.

C4S DB: Comprehensive Collection and Comparison for ChIP-Seq Database

Combining multiple binding profiles, such as transcription factors and histone modifications, is a crucial step in revealing the functions of complex biological systems. Although a massive amount of chromatin immunoprecipitation followed by sequencing (ChIP-seq) data is available, existing ChIP-seq databases or repositories focus on individual experiments, and it is difficult to elucidate orchestrated regulation by DNA-binding elements. We developed the Comprehensive Collection and Comparison for ChIP-Seq Database (C4S DB) to provide researchers with insights into the combination of DNA binding elements based on quality-assessed public ChIP-seq data. The C4S DB is based on > 16,000 human ChIP-seq experiments and provides two main web interfaces to discover the relationships between ChIP-seq data. "Gene browser" illustrates the landscape of distributions of binding elements around a specified gene, and "global similarity," a hierarchical clustering heatmap based on a similarity between two ChIP-seq experiments, gives an overview of genome-wide relations of regulatory elements. These functions promote the identification or evaluation of both gene-specific and genome-wide colocalization or mutually exclusive localization. Modern web technologies allow users to search for and aggregate large-scale experimental data through interactive web interfaces with quick responses. The C4S DB is available at https://c4s.site.

Circlehunter: a tool to identify extrachromosomal circular DNA from ATAC-Seq data

In cancer, extrachromosomal circular DNA (ecDNA), or megabase-pair amplified circular DNA, plays an essential role in intercellular heterogeneity and tumor cell revolution because of its non-Mendelian inheritance. We developed circlehunter ( https://github.com/suda-huanglab/circlehunter ), a tool for identifying ecDNA from ATAC-Seq data using the enhanced chromatin accessibility of ecDNA. Using simulated data, we showed that circlehunter has an F1 score of 0.93 at 30× local depth and read lengths as short as 35 bp. Based on 1312 ecDNAs predicted from 94 publicly available datasets of ATAC-Seq assays, we found 37 oncogenes contained in these ecDNAs with amplification characteristics. In small cell lung cancer cell lines, ecDNA containing MYC leads to amplification of MYC and cis-regulates the expression of NEUROD1, resulting in an expression pattern consistent with the NEUROD1 high expression subtype and sensitive to Aurora kinase inhibitors. This showcases that circlehunter could serve as a valuable pipeline for the investigation of tumorigenesis.

Optimal HSF1 activation in response to acute cold stress in BAT requires nuclear TXNIP

While TXNIP (thioredoxin interacting protein) in the plasma membrane and vesicular location is known to negatively regulate cellular glucose uptake by facilitating glucose transporter endocytosis, the function of TXNIP in the nucleus is far less understood. Herein, we sought to determine the function of nuclear TXNIP in vivo, using a new HA-tagged TXNIP knock-in mouse model. We observed that TXNIP can be found in the nucleus of a variety of cells from different tissues including hepatocytes (liver), enterocytes (small intestine), exocrine cells (pancreas), and brown adipocytes (BAT). Further investigations into the role of nuclear TXNIP in BAT revealed that cold stress rapidly and transiently activated HSF1 (heat shock factor 1). HSF1 interaction with TXNIP during its activation is required for optimal HSF1 directed cold shock response in BAT.

The size and composition of haplotype reference panels impact the accuracy of imputation from low-pass sequencing in cattle

BACKGROUND

Low-pass sequencing followed by sequence variant genotype imputation is an alternative to the routine microarray-based genotyping in cattle. However, the impact of haplotype reference panels and their interplay with the coverage of low-pass whole-genome sequencing data have not been sufficiently explored in typical livestock settings where only a small number of reference samples is available.

METHODS

Sequence variant genotyping accuracy was compared between two variant callers, GATK and DeepVariant, in 50 Brown Swiss cattle with sequencing coverages ranging from 4- to 63-fold. Haplotype reference panels of varying sizes and composition were built with DeepVariant based on 501 individuals from nine breeds. High-coverage sequence data for 24 Brown Swiss cattle were downsampled to between 0.01- and 4-fold to mimic low-pass sequencing. GLIMPSE was used to infer sequence variant genotypes from the low-pass sequencing data using different haplotype reference panels. The accuracy of the sequence variant genotypes that were inferred from low-pass sequencing data was compared with sequence variant genotypes called from high-coverage data.

RESULTS

DeepVariant was used to establish bovine haplotype reference panels because it outperformed GATK in all evaluations. Within-breed haplotype reference panels were more accurate and efficient to impute sequence variant genotypes from low-pass sequencing than equally-sized multibreed haplotype reference panels for all target sample coverages and allele frequencies. F1 scores greater than 0.9, which indicate high harmonic means of recall and precision of called genotypes, were achieved with 0.25-fold sequencing coverage when large breed-specific haplotype reference panels (n = 150) were used. In absence of such large within-breed haplotype panels, variant genotyping accuracy from low-pass sequencing could be increased either by adding non-related samples to the haplotype reference panel or by increasing the coverage of the low-pass sequencing data. Sequence variant genotyping from low-pass sequencing was substantially less accurate when the reference panel lacked individuals from the target breed.

CONCLUSIONS

Variant genotyping is more accurate with DeepVariant than GATK. DeepVariant is therefore suitable to establish bovine haplotype reference panels. Medium-sized breed-specific haplotype reference panels and large multibreed haplotype reference panels enable accurate imputation of low-pass sequencing data in a typical cattle breed.

Genetically defined individual reference ranges for tryptase limit unnecessary procedures and unmask myeloid neoplasms

Serum tryptase is a biomarker used to aid in the identification of certain myeloid neoplasms, most notably systemic mastocytosis, where basal serum tryptase (BST) levels >20 ng/mL are a minor criterion for diagnosis. Although clonal myeloid neoplasms are rare, the common cause for elevated BST levels is the genetic trait hereditary α-tryptasemia (HαT) caused by increased germline TPSAB1 copy number. To date, the precise structural variation and mechanism(s) underlying elevated BST in HαT and the general clinical utility of tryptase genotyping, remain undefined. Through cloning, long-read sequencing, and assembling of the human tryptase locus from an individual with HαT, and validating our findings in vitro and in silico, we demonstrate that BST elevations arise from overexpression of replicated TPSAB1 loci encoding canonical α-tryptase protein owing to coinheritance of a linked overactive promoter element. Modeling BST levels based on TPSAB1 replication number, we generate new individualized clinical reference values for the upper limit of normal. Using this personalized laboratory medicine approach, we demonstrate the clinical utility of tryptase genotyping, finding that in the absence of HαT, BST levels >11.4 ng/mL frequently identify indolent clonal mast cell disease. Moreover, substantial BST elevations (eg, >100 ng/mL), which would ordinarily prompt bone marrow biopsy, can result from TPSAB1 replications alone and thus be within normal limits for certain individuals with HαT.

Catch me if you can: capturing microbial community transformation by extracellular DNA using Hi-C sequencing

The transformation of environmental microorganisms by extracellular DNA is an overlooked mechanism of horizontal gene transfer and evolution. It initiates the acquisition of exogenous genes and propagates antimicrobial resistance alongside vertical and conjugative transfers. We combined mixed-culture biotechnology and Hi-C sequencing to elucidate the transformation of wastewater microorganisms with a synthetic plasmid encoding GFP and kanamycin resistance genes, in the mixed culture of chemostats exposed to kanamycin at concentrations representing wastewater, gut and polluted environments (0.01-2.5-50-100 mg L^-1). We found that the phylogenetically distant Gram-negative Runella (102 Hi-C links), Bosea (35), Gemmobacter (33) and Zoogloea (24) spp., and Gram-positive Microbacterium sp. (90) were transformed by the foreign plasmid, under high antibiotic exposure (50 mg L^-1). In addition, the antibiotic pressure shifted the origin of aminoglycoside resistance genes from genomic DNA to mobile genetic elements on plasmids accumulating in microorganisms. These results reveal the power of Hi-C sequencing to catch and surveil the transfer of xenogenetic elements inside microbiomes.

LKB1 controls inflammatory potential through CRTC2-dependent histone acetylation

Deregulated inflammation is a critical feature driving the progression of tumors harboring mutations in the liver kinase B1 (LKB1), yet the mechanisms linking LKB1 mutations to deregulated inflammation remain undefined. Here, we identify deregulated signaling by CREB-regulated transcription coactivator 2 (CRTC2) as an epigenetic driver of inflammatory potential downstream of LKB1 loss. We demonstrate that LKB1 mutations sensitize both transformed and non-transformed cells to diverse inflammatory stimuli, promoting heightened cytokine and chemokine production. LKB1 loss triggers elevated CRTC2-CREB signaling downstream of the salt-inducible kinases (SIKs), increasing inflammatory gene expression in LKB1-deficient cells. Mechanistically, CRTC2 cooperates with the histone acetyltransferases CBP/p300 to deposit histone acetylation marks associated with active transcription (i.e., H3K27ac) at inflammatory gene loci, promoting cytokine expression. Together, our data reveal a previously undefined anti-inflammatory program, regulated by LKB1 and reinforced through CRTC2-dependent histone modification signaling, that links metabolic and epigenetic states to cell-intrinsic inflammatory potential.

Regulatory de novo mutations underlying intellectual disability

The genetic aetiology of a major fraction of patients with intellectual disability (ID) remains unknown. De novo mutations (DNMs) in protein-coding genes explain up to 40% of cases, but the potential role of regulatory DNMs is still poorly understood. We sequenced 63 whole genomes from 21 ID probands and their unaffected parents. In addition, we analysed 30 previously sequenced genomes from exome-negative ID probands. We found that regulatory DNMs were selectively enriched in fetal brain-specific enhancers as compared with adult brain enhancers. DNM-containing enhancers were associated with genes that show preferential expression in the prefrontal cortex. Furthermore, we identified recurrently mutated enhancer clusters that regulate genes involved in nervous system development (CSMD1, OLFM1, and POU3F3). Most of the DNMs from ID probands showed allele-specific enhancer activity when tested using luciferase assay. Using CRISPR-mediated mutation and editing of epigenomic marks, we show that DNMs at regulatory elements affect the expression of putative target genes. Our results, therefore, provide new evidence to indicate that DNMs in fetal brain-specific enhancers play an essential role in the aetiology of ID.

Widespread somatic L1 retrotransposition in normal colorectal epithelium

Throughout an individual's lifetime, genomic alterations accumulate in somatic cells1-11. However, the mutational landscape induced by retrotransposition of long interspersed nuclear element-1 (L1), a widespread mobile element in the human genome12-14, is poorly understood in normal cells. Here we explored the whole-genome sequences of 899 single-cell clones established from three different cell types collected from 28 individuals. We identified 1,708 somatic L1 retrotransposition events that were enriched in colorectal epithelium and showed a positive relationship with age. Fingerprinting of source elements showed 34 retrotransposition-competent L1s. Multidimensional analysis demonstrated that (1) somatic L1 retrotranspositions occur from early embryogenesis at a substantial rate, (2) epigenetic on/off of a source element is preferentially determined in the early organogenesis stage, (3) retrotransposition-competent L1s with a lower population allele frequency have higher retrotransposition activity and (4) only a small fraction of L1 transcripts in the cytoplasm are finally retrotransposed in somatic cells. Analysis of matched cancers further suggested that somatic L1 retrotransposition rate is substantially increased during colorectal tumourigenesis. In summary, this study illustrates L1 retrotransposition-induced somatic mosaicism in normal cells and provides insights into the genomic and epigenomic regulation of transposable elements over the human lifetime.

High-resolution structural variants catalogue in a large-scale whole genome sequenced bovine family cohort data

BACKGROUND

Structural variants (SVs) are chromosomal segments that differ between genomes, such as deletions, duplications, insertions, inversions and translocations. The genomics revolution enabled the discovery of sub-microscopic SVs via array and whole-genome sequencing (WGS) data, paving the way to unravel the functional impact of SVs. Recent human expression QTL mapping studies demonstrated that SVs play a disproportionally large role in altering gene expression, underlining the importance of including SVs in genetic analyses. Therefore, this study aimed to generate and explore a high-quality bovine SV catalogue exploiting a unique cattle family cohort data (total 266 samples, forming 127 trios).

RESULTS

We curated 13,731 SVs segregating in the population, consisting of 12,201 deletions, 1,509 duplications, and 21 multi-allelic CNVs (> 50-bp). Of these, we validated a subset of copy number variants (CNVs) utilising a direct genotyping approach in an independent cohort, indicating that at least 62% of the CNVs are true variants, segregating in the population. Among gene-disrupting SVs, we prioritised two likely high impact duplications, encompassing ORM1 and POPDC3 genes, respectively. Liver expression QTL mapping results revealed that these duplications are likely causing altered gene expression, confirming the functional importance of SVs. Although most of the accurately genotyped CNVs are tagged by single nucleotide polymorphisms (SNPs) ascertained in WGS data, most CNVs were not captured by individual SNPs obtained from a 50K genotyping array.

CONCLUSION

We generated a high-quality SV catalogue exploiting unique whole genome sequenced bovine family cohort data. Two high impact duplications upregulating the ORM1 and POPDC3 are putative candidates for postpartum feed intake and hoof health traits, thus warranting further investigation. Generally, CNVs were in low LD with SNPs on the 50K array. Hence, it remains crucial to incorporate CNVs via means other than tagging SNPs, such as investigation of tagging haplotypes, direct imputation of CNVs, or direct genotyping as done in the current study. The SV catalogue and the custom genotyping array generated in the current study will serve as valuable resources accelerating utilisation of full spectrum of genetic variants in bovine genomes.

RNA-Seq-Pop: Exploiting the sequence in RNA sequencing-A Snakemake workflow reveals patterns of insecticide resistance in the malaria vector Anopheles gambiae

We provide a reproducible and scalable Snakemake workflow, called RNA-Seq-Pop, which provides end-to-end analysis of RNA sequencing data sets. The workflow allows the user to perform quality control, perform differential expression analyses and call genomic variants. Additional options include the calculation of allele frequencies of variants of interest, summaries of genetic variation and population structure, and genome-wide selection scans, together with clear visualizations. RNA-Seq-Pop is applicable to any organism, and we demonstrate the utility of the workflow by investigating pyrethroid resistance in selected strains of the major malaria mosquito, Anopheles gambiae. The workflow provides additional modules specifically for An. gambiae, including estimating recent ancestry and determining the karyotype of common chromosomal inversions. The Busia laboratory colony used for selections was collected in Busia, Uganda, in November 2018. We performed a comparative analysis of three groups: a parental G24 Busia strain; its deltamethrin-selected G28 offspring; and the susceptible reference strain Kisumu. Measures of genetic diversity reveal patterns consistent with that of laboratory colonization and selection, with the parental Busia strain exhibiting the highest nucleotide diversity, followed by the selected Busia offspring, and finally, Kisumu. Differential expression and variant analyses reveal that the selected Busia colony exhibits a number of distinct mechanisms of pyrethroid resistance, including the Vgsc-995S target-site mutation, upregulation of SAP genes, P450s and a cluster of carboxylesterases. During deltamethrin selections, the 2La chromosomal inversion rose in frequency (from 33% to 86%), supporting a previous link with pyrethroid resistance. RNA-Seq-Pop is hosted at: github.com/sanjaynagi/rna-seq-pop. We anticipate that the workflow will provide a useful tool to facilitate reproducible, transcriptomic studies in An. gambiae and other taxa.

Incomplete dosage balance and dosage compensation in the ZZ/ZW Gila monster (Heloderma suspectum) revealed by de novo genome assembly

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z-chromosome scaffolds. These putative Z-chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Swi/Snf Chromatin Remodeling Regulates Transcriptional Interference and Gene Repression

Alternative transcription start sites can affect transcript isoform diversity and translation levels. In a recently described form of gene regulation, coordinated transcriptional and translational interference results in transcript isoform-dependent changes in protein expression. Specifically, a long undecoded transcript isoform (LUTI) is transcribed from a gene-distal promoter, interfering with expression of the gene-proximal promoter. While transcriptional and chromatin features associated with LUTI expression have been described, the mechanism underlying LUTI-based transcriptional interference is not well understood. Using an unbiased genetic approach followed by integrated genomic analysis, we uncovered that the Swi/Snf chromatin remodeling complex is required for co-transcriptional nucleosome remodeling that leads to LUTI-based repression. We identified genes with tandem promoters that rely on Swi/Snf function for transcriptional interference during protein folding stress, including LUTI-regulated genes. To our knowledge, this study is the first to observe Swi/Snf's direct involvement in gene repression via a cis transcriptional interference mechanism.

Chromosome-level genome assembly and population genomic resource to accelerate orphan crop lablab breeding

Under-utilised orphan crops hold the key to diversified and climate-resilient food systems. Here, we report on orphan crop genomics using the case of Lablab purpureus (L.) Sweet (lablab) - a legume native to Africa and cultivated throughout the tropics for food and forage. Our Africa-led plant genome collaboration produces a high-quality chromosome-scale assembly of the lablab genome. Our assembly highlights the genome organisation of the trypsin inhibitor genes - an important anti-nutritional factor in lablab. We also re-sequence cultivated and wild lablab accessions from Africa confirming two domestication events. Finally, we examine the genetic and phenotypic diversity in a comprehensive lablab germplasm collection and identify genomic loci underlying variation of important agronomic traits in lablab. The genomic data generated here provide a valuable resource for lablab improvement. Our inclusive collaborative approach also presents an example that can be explored by other researchers sequencing indigenous crops, particularly from low and middle-income countries (LMIC).

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts the influence of few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that a few large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling, where some (but not all) of these QTL loci were shared across all populations, and some (but not all) of these loci showed signatures of recent natural selection in the corresponding wild population. Some strong candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

Chromosome-Level Genome Assembly of the Blue Mussel Mytilus chilensis Reveals Molecular Signatures Facing the Marine Environment

The blue mussel Mytilus chilensis is an endemic and key socioeconomic species inhabiting the southern coast of Chile. This bivalve species supports a booming aquaculture industry, which entirely relies on artificially collected seeds from natural beds that are translocated to diverse physical-chemical ocean farming conditions. Furthermore, mussel production is threatened by a broad range of microorganisms, pollution, and environmental stressors that eventually impact its survival and growth. Herein, understanding the genomic basis of the local adaption is pivotal to developing sustainable shellfish aquaculture. We present a high-quality reference genome of M. chilensis, which is the first chromosome-level genome for a Mytilidae member in South America. The assembled genome size was 1.93 Gb, with a contig N50 of 134 Mb. Through Hi-C proximity ligation, 11,868 contigs were clustered, ordered, and assembled into 14 chromosomes in congruence with the karyological evidence. The M. chilensis genome comprises 34,530 genes and 4795 non-coding RNAs. A total of 57% of the genome contains repetitive sequences with predominancy of LTR-retrotransposons and unknown elements. Comparative genome analysis of M. chilensis and M. coruscus was conducted, revealing genic rearrangements distributed into the whole genome. Notably, transposable Steamer-like elements associated with horizontal transmissible cancer were explored in reference genomes, suggesting putative relationships at the chromosome level in Bivalvia. Genome expression analysis was also conducted, showing putative genomic differences between two ecologically different mussel populations. The evidence suggests that local genome adaptation and physiological plasticity can be analyzed to develop sustainable mussel production. The genome of M. chilensis provides pivotal molecular knowledge for the Mytilus complex.

Synthetic regulatory genomics uncovers enhancer context dependence at the Sox2 locus

Sox2 expression in mouse embryonic stem cells (mESCs) depends on a distal cluster of DNase I hypersensitive sites (DHSs), but their individual contributions and degree of interdependence remain a mystery. We analyzed the endogenous Sox2 locus using Big-IN to scarlessly integrate large DNA payloads incorporating deletions, rearrangements, and inversions affecting single or multiple DHSs, as well as surgical alterations to transcription factor (TF) recognition sequences. Multiple mESC clones were derived for each payload, sequence-verified, and analyzed for Sox2 expression. We found that two DHSs comprising a handful of key TF recognition sequences were each sufficient for long-range activation of Sox2 expression. By contrast, three nearby DHSs were entirely context dependent, showing no activity alone but dramatically augmenting the activity of the autonomous DHSs. Our results highlight the role of context in modulating genomic regulatory element function, and our synthetic regulatory genomics approach provides a roadmap for the dissection of other genomic loci.

Multiomic analysis of malignant pleural mesothelioma identifies molecular axes and specialized tumor profiles driving intertumor heterogeneity

Malignant pleural mesothelioma (MPM) is an aggressive cancer with rising incidence and challenging clinical management. Through a large series of whole-genome sequencing data, integrated with transcriptomic and epigenomic data using multiomics factor analysis, we demonstrate that the current World Health Organization classification only accounts for up to 10% of interpatient molecular differences. Instead, the MESOMICS project paves the way for a morphomolecular classification of MPM based on four dimensions: ploidy, tumor cell morphology, adaptive immune response and CpG island methylator profile. We show that these four dimensions are complementary, capture major interpatient molecular differences and are delimited by extreme phenotypes that-in the case of the interdependent tumor cell morphology and adapted immune response-reflect tumor specialization. These findings unearth the interplay between MPM functional biology and its genomic history, and provide insights into the variations observed in the clinical behavior of patients with MPM.

Analyzing genomic and epigenetic profiles in single cells by hybrid transposase (scGET-seq)

scGET-seq simultaneously profiles euchromatin and heterochromatin. scGET-seq exploits the concurrent action of transposase Tn5 and its hybrid form TnH, which targets H3K9me3 domains. Here we present a step-by-step protocol to profile single cells by scGET-seq using a 10× Chromium Controller. We describe steps for transposomes preparation and validation. We detail nuclei preparation and transposition, followed by encapsulation, library preparation, sequencing, and data analysis. For complete details on the use and execution of this protocol, please refer to Tedesco et al. (2022)1 and de Pretis and Cittaro (2022).2.

Detailed molecular and epigenetic characterization of the pig IPEC-J2 and chicken SL-29 cell lines

The pig IPEC-J2 and chicken SL-29 cell lines are of interest because of their untransformed nature and wide use in functional studies. Molecular characterization of these cell lines is important to gain insight into possible molecular aberrations. The aim of this paper is to provide a molecular and epigenetic characterization of the IPEC-J2 and SL-29 cell lines, a cell-line reference for the FAANG community, and future biomedical research. Whole genome sequencing, gene expression, DNA methylation, chromatin accessibility, and ChIP-seq of four histone marks (H3K4me1, H3K4me3, H3K27ac, H3K27me3) and an insulator (CTCF) are used to achieve these aims. Heteroploidy (aneuploidy) of various chromosomes was observed from whole genome sequencing analysis in both cell lines. Furthermore, higher gene expression for genes located on chromosomes with aneuploidy in comparison to diploid chromosomes was observed. Regulatory complexity of gene expression, DNA methylation, and chromatin accessibility was investigated through an integrative approach.

Multi-omics analysis uncovers tumor ecosystem dynamics during neoadjuvant toripalimab plus nab-paclitaxel and S-1 for esophageal squamous cell carcinoma: a single-center, open-label, single-arm phase 2 trial

BACKGROUND

Immune checkpoint inhibitors combined with chemotherapy as a neoadjuvant therapy have been applied to the treatment of esophageal squamous cell carcinoma (ESCC). However, the optimal regimen needs to be further explored, particularly for older patients, and the mechanisms by which the immune checkpoint inhibitor combined with chemotherapy modulates the evolution of ESCC are unknown.

METHODS

In this single-arm phase 2 trial, patients with resectable (stage II/III/IV without metastasis) ESCC were enrolled and received nanoparticle albumin-bound (nab) paclitaxel for two cycles and oral S-1 for 2 weeks, combined with intravenous toripalimab for two cycles before surgery. Combination postoperative adjuvant therapy was administered. The primary outcome was the major pathological response (MPR). Secondary outcomes included pathological complete response (pCR), overall response rate (ORR), disease control rate (DCR), disease-free survival (DFS), overall survival (OS), improvement in Stooler's dysphagia score and degree of daily living ability (dADL). Biopsies and plasma pre- and post-neoadjuvant therapy were performed using whole-exome sequencing, transcriptome sequencing, immunohistochemistry (IHC) for PD-L1, multiplex immunofluorescence (mIF) and proximity extension assay technology (PEA) for 92 proteins.

FINDINGS

From November 2019 to July 2021, 60 patients were enrolled. After neoadjuvant therapy, R0 resection was achieved in 55 (98.21%) patients. MPR was identified in 27 patients (49.09%), and 16 patients (29.09%) achieved pCR. Patients with PR, SD and PD were 37 (61.67%), 21 (35.00%) and 2 (3.33%), respectively. The overall staging, Stooler dysphagia scores and dADL were significantly decreased after treatment. 11 patients (18.3%) experienced grade ≥3 AEs. Compared to PD-L1-Low patients, PD-L1-High patients had a significantly higher ratio of PR. During therapy, the tumor mutation burden (TMB) and tumor neoantigen burden (TNB) were significantly decreased in patients with PR. Differential clonal evolution within tumors was demonstrated by analysis of intratumoral heterogeneity. Transcriptome analyses revealed that the infiltration of CD4+ T lymphocytes at baseline was associated with clinical outcome. During therapy, CD8+ T cells and CD4+ T cells were increased in all patients; however, exhausted cells, nTregs and iTregs were significantly increased in patients with non-MPR. Protein analyses revealed that the levels of IFN-γ, Gal.1 and LAMP3 can predict the clinical benefit. In addition, the expression of CD83, TNFRSF4, TNFSF14, VEGFR2, ADA, ARG1, and HO-1 was associated with serious AEs. More importantly, the integration of CD4+ T cells with plasma protein of IFN-γ, Gal.1 or LAMP3 could further distinguish responders from non-responders.

INTERPRETATION

In this study, neoadjuvant therapy with toripalimab, nab-paclitaxel and S-1 was less toxic and showed promising antitumor activity in patients with resectable ESCC. Changes in the genome, transcriptome, PD-L1 expression and serum proteins were comprehensively analyzed and correlated with clinical outcomes, which provides insight into the mechanism of action of toripalimab combined with nab-paclitaxel and S-1 in patients with ESCC.

FUNDING

This study was funded by Major projects of the ministry of science and technology of the 13th five-year plan of China [grant number: 2018ZX09201013].

Triap1 upregulation promotes escape from mitotic-slippage-induced G1 arrest

Drugs targeting microtubules rely on the mitotic checkpoint to arrest cell proliferation. The prolonged mitotic arrest induced by such drugs is followed by a G1 arrest. Here, we follow for several weeks the fate of G1-arrested human cells after treatment with nocodazole. We find that a small fraction of cells escapes from the arrest and resumes proliferation. These escaping cells experience reduced DNA damage and p21 activation. Cells surviving treatment are enriched for anti-apoptotic proteins, including Triap1. Increasing Triap1 levels allows cells to survive the first treatment with reduced DNA damage and lower levels of p21; accordingly, decreasing Triap1 re-sensitizes cells to nocodazole. We show that Triap1 upregulation leads to the retention of cytochrome c in the mitochondria, opposing the partial activation of caspases caused by nocodazole. In summary, our results point to a potential role of Triap1 upregulation in the emergence of resistance to drugs that induce prolonged mitotic arrest.

The Driverless Triple-Wild-Type (BRAF, RAS, KIT) Cutaneous Melanoma: Whole Genome Sequencing Discoveries

The genetic makeup of the triple-wild-type melanoma (BRAF, NRAS and NF1) has been known for some time, but those studies grouped together rare histopathological versions with common ones, as well as mucosal and even uveal ones. Here we used whole genome sequencing to genetically characterize the triple-wild-type melanoma (TWM), termed here as BRAF, RAS and KIT wild type (the most frequent oncogenic drivers of skin melanoma), using the most common histological forms and excluding rare ones. All these tumors except one were clearly induced by UV based on the mutational signature. The tumor mutational burden was low in TWM, except in the NF1 mutant forms, and a relatively high frequency of elevated LOH scores suggested frequent homologue recombination deficiency, but this was only confirmed by the mutation signature in one case. Furthermore, all these TWMs were microsatellite-stabile. In this driverless setting, we revealed rare oncogenic drivers known from melanoma or other cancer types and identified rare actionable tyrosine kinase mutations in NTRK1, RET and VEGFR1. Mutations of TWM identified genes involved in antitumor immunity (negative and positive predictors of immunotherapy), Ca++ and BMP signaling. The two regressed melanomas of this cohort shared a 17-gene mutation signature, containing genes involved in antitumor immunity and several cell surface receptors. Even with this comprehensive genomic approach, a few cases remained driverless, suggesting that unrecognized drivers are hiding among passenger mutations.