Genomic Classification and Individualized Prognosis in Multiple Myeloma

PURPOSE

Outcomes for patients with newly diagnosed multiple myeloma (NDMM) are heterogenous, with overall survival (OS) ranging from months to over 10 years.

METHODS

To decipher and predict the molecular and clinical heterogeneity of NDMM, we assembled a series of 1,933 patients with available clinical, genomic, and therapeutic data.

RESULTS

Leveraging a comprehensive catalog of genomic drivers, we identified 12 groups, expanding on previous gene expression-based molecular classifications. To build a model predicting individualized risk in NDMM (IRMMa), we integrated clinical, genomic, and treatment variables. To correct for time-dependent variables, including high-dose melphalan followed by autologous stem-cell transplantation (HDM-ASCT), and maintenance therapy, a multi-state model was designed. The IRMMa model accuracy was significantly higher than all comparator prognostic models, with a c-index for OS of 0.726, compared with International Staging System (ISS; 0.61), revised-ISS (0.572), and R2-ISS (0.625). Integral to model accuracy was 20 genomic features, including 1q21 gain/amp, del 1p, TP53 loss, NSD2 translocations, APOBEC mutational signatures, and copy-number signatures (reflecting the complex structural variant chromothripsis). IRMMa accuracy and superiority compared with other prognostic models were validated on 256 patients enrolled in the GMMG-HD6 (ClinicalTrials.gov identifier: NCT02495922) clinical trial. Individualized patient risks were significantly affected across the 12 genomic groups by different treatment strategies (ie, treatment variance), which was used to identify patients for whom HDM-ASCT is particularly effective versus patients for whom the impact is limited.

CONCLUSION

Integrating clinical, demographic, genomic, and therapeutic data, to our knowledge, we have developed the first individualized risk-prediction model enabling personally tailored therapeutic decisions for patients with NDMM.

Genomic and Immunophenotypic Landscape of Acquired Resistance to PD-(L)1 Blockade in Non-Small-Cell Lung Cancer

PURPOSE

Although immune checkpoint inhibitors (ICI) have extended survival in patients with non-small-cell lung cancer (NSCLC), acquired resistance (AR) to ICI frequently develops after an initial benefit. However, the mechanisms of AR to ICI in NSCLC are largely unknown.

METHODS

Comprehensive tumor genomic profiling, machine learning-based assessment of tumor-infiltrating lymphocytes, multiplexed immunofluorescence, and/or HLA-I immunohistochemistry (IHC) were performed on matched pre- and post-ICI tumor biopsies from patients with NSCLC treated with ICI at the Dana-Farber Cancer Institute who developed AR to ICI. Two additional cohorts of patients with intervening chemotherapy or targeted therapies between biopsies were included as controls.

RESULTS

We performed comprehensive genomic profiling and immunophenotypic characterization on samples from 82 patients with NSCLC and matched pre- and post-ICI biopsies and compared findings with a control cohort of patients with non-ICI intervening therapies between biopsies (chemotherapy, N = 32; targeted therapies, N = 89; both, N = 17). Putative resistance mutations were identified in 27.8% of immunotherapy-treated cases and included acquired loss-of-function mutations in STK11, B2M, APC, MTOR, KEAP1, and JAK1/2; these acquired alterations were not observed in the control groups. Immunophenotyping of matched pre- and post-ICI samples demonstrated significant decreases in intratumoral lymphocytes, CD3e+ and CD8a+ T cells, and PD-L1-PD1 engagement, as well as increased distance between tumor cells and CD8+PD-1+ T cells. There was a significant decrease in HLA class I expression in the immunotherapy cohort at the time of AR compared with the chemotherapy (P = .005) and the targeted therapy (P = .01) cohorts.

CONCLUSION

These findings highlight the genomic and immunophenotypic heterogeneity of ICI resistance in NSCLC, which will need to be considered when developing novel therapeutic strategies aimed at overcoming resistance.

Multiomics data integration, limitations, and prospects to reveal the metabolic activity of the coral holobiont

Since their radiation in the Middle Triassic period ∼240 million years ago, stony corals have survived past climate fluctuations and five mass extinctions. Their long-term survival underscores the inherent resilience of corals, particularly when considering the nutrient-poor marine environments in which they have thrived. However, coral bleaching has emerged as a global threat to coral survival, requiring rapid advancements in coral research to understand holobiont stress responses and allow for interventions before extensive bleaching occurs. This review encompasses the potential, as well as the limits, of multiomics data applications when applied to the coral holobiont. Synopses for how different omics tools have been applied to date and their current restrictions are discussed, in addition to ways these restrictions may be overcome, such as recruiting new technology to studies, utilizing novel bioinformatics approaches, and generally integrating omics data. Lastly, this review presents considerations for the design of holobiont multiomics studies to support lab-to-field advancements of coral stress marker monitoring systems. Although much of the bleaching mechanism has eluded investigation to date, multiomic studies have already produced key findings regarding the holobiont's stress response, and have the potential to advance the field further.

[Clinical and genetic analysis of two pedigrees affected with Carnitine-acylcarnitine translocase deficiency due to variant of SLC25A20 gene]

OBJECTIVE

To analyze the clinical phenotype and genotypes of two children with Carnitine-acylcarnitine translocase deficiency (CACTD).

METHODS

Two children diagnosed with CACTD at the Gansu Provincial Maternal and Child Health Care Hospital respectively on January 3 and November 19, 2018 were selected as the study subjects. Trio-whole exome sequencing (trio-WES) was carried out, and candidate variants were validated through Sanger sequencing and pathogenicity analysis.

RESULTS

Both children were males and had manifested mainly with hypoglycemia. Trio-WES and Sanger sequencing showed that child 1 had harbored compound heterozygous variants of the SLC25A20 gene, namely c.49G>C (p.Gly17Arg) and c.106-2A>G, which were inherited from his father and mother, respectively. Child 2 had harbored homozygous c.199-10T>G variants of the SLC25A20 gene, which were inherited from both of his parents. Among these, the c.106-2A>G and c.49G>C variants were unreported previously. Based on the guidelines from the American College of Medical Genetics and Genomics (ACMG), the c.49G>C (p.Gly17Arg), c.106-2A>G, and c.199-10T>G variants were classified as likely pathogenic (PM2_supporting+PP3+PM3_strong+PP4), pathogenic (PVS1+PM2_supporting+PM5+PP3), and pathogenic (PVS1+PM2_supporting+PP3+PP5), respectively.

CONCLUSION

Combined with their clinical phenotype and genetic analysis, both children were diagnosed with CACTD. Above finding has provided a basis for their treatment as well as genetic counseling and prenatal diagnosis for their families.

The race-based stress reduction intervention (RiSE) study on African American women in NYC and Chicago: Design and methods for complex genomic analysis

