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

Senescent characteristics of human corneal endothelial cells upon ultraviolet-A exposure

PURPOSE

The objective of this study was to investigate the senescent phenotypes of human corneal endothelial cells (hCEnCs) upon treatment with ultraviolet (UV)-A.

METHODS

We assessed cell morphology, senescence-associated β-galactosidase (SA-β-gal) activity, cell proliferation and expression of senescence markers (p16 and p21) in hCEnCs exposed to UV-A radiation, and senescent hCEnCs induced by ionizing radiation (IR) were used as positive controls. We performed RNA sequencing and proteomics analyses to compare gene and protein expression profiles between UV-A- and IR-induced senescent hCEnCs, and we also compared the results to non-senescent hCEnCs.

RESULTS

Cells exposed to 5 J/cm2 of UV-A or to IR exhibited typical senescent phenotypes, including enlargement, increased SA-β-gal activity, decreased cell proliferation and elevated expression of p16 and p21. RNA-Seq analysis revealed that 83.9% of the genes significantly upregulated and 82.6% of the genes significantly downregulated in UV-A-induced senescent hCEnCs overlapped with the genes regulated in IR-induced senescent hCEnCs. Proteomics also revealed that 93.8% of the proteins significantly upregulated in UV-A-induced senescent hCEnCs overlapped with those induced by IR. In proteomics analyses, senescent hCEnCs induced by UV-A exhibited elevated expression levels of several factors part of the senescence-associated secretory phenotype.

CONCLUSIONS

In this study, where senescence was induced by UV-A, a more physiological stress for hCEnCs compared to IR, we determined that UV-A modulated the expression of many genes and proteins typically altered upon IR treatment, a more conventional method of senescence induction, even though UV-A also modulated specific pathways unrelated to IR.

Multi-omic features of oesophageal adenocarcinoma in patients treated with preoperative neoadjuvant therapy

Oesophageal adenocarcinoma is a poor prognosis cancer and the molecular features underpinning response to treatment remain unclear. We investigate whole genome, transcriptomic and methylation data from 115 oesophageal adenocarcinoma patients mostly from the DOCTOR phase II clinical trial (Australian New Zealand Clinical Trials Registry-ACTRN12609000665235), with exploratory analysis pre-specified in the study protocol of the trial. We report genomic features associated with poorer overall survival, such as the APOBEC mutational and RS3-like rearrangement signatures. We also show that positron emission tomography non-responders have more sub-clonal genomic copy number alterations. Transcriptomic analysis categorises patients into four immune clusters correlated with survival. The immune suppressed cluster is associated with worse survival, enriched with myeloid-derived cells, and an epithelial-mesenchymal transition signature. The immune hot cluster is associated with better survival, enriched with lymphocytes, myeloid-derived cells, and an immune signature including CCL5, CD8A, and NKG7. The immune clusters highlight patients who may respond to immunotherapy and thus may guide future clinical trials.

Efficient in vitro and in vivo RNA editing via recruitment of endogenous ADARs using circular guide RNAs

Recruiting endogenous adenosine deaminases using exogenous guide RNAs to edit cellular RNAs is a promising therapeutic strategy, but editing efficiency and durability remain low using current guide RNA designs. In this study, we engineered circular ADAR-recruiting guide RNAs (cadRNAs) to enable more efficient programmable adenosine-to-inosine RNA editing without requiring co-delivery of any exogenous proteins. Using these cadRNAs, we observed robust and durable RNA editing across multiple sites and cell lines, in both untranslated and coding regions of RNAs, and high transcriptome-wide specificity. Additionally, we increased transcript-level specificity for the target adenosine by incorporating interspersed loops in the antisense domains, reducing bystander editing. In vivo delivery of cadRNAs via adeno-associated viruses enabled 53% RNA editing of the mPCSK9 transcript in C57BL/6J mice livers and 12% UAG-to-UGG RNA correction of the amber nonsense mutation in the IDUA-W392X mouse model of mucopolysaccharidosis type I-Hurler syndrome. cadRNAs enable efficient programmable RNA editing in vivo with diverse protein modulation and gene therapeutic applications.

Therapeutic homology-independent targeted integration in retina and liver

Challenges to the widespread application of gene therapy with adeno-associated viral (AAV) vectors include dominant conditions due to gain-of-function mutations which require allele-specific knockout, as well as long-term transgene expression from proliferating tissues, which is hampered by AAV DNA episomal status. To overcome these challenges, we used CRISPR/Cas9-mediated homology-independent targeted integration (HITI) in retina and liver as paradigmatic target tissues. We show that AAV-HITI targets photoreceptors of both mouse and pig retina, and this results in significant improvements to retinal morphology and function in mice with autosomal dominant retinitis pigmentosa. In addition, we show that neonatal systemic AAV-HITI delivery achieves stable liver transgene expression and phenotypic improvement in a mouse model of a severe lysosomal storage disease. We also show that HITI applications predominantly result in on-target editing. These results lay the groundwork for the application of AAV-HITI for the treatment of diseases affecting various organs.

Inter-species Metabolic Interactions in an In-vitro Minimal Human Gut Microbiome of Core Bacteria

Knowledge of the functional roles and interspecies interactions are crucial for improving our understanding of the human intestinal microbiome in health and disease. However, the complexity of the human intestinal microbiome and technical challenges in investigating it pose major challenges. In this proof-of-concept study, we rationally designed, assembled and experimentally tested a synthetic Diet-based Minimal Microbiome (Db-MM) consisting of ten core intestinal bacterial species that together are capable of efficiently converting dietary fibres into short chain fatty acids (SCFAs). Despite their genomic potential for metabolic competition, all ten bacteria coexisted during growth on a mixture of dietary fibres, including pectin, inulin, xylan, cellobiose and starch. By integrated analyses of metabolite production, community composition and metatranscriptomics-based gene expression data, we identified interspecies metabolic interactions leading to production of key SCFAs such as butyrate and propionate. While public goods, such as sugars liberated from colonic fibres, are harvested by non-degraders, some species thrive by cross-feeding on energetically challenging substrates, including the butyrogenic conversion of acetate and lactate. Using a reductionist approach in an in-vitro system combined with functional measurements, our study provides key insights into the complex interspecies metabolic interactions between core intestinal bacterial species.

Global population genomic signature of Spodoptera frugiperda (fall armyworm) supports complex introduction events across the Old World

Native to the Americas, the invasive Spodoptera frugiperda (fall armyworm; FAW) was reported in West Africa in 2016, followed by its chronological detection across the Old World and the hypothesis of an eastward Asia expansion. We explored population genomic signatures of American and Old World FAW and identified 12 maternal mitochondrial DNA genome lineages across the invasive range. 870 high-quality nuclear single nucleotide polymorphic DNA markers identified five distinct New World population clusters, broadly reflecting FAW native geographical ranges and the absence of host-plant preferences. We identified unique admixed Old World populations, and admixed and non-admixed Asian FAW individuals, all of which suggested multiple introductions underpinning the pest’s global spread. Directional gene flow from the East into eastern Africa was also detected, in contrast to the west-to-east spread hypothesis. Our study demonstrated the potential of population genomic approaches via international partnership to address global emerging pest threats and biosecurity challenges.

This population genomics study identifies the complex multiple introduction history of Spodoptera frugiperda (fall armyworm) from the Americas, into Africa and Asia. This provides new insight into the ‘east-to-west’ directionality of gene flow, and suggests ample genomic exchange at the nuclear level.

Identification of neoantigens for individualized therapeutic cancer vaccines

