Non-invasive detection of lymphoma with circulating tumor DNA features and protein tumor markers.

Background: According to GLOBOCAN 2020, lymphoma ranked as the 9th most common cancer and the 12th leading cause of cancer-related deaths worldwide. Traditional diagnostic methods rely on the invasive excisional lymph node biopsy, which is an invasive approach with some limitations. Most lymphoma patients are diagnosed at an advanced stage since they are asymptomatic at the beginning, which has significantly impacted treatment efficacy and prognosis of the disease. Method: This study assessed the performance and utility of a newly developed blood-based assay (SeekInCare) for lymphoma early detection. SeekInCare utilized protein tumor markers and a comprehensive set of cancer-associated genomic features, including copy number aberration (CNA), fragment size (FS), end motif, and lymphoma-related virus, which were profiled by shallow WGS of cfDNA. Results: Protein marker CA125 could be used for lymphoma detection independent of gender, and the sensitivity was 27.8% at specificity of 98.0%. After integrating these multi-dimensional features, 77.8% sensitivity was achieved at specificity of 98.0%, while its NPV and PPV were both more than 92% for lymphoma detection. The sensitivity of early-stage (I-II) lymphoma was up to 51.3% (47.4% and 55.0% for stage I and II respectively). After 2 cycles of treatment, the molecular response of SeekInCare was correlated with the clinical outcome. Conclusion: In summary, a blood-based assay can be an alternative to detect lymphoma with adequate performance. This approach becomes particularly valuable in cases where obtaining tissue biopsy is difficult to obtain or inconclusive.

Circulating cell-free DNA fragmentation is a stepwise and conserved process linked to apoptosis.

BACKGROUND: Circulating cell-free DNA (cfDNA) is a pool of short DNA fragments mainly released from apoptotic hematopoietic cells. Nevertheless, the precise physiological process governing the DNA fragmentation and molecular profile of cfDNA remains obscure. To dissect the DNA fragmentation process, we use a human leukemia cell line HL60 undergoing apoptosis to analyze the size distribution of DNA fragments by shallow whole-genome sequencing (sWGS). Meanwhile, we also scrutinize the size profile of plasma cfDNA in 901 healthy human subjects and 38 dogs, as well as 438 patients with six common cancer types by sWGS. RESULTS: Distinct size distribution profiles were observed in the HL60 cell pellet and supernatant, suggesting fragmentation is a stepwise process. Meanwhile, C-end preference was seen in both intracellular and extracellular cfDNA fragments. Moreover, the cfDNA profiles are characteristic and conserved across mammals. Compared with healthy subjects, distinct cfDNA profiles with a higher proportion of short fragments and lower C-end preference were found in cancer patients. CONCLUSIONS: Our study provides new insight into fragmentomics of circulating cfDNA processing, which will be useful for early diagnosis of cancer and surveillance during cancer progression.

Whole-genome sequencing as an alternative to analyze copy number abnormalities in acute myeloid leukemia and myelodysplastic syndrome.

Copy number aberrations (CNA) are the core determinants for diagnosis, risk stratification and prognosis in acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS). In this study, a shallow whole-genome sequencing-based assay, LeukoPrint, was utilized to depict genomic CNA profiles from the bone marrow of 137 newly diagnosed AML/MDS patients. It demonstrated 98.1% concordance of CNA profiles with cytogenetics and/or fluorescence in situ hybridization (FISH). It is advantageous in detecting CNAs of short segments (1 Mb) and from samples with low leukemic cell content, more accurate for describing complex karyotypes and less confounded by subjective bias. LeukoPrint improved the overall diagnostic yield by redefining the risk categories for 16 patients by presenting new information. In summary, LeukoPrint provided an automated, convenient, and cost-effective approach to describe genomic CNA profiles. It brought greater diagnostic yield and risk stratification information by incorporating into the routine cytogenetics based on the CNA-related criteria of standard ELN/IPSS-R guidelines.

Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent 1-year global expansion.

A novel coronavirus, SARS-CoV-2, has caused over 274 million cases and over 5.3 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of the cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to the early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e., M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic 1 year later. By analyzing 4013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of the initial phase of global ignition over a time span of 14 weeks. Eleven major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted of half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implies that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261 323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e., including all the publicly available viral genomes) to recapitulate our findings over 1-year time span. By December 25, 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the five variants of concern (VOC) are the descendants of WE1 type. This study demonstrates that viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time.

Comprehensive Analysis of Respiratory Burst Oxidase Homologs (Rboh) Gene Family and Function of GbRboh5/18 on Verticillium Wilt Resistance in Gossypium barbadense.