RiSE study aims to evaluate a race-based stress-reduction intervention as an effective strategy to improve coping and decrease stress-related symptoms, inflammatory burden, and modify DNA methylation of stress response-related genes in older AA women. This article will describe genomic analytic methods to be utilized in this longitudinal, randomized clinical trial of older adult AA women in Chicago and NYC that examines the effect of the RiSE intervention on DNAm pre- and post-intervention, and its overall influence on inflammatory burden. Salivary DNAm will be measured at baseline and 6 months following the intervention, using the Oragene-DNA kit. Measures of perceived stress, depressive symptoms, fatigue, sleep, inflammatory burden, and coping strategies will be assessed at 4 time points including at baseline, 4 weeks, 8 weeks, and 6 months. Genomic data analysis will include the use of pre-processed and quality-controlled methylation data expressed as beta (β) values. Association analyses will be performed to detect differentially methylated sites on the targeted candidate genes between the intervention and non-intervention groups using the Δβ (changes in methylation) with adjustment for age, health behaviors, early life adversity, hybridization batch, and top principal components of the probes as covariates. To account for multiple testing, we will use FDR adjustment with a corrected p-value of <0.05 regarded as statistically significant. To assess the relationship between inflammatory burden and Δβ among the study samples, we will repeat association analyses with the inclusion of individual inflammation protein measures. ANCOVA will be used because it is more statistically powerful to detect differences.

Strengthening Bordetella pertussis genomic surveillance by direct sequencing of residual positive specimens

Whole-genome sequencing (WGS) of microbial pathogens recovered from patients with infectious disease facilitates high-resolution strain characterization and molecular epidemiology. However, increasing reliance on culture-independent methods to diagnose infectious diseases has resulted in few isolates available for WGS. Here, we report a novel culture-independent approach to genome characterization of Bordetella pertussis, the causative agent of pertussis and a paradigm for insufficient genomic surveillance due to limited culture of clinical isolates. Sequencing libraries constructed directly from residual pertussis-positive diagnostic nasopharyngeal specimens were hybridized with biotinylated RNA "baits" targeting B. pertussis fragments within complex mixtures that contained high concentrations of host and microbial background DNA. Recovery of B. pertussis genome sequence data was evaluated with mock and pooled negative clinical specimens spiked with reducing concentrations of either purified DNA or inactivated cells. Targeted enrichment increased the yield of B. pertussis sequencing reads up to 90% while simultaneously decreasing host reads to less than 10%. Filtered sequencing reads provided sufficient genome coverage to perform characterization via whole-genome single nucleotide polymorphisms and whole-genome multilocus sequencing typing. Moreover, these data were concordant with sequenced isolates recovered from the same specimens such that phylogenetic reconstructions from either consistently clustered the same putatively linked cases. The optimized protocol is suitable for nasopharyngeal specimens with diagnostic IS481 Ct < 35 and >10 ng DNA. Routine implementation of these methods could strengthen surveillance and study of pertussis resurgence by capturing additional cases with genomic characterization.

The genomic evolutionary dynamics and global circulation patterns of respiratory syncytial virus

Respiratory syncytial virus (RSV) is a leading cause of acute lower respiratory tract infection in young children and the second leading cause of infant death worldwide. While global circulation has been extensively studied for respiratory viruses such as seasonal influenza, and more recently also in great detail for SARS-CoV-2, a lack of global multi-annual sampling of complete RSV genomes limits our understanding of RSV molecular epidemiology. Here, we capitalise on the genomic surveillance by the INFORM-RSV study and apply phylodynamic approaches to uncover how selection and neutral epidemiological processes shape RSV diversity. Using complete viral genome sequences, we show similar patterns of site-specific diversifying selection among RSVA and RSVB and recover the imprint of non-neutral epidemic processes on their genealogies. Using a phylogeographic approach, we provide evidence for air travel governing the global patterns of RSVA and RSVB spread, which results in a considerable degree of phylogenetic mixing across countries. Our findings highlight the potential of systematic global RSV genomic surveillance for transforming our understanding of global RSV spread.

Clinical manifestations, antimicrobial resistance and genomic feature analysis of multidrug-resistant Elizabethkingia strains

BACKGROUND

Elizabethkingia is emerging as an opportunistic pathogen in humans. The aim of this study was to investigate the clinical epidemiology, antimicrobial susceptibility, virulence factors, and genome features of Elizabethkingia spp.

METHODS

Clinical data from 71 patients who were diagnosed with Elizabethkingia-induced pneumonia and bacteremia between August 2019 and September 2021 were analyzed. Whole-genome sequencing was performed on seven isolates, and the results were compared with a dataset of 83 available Elizabethkingia genomes. Genomic features, Kyoto Encyclopedia of Genes and Genomes (KEGG) results and clusters of orthologous groups (COGs) were analyzed.

RESULTS

The mean age of the patients was 56.9 ± 20.7 years, and the in-hospital mortality rate was 29.6% (21/71). Elizabethkingia strains were obtained mainly from intensive care units (36.6%, 26/71) and emergency departments (32.4%, 23/71). The majority of the strains were isolated from respiratory tract specimens (85.9%, 61/71). All patients had a history of broad-spectrum antimicrobial exposure. Hospitalization for invasive mechanical ventilation or catheter insertion was found to be a risk factor for infection. The isolates displayed a high rate of resistance to cephalosporins and carbapenems, but all were susceptible to minocycline and colistin. Genomic analysis identified five β-lactamase genes (blaGOB, blaBlaB, blaCME, blaOXA, and blaTEM) responsible for β-lactam resistance and virulence genes involved in stress adaptation (ureB/G, katA/B, and clpP), adherence (groEL, tufA, and htpB) and immune modulation (gmd, tviB, cps4J, wbtIL, cap8E/D/G, and rfbC). Functional analysis of the COGs revealed that "metabolism" constituted the largest category within the core genome, while "information storage and processing" was predominant in both the accessory and unique genomes. The unique genes in our 7 strains were mostly enriched in KEGG pathways related to microRNAs in cancer, drug resistance (β-lactam and vancomycin), ABC transporters, biological metabolism and biosynthesis, and nucleotide excision repair mechanisms.

CONCLUSION

The Elizabethkingia genus exhibits multidrug resistance and carries carbapenemase genes. This study presents a comparative genomic analysis of Elizabethkingia, providing knowledge that facilitates a better understanding of this microorganism.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Genetic linkage mapping and quantitative trait locus (QTL) analysis of sweet basil (Ocimum basilicum L.) to identify genomic regions associated with cold tolerance and major volatiles

