Empowering the biomedical research community: Innovative SAS deployment on the All of Us Researcher Workbench.

OBJECTIVES: The All of Us Research Program is a precision medicine initiative aimed at establishing a vast, diverse biomedical database accessible through a cloud-based data analysis platform, the Researcher Workbench (RW). Our goal was to empower the research community by co-designing the implementation of SAS in the RW alongside researchers to enable broader use of All of Us data. MATERIALS AND METHODS: Researchers from various fields and with different SAS experience levels participated in co-designing the SAS implementation through user experience interviews. RESULTS: Feedback and lessons learned from user testing informed the final design of the SAS application. DISCUSSION: The co-design approach is critical for reducing technical barriers, broadening All of Us data use, and enhancing the user experience for data analysis on the RW. CONCLUSION: Our co-design approach successfully tailored the implementation of the SAS application to researchers' needs. This approach may inform future software implementations on the RW.

HPV18 E6/E7 activates Ca2+ influx to promote the malignant progression of cervical cancer by inhibiting Ca2+ binding protein 1 expression.

Mounting studies have shown that the oncoproteins E6 and E7 encoded by the human papillomavirus (HPV) genome are essential in HPV-induced cervical cancer (CC). Ca2+ binding protein 1 (CABP1), a downstream target of HPV18-positive HeLa cells that interferes with E6/E7 expression, was identified through screening the GEO Database (GSE6926). It was confirmed to be down-regulated in CC through TCGA prediction and in vitro detection. Subsequent in vitro experiments revealed that knocking down E6/E7 inhibited cell proliferation, migration, and invasion, whereas knocking down CABP1 promoted these processes. Simultaneously knocking down CABP1 reversed these effects. Additionally, the results were validated in vivo. Previous studies have indicated that CABP1 can regulate Ca2+ channels, influencing Ca2+ influx and tumor progression. In this study, it was observed that knocking down CABP1 enhanced Ca2+ inflow, as demonstrated by flow cytometry and confocal microscopy. Knocking down E6/E7 inhibited these processes, whereas simultaneously knocking down E6/E7 and CABP1 restored the inhibitory effect of knocking down E6/E7 on Ca2+ inflow. To further elucidate that E6/E7 promotes CC progression by inhibiting CABP1 expression and activating Ca2+ influx, BAPTA/AM treatment was administered during CABP1 knockdown. It was discovered that Ca2+ chelation could reverse the effect of CABP1 knockdown on CC cells. In conclusion, our results offer a novel target for the diagnosis and treatment of HPV-induced CC.

Integrative genomics reveals the polygenic basis of seedlessness in grapevine.

Seedlessness is a crucial quality trait in table grape (Vitis vinifera L.) breeding. However, the development of seeds involved intricate regulations, and the polygenic basis of seed abortion remains unclear. Here, we combine comparative genomics, population genetics, quantitative genetics, and integrative genomics to unravel the evolution and polygenic basis of seedlessness in grapes. We generated the haplotype-resolved genomes for two seedless grape cultivars, "Thompson Seedless" (TS, syn. "Sultania") and "Black Monukka" (BM). Comparative genomics identified a ∼4.25 Mb hemizygous inversion on Chr10 specific in seedless cultivars, with seedless-associated genes VvTT16 and VvSUS2 located at breakpoints. Population genomic analyses of 548 grapevine accessions revealed two distinct clusters of seedless cultivars, and the identity-by-descent (IBD) results indicated that the origin of the seedlessness trait could be traced back to "Sultania." Introgression, rather than convergent selection, shaped the evolutionary history of seedlessness in grape improvement. Genome-wide association study (GWAS) analysis identified 110 quantitative trait loci (QTLs) associated with 634 candidate genes, including previously unidentified candidate genes, such as three 11S GLOBULIN SEED STORAGE PROTEIN and two CYTOCHROME P450 genes, and well-known genes like VviAGL11. Integrative genomic analyses resulted in 339 core candidate genes categorized into 13 functional categories related to seed development. Machine learning-based genomic selection achieved a remarkable prediction accuracy of 97% for seedlessness in grapevines. Our findings highlight the polygenic nature of seedlessness and provide candidate genes for molecular genetics and an effective prediction for seedlessness in grape genomic breeding.

Targeting TANK-binding kinase 1 attenuates painful diabetic neuropathy via inhibiting microglia pyroptosis.