Somatic mutations in cancer cells can generate tumour-specific neoepitopes, which are recognized by autologous T cells in the host. As neoepitopes are not subject to central immune tolerance and are not expressed in healthy tissues, they are attractive targets for therapeutic cancer vaccines. Because the vast majority of cancer mutations are unique to the individual patient, harnessing the full potential of this rich source of targets requires individualized treatment approaches. Many computational algorithms and machine-learning tools have been developed to identify mutations in sequence data, to prioritize those that are more likely to be recognized by T cells and to design tailored vaccines for every patient. In this Review, we fill the gaps between the understanding of basic mechanisms of T cell recognition of neoantigens and the computational approaches for discovery of somatic mutations and neoantigen prediction for cancer immunotherapy. We present a new classification of neoantigens, distinguishing between guarding, restrained and ignored neoantigens, based on how they confer proficient antitumour immunity in a given clinical context. Such context-based differentiation will contribute to a framework that connects neoantigen biology to the clinical setting and medical peculiarities of cancer, and will enable future neoantigen-based therapies to provide greater clinical benefit.

Exome sequencing of individuals with Huntington's disease implicates FAN1 nuclease activity in slowing CAG expansion and disease onset

The age at onset of motor symptoms in Huntington's disease (HD) is driven by HTT CAG repeat length but modified by other genes. In this study, we used exome sequencing of 683 patients with HD with extremes of onset or phenotype relative to CAG length to identify rare variants associated with clinical effect. We discovered damaging coding variants in candidate modifier genes identified in previous genome-wide association studies associated with altered HD onset or severity. Variants in FAN1 clustered in its DNA-binding and nuclease domains and were associated predominantly with earlier-onset HD. Nuclease activities of purified variants in vitro correlated with residual age at motor onset of HD. Mutating endogenous FAN1 to a nuclease-inactive form in an induced pluripotent stem cell model of HD led to rates of CAG expansion similar to those observed with complete FAN1 knockout. Together, these data implicate FAN1 nuclease activity in slowing somatic repeat expansion and hence onset of HD.

Fast and accurate metagenotyping of the human gut microbiome with GT-Pro

Single nucleotide polymorphisms (SNPs) in metagenomics are used to quantify population structure, track strains and identify genetic determinants of microbial phenotypes. However, existing alignment-based approaches for metagenomic SNP detection require high-performance computing and enough read coverage to distinguish SNPs from sequencing errors. To address these issues, we developed the GenoTyper for Prokaryotes (GT-Pro), a suite of methods to catalog SNPs from genomes and use unique k-mers to rapidly genotype these SNPs from metagenomes. Compared to methods that use read alignment, GT-Pro is more accurate and two orders of magnitude faster. Using high-quality genomes, we constructed a catalog of 104 million SNPs in 909 human gut species and used unique k-mers targeting this catalog to characterize the global population structure of gut microbes from 7,459 samples. GT-Pro enables fast and memory-efficient metagenotyping of millions of SNPs on a personal computer.

Reprogramming of fibroblasts into expandable cardiovascular progenitor cells via small molecules in xeno-free conditions

A major hurdle in cardiac cell therapy is the lack of a bona fide autologous stem-cell type that can be expanded long-term and has authentic cardiovascular differentiation potential. Here we report that a proliferative cell population with robust cardiovascular differentiation potential can be generated from mouse or human fibroblasts via a combination of six small molecules. These chemically induced cardiovascular progenitor cells (ciCPCs) self-renew long-term in fully chemically defined and xeno-free conditions, with faithful preservation of the CPC phenotype and of cardiovascular differentiation capacity in vitro and in vivo. Transplantation of ciCPCs into infarcted mouse hearts improved animal survival and cardiac function up to 13 weeks post-infarction. Mechanistically, activated fibroblasts revert to a plastic state permissive to cardiogenic signals, enabling their reprogramming into ciCPCs. Expanded autologous cardiovascular cells may find uses in drug discovery, disease modelling and cardiac cell therapy.

Droplet-microfluidics-assisted sequencing of HIV proviruses and their integration sites in cells from people on antiretroviral therapy

The human immunodeficiency virus (HIV) integrates its genome in that of infected cells and may enter an inactive state of reversible latency that cannot be targeted using antiretroviral therapy. The resulting HIV DNA is termed a provirus. Sequencing individual proviruses with the adjacent human cellular junctions may elucidate mechanisms of infected cell persistence in humans. Here, we introduce a high throughput microfluidic assay where droplet-based whole genome amplification of the HIV DNA in its native context is followed by a polymerase chain reaction to tag droplets containing proviruses for sequencing, resulting in the assembly of full-length viral genomes connected to their contiguous HIV-human DNA junctions, regardless of the 150 million-fold higher amount of human DNA present in the background. We analyzed infected cells from patients with HIV receiving suppressive antiretroviral therapy, resulting in the detection and sequencing of paired proviral genomes and integration sites, 90% of which weren’t recovered by commonly used nested PCR methods. The sequencing of individual proviral genomes with their integration sites could improve the genetic analysis of persistent HIV-infected cell reservoirs.

Resistomes and microbiome of meat trimmings and colon content from culled cows raised in conventional and organic production systems

Background

The potential to distribute bacteria resistant to antimicrobial drugs in the meat supply is a public health concern. Market cows make up a fifth of the U.S. beef produced but little is known about the entire population of bacteria (the microbiome) and entirety of all resistance genes (the resistome) that are found in this population. The objective of this study was to characterize and compare the resistomes and microbiome of beef, dairy, and organic dairy market cows at slaughter.

Methods

Fifty-four (N = 54) composite samples of both colon content and meat trimmings rinsate samples were collected over six visits to two harvest facilities from cows raised in three different production systems: conventional beef, conventional dairy, and organic dairy (n = 3 samples per visit per production system). Metagenomic DNA obtained from samples were analyzed using target-enriched sequencing (resistome) and 16S rRNA gene sequencing (microbiome).

Results

All colon content samples had at least one identifiable antimicrobial resistance gene (ARG), while 21 of the 54 meat trimmings samples harbored at least one identifiable ARGs. Tetracycline ARGs were the most abundant class in both colon content and carcass meat trimmings. The resistome found on carcass meat trimmings was not significantly different by production system (P = 0.84, R² = 0.00) or harvest facility (P = 0.10, R² = 0.09). However, the resistome of colon content differed (P = 0.01; R² = 0.05) among production systems, but not among the harvest facilities (P = 0.41; R² = 0.00). Amplicon sequencing revealed differences (P < 0.05) in microbial populations in both meat trimmings and colon content between harvest facilities but not production systems (P > 0.05).

Conclusions

These data provide a baseline characterization of an important segment of the beef industry and highlight the effect that the production system where cattle are raised and the harvest facilities where an animal is processed can impact associated microbiome and resistomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-022-00166-z.

The genome of an apodid holothuroid (Chiridota heheva) provides insights into its adaptation to a deep-sea reducing environment