Chilling sensitivity is one of the greatest challenges affecting the marketability and profitability of sweet basil (Ocimum basilicum L.) in the US and worldwide. Currently, there are no sweet basils commercially available with significant chilling tolerance and traditional aroma profiles. This study was conducted to identify quantitative trait loci (QTLs) responsible for chilling tolerance and aroma compounds in a biparental mapping population, including the Rutgers advanced breeding line that served as a chilling tolerant parent, 'CB15', the chilling sensitive parent, 'Rutgers Obsession DMR' and 200 F2 individuals. Chilling tolerance was assessed by percent necrosis using machine learning and aroma profiling was evaluated using gas chromatography (GC) mass spectrometry (MS). Single nucleotide polymorphism (SNP) markers were generated from genomic sequences derived from double digestion restriction-site associated DNA sequencing (ddRADseq) and converted to genotype data using a reference genome alignment. A genetic linkage map was constructed and five statistically significant QTLs were identified in response to chilling temperatures with possible interactions between QTLs. The QTL on LG24 (qCH24) demonstrated the largest effect for chilling response and was significant in all three replicates. No QTLs were identified for linalool, as the population did not segregate sufficiently to detect this trait. Two significant QTLs were identified for estragole (also known as methyl chavicol) with only qEST1 on LG1 being significant in the multiple-QTL model (MQM). QEUC26 was identified as a significant QTL for eucalyptol (also known as 1,8-cineole) on LG26. These QTLs may represent key mechanisms for chilling tolerance and aroma in basil, providing critical knowledge for future investigation of these phenotypic traits and molecular breeding.

Habitat escalated adaptive therapy (HEAT): a phase 2 trial utilizing radiomic habitat-directed and genomic-adjusted radiation dose (GARD) optimization for high-grade soft tissue sarcoma

BACKGROUND

Soft tissue sarcomas (STS), have significant inter- and intra-tumoral heterogeneity, with poor response to standard neoadjuvant radiotherapy (RT). Achieving a favorable pathologic response (FPR ≥ 95%) from RT is associated with improved patient outcome. Genomic adjusted radiation dose (GARD), a radiation-specific metric that quantifies the expected RT treatment effect as a function of tumor dose and genomics, proposed that STS is significantly underdosed. STS have significant radiomic heterogeneity, where radiomic habitats can delineate regions of intra-tumoral hypoxia and radioresistance. We designed a novel clinical trial, Habitat Escalated Adaptive Therapy (HEAT), utilizing radiomic habitats to identify areas of radioresistance within the tumor and targeting them with GARD-optimized doses, to improve FPR in high-grade STS.

METHODS

Phase 2 non-randomized single-arm clinical trial includes non-metastatic, resectable high-grade STS patients. Pre-treatment multiparametric MRIs (mpMRI) delineate three distinct intra-tumoral habitats based on apparent diffusion coefficient (ADC) and dynamic contrast enhanced (DCE) sequences. GARD estimates that simultaneous integrated boost (SIB) doses of 70 and 60 Gy in 25 fractions to the highest and intermediate radioresistant habitats, while the remaining volume receives standard 50 Gy, would lead to a > 3 fold FPR increase to 24%. Pre-treatment CT guided biopsies of each habitat along with clip placement will be performed for pathologic evaluation, future genomic studies, and response assessment. An mpMRI taken between weeks two and three of treatment will be used for biological plan adaptation to account for tumor response, in addition to an mpMRI after the completion of radiotherapy in addition to pathologic response, toxicity, radiomic response, disease control, and survival will be evaluated as secondary endpoints. Furthermore, liquid biopsy will be performed with mpMRI for future ancillary studies.

DISCUSSION

This is the first clinical trial to test a novel genomic-based RT dose optimization (GARD) and to utilize radiomic habitats to identify and target radioresistance regions, as a strategy to improve the outcome of RT-treated STS patients. Its success could usher in a new phase in radiation oncology, integrating genomic and radiomic insights into clinical practice and trial designs, and may reveal new radiomic and genomic biomarkers, refining personalized treatment strategies for STS.

TRIAL REGISTRATION

NCT05301283.

TRIAL STATUS

The trial started recruitment on March 17, 2022.

A fast comparative genome browser for diverse bacteria and archaea

Genome sequencing has revealed an incredible diversity of bacteria and archaea, but there are no fast and convenient tools for browsing across these genomes. It is cumbersome to view the prevalence of homologs for a protein of interest, or the gene neighborhoods of those homologs, across the diversity of the prokaryotes. We developed a web-based tool, fast.genomics, that uses two strategies to support fast browsing across the diversity of prokaryotes. First, the database of genomes is split up. The main database contains one representative from each of the 6,377 genera that have a high-quality genome, and additional databases for each taxonomic order contain up to 10 representatives of each species. Second, homologs of proteins of interest are identified quickly by using accelerated searches, usually in a few seconds. Once homologs are identified, fast.genomics can quickly show their prevalence across taxa, view their neighboring genes, or compare the prevalence of two different proteins. Fast.genomics is available at https://fast.genomics.lbl.gov.

OASIS: An interpretable, finite-sample valid alternative to Pearson's X2 for scientific discovery

Contingency tables, data represented as counts matrices, are ubiquitous across quantitative research and data-science applications. Existing statistical tests are insufficient however, as none are simultaneously computationally efficient and statistically valid for a finite number of observations. In this work, motivated by a recent application in reference-free genomic inference [K. Chaung et al., Cell 186, 5440-5456 (2023)], we develop Optimized Adaptive Statistic for Inferring Structure (OASIS), a family of statistical tests for contingency tables. OASIS constructs a test statistic which is linear in the normalized data matrix, providing closed-form P-value bounds through classical concentration inequalities. In the process, OASIS provides a decomposition of the table, lending interpretability to its rejection of the null. We derive the asymptotic distribution of the OASIS test statistic, showing that these finite-sample bounds correctly characterize the test statistic's P-value up to a variance term. Experiments on genomic sequencing data highlight the power and interpretability of OASIS. Using OASIS, we develop a method that can detect SARS-CoV-2 and Mycobacterium tuberculosis strains de novo, which existing approaches cannot achieve. We demonstrate in simulations that OASIS is robust to overdispersion, a common feature in genomic data like single-cell RNA sequencing, where under accepted noise models OASIS provides good control of the false discovery rate, while Pearson's [Formula: see text] consistently rejects the null. Additionally, we show in simulations that OASIS is more powerful than Pearson's [Formula: see text] in certain regimes, including for some important two group alternatives, which we corroborate with approximate power calculations.

All of Us Program Illustrates What Less-Diverse Genomics Studies Might Be Missing

This Medical News story discusses new findings from genomics studies involving diverse cohorts.

A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network

BACKGROUND

The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice.

METHODS

FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription.

RESULTS

FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency.

CONCLUSIONS

Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

Plant genome information facilitates plant functional genomics

MAIN CONCLUSION

In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