Respiratory burst oxidase homologs (Rbohs) play a predominant role in reactive oxygen species (ROS) production, which is crucial in plant growth, differentiation, as well as their responses to biotic and abiotic stresses. To date, however, there is little knowledge about the function of cotton Rboh genes. Here, we identified a total of 87 Rbohs from five sequenced Gossypium species (the diploids Gossypium arboreum, Gossypium raimondii, and Gossypium australe, and the allotetraploids Gossypium hirsutum and Gossypium barbadense) via BLAST searching their genomes. Phylogenetic analysis of the putative 87 cotton Rbohs revealed that they were divided into seven clades. All members within the same clade are generally similar to each other in terms of gene structure and conserved domain arrangement. In G. barbadense, the expression levels of GbRbohs in the CladeD were induced in response to a fungal pathogen and to hormones (i.e., jasmonic acid and abscisic acid), based upon which the main functional member in CladeD was discerned to be GbRboh5/18. Further functional and physiological analyses showed that the knock-down of GbRboh5/18 expression attenuates plant resistance to Verticillium dahliae infection. Combined with the molecular and biochemical analyses, we found less ROS accumulation in GbRboh5/18-VIGS plants than in control plants after V. dahliae infection. Overexpression of GbRboh5/18 in G. barbadense resulted in more ROS accumulation than in control plants. These results suggest that GbRboh5/18 enhances the cotton plants' resistance against V. dahliae by elevating the levels of ROS accumulation. By integrating phylogenetic, molecular, and biochemical approaches, this comprehensive study provides a detailed overview of the number, phylogeny, and evolution of the Rboh gene family from five sequenced Gossypium species, as well as elucidating the function of GbRboh5/18 for plant resistance against V. dahliae. This study sheds fresh light on the molecular evolutionary properties and function of Rboh genes in cotton, and provides a reference for improving cotton's responses to the pathogen V. dahliae.

Genome sequencing of the Australian wild diploid species Gossypium australe highlights disease resistance and delayed gland morphogenesis.

The diploid wild cotton species Gossypium australe possesses excellent traits including resistance to disease and delayed gland morphogenesis, and has been successfully used for distant breeding programmes to incorporate disease resistance traits into domesticated cotton. Here, we sequenced the G. australe genome by integrating PacBio, Illumina short read, BioNano (DLS) and Hi-C technologies, and acquired a high-quality reference genome with a contig N50 of 1.83 Mb and a scaffold N50 of 143.60 Mb. We found that 73.5% of the G. australe genome is composed of various repeat sequences, differing from those of G. arboreum (85.39%), G. hirsutum (69.86%) and G. barbadense (69.83%). The G. australe genome showed closer collinear relationships with the genome of G. arboreum than G. raimondii and has undergone less extensive genome reorganization than the G. arboreum genome. Selection signature and transcriptomics analyses implicated multiple genes in disease resistance responses, including GauCCD7 and GauCBP1, and experiments revealed induction of both genes by Verticillium dahliae and by the plant hormones strigolactone (GR24), salicylic acid (SA) and methyl jasmonate (MeJA). Experiments using a Verticillium-resistant domesticated G. barbadense cultivar confirmed that knockdown of the homologues of these genes caused a significant reduction in resistance against Verticillium dahliae. Moreover, knockdown of a newly identified gland-associated gene GauGRAS1 caused a glandless phenotype in partial tissues using G. australe. The G. australe genome represents a valuable resource for cotton research and distant relative breeding as well as for understanding the evolutionary history of crop genomes.

Characterization, Expression, and Functional Analysis of a Novel NAC Gene Associated with Resistance to Verticillium Wilt and Abiotic Stress in Cotton.

Elucidating the mechanism of resistance to biotic and abiotic stress is of great importance in cotton. In this study, a gene containing the NAC domain, designated GbNAC1, was identified from Gossypium barbadense L. Homologous sequence alignment indicated that GbNAC1 belongs to the TERN subgroup. GbNAC1 protein localized to the cell nucleus. GbNAC1 was expressed in roots, stems, and leaves, and was especially highly expressed in vascular bundles. Functional analysis showed that cotton resistance to Verticillium wilt was reduced when the GbNAC1 gene was silenced using the virus-induced gene silencing (VIGS) method. GbNAC1-overexpressing Arabidopsis showed enhanced resistance to Verticillium dahliae compared to wild-type. Thus, GbNAC1 is involved in the positive regulation of resistance to Verticillium wilt. In addition, analysis of GbNAC1-overexpressing Arabidopsis under different stress treatments indicated that it is involved in plant growth, development, and response to various abiotic stresses (ABA, mannitol, and NaCl). This suggests that GbNAC1 plays an important role in resistance to biotic and abiotic stresses in cotton. This study provides a foundation for further study of the function of NAC genes in cotton and other plants.

To Be a Flower or Fruiting Branch: Insights Revealed by mRNA and Small RNA Transcriptomes from Different Cotton Developmental Stages.

The architecture of the cotton plant, including fruit branch formation and flowering pattern, is the most important characteristic that directly influences light exploitation, yield and cost of planting. Nulliplex branch is a useful phenotype to study cotton architecture. We used RNA sequencing to obtain mRNA and miRNA profiles from nulliplex- and normal-branch cotton at three developmental stages. The differentially expressed genes (DEGs) and miRNAs were identified that preferentially/specifically expressed in the pre-squaring stage, which is a key stage controlling the transition from vegetative to reproductive growth. The DEGs identified were primarily enriched in RNA, protein, and signalling categories in Gossypium barbadense and Gossypium hirsutum. Interestingly, during the pre-squaring stage, the DEGs were predominantly enriched in transcription factors in both G. barbadense and G. hirsutum, and these transcription factors were mainly involved in branching and flowering. Related miRNAs were also identified. The results showed that fruit branching in cotton is controlled by molecular pathways similar to those in Arabidopsis and that multiple regulated pathways may affect the development of floral buds. Our study showed that the development of fruit branches is closely related to flowering induction and provides insight into the molecular mechanisms of branch and flower development in cotton.