BACKGROUND: Painful diabetic neuropathy (PDN) is closely linked to inflammation, which has been demonstrated to be associated with pyroptosis. Emerging evidence has implicated TANK-binding kinase 1 (TBK1) in various inflammatory diseases. However, it remains unknown whether activated TBK1 causes hyperalgesia via pyroptosis. METHODS: PDN mice model of type 1 or type 2 diabetic was induced by C57BL/6J or BKS-DB mice with Lepr gene mutation. For type 2 diabetes PDN model, TBK1-siRNA, Caspase-1 inhibitor Ac-YVAD-cmk or TBK1 inhibitor amlexanox (AMX) were delivered by intrathecal injection or intragastric administration. The pain threshold and plantar skin blood perfusion were evaluated through animal experiments. The assessments of spinal cord, dorsal root ganglion, sciatic nerve, plantar skin and serum included western blotting, immunofluorescence, ELISA, and transmission electron microscopy. RESULTS: In the PDN mouse model, we found that TBK1 was significantly activated in the spinal dorsal horn (SDH) and mainly located in microglia, and intrathecal injection of chemically modified TBK1-siRNA could improve hyperalgesia. Herein, we described the mechanism that TBK1 could activate the noncanonical nuclear factor κB (NF-κB) pathway, mediate the activation of NLRP3 inflammasome, trigger microglia pyroptosis, and ultimately induce PDN, which could be reversed following TBK1-siRNA injection. We also found that systemic administration of AMX, a TBK1 inhibitor, could effectively improve peripheral nerve injury. These results revealed the key role of TBK1 in PDN and that TBK1 inhibitor AMX could be a potential strategy for treating PDN. CONCLUSIONS: Our findings revealed a novel causal role of TBK1 in pathogenesis of PDN, which raises the possibility of applying amlexanox to selectively target TBK1 as a potential therapeutic strategy for PDN.

Identification of immune-associated biomarkers of diabetes nephropathy tubulointerstitial injury based on machine learning: a bioinformatics multi-chip integrated analysis.

BACKGROUND: Diabetic nephropathy (DN) is a major microvascular complication of diabetes and has become the leading cause of end-stage renal disease worldwide. A considerable number of DN patients have experienced irreversible end-stage renal disease progression due to the inability to diagnose the disease early. Therefore, reliable biomarkers that are helpful for early diagnosis and treatment are identified. The migration of immune cells to the kidney is considered to be a key step in the progression of DN-related vascular injury. Therefore, finding markers in this process may be more helpful for the early diagnosis and progression prediction of DN. METHODS: The gene chip data were retrieved from the GEO database using the search term ' diabetic nephropathy '. The ' limma ' software package was used to identify differentially expressed genes (DEGs) between DN and control samples. Gene set enrichment analysis (GSEA) was performed on genes obtained from the molecular characteristic database (MSigDB. The R package 'WGCNA' was used to identify gene modules associated with tubulointerstitial injury in DN, and it was crossed with immune-related DEGs to identify target genes. Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed on differentially expressed genes using the 'ClusterProfiler' software package in R. Three methods, least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE) and random forest (RF), were used to select immune-related biomarkers for diagnosis. We retrieved the tubulointerstitial dataset from the Nephroseq database to construct an external validation dataset. Unsupervised clustering analysis of the expression levels of immune-related biomarkers was performed using the 'ConsensusClusterPlus 'R software package. The urine of patients who visited Dongzhimen Hospital of Beijing University of Chinese Medicine from September 2021 to March 2023 was collected, and Elisa was used to detect the mRNA expression level of immune-related biomarkers in urine. Pearson correlation analysis was used to detect the effect of immune-related biomarker expression on renal function in DN patients. RESULTS: Four microarray datasets from the GEO database are included in the analysis : GSE30122, GSE47185, GSE99340 and GSE104954. These datasets included 63 DN patients and 55 healthy controls. A total of 9415 genes were detected in the data set. We found 153 differentially expressed immune-related genes, of which 112 genes were up-regulated, 41 genes were down-regulated, and 119 overlapping genes were identified. GO analysis showed that they were involved in various biological processes including leukocyte-mediated immunity. KEGG analysis showed that these target genes were mainly involved in the formation of phagosomes in Staphylococcus aureus infection. Among these 119 overlapping genes, machine learning results identified AGR2, CCR2, CEBPD, CISH, CX3CR1, DEFB1 and FSTL1 as potential tubulointerstitial immune-related biomarkers. External validation suggested that the above markers showed diagnostic efficacy in distinguishing DN patients from healthy controls. Clinical studies have shown that the expression of AGR2, CX3CR1 and FSTL1 in urine samples of DN patients is negatively correlated with GFR, the expression of CX3CR1 and FSTL1 in urine samples of DN is positively correlated with serum creatinine, while the expression of DEFB1 in urine samples of DN is negatively correlated with serum creatinine. In addition, the expression of CX3CR1 in DN urine samples was positively correlated with proteinuria, while the expression of DEFB1 in DN urine samples was negatively correlated with proteinuria. Finally, according to the level of proteinuria, DN patients were divided into nephrotic proteinuria group (n = 24) and subrenal proteinuria group. There were significant differences in urinary AGR2, CCR2 and DEFB1 between the two groups by unpaired t test (P < 0.05). CONCLUSIONS: Our study provides new insights into the role of immune-related biomarkers in DN tubulointerstitial injury and provides potential targets for early diagnosis and treatment of DN patients. Seven different genes ( AGR2, CCR2, CEBPD, CISH, CX3CR1, DEFB1, FSTL1 ), as promising sensitive biomarkers, may affect the progression of DN by regulating immune inflammatory response. However, further comprehensive studies are needed to fully understand their exact molecular mechanisms and functional pathways in DN.