Emergence and genomic chion of Proteus mirabilis harboring blaNDM-1 in Korean companion dogs

Proteus mirabilis is a commensal bacterium dwelling in the gastrointestinal (GI) tract of humans and animals. Although New Delhi metallo-β-lactamase 1 (NDM-1) producing P. mirabilis is emerging as a threat, its epidemiology in our society remains largely unknown. LHPm1, the first P. mirabilis isolate harboring NDM-1, was detected from a companion dog that resides with a human owner. The whole-genome study revealed 20 different antimicrobial resistance (AMR) genes against various classes of antimicrobial agents, which corresponded to the MIC results. Genomic regions, including MDR genes, were identified with multiple variations and visualized in a comparative manner. In the whole-genome epidemiological analysis, multiple phylogroups were identified, revealing the genetic relationship of LHPm1 with other P. mirabilis strains carrying various AMR genes. These genetic findings offer comprehensive insights into NDM-1-producing P. mirabilis, underscoring the need for urgent control measures and surveillance programs using a "one health approach".

Molecular subtyping improves breast cancer diagnosis in the Copenhagen Breast Cancer Genomics Study

BACKGROUNDIntrinsic molecular subtypes define distinct biological breast cancers and can be used to further improve diagnosis and risk allocation.METHODSThe Copenhagen Breast Cancer Genomics Study (CBCGS) prospectively included women diagnosed with breast cancer at Rigshospitalet from 2014 to 2021. Eligible patients were females with a primary invasive breast cancer (T1c, if N0M0; otherwise, any T, any N, or any M stage) and no prior malignancy. All patients underwent molecular profiling with the CIT256 and PAM50 molecular profile.RESULTSIn the study period, 2,816 patients were included in the CBCGS. Molecular subtyping showed an increase in nonluminal (molecular-apocrine, luminal C, and Basal-like) as compared with luminal (luminal A, luminal B, and Normal-like) subtypes with increasing stage from I to IV. Across all stages, we found a significant difference in survival among subtypes; 91% of patients with LumA were alive at 5 years compared with 91% for LumB, 84% for LumC, 82% for mApo, and 80% for Basal-like. We identified 442 tumors (16%) that were discordant in subtype between CIT256 and IHC. Discordant subtype proved to be a risk factor of death among patients with IHC luminal breast cancer (hazard ratio [HR], 2.08; 95% CI, 1.51-2.86) in a multivariable Cox regression analysis. Discordance occurred more often among patients with N3, stage IV, or grade III disease.CONCLUSIONOur findings indicate that molecular subtypes are a predominant classification for survival. Assessment is particularly crucial for patients with IHC luminal breast cancer with known high-risk factors, since they are at an increased risk of harboring an aggressive molecular subtype.

The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

Genomic and phenotypic signatures provide insights into the wide adaptation of a global plant invader

Invasive alien species are primary drivers of biodiversity loss and species extinction. Smooth cordgrass (Spartina alterniflora) is one of the most aggressive invasive plants in coastal ecosystems around the world. However, the genomic bases and evolutionary mechanisms underlying its invasion success have remained largely unknown. Here, we assembled a chromosome-level reference genome and performed phenotypic and population genomic analyses between native US and introduced Chinese populations. Our phenotypic comparisons showed that introduced Chinese populations have evolved competitive traits, such as early flowering time and greater plant biomass, during secondary introductions along China's coast. Population genomic and transcriptomic inferences revealed distinct evolutionary trajectories of low- and high-latitude Chinese populations. In particular, genetic mixture among different source populations, together with independent natural selection acting on distinct target genes, may have resulted in high genome dynamics of the introduced Chinese populations. Our study provides novel phenotypic and genomic evidence showing how smooth cordgrass rapidly adapts to variable environmental conditions in its introduced ranges. Moreover, candidate genes related to flowering time, fast growth, and stress tolerance (i.e., salinity and submergence) provide valuable genetic resources for future improvement of cereal crops.

A chromosomal-scale genome assembly of modern cultivated hybrid sugarcane provides insights into origination and evolution

Sugarcane is a vital crop with significant economic and industrial value. However, the cultivated sugarcane's ultra-complex genome still needs to be resolved due to its high ploidy and extensive recombination between the two subgenomes. Here, we generate a chromosomal-scale, haplotype-resolved genome assembly for a hybrid sugarcane cultivar ZZ1. This assembly contains 10.4 Gb genomic sequences and 68,509 annotated genes with defined alleles in two sub-genomes distributed in 99 original and 15 recombined chromosomes. RNA-seq data analysis shows that sugar accumulation-associated gene families have been primarily expanded from the ZZSO subgenome. However, genes responding to pokkah boeng disease susceptibility have been derived dominantly from the ZZSS subgenome. The region harboring the possible smut resistance genes has expanded significantly. Among them, the expansion of WAK and FLS2 families is proposed to have occurred during the breeding of ZZ1. Our findings provide insights into the complex genome of hybrid sugarcane cultivars and pave the way for future genomics and molecular breeding studies in sugarcane.

HEMU: An integrated comparative genomics database and analysis platform for Andropogoneae grasses

This study reports an online database and analysis platform HEMU, which integrates 75 genome assemblies from 20 unique species, large amounts of multi-omics data, and six sophisticated analysis toolkits. HEMU will facilitate comparative genomics analysis within the tribe Andropogoneae.

Mo than meets the eye: genomic insights into molybdoenzyme diversity of Seleniivibrio woodruffii strain S4T

Seleniivibrio woodruffii strain S4T is an obligate anaerobe belonging to the phylum Deferribacterota. It was isolated for its ability to respire selenate and was also found to respire arsenate. The high-quality draft genome of this bacterium is 2.9 Mbp, has a G+C content of 48%, 2762 predicted genes of which 2709 are protein-coding, and 53 RNA genes. An analysis of the genome focusing on the genes encoding for molybdenum-containing enzymes (molybdoenzymes) uncovered a remarkable number of genes encoding for members of the dimethylsulfoxide reductase family of proteins (DMSOR), including putative reductases for selenate and arsenate respiration, as well as genes for nitrogen fixation. Respiratory molybdoenzymes catalyze redox reactions that transfer electrons to a variety of substrates that can act as terminal electron acceptors for energy generation. Seleniivibrio woodruffii strain S4T also has essential genes for molybdate transporters and the biosynthesis of the molybdopterin guanine dinucleotide cofactors characteristic of the active centers of DMSORs. Phylogenetic analysis revealed candidate respiratory DMSORs spanning nine subfamilies encoded within the genome. Our analysis revealed the untapped potential of this interesting microorganism and expanded our knowledge of molybdoenzyme co-occurrence.

Genetic interactions reveal distinct biological and therapeutic implications in breast cancer

