Uncovering ghost introgression through genomic analysis of a distinct eastern Asian hickory species.

Ghost introgression, or the transfer of genetic material from extinct or unsampled lineages to sampled species, has attracted much attention. However, conclusive evidence for ghost introgression, especially in plant species, remains scarce. Here, we newly assembled chromosome-level genomes for both Carya sinensis and Carya cathayensis, and additionally re-sequenced the whole genomes of 43 C. sinensis individuals as well as 11 individuals representing 11 diploid hickory species. These genomic datasets were used to investigate the reticulation and bifurcation patterns within the genus Carya (Juglandaceae), with a particular focus on the beaked hickory C. sinensis. By combining the D-statistic and BPP methods, we obtained compelling evidence that supports the occurrence of ghost introgression in C. sinensis from an extinct ancestral hickory lineage. This conclusion was reinforced through the phylogenetic network analysis and a genome scan method VolcanoFinder, the latter of which can detect signatures of adaptive introgression from unknown donors. Our results not only dispel certain misconceptions about the phylogenetic history of C. sinensis but also further refine our understanding of Carya's biogeography via divergence estimates. Moreover, the successful integration of the D-statistic and BPP methods demonstrates their efficacy in facilitating a more precise identification of introgression types.

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information.

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

Genomic insights into biased allele loss and increased gene numbers after genome duplication in autotetraploid Cyclocarya paliurus.

BACKGROUND: Autopolyploidy is a valuable model for studying whole-genome duplication (WGD) without hybridization, yet little is known about the genomic structural and functional changes that occur in autopolyploids after WGD. Cyclocarya paliurus (Juglandaceae) is a natural diploid-autotetraploid species. We generated an allele-aware autotetraploid genome, a chimeric chromosome-level diploid genome, and whole-genome resequencing data for 106 autotetraploid individuals at an average depth of 60 × per individual, along with 12 diploid individuals at an average depth of 90 × per individual. RESULTS: Autotetraploid C. paliurus had 64 chromosomes clustered into 16 homologous groups, and the majority of homologous chromosomes demonstrated similar chromosome length, gene numbers, and expression. The regions of synteny, structural variation and nonalignment to the diploid genome accounted for 81.3%, 8.8% and 9.9% of the autotetraploid genome, respectively. Our analyses identified 20,626 genes (69.18%) with four alleles and 9191 genes (30.82%) with one, two, or three alleles, suggesting post-polyploid allelic loss. Genes with allelic loss were found to occur more often in proximity to or within structural variations and exhibited a marked overlap with transposable elements. Additionally, such genes showed a reduced tendency to interact with other genes. We also found 102 genes with more than four copies in the autotetraploid genome, and their expression levels were significantly higher than their diploid counterparts. These genes were enriched in enzymes involved in stress response and plant defense, potentially contributing to the evolutionary success of autotetraploids. Our population genomic analyses suggested a single origin of autotetraploids and recent divergence (~ 0.57 Mya) from diploids, with minimal interploidy admixture. CONCLUSIONS: Our results indicate the potential for genomic and functional reorganization, which may contribute to evolutionary success in autotetraploid C. paliurus.

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species.

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae).

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes, and its only congeneric species, P. strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics, and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Host insect specificity and interspecific competition drive parasitoid diversification in a plant-insect community.

Ecological interactions among plants, insect herbivores, and parasitoids are pervasive in nature and play important roles in community assembling, but the codiversification of tri-trophic interactions has received less attention. Here we compare pairwise codiversification patterns between a set of 22 fig species, their herbivorous pollinating and galling wasps, and their parasitoids. The parasitoid phylogeny showed significant congruence and more cospeciation events with host insects phylogeny than with host plants. These results suggest that parasitoid phylogeny and speciation is more closely related to their host insects than to their host plants. The pollinating wasps hosted more parasitoid species than gallers and indicated a more intense interspecific competition among parasitoids associated with pollinators. Closer matching and fewer evolutionary host shifts were found between parasitoids and galler hosts than between parasitoids and pollinator hosts. These results suggest that interspecific competition among parasitoids, rather than resource availability of host wasps, is the main driver of the codiversification pattern in this community. Therefore, our study highlights the important role of interspecific competition among high trophic level insects in plant-insect tri-trophic community assembling.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae).

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes.

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