Cold seeps and hydrothermal vents are deep-sea reducing environments that are characterized by lacking oxygen and photosynthesis-derived nutrients. Most animals acquire nutrition in cold seeps or hydrothermal vents by maintaining epi- or endosymbiotic relationship with chemoautotrophic microorganisms. Although several seep- and vent-dwelling animals hosting symbiotic microbes have been well-studied, the genomic basis of adaptation to deep-sea reducing environment in nonsymbiotic animals is still lacking. Here, we report a high-quality genome of Chiridota heheva Pawson & Vance, 2004, which thrives by extracting organic components from sediment detritus and suspended material, as a reference for nonsymbiotic animal's adaptation to deep-sea reducing environments. The expansion of the aerolysin-like protein family in C. heheva compared with other echinoderms might be involved in the disintegration of microbes during digestion. Moreover, several hypoxia-related genes (Pyruvate Kinase M2, PKM2; Phospholysine Phosphohistidine Inorganic Pyrophosphate Phosphatase, LHPP; Poly(A)-specific Ribonuclease Subunit PAN2, PAN2; and Ribosomal RNA Processing 9, RRP9) were subject to positive selection in the genome of C. heheva, which contributes to their adaptation to hypoxic environments.

Chromosome-scale and haplotype-resolved genome assembly of a tetraploid potato cultivar

Potato is the most widely produced tuber crop worldwide. However, reconstructing the four haplotypes of its autotetraploid genome remained an unsolved challenge. Here, we report the 3.1 Gb haplotype-resolved (at 99.6% precision), chromosome-scale assembly of the potato cultivar ‘Otava’ based on high-quality long reads, single-cell sequencing of 717 pollen genomes and Hi-C data. Unexpectedly, ~50% of the genome was identical-by-descent due to recent inbreeding, which was contrasted by highly abundant structural rearrangements involving ~20% of the genome. Among 38,214 genes, only 54% were present in all four haplotypes with an average of 3.2 copies per gene. Taking the leaf transcriptome as an example, 11% of the genes were differently expressed in at least one haplotype, where 25% of them were likely regulated through allele-specific DNA methylation. Our work sheds light on the recent breeding history of potato, the functional organization of its tetraploid genome and has the potential to strengthen the future of genomics-assisted breeding.

Haplotype-resolved genome assembly of the tetraploid potato cultivar ‘Otava’ sheds light on functional organization of the tetraploid genome and provides the potential for genomics-assisted breeding.

Complex hereditary peripheral neuropathies caused by novel variants in mitochondrial-related nuclear genes

Mitochondrial disorders are a group of clinically and genetically heterogeneous multisystem disorders and peripheral neuropathy is frequently described in the context of mutations in mitochondrial-related nuclear genes. This study aimed to identify the causative mutations in mitochondrial-related nuclear genes in suspected hereditary peripheral neuropathy patients. We enrolled a large Japanese cohort of clinically suspected hereditary peripheral neuropathy patients who were mutation negative in the prescreening of the known Charcot–Marie–Tooth disease-causing genes. We performed whole-exome sequencing on 247 patients with autosomal recessive or sporadic inheritance for further analysis of 167 mitochondrial-related nuclear genes. We detected novel bi-allelic likely pathogenic/pathogenic variants in four patients, from four mitochondrial-related nuclear genes: pyruvate dehydrogenase beta-polypeptide (PDHB), mitochondrial poly(A) polymerase (MTPAP), hydroxyacyl-CoA dehydrogenase/3-ketoacyl-CoA thiolase/enoyl-CoA hydratase, beta subunit (HADHB), and succinate-CoA ligase ADP-forming beta subunit (SUCLA2). All these patients showed sensory and motor axonal polyneuropathy, combined with central nervous system or multisystem involvements. The pathological analysis of skeletal muscles revealed mild neurogenic changes without significant mitochondrial abnormalities. Targeted screening of mitochondria-related nuclear genes should be considered for patients with complex hereditary axonal polyneuropathy, accompanied by central nervous system dysfunctions, or with unexplainable multisystem disorders.

Transcriptomic analysis of frontotemporal lobar degeneration with TDP-43 pathology reveals cellular alterations across multiple brain regions

Frontotemporal lobar degeneration (FTLD) is a group of heterogeneous neurodegenerative disorders affecting the frontal and temporal lobes of the brain. Nuclear loss and cytoplasmic aggregation of the RNA-binding protein TDP-43 represents the major FTLD pathology, known as FTLD-TDP. To date, there is no effective treatment for FTLD-TDP due to an incomplete understanding of the molecular mechanisms underlying disease development. Here we compared postmortem tissue RNA-seq transcriptomes from the frontal cortex, temporal cortex, and cerebellum between 28 controls and 30 FTLD-TDP patients to profile changes in cell-type composition, gene expression and transcript usage. We observed downregulation of neuronal markers in all three regions of the brain, accompanied by upregulation of microglia, astrocytes, and oligodendrocytes, as well as endothelial cells and pericytes, suggesting shifts in both immune activation and within the vasculature. We validate our estimates of neuronal loss using neuropathological atrophy scores and show that neuronal loss in the cortex can be mainly attributed to excitatory neurons, and that increases in microglial and endothelial cell expression are highly correlated with neuronal loss. All our analyses identified a strong involvement of the cerebellum in the neurodegenerative process of FTLD-TDP. Altogether, our data provides a detailed landscape of gene expression alterations to help unravel relevant disease mechanisms in FTLD.

Expressed prognostic biomarkers for primary prostate cancer independent of multifocality and transcriptome heterogeneity

The majority of prostate cancer patients are diagnosed with multiple primary malignant foci. The distinct foci are exceptionally heterogeneous with regard to DNA mutations, but whether this is recapitulated at the transcriptome level remains unknown. In this study, inter- and intrafocal heterogeneity has been assessed by whole-transcriptome sequencing of 87 tissue samples from 23 patients with localized prostate cancer treated with radical prostatectomy. From each patient, multiple samples were taken from one or more malignant foci, in addition to one sample from benign prostate tissue. Transcriptomic profiles of different malignant foci from the same patient showed a similar level of heterogeneity as tumors from different patients. This applies to expression of genes, fusion genes, and somatic mutations. Within-patient pair-wise analyses identified expression patterns linked to ETS status and extraprostatic extension. A set of 62 genes were found with low intrapatient heterogeneity and high interpatient heterogeneity, retaining stable expression profiles across foci within the same patient. Among these, 16 genes are associated with biochemical recurrence in a separately published study and are therefore nominated as biomarkers with prognostic value regardless of which malignant focus is sampled. In conclusion, an extensive heterogeneity in multifocal prostate cancer is confirmed at the gene expression level. Diagnostic biomarkers were identified for ETS positive samples and samples from extraprostatic extensions. Finally, prognostic biomarkers independent of multifocal heterogeneity were found.

Optimized analysis parameters of variant calling for whole genome-based phylogeny of Mycobacteroides abscessus

Whole-genome sequence (WGS) analysis provides the best resolution for reconstructing bacterial phylogeny. However, the resulting tree could vary according to parameters used in the WGS pipeline, making it difficult to compare results across multiple studies. This study compares effects on phylogenies when applying different parameter stringencies. We used as the study model to optimize parameters strains of Mycobacteroides abscessus serially isolated at various intervals, isolates known to represent persistent infection (PI) cases or re-infection (RI) cases and isolates from different subspecies. Un-optimized parameters with low stringency provided an excessive number of SNPs (823) compared to the optimized setting (3 SNPs) between paired strains isolated 1 day apart from PI cases, discordant tree topology and misclassification of subspecies and of instances of RI. We demonstrated that using high-quality variants provides more accuracy for recognizing serial isolates of the same clone versus different clones and for phylogenetic analysis of M. abscessus. Our approach might be used as a model for analyses requiring phylogenetic reconstruction of other bacteria.

Genome binning of viral entities from bulk metagenomics data