Phased telomere-to-telomere reference genome and pangenome reveal an expansion of resistance genes during apple domestication.

The cultivated apple (Malus domestica Borkh.) is a cross-pollinated perennial fruit tree of great economic importance. Earlier versions of apple reference genomes were unphased, fragmented, and lacked comprehensive insights into the apple's highly heterozygous genome, which impeded advances in genetic studies and breeding programs. In this study, we assembled a haplotype-resolved telomere-to-telomere (T2T) reference genome for the diploid apple cultivar Golden Delicious. Subsequently, we constructed a pangenome based on 12 assemblies from wild and cultivated species to investigate the dynamic changes of functional genes. Our results revealed the gene gain and loss events during apple domestication. Compared with cultivated species, more gene families in wild species were significantly enriched in oxidative phosphorylation, pentose metabolic process, responses to salt, and abscisic acid biosynthesis process. Our analyses also demonstrated a higher prevalence of different types of resistance gene analogs (RGAs) in cultivars than their wild relatives, partially attributed to segmental and tandem duplication events in certain RGAs classes. Structural variations, mainly deletions and insertions, have affected the presence and absence of TIR-NB-ARC-LRR, NB-ARC-LRR, and CC-NB-ARC-LRR genes. Additionally, hybridization/introgression from wild species has also contributed to the expansion of resistance genes in domesticated apples. Our haplotype-resolved T2T genome and pangenome provide important resources for genetic studies of apples, emphasizing the need to study the evolutionary mechanisms of resistance genes in apple breeding.

Evolutionary genomics of climatic adaptation and resilience to climate change in alfalfa.

Given the escalating impact of climate change on agriculture and food security, gaining insights into the evolutionary dynamics of climatic adaptation and uncovering climate-adapted variation can empower the breeding of climate-resilient crops to face future climate change. Alfalfa (Medicago sativa subsp. sativa), the queen of forages, shows remarkable adaptability across diverse global environments, making it an excellent model for investigating species responses to climate change. In this study, we performed population genomic analyses using genome resequencing data from 702 accessions of 24 Medicago species to unravel alfalfa's climatic adaptation and genetic susceptibility to future climate change. We found that interspecific genetic exchange has contributed to the gene pool of alfalfa, particularly enriching defense and stress-response genes. Intersubspecific introgression between M. sativa subsp. falcata (subsp. falcata) and alfalfa not only aids alfalfa's climatic adaptation but also introduces genetic burden. A total of 1671 genes were associated with climatic adaptation, and 5.7% of them were introgressions from subsp. falcata. By integrating climate-associated variants and climate data, we identified populations that are vulnerable to future climate change, particularly in higher latitudes of the Northern Hemisphere. These findings serve as a clarion call for targeted conservation initiatives and breeding efforts. We also identified pre-adaptive populations that demonstrate heightened resilience to climate fluctuations, illuminating a pathway for future breeding strategies. Collectively, this study enhances our understanding about the local adaptation mechanisms of alfalfa and facilitates the breeding of climate-resilient alfalfa cultivars, contributing to effective agricultural strategies for facing future climate change.

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape.

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Treatment of diabetic foot with moxibustion: Clinical evidence from meta-analysis.