Impact of Ghost Introgression on Coalescent-Based Species Tree Inference and Estimation of Divergence Time.

The species studied in any evolutionary investigation generally constitute a small proportion of all the species currently existing or that have gone extinct. It is therefore likely that introgression, which is widespread across the tree of life, involves "ghosts," that is, unsampled, unknown, or extinct lineages. However, the impact of ghost introgression on estimations of species trees has rarely been studied and is poorly understood. Here, we use mathematical analysis and simulations to examine the robustness of species tree methods based on the multispecies coalescent model to introgression from a ghost or extant lineage. We found that many results originally obtained for introgression between extant species can easily be extended to ghost introgression, such as the strongly interactive effects of incomplete lineage sorting (ILS) and introgression on the occurrence of anomalous gene trees. The relative performance of the summary species tree method (ASTRAL) and the full-likelihood method (*BEAST) varies under different introgression scenarios, with the former being more robust to gene flow between nonsister species, whereas the latter performing better under certain conditions of ghost introgression. When an outgroup ghost (defined as a lineage that diverged before the most basal species under investigation) acts as the donor of the introgressed genes, the time of root divergence among the investigated species generally was overestimated, whereas ingroup introgression, as commonly perceived, can only lead to underestimation. In many cases of ingroup introgression that may or may not involve ghost lineages, the stronger the ILS, the higher the accuracy achieved in estimating the time of root divergence, although the topology of the species tree is more prone to be biased by the effect of introgression. [Anomalous gene trees; divergence time; ghost introgression; multispecies coalescent; simulation; species tree.].

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs.

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost, a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.

Population-genomic analyses reveal bottlenecks and asymmetric introgression from Persian into iron walnut during domestication.

BACKGROUND: Persian walnut, Juglans regia, occurs naturally from Greece to western China, while its closest relative, the iron walnut, Juglans sigillata, is endemic in southwest China; both species are cultivated for their nuts and wood. Here, we infer their demographic histories and the time and direction of possible hybridization and introgression between them. RESULTS: We use whole-genome resequencing data, different population-genetic approaches (PSMC and GONE), and isolation-with-migration models (IMa3) on individuals from Europe, Iran, Kazakhstan, Pakistan, and China. IMa3 analyses indicate that the two species diverged from each other by 0.85 million years ago, with unidirectional gene flow from eastern J. regia and its ancestor into J. sigillata, including the shell-thickness gene. Within J. regia, a western group, located from Europe to Iran, and an eastern group with individuals from northern China, experienced dramatically declining population sizes about 80 generations ago (roughly 2400 to 4000 years), followed by an expansion at about 40 generations, while J. sigillata had a constant population size from about 100 to 20 generations ago, followed by a rapid decline. CONCLUSIONS: Both J. regia and J. sigillata appear to have suffered sudden population declines during their domestication, suggesting that the bottleneck scenario of plant domestication may well apply in at least some perennial crop species. Introgression from introduced J. regia appears to have played a role in the domestication of J. sigillata.

Habitat effects on reproductive phenotype, pollinator behavior, fecundity, and mating outcomes of a bumble bee-pollinated herb.