Despite the accelerating number of uncultivated virus sequences discovered in metagenomics and their apparent importance for health and disease, the human gut virome and its interactions with bacteria in the gastrointestinal tract are not well understood. This is partly due to a paucity of whole-virome datasets and limitations in current approaches for identifying viral sequences in metagenomics data. Here, combining a deep-learning based metagenomics binning algorithm with paired metagenome and metavirome datasets, we develop Phages from Metagenomics Binning (PHAMB), an approach that allows the binning of thousands of viral genomes directly from bulk metagenomics data, while simultaneously enabling clustering of viral genomes into accurate taxonomic viral populations. When applied on the Human Microbiome Project 2 (HMP2) dataset, PHAMB recovered 6,077 high-quality genomes from 1,024 viral populations, and identified viral-microbial host interactions. PHAMB can be advantageously applied to existing and future metagenomes to illuminate viral ecological dynamics with other microbiome constituents.

A comparison of experimental assays and analytical methods for genome-wide identification of active enhancers

Mounting evidence supports the idea that transcriptional patterns serve as more specific identifiers of active enhancers than histone marks; however, the optimal strategy to identify active enhancers both experimentally and computationally has not been determined. Here, we compared 13 genome-wide RNA sequencing assays in K562 cells and showed that the nuclear run-on followed by cap-selection assay (GRO/PRO-cap) has advantages in eRNA detection and active enhancer identification. We also introduced a tool, Peak Identifier for Nascent Transcript Starts (PINTS), to identify active promoters and enhancers genome-wide and pinpoint the precise location of the 5′ transcription start sites. Finally, we compiled a comprehensive enhancer candidate compendium based on the detected eRNA TSSs available in 120 cell and tissue types that can be accessed at https://pints.yulab.org. With the knowledge of the best available assays and pipelines, this large-scale annotation of candidate enhancers will pave the way for selection and characterization of their functions in a time- and labor-efficient manner in the future.

Integrated single-cell multiomics analysis reveals novel candidate markers for prognosis in human pancreatic ductal adenocarcinoma

The epigenomic abnormality of pancreatic ductal adenocarcinoma (PDAC) has rarely been investigated due to its strong heterogeneity. Here, we used single-cell multiomics sequencing to simultaneously analyze the DNA methylome, chromatin accessibility and transcriptome in individual tumor cells of PDAC patients. We identified normal epithelial cells in the tumor lesion, which have euploid genomes, normal patterns of DNA methylation, and chromatin accessibility. Using all these normal epithelial cells as controls, we determined that DNA demethylation in the cancer genome was strongly enriched in heterochromatin regions but depleted in euchromatin regions. There were stronger negative correlations between RNA expression and promoter DNA methylation in cancer cells compared to those in normal epithelial cells. Through in-depth integrated analyses, a set of novel candidate biomarkers were identified, including ZNF667 and ZNF667-AS1, whose expressions were linked to a better prognosis for PDAC patients by affecting the proliferation of cancer cells. Our work systematically revealed the critical epigenomic features of cancer cells in PDAC patients at the single-cell level.

Osmotic stress induces long-term biofilm survival in Liberibacter crescens

Citrus greening, also known as Huanglongbing (HLB), is a devastating citrus plant disease caused predominantly by Liberibacter asiaticus. While nearly all Liberibacter species remain uncultured, here we used the culturable L. crescens BT-1 as a model to examine physiological changes in response to the variable osmotic conditions and nutrient availability encountered within the citrus host. Similarly, physiological responses to changes in growth temperature and dimethyl sulfoxide concentrations were also examined, due to their use in many of the currently employed therapies to control the spread of HLB. Sublethal heat stress was found to induce the expression of genes related to tryptophan biosynthesis, while repressing the expression of ribosomal proteins. Osmotic stress induces expression of transcriptional regulators involved in expression of extracellular structures, while repressing the biosynthesis of fatty acids and aromatic amino acids. The effects of osmotic stress were further evaluated by quantifying biofilm formation of L. crescens in presence of increasing sucrose concentrations at different stages of biofilm formation, where sucrose-induced osmotic stress delayed initial cell attachment while enhancing long-term biofilm viability. Our findings revealed that exposure to osmotic stress is a significant contributing factor to the long term survival of L. crescens and, possibly, to the pathogenicity of other Liberibacter species.

Whole genome re-sequencing reveals the genetic diversity and evolutionary patterns of Eucommia ulmoides

Eucommia ulmoides (E. ulmoides) is a deciduous perennial tree belonging to the order Garryales, and is known as "living fossil" plant, along with ginkgo (Ginkgo biloba), metaspaca (Metasequoia glyptostroboides) and dove tree (Davidia involucrata Baill). However, the genetic diversity and population structure of E. ulmoides are still  ambiguous nowdays. In this study, we re-sequenced the genomes of 12 E. ulmoides accessions from different major climatic geography regions in China to elucidate the genetic diversity, population structure and evolutionary pattern. By integration of phylogenetic analysis, principal component analysis and population structure analysis based on a number of high-quality SNPs, a total of 12 E. ulmoides accessions were clustered into four different groups. This result is consistent with their geographical location except for group samples from Shanghai and Hunan province. E. ulmoides accessions from Hunan province exhibited a closer genetic relationship with E. ulmoides accessions from Shanghai in China compared with other regions, which is also supported by the result of population structure analyses. Genetic diversity analysis further revealed that E. ulmoides samples in Shanghai and Hunan province were with higher genetic diversity than those in other regions in this study. In addition, we treated the E. ulmoides materials from Shanghai and Hunan province as group A, and the other materials from other places as group B, and then analyzed the evolutionary pattern of E. ulmoides. The result showed the significant differentiation (Fst = 0.1545) between group A and group B. Some candidate highly divergent genome regions were identified in group A by selective sweep analyses, and the function analysis of candidate genes in these regions showed that biological regulation processes could be correlated with the Eu-rubber biosynthesis. Notably, nine genes were identified from selective sweep regions. They were involved in the Eu-rubber biosynthesis and expressed in rubber containing tissues. The genetic diversity research and evolution model of E. ulmoides were preliminarily explored in this study, which laid the foundation for the protection of germplasm resources and the development and utilization of multipurpose germplasm resources in the future.

Long-read sequencing of the zebrafish genome reorganizes genomic architecture

Background

Nanopore sequencing technology has revolutionized the field of genome biology with its ability to generate extra-long reads that can resolve regions of the genome that were previously inaccessible to short-read sequencing platforms. Over 50% of the zebrafish genome consists of difficult to map, highly repetitive, low complexity elements that pose inherent problems for short-read sequencers and assemblers.

Results

We used long-read nanopore sequencing to generate a de novo assembly of the zebrafish genome and compared our assembly to the current reference genome, GRCz11. The new assembly identified 1697 novel insertions and deletions over one kilobase in length and placed 106 previously unlocalized scaffolds. We also discovered additional sites of retrotransposon integration previously unreported in GRCz11 and observed the expression of these transposable elements in adult zebrafish under physiologic conditions, implying they have active mobility in the zebrafish genome and contribute to the ever-changing genomic landscape.

Conclusions

We used nanopore sequencing to improve upon and resolve the issues plaguing the current zebrafish reference assembly, GRCz11. Zebrafish is a prominent model of human disease, and our corrected assembly will be useful for studies relying on interspecies comparisons and precise linkage of genetic events to disease phenotypes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08349-3.

PΨFinder: a practical tool for the identification and visualization of novel pseudogenes in DNA sequencing data

BACKGROUND

Processed pseudogenes (PΨgs) are disabled gene copies that are transcribed and may affect expression of paralogous genes. Moreover, their insertion in the genome can disrupt the structure or the regulatory region of a gene, affecting its expression level. These events have been identified as occurring mutations during cancer development, thus being able to identify PΨgs and their location will improve their impact on diagnostic testing, not only in cancer but also in inherited disorders.