Co-occurrence and mutual exclusivity of genomic alterations may reflect the existence of genetic interactions, potentially shaping distinct biological phenotypes and impacting therapeutic response in breast cancer. However, our understanding of them remains limited. Herein, we investigate a large-scale multi-omics cohort (n = 873) and a real-world clinical sequencing cohort (n = 4,405) including several clinical trials with detailed treatment outcomes and perform functional validation in patient-derived organoids, tumor fragments, and in vivo models. Through this comprehensive approach, we construct a network comprising co-alterations and mutually exclusive events and characterize their therapeutic potential and underlying biological basis. Notably, we identify associations between TP53mut-AURKAamp and endocrine therapy resistance, germline BRCA1mut-MYCamp and improved sensitivity to PARP inhibitors, and TP53mut-MYBamp and immunotherapy resistance. Furthermore, we reveal that precision treatment strategies informed by co-alterations hold promise to improve patient outcomes. Our study highlights the significance of genetic interactions in guiding genome-informed treatment decisions beyond single driver alterations.

Failure of artemether-lumefantrine therapy in travellers returning to Belgium with Plasmodium falciparum malaria: an observational case series with genomic analysis

BACKGROUND

Failure of artemisinin-based combination therapy is increasingly reported in patients with Plasmodium falciparum malaria in sub-Saharan Africa. We aimed to describe the clinical and genomic characteristics of recent cases of P. falciparum malaria failing artemether-lumefantrine in Belgium.

METHODS

Travel-related cases of malaria confirmed at the national reference laboratory of the Institute of Tropical Medicine, Antwerp, Belgium, were reviewed. All cases for which attending clinicians reported persistence (beyond Day 3 post-treatment initiation, i.e. early failure) or recrudescence (from Day 7 to 42, i.e. late failure) of P. falciparum parasites despite adequate drug intake were analysed. Both initial and persistent/recurrent samples were submitted to next generation sequencing to investigate resistance-conferring mutations.

RESULTS

From July 2022 to June 2023, eight P. falciparum cases of failure with artemether-lumefantrine therapy were reported (early failure = 1; late failure = 7). All travellers were returning from sub-Saharan Africa, most (6/8) after a trip to visit friends and relatives. PfKelch13 (PF3D7_1343700) mutations associated with resistance to artemisinin were found in two travellers returning from East Africa, including the validated marker R561H in the patient with early failure and the candidate marker A675V in a patient with late failure. Additional mutations were detected that could contribute to decreased susceptibility to artemisinin in another three cases, lumefantrine in six cases and proguanil in all eight participants. Various regimens were used to treat the persistent/recrudescent cases, with favourable outcome.

CONCLUSION

Within a 12-month period, we investigated eight travellers returning from sub-Saharan Africa with P. falciparum malaria and in whom artemether-lumefantrine failure was documented. Mutations conferring resistance to antimalarials were found in all analysed blood samples, especially against lumefantrine and proguanil, but also artemisinin. There is a pressing need for systematic genomic surveillance of resistance to antimalarials in international travellers with P. falciparum malaria, especially those experiencing treatment failure.

Genomic analyses correspond with deep persistence of peoples of Blackfoot Confederacy from glacial times

Mutually beneficial partnerships between genomics researchers and North American Indigenous Nations are rare yet becoming more common. Here, we present one such partnership that provides insight into the peopling of the Americas and furnishes another line of evidence that can be used to further treaty and Indigenous rights. We show that the genomics of sampled individuals from the Blackfoot Confederacy belong to a previously undescribed ancient lineage that diverged from other genomic lineages in the Americas in Late Pleistocene times. Using multiple complementary forms of knowledge, we provide a scenario for Blackfoot population history that fits with oral tradition and provides a plausible model for the evolutionary process of the peopling of the Americas.

A chromosome-level genome assembly of the spider mite Tetranychus piercei McGregor

Despite the rapid advances in sequencing technology, limited genomic resources are currently available for phytophagous spider mites, which include many important agricultural pests. One of these pests is Tetranychus piercei (McGregor), a serious banana pest in East Asia exhibiting remarkable tolerance to high temperature. In this study, we assembled a high-quality genome of T. piercei using a combination of PacBio long reads and Illumina short reads sequencing. With the assistance of chromatin conformation capture technology, 99.9% of the contigs were anchored into three pseudochromosomes with a total size of 86.02 Mb. Repetitive elements, accounting for 14.16% of this genome (12.20 Mb), are predominantly composed of long-terminal repeats (30.7%). By combining evidence of ab initio prediction, transcripts, and homologous proteins, we annotated 11,881 protein-coding genes. Both the genome and proteins have high BUSCO completeness scores (>94%). This high-quality genome, along with reliable annotation, provides a valuable resource for investigating the high-temperature tolerance of this species and exploring the genomic basis that underlies the host range evolution of spider mites.

A rapid approach for sex assignment by RAD-seq using a reference genome

Sex identification is a common objective in molecular ecology. While many vertebrates display sexual dimorphism, determining the sex can be challenging in certain situations, such as species lacking clear sex-related phenotypic characteristics or in studies using non-invasive methods. In these cases, DNA analyses serve as valuable tools not only for sex determination but also for validating sex assignment based on phenotypic traits. In this study, we developed a bioinformatic framework for sex assignment using genomic data obtained through GBS, and having an available closely related genome assembled at the chromosome level. Our method consists of two ad hoc indexes that rely on the different properties of the mammalian heteromorphic sex chromosomes. For this purpose, we mapped RAD-seq loci to a reference genome and then obtained missingness and coverage depth values for the autosomes and X and Y chromosomes of each individual. Our methodology successfully determined the sex of 165 fur seals that had been phenotypically sexed in a previous study and 40 sea lions sampled in a non-invasive way. Additionally, we evaluated the accuracy of each index in sequences with varying average coverage depths, with Index Y proving greater reliability and robustness in assigning sex to individuals with low-depth coverage. We believe that the approach presented here can be extended to any animal taxa with known heteromorphic XY/ZW sex chromosome systems and that it can tolerate various qualities of GBS sequencing data.

A chromosome-scale fishing cat reference genome for the evaluation of potential germline risk variants

The fishing cat, Prionailurus viverrinus, faces a population decline, increasing the importance of maintaining healthy zoo populations. Unfortunately, zoo-managed individuals currently face a high prevalence of transitional cell carcinoma (TCC), a form of bladder cancer. To investigate the genetics of inherited diseases among captive fishing cats, we present a chromosome-scale assembly, generate the pedigree of the zoo-managed population, reaffirm the close genetic relationship with the Asian leopard cat (Prionailurus bengalensis), and identify 7.4 million single nucleotide variants (SNVs) and 23,432 structural variants (SVs) from whole genome sequencing (WGS) data of healthy and TCC cats. Only BRCA2 was found to have a high recurrent number of missense mutations in fishing cats diagnosed with TCC when compared to inherited human cancer risk variants. These new fishing cat genomic resources will aid conservation efforts to improve their genetic fitness and enhance the comparative study of feline genomes.