PREMISE: Fecundity and mating outcomes commonly differ among plant populations occupying contrasting environments. If self-pollination occurs primarily among flowers within plants, contrasting reproductive outcomes among populations must reflect environmental effects on plant-pollinator interactions. Specifically, local conditions could affect features of plant phenotypes that influence pollinator behavior, in turn modifying plant reproductive outcomes. METHODS: We compared phenotypes, pollinator abundance and behavior, and female fecundity and mating in two meadow populations and two forest populations of Aconitum kusnezoffii within 3 km of each other. Mating outcomes were assessed using microsatellites. RESULTS: Meadow plants generally produced more, shorter ramets with more, larger flowers, but less nectar per flower than forest plants. These differences likely largely represent phenotypic plasticity. Individual bumble bees visited more flowers on forest plants, likely because the more abundant bees in the meadows depleted nectar availability, as indicated by briefer visits to individual flowers. Despite similar fruit set in both habitats, forest plants set more seeds per fruit. Nevertheless, meadow plants produced more seeds overall, owing to sevenfold greater flower production. Consistent with individual bees visiting fewer flowers on meadow plants, more of their seeds were outcrossed. However, the outcrossed seeds of forest plants included more male mates. CONCLUSIONS: Reproductive outcomes can vary among populations of animal-pollinated plants as a result of differences in the availability of effective pollinators and environmental effects on plant phenotypes, and their functional consequences for pollinator behavior that governs pollen dispersal.

The rediscovery of Caryapoilanei (Juglandaceae) after 63 years reveals a new record from China.

Despite having been first published in 1941, Caryapoilanei (A.Chev.) J.-F.Leroy is only known from three collections in Vietnam, Laos and Thailand. It has not been recollected since then and was long suspected to have become extinct through repeated deforestation events. Here, we report the rediscovery, and meanwhile the first new record in China, of this extremely rare gigantic hickory species at Yunnan province 63 years after its last collection in 1958. Three small patchy subpopulations were found with a total of about 50 adult trees having diameter at breast height (DBH) larger than 60 cm, together with some seedlings and saplings, but the fruit set was low. Based on new and fresh material, we present a revised morphological description of C.poilanei, and an updated distribution map for the species. In addition, we also provide a key for the hickories in China. Lastly, we suggest C.poilanei should be listed as a Critically Endangered (CR) species according to the latest IUCN Red List Categories and Criteria.

Dead-End Hybridization in Walnut Trees Revealed by Large-Scale Genomic Sequence Data.

Although hybridization plays a large role in speciation, some unknown fraction of hybrid individuals never reproduces, instead remaining as genetic dead-ends. We investigated a morphologically distinct and culturally important Chinese walnut, Juglans hopeiensis, suspected to have arisen from hybridization of Persian walnut (J. regia) with Asian butternuts (J. cathayensis, J. mandshurica, and hybrids between J. cathayensis and J. mandshurica). Based on 151 whole-genome sequences of the relevant taxa, we discovered that all J. hopeiensis individuals are first-generation hybrids, with the time for the onset of gene flow estimated as 370,000 years, implying both strong postzygotic barriers and the presence of J. regia in China by that time. Six inversion regions enriched for genes associated with pollen germination and pollen tube growth may be involved in the postzygotic barriers that prevent sexual reproduction in the hybrids. Despite its long-recurrent origination and distinct traits, J. hopeiensis does not appear on the way to speciation.

Different rates of pollen and seed gene flow cause branch-length and geographic cytonuclear discordance within Asian butternuts.

Topological cytonuclear discordance is commonly observed in plant phylogenetic and phylogeographic studies, yet few studies have attempted to detect two other forms of cytonuclear discordance (branch length and geographical) and to uncover the causes of the discordance. We used the whole nuclear and chloroplast genome data from 80 individual Asian butternuts to reveal the pattern and processes of cytonuclear discordance. Our findings indicate that the chloroplast genome had substantially deeper divergence (branch-length discordance) and a steeper cline in the contact zone (geographic discordance) compared with the nuclear genome. After various hypothesis have been tested, the results suggest that incomplete lineage sorting, positive selection and cytonuclear incompatibility are probably insufficient to explain this pattern. However, isolation-by-distance analysis and gene flow estimation point to a much higher level of gene flow by pollen compared with by seeds, which may have slowed down lineage divergence and mediated wider contact for nuclear genome compared with the chloroplast genome. Altogether, this study highlights a critical role of sex-biased dispersal in causing discordance between the nuclear and plastid genome of Asian butternuts. Given its ubiquity among plants, asymmetric gene flow should be given a high priority in future studies of cytonuclear discordance.