To assess the efficacy of moxibustion for diabetic foot, and compile the findings of randomised clinical trials. China National Knowledge Infrastructure Database, Medicine, WanFang Database, Embase, Chinese Scientific Journal Database and Web of Science were from the establishment to January, 2024 were searched. Randomised controlled trials, which evaluated the effects of moxibustion were included. A total of 12 randomised controlled trials involving 1196 patients were included. According to the pooled results of this meta-analysis, effective rate (relative risk 1.16, 95% confidence intervals, CI [1.11, 1.22]), healing time (mean difference [MD] -6.27, 95% CI [-8.68, -3.86]), wound area (MD 3.46, 95% CI [0.84, 6.09]), and ankle brachial index (MD 0.14, 95% CI [0.03, 0.24]) were statistically significant compared to the control group. This study suggests that moxibustion treatment has the potential for improving symptoms of diabetic foot. However, future in-depth research on the benefits and harms of moxibustion for the diabetic foot is needed before it can be accepted as an evidence-based treatment.

Population comparative genomics discovers gene gain and loss during grapevine domestication.

Plant domestication are evolutionary experiments conducted by early farmers since thousands years ago, during which the crop wild progenitors are artificially selected for desired agronomic traits along with dramatic genomic variation in the course of moderate to severe bottlenecks. However, previous investigations are mainly focused on small-effect variants, while changes in gene contents are rarely investigated due to the lack of population-level assemblies for both the crop and its wild relatives. Here, we applied comparative genomic analyses to discover gene gain and loss during grapevine domestication using long-read assemblies of representative population samples for both domesticated grapevines (V. vinifera ssp. vinifera) and their wild progenitors (V. vinifera ssp. sylvestris). Only ∼7% of gene families were shared by 16 Vitis genomes while ∼8% of gene families were specific to each accession, suggesting dramatic variations of gene contents in grapevine genomes. Compared to wild progenitors, the domesticated accessions exhibited an increased presence of genes associated with asexual reproduction, while the wild progenitors showcased a higher abundance of genes related to pollination, revealing the transition from sexual reproduction to clonal propagation during domestication processes. Moreover, the domesticated accessions harbored fewer disease-resistance genes than wild progenitors. The SVs occurred frequently in aroma and disease-resistance related genes between domesticated grapevines and wild progenitors, indicating the rapid diversification of these genes during domestication. Our study provides insights and resources for biological studies and breeding programs in grapevine.

The protective effects of sevoflurane on subarachnoid hemorrhage.

Sevoflurane has become an important volatile anesthetic in clinic and has been widely studied in recent years. Numerous studies have demonstrated the efficacy of sevoflurane in safeguarding against brain damage across various domains. For example, it has played a neuroprotective role in subarachnoid hemorrhage (SAH), traumatic brain injury, and ischemia/reperfusion injury. The ensuing critique will focus on the significance of sevoflurane in experimental SAH and shed light on the underlying mechanisms. The findings of the current investigation demonstrate that sevoflurane possesses neuroprotective capabilities and clarify that it effectively attenuates secondary damage resulting from SAH through anti-inflammatory and anti-apoptotic pathways. More specifically, sevoflurane is observed to mitigate arterial vasospasm, diminish microvascular thrombosis, and alleviate cerebral edema. In light of these discoveries, we maintain that sevoflurane exhibits significant promise in the management of SAH, and it merits additional investigation to facilitate its prompt clinical implementation. Therefore, a thorough understanding of the neuroprotective properties of sevoflurane is beneficial to exploring novel therapeutic solutions for SAH and providing clinicians with alternative treatment modalities.

CAFs-derived exosomes promote the development of cervical cancer by regulating miR-18a-5p-TMEM170B signaling axis.