The chromosomal-scale genome sequencing and assembly of Athetis lepigone

Athetis lepigone is an emerging highly polyphagous insect pest reported to cause crop damage in several European and Asian countries. However, our understanding of its genetic adaptation mechanisms has been limited due to lack of high-quality genetic resources. In this study, we present a chromosomal-level genome of A. lepigone, representing the first species in the genus of Athetis. We employed PacBio long-read sequencing and Hi-C technologies to generate 612.49 Mb genome assembly which contains 42.43% repeat sequences with a scaffold N50 of 20.9 Mb. The contigs were successfully clustered into 31 chromosomal-size scaffolds with 37% GC content. BUSCO assessment revealed a genome completeness of 97.4% with 96.3 identified as core Arthropoda single copy orthologs. Among the 17,322 genes that were predicted, 15,965 genes were functionally annotated, representing a coverage of 92.17%. Furthermore, we revealed 106 P450, 37 GST, 27 UGT, and 74 COE gene families in the genome of A. lepigone. This genome provides a significant and invaluable genomic resource for further research across the entire genus of Athetis.

Complete plastid genome of Iris orchioides and comparative analysis with 19 Iris plastomes

Iris is a cosmopolitan genus comprising approximately 280 species distributed throughout the Northern Hemisphere. Although Iris is the most diverse group in the Iridaceae, the number of taxa is debatable owing to various taxonomic issues. Plastid genomes have been widely used for phylogenetic research in plants; however, only limited number of plastid DNA markers are available for phylogenetic study of the Iris. To understand the genomic features of plastids within the genus, including its structural and genetic variation, we newly sequenced and analyzed the complete plastid genome of I. orchioides and compared it with those of 19 other Iris taxa. Potential plastid markers for phylogenetic research were identified by computing the sequence divergence and phylogenetic informativeness. We then tested the utility of the markers with the phylogenies inferred from the markers and whole-plastome data. The average size of the plastid genome was 152,926 bp, and the overall genomic content and organization were nearly identical among the 20 Iris taxa, except for minor variations in the inverted repeats. We identified 10 highly informative regions (matK, ndhF, rpoC2, ycf1, ycf2, rps15-ycf, rpoB-trnC, petA-psbJ, ndhG-ndhI and psbK-trnQ) and inferred a phylogeny from each region individually, as well as from their concatenated data. Remarkably, the phylogeny reconstructed from the concatenated data comprising three selected regions (rpoC2, ycf1 and ycf2) exhibited the highest congruence with the phylogeny derived from the entire plastome dataset. The result suggests that this subset of data could serve as a viable alternative to the complete plastome data, especially for molecular diagnoses among closely related Iris taxa, and at a lower cost.

Adaptive functions of structural variants in human brain development

Quantifying the structural variants (SVs) in nonhuman primates could provide a niche to clarify the genetic backgrounds underlying human-specific traits, but such resource is largely lacking. Here, we report an accurate SV map in a population of 562 rhesus macaques, verified by in-house benchmarks of eight macaque genomes with long-read sequencing and another one with genome assembly. This map indicates stronger selective constrains on inversions at regulatory regions, suggesting a strategy for prioritizing them with the most important functions. Accordingly, we identified 75 human-specific inversions and prioritized them. The top-ranked inversions have substantially shaped the human transcriptome, through their dual effects of reconfiguring the ancestral genomic architecture and introducing regional mutation hotspots at the inverted regions. As a proof of concept, we linked APCDD1, located on one of these inversions and down-regulated specifically in humans, to neuronal maturation and cognitive ability. We thus highlight inversions in shaping the human uniqueness in brain development.

Genomic characterization of a novel Hepatovirus identified in Maranhão state, Brazil

Bats are efficient reservoirs of a number of viruses with zoonotic potential, and are involved directly in the transmission cycle of many zoonoses. In the present study, which is part of a larger project that is documenting the viromes of the bat species found in the Mid-North states of Maranhão and Piauí, we analyzed 16 pooled samples obtained from four species of bat of the genus Artibeus-Artibeus obscurus, Artibeus cinereus, Artibeus lituratus and Artibeus planirostris. We describe and identify a Hepatovirus, denominated Hepatovirus H isolate sotense, which was found in a pool of internal organs (liver and lungs) extracted from a specimen of A. planirostris, a frugivorous bat, collected in the Cerrado biome of Maranhão state. This material was analyzed using new generation sequencing, which produced a contig of 7390 nucleotides and presented a degree of identity with a number of existing Hepatovirus sequences available for bats (amino acid identity of 61.5% with Bat hepatovirus C of Miniopterus cf. manavi, 66.6% with Bat hepatovirus G of Coleura afra, 67.4% with Hepatovirus G2 of Rhinolophus landeri, and 75.3% with Hepatovirus H2 of Rhinolophus landeri). The analysis of the functional domains of this contig confirmed a pattern consistent with the characteristics of the genus Hepatovirus (Picornaviridae). In the phylogenetic tree with several other Hepatovirus species, this genome also grouped in a monophyletic clade with Hepatovirus H (HepV-H1; HepV-H2, and HepV-H3) albeit on an external branch, which suggests that it may be a distinct genotype within this species. This is the first isolate of Hepatovirus H identified in bats from South America, and represents an important discovery, given that most studies of viruses associated with bats in the state of Maranhão have focused on the family Rhabdoviridae.

Targeted phasing of 2-200 kilobase DNA fragments with a short-read sequencer and a single-tube linked-read library method

In the human genome, heterozygous sites refer to genomic positions with a different allele or nucleotide variant on the maternal and paternal chromosomes. Resolving these allelic differences by chromosomal copy, also known as phasing, is achievable on a short-read sequencer when using a library preparation method that captures long-range genomic information. TELL-Seq is a library preparation that captures long-range genomic information with the aid of molecular identifiers (barcodes). The same barcode is used to tag the reads derived from the same long DNA fragment within a range of up to 200 kilobases (kb), generating linked-reads. This strategy can be used to phase an entire genome. Here, we introduce a TELL-Seq protocol developed for targeted applications, enabling the phasing of enriched loci of varying sizes, purity levels, and heterozygosity. To validate this protocol, we phased 2-200 kb loci enriched with different methods: CRISPR/Cas9-mediated excision coupled with pulse-field electrophoresis for the longest fragments, CRISPR/Cas9-mediated protection from exonuclease digestion for mid-size fragments, and long PCR for the shortest fragments. All selected loci have known clinical relevance: BRCA1, BRCA2, MLH1, MSH2, MSH6, APC, PMS2, SCN5A-SCN10A, and PKI3CA. Collectively, the analyses show that TELL-Seq can accurately phase 2-200 kb targets using a short-read sequencer.

Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle

The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from -0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075-35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.

Developmental trajectories and cooperating genomic events define molecular subtypes of BCR::ABL1-positive ALL

ABSTRACT

Distinct diagnostic entities within BCR::ABL1-positive acute lymphoblastic leukemia (ALL) are currently defined by the International Consensus Classification of myeloid neoplasms and acute leukemias (ICC): "lymphoid only", with BCR::ABL1 observed exclusively in lymphatic precursors, vs "multilineage", where BCR::ABL1 is also present in other hematopoietic lineages. Here, we analyzed transcriptomes of 327 BCR::ABL1-positive patients with ALL (age, 2-84 years; median, 46 years) and identified 2 main gene expression clusters reproducible across 4 independent patient cohorts. Fluorescence in situ hybridization analysis of fluorescence-activated cell-sorted hematopoietic compartments showed distinct BCR::ABL1 involvement in myeloid cells for these clusters (n = 18/18 vs n = 3/16 patients; P < .001), indicating that a multilineage or lymphoid BCR::ABL1 subtype can be inferred from gene expression. Further subclusters grouped samples according to cooperating genomic events (multilineage: HBS1L deletion or monosomy 7; lymphoid: IKZF1-/- or CDKN2A/PAX5 deletions/hyperdiploidy). A novel HSB1L transcript was highly specific for BCR::ABL1 multilineage cases independent of HBS1L genomic aberrations. Treatment on current German Multicenter Study Group for Adult ALL (GMALL) protocols resulted in comparable disease-free survival (DFS) for multilineage vs lymphoid cluster patients (3-year DFS: 70% vs 61%; P = .530; n = 91). However, the IKZF1-/- enriched lymphoid subcluster was associated with inferior DFS, whereas hyperdiploid cases showed a superior outcome. Thus, gene expression clusters define underlying developmental trajectories and distinct patterns of cooperating events in BCR::ABL1-positive ALL with prognostic relevance.

From Base Pairs to City Squares: Comprehensive Precision Oncology for the Future

SUMMARY

An increased understanding of the role of the social determinants of health in cancer prevention, cancer care, and outcomes can lead to their integration into genetics and genomics as well as informing interventions and clinical trials, creating a comprehensive precision oncology framework.

Integrated meta-analysis and transcriptomics pinpoint genomic loci and novel candidate genes associated with submergence tolerance in rice

BACKGROUND

Due to rising costs, water shortages, and labour shortages, farmers across the globe now prefer a direct seeding approach. However, submergence stress remains a major bottleneck limiting the success of this approach in rice cultivation. The merger of accumulated rice genetic resources provides an opportunity to detect key genomic loci and candidate genes that influence the flooding tolerance of rice.

RESULTS

In the present study, a whole-genome meta-analysis was conducted on 120 quantitative trait loci (QTL) obtained from 16 independent QTL studies reported from 2004 to 2023. These QTL were confined to 18 meta-QTL (MQTL), and ten MQTL were successfully validated by independent genome-wide association studies from diverse natural populations. The mean confidence interval (CI) of the identified MQTL was 3.44 times narrower than the mean CI of the initial QTL. Moreover, four core MQTL loci with genetic distance less than 2 cM were obtained. By combining differentially expressed genes (DEG) from two transcriptome datasets with 858 candidate genes identified in the core MQTL regions, we found 38 common differentially expressed candidate genes (DECGs). In silico expression analysis of these DECGs led to the identification of 21 genes with high expression in embryo and coleoptile under submerged conditions. These DECGs encode proteins with known functions involved in submergence tolerance including WRKY, F-box, zinc fingers, glycosyltransferase, protein kinase, cytochrome P450, PP2C, hypoxia-responsive family, and DUF domain. By haplotype analysis, the 21 DECGs demonstrated distinct genetic differentiation and substantial genetic distance mainly between indica and japonica subspecies. Further, the MQTL7.1 was successfully validated using flanked marker S2329 on a set of genotypes with phenotypic variation.

CONCLUSION

This study provides a new perspective on understanding the genetic basis of submergence tolerance in rice. The identified MQTL and novel candidate genes lay the foundation for marker-assisted breeding/engineering of flooding-tolerant cultivars conducive to direct seeding.

Multienvironment genomic prediction in tetraploid potato

Multienvironment genomic prediction was applied to tetraploid potato using 147 potato varieties, tested for 2 years, in 3 locations representative of 3 distinct regions in Europe. Different prediction scenarios were investigated to help breeders predict genotypic performance in the regions from one year to the next, for genotypes that were tested this year (scenario 1), as well as new genotypes (scenario 3). In scenario 2, we predicted new genotypes for any one of the 6 trials, using all the information that is available. The choice of prediction model required assessment of the variance-covariance matrix in a mixed model that takes into account heterogeneity of genetic variances and correlations. This was done for each analyzed trait (tuber weight, tuber length, and dry matter) where examples of both limited and higher degrees of heterogeneity was observed. This explains why dry matter did not need complex multienvironment modeling to combine environments and increase prediction ability, while prediction in tuber weight, improved only when models were flexible enough to capture the heterogeneous variances and covariances between environments. We also found that the prediction abilities in a target trial condition decreased, if trials with a low genetic correlation to the target were included when training the model. Genomic prediction in tetraploid potato can work once there is clarity about the prediction scenario, a suitable training set is created, and a multienvironment prediction model is chosen based on the patterns of G×E indicated by the genetic variances and covariances.

Enhancing grapevine breeding efficiency through genomic prediction and selection index

Grapevine (Vitis vinifera) breeding reaches a critical point. New cultivars are released every year with resistance to powdery and downy mildews. However, the traditional process remains time-consuming, taking 20-25 years, and demands the evaluation of new traits to enhance grapevine adaptation to climate change. Until now, the selection process has relied on phenotypic data and a limited number of molecular markers for simple genetic traits such as resistance to pathogens, without a clearly defined ideotype, and was carried out on a large scale. To accelerate the breeding process and address these challenges, we investigated the use of genomic prediction, a methodology using molecular markers to predict genotypic values. In our study, we focused on 2 existing grapevine breeding programs: Rosé wine and Cognac production. In these programs, several families were created through crosses of emblematic and interspecific resistant varieties to powdery and downy mildews. Thirty traits were evaluated for each program, using 2 genomic prediction methods: Genomic Best Linear Unbiased Predictor and Least Absolute Shrinkage Selection Operator. The results revealed substantial variability in predictive abilities across traits, ranging from 0 to 0.9. These discrepancies could be attributed to factors such as trait heritability and trait characteristics. Moreover, we explored the potential of across-population genomic prediction by leveraging other grapevine populations as training sets. Integrating genomic prediction allowed us to identify superior individuals for each program, using multivariate selection index method. The ideotype for each breeding program was defined collaboratively with representatives from the wine-growing sector.