TcpC inhibits neutrophil extracellular trap formation by enhancing ubiquitination mediated degradation of peptidylarginine deiminase 4.

TcpC is a multifunctional virulence factor of uropathogenic E. coli (UPEC). Neutrophil extracellular trap formation (NETosis) is a crucial anti-infection mechanism of neutrophils. Here we show the influence of TcpC on NETosis and related mechanisms. We show NETosis in the context of a pyelonephritis mouse model induced by TcpC-secreting wild-type E. coli CFT073 (CFT073wt) and LPS-induced in vitro NETosis with CFT073wt or recombinant TcpC (rTcpC)-treated neutrophils are inhibited. rTcpC enters neutrophils through caveolin-mediated endocytosis and inhibits LPS-induced production of ROS, proinflammatory cytokines and protein but not mRNA levels of peptidylarginine deiminase 4 (PAD4). rTcpC treatment enhances PAD4 ubiquitination and accumulation in proteasomes. Moreover, in vitro ubiquitination kit analyses show that TcpC is a PAD4-targetd E3 ubiquitin-ligase. These data suggest that TcpC inhibits NETosis primarily by serving as an E3 ligase that promotes degradation of PAD4. Our findings provide a novel mechanism underlying TcpC-mediated innate immune evasion.

Mechanisms of cancer stem cell senescence: Current understanding and future perspectives.

Cancer stem cells (CSCs) are a small population of heterogeneous tumor cells with the capacity of self-renewal and aberrant differentiation for immortality and divergent lineages of cancer cells. In contrast to bulky tumor cells, CSCs remain less differentiated and resistant to therapy even when targeted with tissue-specific antigenic markers. This makes CSCs responsible for not only tumor initiation, development, but also tumor recurrence. Emerging evidence suggests that CSCs can undergo cell senescence, a non-proliferative state of cells in response to stress. While cell senescence attenuates tumor cell proliferation, it is commonly regarded as a tumor suppressive mechanism. However, mounting research indicates that CSC senescence also provides these cells with the capacity to evade cytotoxic effects from cancer therapy, exacerbating cancer relapse and metastasis. Recent studies demonstrate that senescence drives reprogramming of cancer cell toward stemness and promotes CSC generation. In this review, we highlight the origin, heterogeneity and senescence regulatory mechanisms of CSCs, the complex relationship between CSC senescence and tumor therapy, and the recent beneficial effects of senotherapy on eliminating senescent tumor cells.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers.

BACKGROUND: Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM. RESULTS: This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (pGSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM. CONCLUSIONS: This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Reducing metal/graphene contact resistance via N, N-dimethylacetamide-assisted clean fabrication process.

Contact resistance (RC) is of great importance for radio frequency (RF) applications of graphene, especially graphene field effect transistors (FETs) with short channel. FETs and transmission line model test structures based on chemical vapor deposition grown graphene are fabricated. The effects of employing traditional lithography solvent (Acetone) and strong solvents for photo resist, such as N, N-Dimethylacetamide (ZDMAC) and N-Methyl pyrrolidone (NMP), are systematically investigated. It was found that ZDMAC and NMP have more proficiency than acetone to remove the photo-resist residues and contaminations attached on graphene surface, enabling clean surface of graphene. However, strong solvents are found to destroy the lattice structure of graphene channel and induce defects in graphene lattice. Clean surface contributes to a significant reduction in theRCbetween graphene channel and metal electrode, and the defects introduced on graphene surface underneath metal electrodes also contribute the reduction ofRC. But defects and deformation of lattice will increase the resistance in graphene channel and lead to the compromise of device performance. To address this problem, a mix wet-chemical approach employing both acetone and ZDMAC was developed in our study to realize a 19.07% reduction ofRC, without an unacceptable mass production of defects.