Mounting studies have showed that tumor microenvironment (TME) is crucial for cervical cancer (CC), and cancer-related fibroblasts (CAFs) play a major role in it. Recently, exosomal miRNAs secreted by CAFs have been found to be potential targets for cancer diagnosis and therapy. In this paper, we aimed to investigate the function of CAFs-mediated exosome miR-18a-5p (CAFs-exo-miR-18a-5p) in CC. First, in combination with bioinformatic data analysis of the GEO database (GSE86100) and RT-qPCR of CC clinical tissue samples and cell lines, miR-18a-5p was discovered to be markedly up-regulated in CC. Next, CAFs-secreted exosomes were isolated and it was found that miR-18a-5p expression was dramatically promoted in CC cell lines when treated with CAFs-exos. The CAFs-exo-miR-18a-5p was then elucidated to stimulate the proliferation and migration and inhibit the apoptosis of CC cells. In order to clarify the underlying mechanism, we further screened the target genes of miR-18a-5p. TMEM170B was selected by bioinformatic data analysis of online databases combined with RT-qPCR of CC clinical tissues and cells. Luciferase reporter gene analysis combined with molecular biology experiments further elucidated that miR-18a-5p suppressed TMEM170B expression in CC. Finally, both cell and animal experiments demonstrated that TMEM170B over-expression attenuated the oncogenic effect of CAFs-exo-miR-18a-5p. In conclusion, our study indicates that CAFs-mediated exosome miR-18a-5p promotes the initiation and development of CC by suppressing TMEM170B signaling axis, which provides a possible direction for the diagnosis and therapy of CC.

Defective autophagy contributes to bupivacaine-induced aggravation of painful diabetic neuropathy in db/db mice.

Current evidence suggests that hyperactivated or impaired autophagy can lead to neuronal death. The effect of local anesthetics on painful diabetic neuropathy (PDN) and the role of autophagy in the above pathological process remain unclear, warranting further studies. So, PDN models were established by assessing the paw withdrawal threshold (PWT) and paw withdrawal latency (PWL) in leptin gene-mutation (db/db) mice. Wild type (WT) and PDN mice received intrathecal 0.75% bupivacaine or/with intraperitoneal drug treatment (rapamycin or bafilomycin A1). Subsequently, the PWT and PWL were measured to assess hyperalgesia at 6 h, 24 h, 30 h, and 48 h after intrathecal bupivacaine. Also, sensory nerve conduction velocity (SNCV) and motor nerve conduction velocity (MNCV) were measured before and 48 h after intrathecal bupivacaine treatment. The spinal cord tissue of L4-L6 segments and serum were harvested to evaluate the change of autophagy, oxidative stress, oxidative injury, and apoptosis. We found that bupivacaine induced the activation of autophagy but did not affect the pain threshold, SNCV and MNCV in WT mice at predefined time points. Conversely, bupivacaine lowered autophagosome generation and degradation, slowed SNCV and aggravated spinal dorsal horn neuron oxidative injury and hyperalgesia in PDN mice. The autophagy activator (rapamycin) could decrease spinal dorsal horn neuron oxidative injury, alleviate the alterations in SNCV and hyperalgesia in bupivacaine-treated PDN mice. Meanwhile, the autophagy inhibitor (bafilomycin A1) could exacerbate spinal dorsal horn neuron oxidative injury, the alterations in SNCV and hyperalgesia in bupivacaine-treated PDN mice. Our results showed that bupivacaine could induce defective autophagy, slowed SNCV and aggravate spinal dorsal horn neuron oxidative injury and hyperalgesia in PDN mice. Restoring autophagy may represent a potential therapeutic approach against nerve injury in PDN patients with local anesthesia and analgesia.

The haplotype-resolved T2T genome of teinturier cultivar Yan73 reveals the genetic basis of anthocyanin biosynthesis in grapes.

