Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species.

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Association of PD-1/PD-L1 expression and Epstein--Barr virus infection in patients with invasive breast cancer.

PURPOSE: Causative factors of breast cancer include infections, such as Epstein-Barr virus (EBV) infection. The aim of this study was to analyze the clinicopathological features of EBV-positive (IBC) and determine if EBV affects programmed cell death receptor 1 (PD-1)/PD ligand 1 (PD-L1) expression in IBC, similar to other EBV-infected tumors with PD-L1/PD-1 expression. METHODS: We collected 140 samples of IBC tissues and 25 samples of adjacent tissues. All patients were followed-up by telephone from the day of surgery to December 2020. Chromogenic in-situ hybridization was performed to evaluate EBV-encoded RNA (EBER). Immunohistochemistry was performed to evaluate PD-L1 and PD-1 expressions. The correlation between PD1/PDL1 expression and clinicopathological features was also analyzed. RESULTS: EBER was detected in 57 of 140 (40.7%) IBC tissues and not detected in any adjacent tissue (P < 0.05). Clinicopathologic features of patients were consistent with EBV-associated IBC. EBV infection was correlated with the mass size, menopausal status, axillary lymph node metastasis, vascular invasion, Ki-67 index, clinical stage, and estrogen receptor and progesterone receptor expressions (all P < 0.05), but not with the histological type, invasive ductal carcinoma histological grade, or human epidermal growth factor receptor 2 (HER2) expression (all P > 0.05). The positive rate of PD-1/PD-L1 expression was higher in the EBV-positive group than in the EBV-negative group (P < 0.05). The Kaplan-Meier univariate survival analysis showed that EBV was associated with poor disease-free survival and overall survival in patients with IBC. PD-L1/PD-1 expression could predict a poor prognosis. CONCLUSIONS: In this study, clinicopathologic characteristics of patients were consistent with EBV-infected IBC. Patients with EBV-positive breast cancer were more likely to have elevated PD-1/PDL-1 expression compared to those with EBV-negative breast cancer. This finding could serve as a basis to explore therapeutic targets, particularly immunotherapy, for patients with IBC.

Comparative Genomics and Phylogenomics of East Asian Tulips (Amana, Liliaceae).

The genus Amana Honda (Liliaceae), when it is treated as separate from Tulipa, comprises six perennial herbaceous species that are restricted to China, Japan and the Korean Peninsula. Although all six Amana species have important medicinal and horticultural uses, studies focused on species identification and molecular phylogenetics are few. Here we report the nucleotide sequences of six complete Amana chloroplast (cp) genomes. The cp genomes of Amana range from 150,613 bp to 151,136 bp in length, all including a pair of inverted repeats (25,629-25,859 bp) separated by the large single-copy (81,482-82,218 bp) and small single-copy (17,366-17,465 bp) regions. Each cp genome equivalently contains 112 unique genes consisting of 30 transfer RNA genes, four ribosomal RNA genes, and 78 protein coding genes. Gene content, gene order, AT content, and IR/SC boundary structure are nearly identical among all Amana cp genomes. However, the relative contraction and expansion of the IR/SC borders among the six Amana cp genomes results in length variation among them. Simple sequence repeat (SSR) analyses of these Amana cp genomes indicate that the richest SSRs are A/T mononucleotides. The number of repeats among the six Amana species varies from 54 (A. anhuiensis) to 69 (Amana kuocangshanica) with palindromic (28-35) and forward repeats (23-30) as the most common types. Phylogenomic analyses based on these complete cp genomes and 74 common protein-coding genes strongly support the monophyly of the genus, and a sister relationship between Amana and Erythronium, rather than a shared common ancestor with Tulipa. Nine DNA markers (rps15-ycf1, accD-psaI, petA-psbJ, rpl32-trnL, atpH-atpI, petD-rpoA, trnS-trnG, psbM-trnD, and ycf4-cemA) with number of variable sites greater than 0.9% were identified, and these may be useful for future population genetic and phylogeographic studies of Amana species.