RESULTS

We have implemented PΨFinder (P-psy-finder), a tool that identifies PΨgs, annotates known ones and predicts their insertion site(s) in the genome. The tool screens alignment files and provides user-friendly summary reports and visualizations. To demonstrate its applicability, we scanned 218 DNA samples from patients screened for hereditary colorectal cancer. We detected 423 PΨgs distributed in 96% of the samples, comprising 7 different parent genes. Among these, we confirmed the well-known insertion site of the SMAD4-PΨg within the last intron of the SCAI gene in one sample. While for the ubiquitous CBX3-PΨg, present in 82.6% of the samples, we found it reversed inserted in the second intron of the C15ORF57 gene.

CONCLUSIONS

PΨFinder is a tool that can automatically identify novel PΨgs from DNA sequencing data and determine their location in the genome with high sensitivity (95.92%). It generates high quality figures and tables that facilitate the interpretation of the results and can guide the experimental validation. PΨFinder is a complementary analysis to any mutational screening in the identification of disease-causing mutations within cancer and other diseases.

Establishment of reference standards for multifaceted mosaic variant analysis

Detection of somatic mosaicism in non-proliferative cells is a new challenge in genome research, however, the accuracy of current detection strategies remains uncertain due to the lack of a ground truth. Herein, we sought to present a set of ultra-deep sequenced WES data based on reference standards generated by cell line mixtures, providing a total of 386,613 mosaic single-nucleotide variants (SNVs) and insertion-deletion mutations (INDELs) with variant allele frequencies (VAFs) ranging from 0.5% to 56%, as well as 35,113,417 non-variant and 19,936 germline variant sites as a negative control. The whole reference standard set mimics the cumulative aspect of mosaic variant acquisition such as in the early developmental stage owing to the progressive mixing of cell lines with established genotypes, ultimately unveiling 741 possible inter-sample relationships with respect to variant sharing and asymmetry in VAFs. We expect that our reference data will be essential for optimizing the current use of mosaic variant detection strategies and for developing algorithms to enable future improvements.

Measurement(s)	genotype
Technology Type(s)	DNA sequencing
Factor Type(s)	genotyping
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Environment	cell line

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

An investigation into the role of inherited CEACAM gene family variants and colorectal cancer risk

Objective

This study was designed to determine if CEACAM mutations are associated with inherited risk of colorectal cancer. Recently, protein-truncating mutations in the CEACAM gene family were associated with inherited breast cancer risk. That discovery, along with aberrant expression of CEACAM genes in colorectal cancer tumors and that colorectal cancer and breast cancer share many risk factors, including genetics, inspired our team to search for inherited CEACAM mutations in colorectal cancer cases. Specifically utilizing The Cancer Genome Atlas (TCGA) blood-derived whole-exome sequencing data from the colorectal cancer cohort, rare protein-truncating variants and missense variants were investigated through single variant and aggregation analyses in European American and African American cases and compared to ethnic-matched controls.

Results

A total of 34 and 14 different CEACAM variants were identified in European American and African American colorectal cancer cases, respectively. Nine missense variants were individually associated with risk, two in African Americans and seven in European Americans. No identified protein-truncating variants were associated with CRC risk in either ethnicity. Gene family and gene-specific aggregation analyses did not yield any significant results.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-05907-6.

Radiation exposure elicits a neutrophil-driven response in healthy lung tissue that enhances metastatic colonization

Radiotherapy is one of the most effective approaches to achieve tumor control in cancer patients, although healthy tissue injury due to off-target radiation exposure can occur. In this study, we used a model of acute radiation injury to the lung, in the context of cancer metastasis, to understand the biological link between tissue damage and cancer progression. We exposed healthy mouse lung tissue to radiation before the induction of metastasis and observed a strong enhancement of cancer cell growth. We found that locally activated neutrophils were key drivers of the tumor-supportive preconditioning of the lung microenvironment, governed by enhanced regenerative Notch signaling. Importantly, these tissue perturbations endowed arriving cancer cells with an augmented stemness phenotype. By preventing neutrophil-dependent Notch activation, via blocking degranulation, we were able to significantly offset the radiation-enhanced metastases. This work highlights a pro-tumorigenic activity of neutrophils, which is likely linked to their tissue regenerative functions.

LINE1 are spliced in non-canonical transcript variants to regulate T cell quiescence and exhaustion

How gene expression is controlled to preserve human T cell quiescence is poorly understood. Here we show that non-canonical splicing variants containing long interspersed nuclear element 1 (LINE1) enforce naive CD4+ T cell quiescence. LINE1-containing transcripts are derived from CD4+ T cell-specific genes upregulated during T cell activation. In naive CD4+ T cells, LINE1-containing transcripts are regulated by the transcription factor IRF4 and kept at chromatin by nucleolin; these transcripts act in cis, hampering levels of histone 3 (H3) lysine 36 trimethyl (H3K36me3) and stalling gene expression. T cell activation induces LINE1-containing transcript downregulation by the splicing suppressor PTBP1 and promotes expression of the corresponding protein-coding genes by the elongating factor GTF2F1 through mTORC1. Dysfunctional T cells, exhausted in vitro or tumor-infiltrating lymphocytes (TILs), accumulate LINE1-containing transcripts at chromatin. Remarkably, depletion of LINE1-containing transcripts restores TIL effector function. Our study identifies a role for LINE1 elements in maintaining T cell quiescence and suggests that an abundance of LINE1-containing transcripts is critical for T cell effector function and exhaustion.

KAS-seq: genome-wide sequencing of single-stranded DNA by N3-kethoxal-assisted labeling

Transcription and its dynamics are crucial for gene expression regulation. However, very few methods can directly read out transcriptional activity with low-input material and high temporal resolution. This protocol describes KAS-seq, a robust and sensitive approach for capturing genome-wide single-stranded DNA (ssDNA) profiles using N₃-kethoxal-assisted labeling. We developed N₃-kethoxal, an azido derivative of kethoxal that reacts with deoxyguanosine bases of ssDNA in live cells within 5-10 min at 37 °C, allowing the capture of dynamic changes. Downstream biotinylation of labeled DNA occurs via copper-free click chemistry. Altogether, the KAS-seq procedure involves N₃-kethoxal labeling, DNA isolation, biotinylation, fragmentation, affinity pull-down, library preparation, sequencing and bioinformatics analysis. The pre-library construction labeling and enrichment can be completed in as little as 3-4 h and is applicable to both animal tissue and as few as 1,000 cultured cells. Our recent study shows that ssDNA signals measured by KAS-seq simultaneously reveal the dynamics of transcriptionally engaged RNA polymerase (Pol) II, transcribing enhancers, RNA Pol I and Pol III activities and potentially non-canonical DNA structures with high analytical sensitivity. In addition to the experimental protocol, we also introduce here KAS-pipe, a user-friendly integrative data analysis pipeline for KAS-seq.

Petabase-scale sequence alignment catalyses viral discovery