Teinturier grapes are characterized by the typical accumulation of anthocyanins in grape skin, flesh, and vegetative tissues, endowing them with high utility value in red wine blending and nutrient-enriched foods developing. However, due to the lack of genome information, the mechanism involved in regulating teinturier grape coloring has not yet been elucidated and their genetic utilization research is still insufficient. Here, the cultivar 'Yan73' was used for assembling the telomere-to-telomere (T2T) genome of teinturier grapes by combining the High Fidelity (HiFi)， Hi-C and ultralong Oxford Nanopore Technologies (ONT) reads. Two haplotype genomes were assembled, at the sizes of 501.68 Mb and 493.38 Mb, respectively. In the haplotype 1 genome, the transposable elements (TEs) contained 32.77% of long terminal repeats (LTRs), while in the haplotype 2 genome, 31.53% of LTRs were detected in TEs. Furthermore, obvious inversions were identified in chromosome 18 between the two haplotypes. Transcriptome profiling suggested that the gene expression patterns in 'Cabernet Sauvignon' and 'Yan73' were diverse depending on tissues, developmental stages, and varieties. The transcription program of genes in the anthocyanins biosynthesis pathway between the two cultivars exhibited high similarity in different tissues and developmental stages, whereas the expression levels of numerous genes showed significant differences. Compared with other genes, the expression levels of VvMYBA1 and VvUFGT4 in all samples, VvCHS2 except in young shoots and VvPAL9 except in the E-L23 stage of 'Yan73' were higher than those of 'Cabernet Sauvignon'. Further sequence alignments revealed potential variant gene loci and structure variations of anthocyanins biosynthesis related genes and a 816 bp sequence insertion was found in the promoter of VvMYBA1 of 'Yan73' haplotype 2 genome. The 'Yan73' T2T genome assembly and comparative analysis provided valuable foundations for further revealing the coloring mechanism of teinturier grapes and the genetic improvement of grape coloring traits.

Genomic conservation of crop wild relatives: A case study of citrus.

Conservation of crop wild relatives is critical for plant breeding and food security. The lack of clarity on the genetic factors that lead to endangered status or extinction create difficulties when attempting to develop concrete recommendations for conserving a citrus wild relative: the wild relatives of crops. Here, we evaluate the conservation of wild kumquat (Fortunella hindsii) using genomic, geographical, environmental, and phenotypic data, and forward simulations. Genome resequencing data from 73 accessions from the Fortunella genus were combined to investigate population structure, demography, inbreeding, introgression, and genetic load. Population structure was correlated with reproductive type (i.e., sexual and apomictic) and with a significant differentiation within the sexually reproducing population. The effective population size for one of the sexually reproducing subpopulations has recently declined to ~1,000, resulting in high levels of inbreeding. In particular, we found that 58% of the ecological niche overlapped between wild and cultivated populations and that there was extensive introgression into wild samples from cultivated populations. Interestingly, the introgression pattern and accumulation of genetic load may be influenced by the type of reproduction. In wild apomictic samples, the introgressed regions were primarily heterozygous, and genome-wide deleterious variants were hidden in the heterozygous state. In contrast, wild sexually reproducing samples carried a higher recessive deleterious burden. Furthermore, we also found that sexually reproducing samples were self-incompatible, which prevented the reduction of genetic diversity by selfing. Our population genomic analyses provide specific recommendations for distinct reproductive types and monitoring during conservation. This study highlights the genomic landscape of a wild relative of citrus and provides recommendations for the conservation of crop wild relatives.

Adaptive and maladaptive introgression in grapevine domestication.

Domesticated grapevines spread to Europe around 3,000 years ago. Previous studies have revealed genomic signals of introgression from wild to cultivated grapes in Europe, but the time, mode, genomic pattern, and biological effects of these introgression events have not been investigated. Here, we studied resequencing data from 345 samples spanning the distributional range of wild (Vitis vinifera ssp. sylvestris) and cultivated (V. vinifera ssp. vinifera) grapes. Based on machine learning-based population genetic analyses, we detected evidence for a single domestication of grapevine, followed by continuous gene flow between European wild grapes (EU) and cultivated grapes over the past ~2,000 y, especially from EU to wine grapes. We also inferred that soft-selective sweeps were the dominant signals of artificial selection. Gene pathways associated with the synthesis of aromatic compounds were enriched in regions that were both selected and introgressed, suggesting EU wild grapes were an important resource for improving the flavor of cultivated grapes. Despite the potential benefits of introgression in grape improvement, the introgressed fragments introduced a higher deleterious burden, with most deleterious SNPs and structural variants hidden in a heterozygous state. Cultivated wine grapes have benefited from adaptive introgression with wild grapes, but introgression has also increased the genetic load. In general, our study of beneficial and harmful effects of introgression is critical for genomic breeding of grapevine to take advantage of wild resources.

The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding.