Prospective cohort study of genomic newborn screening: BabyScreen+ pilot study protocol

INTRODUCTION

Newborn bloodspot screening (NBS) is a highly successful public health programme that uses biochemical and other assays to screen for severe but treatable childhood-onset conditions. Introducing genomic sequencing into NBS programmes increases the range of detectable conditions but raises practical and ethical issues. Evidence from prospectively ascertained cohorts is required to guide policy and future implementation. This study aims to develop, implement and evaluate a genomic NBS (gNBS) pilot programme.

METHODS AND ANALYSIS

The BabyScreen+ study will pilot gNBS in three phases. In the preimplementation phase, study materials, including education resources, decision support and data collection tools, will be designed. Focus groups and key informant interviews will also be undertaken to inform delivery of the study and future gNBS programmes. During the implementation phase, we will prospectively recruit birth parents in Victoria, Australia, to screen 1000 newborns for over 600 severe, treatable, childhood-onset conditions. Clinically accredited whole genome sequencing will be performed following standard NBS using the same sample. High chance results will be returned by genetic healthcare professionals, with follow-on genetic and other confirmatory testing and referral to specialist services as required. The postimplementation phase will evaluate the feasibility of gNBS as the primary aim, and assess ethical, implementation, psychosocial and health economic factors to inform future service delivery.

ETHICS AND DISSEMINATION

This project received ethics approval from the Royal Children's Hospital Melbourne Research Ethics Committee: HREC/91500/RCHM-2023, HREC/90929/RCHM-2022 and HREC/91392/RCHM-2022. Findings will be disseminated to policy-makers, and through peer-reviewed journals and conferences.

Distinct genomic contexts predict gene presence-absence variation in different pathotypes of Magnaporthe oryzae

Fungi use the accessory gene content of their pangenomes to adapt to their environments. While gene presence-absence variation contributes to shaping accessory gene reservoirs, the genomic contexts that shape these events remain unclear. Since pangenome studies are typically species-wide and do not analyze different populations separately, it is yet to be uncovered whether presence-absence variation patterns and mechanisms are consistent across populations. Fungal plant pathogens are useful models for studying presence-absence variation because they rely on it to adapt to their hosts, and members of a species often infect distinct hosts. We analyzed gene presence-absence variation in the blast fungus, Magnaporthe oryzae (syn. Pyricularia oryzae), and found that presence-absence variation genes involved in host-pathogen and microbe-microbe interactions may drive the adaptation of the fungus to its environment. We then analyzed genomic and epigenomic features of presence-absence variation and observed that proximity to transposable elements, gene GC content, gene length, expression level in the host, and histone H3K27me3 marks were different between presence-absence variation genes and conserved genes. We used these features to construct a model that was able to predict whether a gene is likely to experience presence-absence variation with high precision (86.06%) and recall (92.88%) in M. oryzae. Finally, we found that presence-absence variation genes in the rice and wheat pathotypes of M. oryzae differed in their number and their genomic context. Our results suggest that genomic and epigenomic features of gene presence-absence variation can be used to better understand and predict fungal pangenome evolution. We also show that substantial intra-species variation can exist in these features.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Whole-genome sequencing and comparative genomic analyses of the medicinal fungus Sanguinoderma infundibulare in Ganodermataceae

Sanguinoderma infundibulare is a newly discovered species of Ganodermataceae known to have high medicinal and ecological values. In this study, the whole-genome sequencing and comparative genomic analyses were conducted to further understand Ganodermataceae's genomic structural and functional characteristics. Using the Illumina NovaSeq and PacBio Sequel platforms, 88 scaffolds were assembled to obtain a 48.99-Mb high-quality genome of S. infundibulare. A total of 14,146 protein-coding genes were annotated in the whole genome, with 98.6% of complete benchmarking universal single-copy orthologs (BUSCO) scores. Comparative genomic analyses were conducted among S. infundibulare, Sanguinoderma rugosum, Ganoderma lucidum, and Ganoderma sinense to determine their intergeneric differences. The 4 species were found to share 4,011 orthogroups, and 24 specific gene families were detected in the genus Sanguinoderma. The gene families associated with carbohydrate esterase in S. infundibulare were significantly abundant, which was reported to be involved in hemicellulose degradation. One specific gene family in Sanguinoderma was annotated with siroheme synthase, which may be related to the typical characteristics of fresh pore surface changing to blood red when bruised. This study enriched the available genome data for the genus Sanguinoderma, elucidated the differences between Ganoderma and Sanguinoderma, and provided insights into the characteristics of the genome structure and function of S. infundibulare.

The Molecular Genetics and Genomics of Plant-Pathogen Interactions

Plants have evolved an intricate immune system to protect themselves from potential pathogens [...].

Shared evolutionary processes shape landscapes of genomic variation in the great apes

For at least the past 5 decades, population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modeling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well-sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations, we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modeling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

Genomic insights into the evolution and adaptation of secondary metabolite gene clusters in fungicolous species Cladobotryum mycophilum ATHUM6906

Mycophilic or fungicolous fungi can be found wherever fungi exist since they are able to colonize other fungi, which occupy a diverse range of habitats. Some fungicolous species cause important diseases on Basidiomycetes, and thus, they are the main reason for the destruction of mushroom cultivations. Nonetheless, despite their ecological significance, their genomic data remain limited. Cladobotryum mycophilum is one of the most aggressive species of the genus, destroying the economically important Agaricus bisporus cultivations. The 40.7 Mb whole genome of the Greek isolate ATHUM6906 is assembled in 16 fragments, including the mitochondrial genome and 2 small circular mitochondrial plasmids, in this study. This genome includes a comprehensive set of 12,282 protein coding, 56 rRNA, and 273 tRNA genes. Transposable elements, CAZymes, and pathogenicity related genes were also examined. The genome of C. mycophilum contained a diverse arsenal of genes involved in secondary metabolism, forming 106 biosynthetic gene clusters, which renders this genome as one of the most BGC abundant among fungicolous species. Comparative analyses were performed for genomes of species of the family Hypocreaceae. Some BGCs identified in C. mycophilum genome exhibited similarities to clusters found in the family Hypocreaceae, suggesting vertical heritage. In contrast, certain BGCs showed a scattered distribution among Hypocreaceae species or were solely found in Cladobotryum genomes. This work provides evidence of extensive BGC losses, horizontal gene transfer events, and formation of novel BGCs during evolution, potentially driven by neutral or even positive selection pressures. These events may increase Cladobotryum fitness under various environmental conditions and potentially during host-fungus interaction.