Public databases contain a planetary collection of nucleic acid sequences, but their systematic exploration has been inhibited by a lack of efficient methods for searching this corpus, which (at the time of writing) exceeds 20 petabases and is growing exponentially1. Here we developed a cloud computing infrastructure, Serratus, to enable ultra-high-throughput sequence alignment at the petabase scale. We searched 5.7 million biologically diverse samples (10.2 petabases) for the hallmark gene RNA-dependent RNA polymerase and identified well over 105 novel RNA viruses, thereby expanding the number of known species by roughly an order of magnitude. We characterized novel viruses related to coronaviruses, hepatitis delta virus and huge phages, respectively, and analysed their environmental reservoirs. To catalyse the ongoing revolution of viral discovery, we established a free and comprehensive database of these data and tools. Expanding the known sequence diversity of viruses can reveal the evolutionary origins of emerging pathogens and improve pathogen surveillance for the anticipation and mitigation of future pandemics.

Immunohistochemical staining patterns of p53 predict the mutational status of TP53 in oral epithelial dysplasia

Next-generation sequencing of oral squamous cell carcinoma (OSCC) has revealed TP53 as the most frequently mutated gene in OSCC mutually exclusive with human papillomavirus infection. Oral epithelial dysplasia (OED) is defined as a precancerous lesion of OSCC by the current World Health Organization (WHO) classification; therefore, it is assumed that TP53 mutations occur in early precancerous conditions such as OED. Here, we conducted an integrated analysis of TP53, including whole coding sequencing of TP53, FISH analysis of the 17p13.1 locus, and immunohistochemical analysis for p53 (p53-IHC), in 40 OED cases. We detected 20 mutations in 16 (40%) OED cases, and four cases, each harbored two mutations. FISH analysis revealed six of 24 cases (25%) had a deletion on 17p13.1, and four cases had concurrent TP53 mutations and 17p13.1 deletion (2-hit). Also, the increased frequency of TP53 mutations in higher degrees of OED implies acquisition of the mutation is a major event toward OSCC. p53-IHC revealed that overall cases could be categorized into four patterns that correlate well with the mutational status of TP53. Especially, two patterns, broad p53 expression type (pattern HI) and p53 null type (pattern LS), strongly correlated with a missense mutation and nonsense mutation, respectively. Furthermore, seven of the 40 cases progressed to SCC, and six of these seven cases presented pattern HI or LS. Therefore, patterns HI and LS have a high risk for malignant transformation if excisional treatment is not performed irrespective of the dysplasia grade. Although the current WHO classification mainly focuses on morphological criteria for the diagnosis of OED, interobserver discrepancy appears in some instances of the OED diagnosis. Our immunohistochemical analysis supports a more accurate pathological diagnosis for OED in cases of low dysplastic changes or of differential diagnosis with non-dysplastic lesions.

The transcriptional regulator HDP1 controls expansion of the inner membrane complex during early sexual differentiation of malaria parasites

Transmission of Plasmodium falciparum and other malaria parasites requires their differentiation from asexual blood stages into gametocytes, the non-replicative sexual stage necessary to infect the mosquito vector. This transition involves changes in gene expression and chromatin reorganization that result in the activation and silencing of stage-specific genes. However, the genomes of malaria parasites have been noted for their limited number of transcriptional and chromatin regulators, and the molecular mediators of these changes remain largely unknown. We recently identified homeodomain protein 1 (HDP1) as a DNA-binding protein, first expressed in gametocytes, that enhances the expression of key genes critical for early sexual differentiation. The discovery of HDP1 marks a new class of transcriptional regulator in malaria parasites outside of the better-characterized ApiAP2 family. Here, using molecular biology, biochemistry and microscopy techniques, we show that HDP1 is essential for gametocyte maturation, facilitating the necessary upregulation of inner membrane complex components during early gametocytogenesis that gives P. falciparum gametocytes their characteristic shape.

Utilizing museomics to trace the complex history and species boundaries in an avian-study system of conservation concern

A taxonomic classification that accurately captures evolutionary history is essential for conservation. Genomics provides powerful tools for delimiting species and understanding their evolutionary relationships. This allows for a more accurate and detailed view on conservation status compared with other, traditionally used, methods. However, from a practical and ethical perspective, gathering sufficient samples for endangered taxa may be difficult. Here, we use museum specimens to trace the evolutionary history and species boundaries in an Asian oriole clade. The endangered silver oriole has long been recognized as a distinct species based on its unique coloration, but a recent study suggested that it might be nested within the maroon oriole-species complex. To evaluate species designation, population connectivity, and the corresponding conservation implications, we assembled a de novo genome and used whole-genome resequencing of historical specimens. Our results show that the silver orioles form a monophyletic lineage within the maroon oriole complex and that maroon and silver forms continued to interbreed after initial divergence, but do not show signs of recent gene flow. Using a genome scan, we identified genes that may form the basis for color divergence and act as reproductive barriers. Taken together, our results confirm the species status of the silver oriole and highlight that taxonomic revision of the maroon forms is urgently needed. Our study demonstrates how genomics and Natural History Collections (NHC) can be utilized to shed light on the taxonomy and evolutionary history of natural populations and how such insights can directly benefit conservation practitioners when assessing wild populations.

Significance of gene variants for the functional biogeography of the near-surface Atlantic Ocean microbiome

Microbial communities are major drivers of global elemental cycles in the oceans due to their high abundance and enormous taxonomic and functional diversity. Recent studies assessed microbial taxonomic and functional biogeography in global oceans but microbial functional biogeography remains poorly studied. Here we show that in the near-surface Atlantic and Southern Ocean between 62°S and 47°N microbial communities exhibit distinct taxonomic and functional adaptations to regional environmental conditions. Richness and diversity showed maxima around 40° latitude and intermediate temperatures, especially in functional genes (KEGG-orthologues, KOs) and gene profiles. A cluster analysis yielded three clusters of KOs but five clusters of genes differing in the abundance of genes involved in nutrient and energy acquisition. Gene profiles showed much higher distance-decay rates than KO and taxonomic profiles. Biotic factors were identified as highly influential in explaining the observed patterns in the functional profiles, whereas temperature and biogeographic province mainly explained the observed taxonomic patterns. Our results thus indicate fine-tuned genetic adaptions of microbial communities to regional biotic and environmental conditions in the Atlantic and Southern Ocean.

The taxonomic and functional diversity of marine microbial communities are shaped by both environmental and biotic factors. Here, the authors investigate the functional biogeography of epipelagic prokaryotic communities along a 13,000-km transect in the Southern and Atlantic Oceans, showing finely tuned genetic adaptations to regional conditions.

Genome resequencing-based high-density genetic map and QTL detection for yield and fiber quality traits in diploid Asiatic cotton (Gossypium arboreum)

Cotton is the most important fiber crop in the world. Asiatic cotton (Gossypium arboreum, genome A₂) is a diploid cotton species producing spinnable fibers and important germplasm for cotton breeding and a significant model for fiber biology. However, the genetic map of Asiatic cotton has been lagging behind tetraploid cottons, as well as other stable crops. This study aimed to construct a high-density SNP genetic map and to map QTLs for important yield and fiber quality traits. Using a recombinant inbred line (RIL) population and genome resequencing technology, we constructed a high-density genetic map that covered 1980.17 cM with an average distance of 0.61 cM between adjacent markers. QTL analysis revealed a total of 297 QTLs for 13 yield and fiber quality traits in three environments, explaining 5.0-37.4% of the phenotypic variance, among which 75 were stably detected in two or three environments. Besides, 47 QTL clusters, comprising 131 QTLs for representative traits, were identified. Our works laid solid foundation for fine mapping and cloning of QTL for yield and fiber quality traits in Asiatic cotton.

Mapping short tandem repeats for liver gene expression traits helps prioritize potential causal variants for complex traits in pigs

BACKGROUND