Grapevine is one of the most economically important crops worldwide. However, the previous versions of the grapevine reference genome tipically consist of thousands of fragments with missing centromeres and telomeres, limiting the accessibility of the repetitive sequences, the centromeric and telomeric regions, and the study of inheritance of important agronomic traits in these regions. Here, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar PN40024 using PacBio HiFi long reads. The T2T reference genome (PN_T2T) is 69 Mb longer with 9018 more genes identified than the 12X.v0 version. We annotated 67% repetitive sequences, 19 centromeres and 36 telomeres, and incorporated gene annotations of previous versions into the PN_T2T assembly. We detected a total of 377 gene clusters, which showed associations with complex traits, such as aroma and disease resistance. Even though PN40024 derives from nine generations of selfing, we still found nine genomic hotspots of heterozygous sites associated with biological processes, such as the oxidation-reduction process and protein phosphorylation. The fully annotated complete reference genome therefore constitutes an important resource for grapevine genetic studies and breeding programs.

Whole-genome re-sequencing, diversity analysis, and stress-resistance analysis of 77 grape rootstock genotypes.

Introduction: Grape rootstocks play critical role in the development of the grape industry over the globe for their higher adaptability to various environments, and the evaluation of their genetic diversity among grape genotypes is necessary to the conservation and utility of genotypes. Methods: To analyze the genetic diversity of grape rootstocks for a better understanding multiple resistance traits, whole-genome re-sequencing of 77 common grape rootstock germplasms was conducted in the present study. Results: About 645 billion genome sequencing data were generated from the 77 grape rootstocks at an average depth of ~15.5×, based on which the phylogenic clusters were generated and the domestication of grapevine rootstocks was explored. The results indicated that the 77 rootstocks originated from five ancestral components. Through phylogenetic, principal components, and identity-by-descent (IBD) analyses, these 77 grape rootstocks were assembled into ten groups. It is noticed that the wild resources of V. amurensis and V. davidii, originating from China and being generally considered to have stronger resistance against biotic and abiotic stresses, were sub-divided from the other populations. Further analysis indicated that a high level of linkage disequilibrium was found among the 77 rootstock genotypes, and a total of 2,805,889 single nucleotide polymorphisms (SNPs) were excavated, GWAS analysis among the grape rootstocks located 631, 13, 9, 2, 810, and 44 SNP loci that were responsible to resistances to phylloxera, root-knot nematodes, salt, drought, cold and waterlogging traits. Discussion: This study generated a significant amount of genomic data from grape rootstocks, thus providing a theoretical basis for further research on the resistance mechanism of grape rootstocks and the breeding of resistant varieties. These findings also reveal that China originated V. amurensis and V. davidii could broaden the genetic background of grapevine rootstocks and be important germplasm used in breeding high stress-resistant grapevine rootstocks.

Structural variation and parallel evolution of apomixis in citrus during domestication and diversification.

Apomixis, or asexual seed formation, is prevalent in Citrinae via a mechanism termed nucellar or adventitious embryony. Here, multiple embryos of a maternal genotype form directly from nucellar cells in the ovule and can outcompete the developing zygotic embryo as they utilize the sexually derived endosperm for growth. Whilst nucellar embryony enables the propagation of clonal plants of maternal genetic constitution, it is also a barrier to effective breeding through hybridization. To address the genetics and evolution of apomixis in Citrinae, a chromosome-level genome of the Hongkong kumquat (Fortunella hindsii) was assembled following a genome-wide variation map including structural variants (SVs) based on 234 Citrinae accessions. This map revealed that hybrid citrus cultivars shelter genome-wide deleterious mutations and SVs into heterozygous states free from recessive selection, which may explain the capability of nucellar embryony in most cultivars during Citrinae diversification. Analyses revealed that parallel evolution may explain the repeated origin of apomixis in different genera of Citrinae. Within Fortunella, we found that apomixis of some varieties originated via introgression. In apomictic Fortunella, the locus associated with apomixis contains the FhRWP gene, encoding an RWP-RK domain-containing protein previously shown to be required for nucellar embryogenesis in Citrus. We found the heterozygous SV in the FhRWP and CitRWP promoters from apomictic Citrus and Fortunella, due to either two or three miniature inverted transposon element (MITE) insertions. A transcription factor, FhARID, encoding an AT-rich interaction domain-containing protein binds to the MITEs in the promoter of apomictic varieties, which facilitates induction of nucellar embryogenesis. This study provides evolutionary genomic and molecular insights into apomixis in Citrinae and has potential ramifications for citrus breeding.