Short tandem repeats (STRs) were recently found to have significant impacts on gene expression and diseases in humans, but their roles on gene expression and complex traits in pigs remain unexplored. This study investigates the effects of STRs on gene expression in liver tissues based on the whole-genome sequences and RNA-Seq data of a discovery cohort of 260 F6 individuals and a validation population of 296 F7 individuals from a heterogeneous population generated from crosses among eight pig breeds.

RESULTS

We identified 5203 and 5868 significantly expression STRs (eSTRs, FDR < 1%) in the F6 and F7 populations, respectively, most of which could be reciprocally validated (π1 = 0.92). The eSTRs explained 27.5% of the cis-heritability of gene expression traits on average. We further identified 235 and 298 fine-mapped STRs through the Bayesian fine-mapping approach in the F6 and F7 pigs, respectively, which were significantly enriched in intron, ATAC peak, compartment A and H3K4me3 regions. We identified 20 fine-mapped STRs located in 100 kb windows upstream and downstream of published complex trait-associated SNPs, which colocalized with epigenetic markers such as H3K27ac and ATAC peaks. These included eSTR of the CLPB, PGLS, PSMD6 and DHDH genes, which are linked with genome-wide association study (GWAS) SNPs for blood-related traits, leg conformation, growth-related traits, and meat quality traits, respectively.

CONCLUSIONS

This study provides insights into the effects of STRs on gene expression traits. The identified eSTRs are valuable resources for prioritizing causal STRs for complex traits in pigs.

The breast pre-cancer atlas illustrates the molecular and micro-environmental diversity of ductal carcinoma in situ

Microenvironmental and molecular factors mediating the progression of Breast Ductal Carcinoma In Situ (DCIS) are not well understood, impeding the development of prevention strategies and the safe testing of treatment de-escalation. We addressed methodological barriers and characterized the mutational, transcriptional, histological, and microenvironmental landscape across 85 multiple microdissected regions from 39 cases. Most somatic alterations, including whole-genome duplications, were clonal, but genetic divergence increased with physical distance. Phenotypic and subtype heterogeneity was frequently associated with underlying genetic heterogeneity and regions with low-risk features preceded those with high-risk features according to the inferred phylogeny. B- and T-lymphocytes spatial analysis identified three immune states, including an epithelial excluded state located preferentially at DCIS regions, and characterized by histological and molecular features of immune escape, independently from molecular subtypes. Such breast pre-cancer atlas with uniquely integrated observations will help scope future expansion studies and build finer models of outcomes and progression risk.

Sequencing and Reconstructing Helminth Mitochondrial Genomes Directly from Genomic Next-Generation Sequencing Data

We present a detailed method for extraction of high-molecular weight genomic DNA suitable for numerous DNA sequencing applications, and a straightforward in silico approach for reconstructing novel mitochondrial (mt) genomes directly from total genomic DNA extracts derived from next-generation sequencing (NGS) data sets. The in silico post-sequencing pipeline described is fast, accurate, and highly efficient, with modest memory requirements that can be performed using a standard desktop computer. The approach is particularly effective for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information currently available and overcomes many of the limitations of traditional strategies. The described methodologies are also applicable for metagenomics sequencing from mixed or pooled samples containing multiple species and subsequent specific assembly of specific mitochondrial genomes.

Advanced genomics identifies growth effectors for proteotoxic ER stress recovery in Arabidopsis thaliana

Adverse environmental and pathophysiological situations can overwhelm the biosynthetic capacity of the endoplasmic reticulum (ER), igniting a potentially lethal condition known as ER stress. ER stress hampers growth and triggers a conserved cytoprotective signaling cascade, the unfolded protein response (UPR) for ER homeostasis. As ER stress subsides, growth is resumed. Despite the pivotal role of the UPR in growth restoration, the underlying mechanisms for growth resumption are yet unknown. To discover these, we undertook a genomics approach in the model plant species Arabidopsis thaliana and mined the gene reprogramming roles of the UPR modulators, basic leucine zipper28 (bZIP28) and bZIP60, in ER stress resolution. Through a network modeling and experimental validation, we identified key genes downstream of the UPR bZIP-transcription factors (bZIP-TFs), and demonstrated their functional roles. Our analyses have set up a critical pipeline for functional gene discovery in ER stress resolution with broad applicability across multicellular eukaryotes.

Ko and Brandizzi use Arabidopsis thaliana to investigate the downstream regulators of two major endoplasmic reticulum (ER) stress-related transcription factors, bZIP60 and bZIP28. Their results provide further insight on how two modulators of the unfolded protein response contribute to growth recovery from ER stress.

Mutational landscape of paired primary and synchronous metastatic lymph node in chemotherapy naive gallbladder cancer

BACKGROUND

Comprehensive genomic analysis of paired primary tumors and their metastatic lesions may provide new insights into the biology of metastatic processes and therefore guide the development of novel strategies for intervention. To date, our knowledge of the genetic divergence and phylogenetic relationships in gallbladder cancer (GBC) is limited.

METHODS

We performed whole exome sequencing for 5 patients with primary tumor, metastatic lymph node (LNM) and corresponding normal tissue. Mutations, mutation signatures and copy number variations were analyzed with state-of-art bioinformatics methods. Phylogenetic tree was also generated to infer metastatic pattern.

RESULTS

Five driver mutations were detected in these patients. Among which, TP53 was the only shared mutation between primary tumor and LNM. Although tumor mutational burden was comparable between primary tumor and LNM, higher mutation burden was observed in LNM of one patient. Copy number variations (CNVs) burden was higher in LNM than their primary tumor. Phylogenetic analysis indicated both linear and parallel progression of metastasis exist in these patients. TP53 mutation and CNVs were homogenously between primary tumor and LNM.

CONCLUSIONS

High consistence of genetic landscape were shown between primary tumor and LNM in GBC. However, heterogenicity still exist between primary tumor and LNM in particular patients in term of driver mutation, TMB and CNV burden. Phylogenetic analysis indicated both Linear and parallel progression of metastasis were exist among these patients.

Mendelian randomization analyses support causal relationships between blood metabolites and the gut microbiome

The gut microbiome has been implicated in a variety of physiological states, but controversy over causality remains unresolved. Here, we performed bidirectional Mendelian randomization analyses on 3,432 Chinese individuals with whole-genome, whole-metagenome, anthropometric and blood metabolic trait data. We identified 58 causal relationships between the gut microbiome and blood metabolites, and replicated 43 of them. Increased relative abundances of fecal Oscillibacter and Alistipes were causally linked to decreased triglyceride concentration. Conversely, blood metabolites such as glutamic acid appeared to decrease fecal Oxalobacter, and members of Proteobacteria were influenced by metabolites such as 5-methyltetrahydrofolic acid, alanine, glutamate and selenium. Two-sample Mendelian randomization with data from Biobank Japan partly corroborated results with triglyceride and with uric acid, and also provided causal support for published fecal bacterial markers for cancer and cardiovascular diseases. This study illustrates the value of human genetic information to help prioritize gut microbial features for mechanistic and clinical studies.

Two divergent haplotypes from a highly heterozygous lychee genome suggest independent domestication events for early and late-maturing cultivars

Lychee is an exotic tropical fruit with a distinct flavor. The genome of cultivar ‘Feizixiao’ was assembled into 15 pseudochromosomes, totaling ~470 Mb. High heterozygosity (2.27%) resulted in two complete haplotypic assemblies. A total of 13,517 allelic genes (42.4%) were differentially expressed in diverse tissues. Analyses of 72 resequenced lychee accessions revealed two independent domestication events. The extremely early maturing cultivars preferentially aligned to one haplotype were domesticated from a wild population in Yunnan, whereas the late-maturing cultivars that mapped mostly to the second haplotype were domesticated independently from a wild population in Hainan. Early maturing cultivars were probably developed in Guangdong via hybridization between extremely early maturing cultivar and late-maturing cultivar individuals. Variable deletions of a 3.7 kb region encompassed by a pair of CONSTANS-like genes probably regulate fruit maturation differences among lychee cultivars. These genomic resources provide insights into the natural history of lychee domestication and will accelerate the improvement of lychee and related crops.

Two divergent haplotypes from a highly heterozygous lychee genome of the cultivar ‘Feizixiao’ and resequencing of 72 lychee accessions provide insights into the genome evolution and domestication history of lychee.

Generating lineage-resolved, complete metagenome-assembled genomes from complex microbial communities

Microbial communities might include distinct lineages of closely related organisms that complicate metagenomic assembly and prevent the generation of complete metagenome-assembled genomes (MAGs). Here we show that deep sequencing using long (HiFi) reads combined with Hi-C binning can address this challenge even for complex microbial communities. Using existing methods, we sequenced the sheep fecal metagenome and identified 428 MAGs with more than 90% completeness, including 44 MAGs in single circular contigs. To resolve closely related strains (lineages), we developed MAGPhase, which separates lineages of related organisms by discriminating variant haplotypes across hundreds of kilobases of genomic sequence. MAGPhase identified 220 lineage-resolved MAGs in our dataset. The ability to resolve closely related microbes in complex microbial communities improves the identification of biosynthetic gene clusters and the precision of assigning mobile genetic elements to host genomes. We identified 1,400 complete and 350 partial biosynthetic gene clusters, most of which are novel, as well as 424 (298) potential host-viral (host-plasmid) associations using Hi-C data.

Decoding gene regulation in the fly brain

The Drosophila brain is a frequently used model in neuroscience. Single-cell transcriptome analysis^1-6, three-dimensional morphological classification⁷ and electron microscopy mapping of the connectome^8,9 have revealed an immense diversity of neuronal and glial cell types that underlie an array of functional and behavioural traits in the fly. The identities of these cell types are controlled by gene regulatory networks (GRNs), involving combinations of transcription factors that bind to genomic enhancers to regulate their target genes. Here, to characterize GRNs at the cell-type level in the fly brain, we profiled the chromatin accessibility of 240,919 single cells spanning 9 developmental timepoints and integrated these data with single-cell transcriptomes. We identify more than 95,000 regulatory regions that are used in different neuronal cell types, of which 70,000 are linked to developmental trajectories involving neurogenesis, reprogramming and maturation. For 40 cell types, uniquely accessible regions were associated with their expressed transcription factors and downstream target genes through a combination of motif discovery, network inference and deep learning, creating enhancer GRNs. The enhancer architectures revealed by DeepFlyBrain lead to a better understanding of neuronal regulatory diversity and can be used to design genetic driver lines for cell types at specific timepoints, facilitating their characterization and manipulation.

A catalogue of 1,167 genomes from the human gut archaeome

The human gut microbiome plays an important role in health, but its archaeal diversity remains largely unexplored. In the present study, we report the analysis of 1,167 nonredundant archaeal genomes (608 high-quality genomes) recovered from human gastrointestinal tract, sampled across 24 countries and rural and urban populations. We identified previously undescribed taxa including 3 genera, 15 species and 52 strains. Based on distinct genomic features, we justify the split of the Methanobrevibacter smithii clade into two separate species, with one represented by the previously undescribed 'Candidatus Methanobrevibacter intestini'. Patterns derived from 28,581 protein clusters showed significant associations with sociodemographic characteristics such as age groups and lifestyle. We additionally show that archaea are characterized by specific genomic and functional adaptations to the host and carry a complex virome. Our work expands our current understanding of the human archaeome and provides a large genome catalogue for future analyses to decipher its impact on human physiology.

Association of rs3027178 polymorphism in the circadian clock gene PER1 with susceptibility to Alzheimer's disease and longevity in an Italian population

Many physiological processes in the human body follow a 24-h circadian rhythm controlled by the circadian clock system. Light, sensed by retina, is the predominant “zeitgeber” able to synchronize the circadian rhythms to the light-dark cycles. Circadian rhythm dysfunction and sleep disorders have been associated with aging and neurodegenerative diseases including mild cognitive impairment (MCI) and Alzheimer’s disease (AD). In the present study, we aimed at investigating the genetic variability of clock genes in AD patients compared to healthy controls from Italy. We also included a group of Italian centenarians, considered as super-controls in association studies given their extreme phenotype of successful aging. We analyzed the exon sequences of eighty-four genes related to circadian rhythms, and the most significant variants identified in this first discovery phase were further assessed in a larger independent cohort of AD patients by matrix assisted laser desorption/ionization-time of flight mass spectrometry. The results identified a significant association between the rs3027178 polymorphism in the PER1 circadian gene with AD, the G allele being protective for AD. Interestingly, rs3027178 showed similar genotypic frequencies among AD patients and centenarians. These results collectively underline the relevance of circadian dysfunction in the predisposition to AD and contribute to the discussion on the role of the relationship between the genetics of age-related diseases and of longevity.

Epigenetic models developed for plains zebras predict age in domestic horses and endangered equids

Effective conservation and management of threatened wildlife populations require an accurate assessment of age structure to estimate demographic trends and population viability. Epigenetic aging models are promising developments because they estimate individual age with high accuracy, accurately predict age in related species, and do not require invasive sampling or intensive long-term studies. Using blood and biopsy samples from known age plains zebras (Equus quagga), we model epigenetic aging using two approaches: the epigenetic clock (EC) and the epigenetic pacemaker (EPM). The plains zebra EC has the potential for broad application within the genus Equus given that five of the seven extant wild species of the genus are threatened. We test the EC's ability to predict age in sister taxa, including two endangered species and the more distantly related domestic horse, demonstrating high accuracy in all cases. By comparing chronological and estimated age in plains zebras, we investigate age acceleration as a proxy of health status. An interaction between chronological age and inbreeding is associated with age acceleration estimated by the EPM, suggesting a cumulative effect of inbreeding on biological aging throughout life.

Temporal transitions in the post-mitotic nervous system of Caenorhabditis elegans

In most animals, the majority of the nervous system is generated and assembled into neuronal circuits during embryonic development¹. However, during juvenile stages, nervous systems still undergo extensive anatomical and functional changes to eventually form a fully mature nervous system by the adult stage^2,3. The molecular changes in post-mitotic neurons across post-embryonic development and the genetic programs that control these temporal transitions are not well understood^4,5. Here, using the model system Caenorhabditis elegans, we comprehensively characterized the distinct functional states (locomotor behaviour) and the corresponding distinct molecular states (transcriptome) of the post-mitotic nervous system across temporal transitions during post-embryonic development. We observed pervasive, neuron-type-specific changes in gene expression, many of which are controlled by the developmental upregulation of the conserved heterochronic microRNA LIN-4 and the subsequent promotion of a mature neuronal transcriptional program through the repression of its target, the transcription factor lin-14. The functional relevance of these molecular transitions are exemplified by a temporally regulated target gene of the LIN-14 transcription factor, nlp-45, a neuropeptide-encoding gene, which we find is required for several distinct temporal transitions in exploratory activity during post-embryonic development. Our study provides insights into regulatory strategies that control neuron-type-specific gene batteries to modulate distinct behavioural states across temporal, sexual and environmental dimensions of post-embryonic development.