The HN protein of Newcastle disease virus induces cell apoptosis through the induction of lysosomal membrane permeabilization.

Lysosomes are acidic organelles that mediate the degradation and recycling of cellular waste materials. Damage to lysosomes can cause lysosomal membrane permeabilization (LMP) and trigger different types of cell death, including apoptosis. Newcastle disease virus (NDV) can naturally infect most birds. Additionally, it serves as a promising oncolytic virus known for its effective infection of tumor cells and induction of intensive apoptotic responses. However, the involvement of lysosomes in NDV-induced apoptosis remains poorly understood. Here, we demonstrate that NDV infection profoundly triggers LMP, leading to the translocation of cathepsin B and D and subsequent mitochondria-dependent apoptosis in various tumor and avian cells. Notably, the released cathepsin B and D exacerbate NDV-induced LMP by inducing the generation of reactive oxygen species. Additionally, we uncover that the viral Hemagglutinin neuraminidase (HN) protein induces the deglycosylation and degradation of lysosome-associated membrane protein 1 (LAMP1) and LAMP2 dependent on its sialidase activity, which finally contributes to NDV-induced LMP and cellular apoptosis. Overall, our findings elucidate the role of LMP in NDV-induced cell apoptosis and provide novel insights into the function of HN during NDV-induced LMP, which provide innovative approaches for the development of NDV-based oncolytic agents.

The NP protein of Newcastle disease virus dictates its oncolytic activity by regulating viral mRNA translation efficiency.

Newcastle disease virus (NDV) has been extensively studied as a promising oncolytic virus for killing tumor cells in vitro and in vivo in clinical trials. However, the viral components that regulate the oncolytic activity of NDV remain incompletely understood. In this study, we systematically compared the replication ability of different NDV genotypes in various tumor cells and identified NP protein determines the oncolytic activity of NDV. On the one hand, NDV strains with phenylalanine (F) at the 450th amino acid position of the NP protein (450th-F-NP) exhibit a loss of oncolytic activity. This phenotype is predominantly associated with genotype VII NDVs. In contrast, the NP protein with a leucine amino acid at this site in other genotypes (450th-L-NP) can facilitate the loading of viral mRNA onto ribosomes more effectively than 450th-F-NP. On the other hand, the NP protein from NDV strains that exhibit strong oncogenicity interacts with eIF4A1 within its 366-489 amino acid region, leading to the inhibition of cellular mRNA translation with a complex 5' UTR structure. Our study provide mechanistic insights into how highly oncolytic NDV strains selectively promote the translation of viral mRNA and will also facilitate the screening of oncolytic strains for oncolytic therapy.

RBM14 inhibits the replication of porcine epidemic diarrhea virus by recruiting p62 to degrade nucleocapsid protein through the activation of autophagy and interferon pathway.

Porcine epidemic diarrhea virus (PEDV) results in PED, which is an infectious intestinal disease with the representative features of diarrhea, vomiting, and dehydration. PEDV infects neonatal piglets, causing high mortality rates. Therefore, elucidating the interaction between the virus and host in preventing and controlling PEDV infection is of immense significance. We found a new antiviral function of the host protein, RNA-binding motif protein 14 (RBM14), which can inhibit PEDV replication via the activation of autophagy and interferon (IFN) signal pathways. We found that RBM14 can recruit cargo receptor p62 to degrade PEDV nucleocapsid (N) protein through the RBM14-p62-autophagosome pathway. Furthermore, RBM14 can also improve the antiviral ability of the hosts through interacting with mitochondrial antiviral signaling protein to induce IFN expression. These results highlight the novel mechanism underlying RBM14-induced viral restriction. This mechanism leads to the degradation of viral N protein via the autophagy pathway and upregulates IFN for inhibiting PEDV replication; thus, offering new ways for preventing and controlling PED.IMPORTANCEPorcine epidemic diarrhea virus (PEDV) is a vital reason for diarrhea in neonatal piglets, which causes high morbidity and mortality rates. There is currently no effective vaccine or drug to treat and prevent infection with the PEDV. During virus infection, the host inhibits virus replication through various antiviral factors, and at the same time, the virus antagonizes the host's antiviral reaction through its own encoded protein, thus completing the process of virus replication. Our study has revealed that the expression of RNA-binding motif protein 14 (RBM14) was downregulated in PEDV infection. We found that RBM14 can recruit cargo receptor p62 to degrade PEDV N protein via the RBM14-p62-autophagosome pathway and interacted with mitochondrial antiviral signaling protein and TRAF3 to activate the interferon signal pathway, resulting in the inhibition of PEDV replication.

Influenza A virus infection activates caspase-8 to enhance innate antiviral immunity by cleaving CYLD and blocking TAK1 and RIG-I deubiquitination.

Caspase-8, an aspartate-specific cysteine protease that primarily functions as an initiator caspase to induce apoptosis, can downregulate innate immunity in part by cleaving RIPK1 and IRF3. However, patients with caspase-8 mutations or deficiency develop immunodeficiency and are prone to viral infections. The molecular mechanism underlying this controversy remains unknown. Whether caspase-8 enhances or suppresses antiviral responses against influenza A virus (IAV) infection remains to be determined. Here, we report that caspase-8 is readily activated in A549 and NL20 cells infected with the H5N1, H5N6, and H1N1 subtypes of IAV. Surprisingly, caspase-8 deficiency and two caspase-8 inhibitors, Z-VAD and Z-IETD, do not enhance but rather downregulate antiviral innate immunity, as evidenced by decreased TBK1, IRF3, IκBα, and p65 phosphorylation, decreased IL-6, IFN-ß, MX1, and ISG15 gene expression; and decreased IFN-ß production but increased virus replication. Mechanistically, caspase-8 cleaves and inactivates CYLD, a tumor suppressor that functions as a deubiquitinase. Caspase-8 inhibition suppresses CYLD cleavage, RIG-I and TAK1 ubiquitination, and innate immune signaling. In contrast, CYLD deficiency enhances IAV-induced RIG-I and TAK1 ubiquitination and innate antiviral immunity. Neither caspase-3 deficiency nor treatment with its inhibitor Z-DEVD affects CYLD cleavage or antiviral innate immunity. Our study provides evidence that caspase-8 activation in two human airway epithelial cell lines does not silence but rather enhances innate immunity by inactivating CYLD.

Influenza virus infection activates TAK1 to suppress RIPK3-independent apoptosis and RIPK1-dependent necroptosis.

Many DNA viruses develop various strategies to inhibit cell death to facilitate their replication. However, whether influenza A virus (IAV), a fast-replicating RNA virus, attenuates cell death remains unknown. Here, we report that IAV infection induces TAK1 phosphorylation in a murine alveolar epithelial cell line (LET1) and a murine fibroblastoma cell line (L929). The TAK1-specific inhibitor 5Z-7-Oxzeneonal (5Z) and TAK1 knockout significantly enhance IAV-induced apoptosis, as evidenced by increased PARP, caspase-8, and caspase-3 cleavage. TAK1 inhibition also increases necroptosis as evidenced by increased RIPK1S166, RIPK3T231/S232, and MLKLS345 phosphorylation. Mechanistically, TAK1 activates IKK, which phosphorylates RIPK1S25 and inhibits its activation. TAK1 also activates p38 and its downstream kinase MK2, which phosphorylates RIPK1S321 but does not affect RIPK1 activation. Further investigation revealed that the RIPK1 inhibitor Nec-1 and RIPK1 knockout abrogate IAV-induced apoptosis and necroptosis; re-expression of wild-type but not kinase-dead (KD)-RIPK1 restores IAV-induced cell death. ZBP1 knockout abrogates IAV-induced cell death, whereas RIPK3 knockout inhibits IAV-induced necroptosis but not apoptosis. 5Z treatment enhances IAV-induced cell death and slightly reduces the inflammatory response in the lungs of H1N1 virus-infected mice and prolongs the survival of IAV-infected mice. Our study provides evidence that IAV activates TAK1 to suppress RIPK1-dependent apoptosis and necroptosis, and that RIPK3 is required for IAV-induced necroptosis but not apoptosis in epithelial cells.

The haemagglutinin-neuraminidase protein of velogenic Newcastle disease virus enhances viral infection through NF-κB-mediated programmed cell death.

The haemagglutinin-neuraminidase (HN) protein, a vital membrane glycoprotein, plays a pivotal role in the pathogenesis of Newcastle disease virus (NDV). Previously, we demonstrated that a mutation in the HN protein is essential for the enhanced virulence of JS/7/05/Ch, a velogenic variant NDV strain originating from the mesogenic vaccine strain Mukteswar. Here, we explored the effects of the HN protein during viral infection in vitro using three viruses: JS/7/05/Ch, Mukteswar, and an HN-replacement chimeric NDV, JS/MukHN. Through microscopic observation, CCK-8, and LDH release assays, we demonstrated that compared with Mukteswar and JS/MukHN, JS/7/05/Ch intensified the cellular damage and mortality attributed to the mutant HN protein. Furthermore, JS/7/05/Ch induced greater levels of apoptosis, as evidenced by the activation of caspase-3/8/9. Moreover, JS/7/05/Ch promoted autophagy, leading to increased autophagosome formation and autophagic flux. Subsequent pharmacological experiments revealed that inhibition of apoptosis and autophagy significantly impacted virus replication and cell viability in the JS/7/05/Ch-infected group, whereas less significant effects were observed in the other two infected groups. Notably, the mutant HN protein enhanced JS/7/05/Ch-induced apoptosis and autophagy by suppressing NF-κB activation, while it mitigated the effects of NF-κB on NDV infection. Overall, our study offers novel insights into the mechanisms underlying the increased virulence of NDV and serves as a reference for the development of vaccines.

Amino acid mutations PB1-V719M and PA-N444D combined with PB2-627K contribute to the pathogenicity of H7N9 in mice.

H7N9 subtype avian influenza viruses (AIVs) cause 1567 human infections and have high mortality, posing a significant threat to public health. Previously, we reported that two avian-derived H7N9 isolates (A/chicken/Eastern China/JTC4/2013 and A/chicken/Eastern China/JTC11/2013) exhibit different pathogenicities in mice. To understand the genetic basis for the differences in virulence, we constructed a series of mutant viruses based on reverse genetics. We found that the PB2-E627K mutation alone was not sufficient to increase the virulence of H7N9 in mice, despite its ability to enhance polymerase activity in mammalian cells. However, combinations with PB1-V719M and/or PA-N444D mutations significantly enhanced H7N9 virulence. Additionally, these combined mutations augmented polymerase activity, thereby intensifying virus replication, inflammatory cytokine expression, and lung injury, ultimately increasing pathogenicity in mice. Overall, this study revealed that virulence in H7N9 is a polygenic trait and identified novel virulence-related residues (PB2-627K combined with PB1-719M and/or PA-444D) in viral ribonucleoprotein (vRNP) complexes. These findings provide new insights into the molecular mechanisms underlying AIV pathogenesis in mammals, with implications for pandemic preparedness and intervention strategies.

Phylogenomic and ecological analyses reveal the spatiotemporal evolution of global pines.

How coniferous forests evolved in the Northern Hemisphere remains largely unknown. Unlike most groups of organisms that generally follow a latitudinal diversity gradient, most conifer species in the Northern Hemisphere are distributed in mountainous areas at middle latitudes. It is of great interest to know whether the midlatitude region has been an evolutionary cradle or museum for conifers and how evolutionary and ecological factors have driven their spatiotemporal evolution. Here, we investigated the macroevolution of Pinus, the largest conifer genus and characteristic of northern temperate coniferous forests, based on nearly complete species sampling. Using 1,662 genes from transcriptome sequences, we reconstructed a robust species phylogeny and reestimated divergence times of global pines. We found that ∼90% of extant pine species originated in the Miocene in sharp contrast to the ancient origin of Pinus, indicating a Neogene rediversification. Surprisingly, species at middle latitudes are much older than those at other latitudes. This finding, coupled with net diversification rate analysis, indicates that the midlatitude region has provided an evolutionary museum for global pines. Analyses of 31 environmental variables, together with a comparison of evolutionary rates of niche and phenotypic traits with a net diversification rate, found that topography played a primary role in pine diversification, and the aridity index was decisive for the niche rate shift. Moreover, fire has forced diversification and adaptive evolution of Pinus Our study highlights the importance of integrating phylogenomic and ecological approaches to address evolution of biological groups at the global scale.

Cyclopeptide Inhibitors Target the N-Terminal Tail of STING and Alleviate Autoinflammation.

Cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) signaling pathway is a crucial component of innate immunity that plays a vital role in protecting against pathogen infections and cellular stress. However, aberrant activation of cGAS-STING pathway is related to inflammatory and autoimmune diseases. Here, we developed cyclopeptide STING inhibitors by cyclizing the N-terminal tail (NTT) of STING. These cyclopeptides selectively inhibited the activation of STING pathway in human or murine cell lines. Mechanistically, the inhibitors directly bound to STING, and subsequently blocked the aggregation and activation of STING. In addition, the optimal inhibitor STi-2 significantly suppressed proinflammatory cytokine production and systemic inflammation in Trex1-/- mice. Overall, our work facilitates the development of specific inhibitors of STING as potential therapies for cGAS-STING associated autoinflammatory diseases.

Phylotranscriptomic and ecological analyses reveal the evolution and morphological adaptation of Abies.

Coniferous forests are under severe threat of the rapid anthropogenic climate warming. Abies (firs), the fourth-largest conifer genus, is a keystone component of the boreal and temperate dark-coniferous forests and harbors a remarkably large number of relict taxa. However, the uncertainty of the phylogenetic and biogeographic history of Abies significantly impedes our prediction of future dynamics and efficient conservation of firs. In this study, using 1,533 nuclear genes generated from transcriptome sequencing and a complete sampling of all widely recognized species, we have successfully reconstructed a robust phylogeny of global firs, in which four clades are strongly supported and all intersectional relationships are resolved, although phylogenetic discordance caused mainly by incomplete lineage sorting and hybridization was detected. Molecular dating and ancestral area reconstruction suggest a Northern Hemisphere high-latitude origin of Abies during the Late Cretaceous, but all extant firs diversified during the Miocene to the Pleistocene, and multiple continental and intercontinental dispersals took place in response to the late Neogene climate cooling and orogenic movements. Notably, four critically endangered firs endemic to subtropical mountains of China, including A. beshanzuensis, A. ziyuanensis, A. fanjingshanensis and A. yuanbaoshanensis from east to west, have different origins and evolutionary histories. Moreover, three hotspots of species richness, including western North America, central Japan, and the Hengduan Mountains, were identified in Abies. Elevation and precipitation, particularly precipitation of the coldest quarter, are the most significant environmental factors driving the global distribution pattern of fir species diversity. Some morphological traits are evolutionarily constrained, and those linked to elevational variation (e.g., purple cone) and cold resistance (e.g., pubescent branch and resinous bud) may have contributed to the diversification of global firs. Our study sheds new light on the spatiotemporal evolution of global firs, which will be of great help to forest management and species conservation in a warming world.

Spatiotemporal Evolution of the Global Species Diversity of Rhododendron.

Evolutionary radiation is a widely recognized mode of species diversification, but its underlying mechanisms have not been unambiguously resolved for species-rich cosmopolitan plant genera. In particular, it remains largely unknown how biological and environmental factors have jointly driven its occurrence in specific regions. Here, we use Rhododendron, the largest genus of woody plants in the Northern Hemisphere, to investigate how geographic and climatic factors, as well as functional traits, worked together to trigger plant evolutionary radiations and shape the global patterns of species richness based on a solid species phylogeny. Using 3,437 orthologous nuclear genes, we reconstructed the first highly supported and dated backbone phylogeny of Rhododendron comprising 200 species that represent all subgenera, sections, and nearly all multispecies subsections, and found that most extant species originated by evolutionary radiations when the genus migrated southward from circumboreal areas to tropical/subtropical mountains, showing rapid increases of both net diversification rate and evolutionary rate of environmental factors in the Miocene. We also found that the geographically uneven diversification of Rhododendron led to a much higher diversity in Asia than in other continents, which was mainly driven by two environmental variables, that is, elevation range and annual precipitation, and were further strengthened by the adaptation of leaf functional traits. Our study provides a good example of integrating phylogenomic and ecological analyses in deciphering the mechanisms of plant evolutionary radiations, and sheds new light on how the intensification of the Asian monsoon has driven evolutionary radiations in large plant genera of the Himalaya-Hengduan Mountains.

Cellular vimentin regulates the infectivity of Newcastle disease virus through targeting of the HN protein.

The haemagglutinin-neuraminidase (HN) protein plays a crucial role in the infectivity and virulence of Newcastle disease virus (NDV). In a previous study, the mutant HN protein was identified as a crucial virulence factor for the velogenic variant NDV strain JS/7/05/Ch, which evolved from the prototypic vaccine strain Mukteswar. Furthermore, macrophages are the main susceptible target cells of NDV. However, the possible involvement of cellular molecules in viral infectivity remains unclear. Herein, we elucidate the crucial role of vimentin, an intermediate filament protein, in regulating NDV infectivity through targeting of the HN protein. Using LC‒MS/MS mass spectrometry and coimmunoprecipitation assays, we identified vimentin as a host protein that differentially interacted with prototypic and mutant HN proteins. Further analysis revealed that the variant NDV strain induced more significant rearrangement of vimentin fibres compared to the prototypic NDV strain and showed an interdependence between vimentin rearrangement and virus replication. Notably, these mutual influences were pronounced in HD11 chicken macrophages. Moreover, vimentin was required for multiple infection processes of the variant NDV strain in HD11 cells, including viral internalization, fusion, and release, while it was not necessary for those of the prototypic NDV strain. Collectively, these findings underscore the pivotal role of vimentin in NDV infection through targeting of the HN protein, providing novel targets for antiviral treatment strategies for NDV.

Reintroduction of highly pathogenic avian influenza A H7N9 virus in southwestern China.

Highly pathogenic (HP) avian influenza A H7N9 virus has emerged in China since 2016. In recent years, it has been most prevalent in northern China. However, several strains of HP H7N9 reappeared in southwestern China (Yunnan Province) in 2021. As a result, we are wondering if these viruses have re-emerged in situ or been reintroduced. Here, we present phylogenetic evidence that the HP H7N9 viruses isolated in Yunnan emigrated from northern to southwestern China in 2020. The northern subregion of China has become a novel epicenter in HP H7N9 dissemination. Meanwhile, a cleavage motif re-emerged due to the T341I mutation, implying a parallel evolution. This cross-region transmission, which originated in non-adjacent provinces and traveled a great geographic distance in an unknown way, indicates that HP H7N9 dissemination did not halt in 2020, even under the shadow of the COVID-19 pandemic. Additional surveillance studies in poultry are required to determine the HP H7N9 virus's geographic distribution and spread.

Monoclonal antibody targeting a novel linear epitope on nucleoprotein confers pan-reactivity to influenza A virus.

Nucleoprotein (NP) functions crucially in the replicative cycle of influenza A virus (IAV) via forming the ribonucleoprotein complex together with PB2, PB1, and PA proteins. As its high conservation, NP ranks one of the hot targets for design of universal diagnostic reagents and antiviral drugs for IAV. Here, we report an anti-NP murine monoclonal antibody (mAb) 5F10 prepared from traditional lymphocyte hybridoma technique with the immunogen of a clade 2.3.4.4 H5N1 subtype avian influenza virus. The specificity of mAb 5F10 to NP protein was confirmed by immunofluorescence assay and western blotting, and the mAb 5F10 could be used in immunoprecipitation and immunohistochemistry assays. Importantly, mAb 5F10 possessed broad-spectrum reactivity against H1~H11 subtypes of avian influenza viruses, including various HA clades of H5Nx subtype. In addition, mAb 5F10 also showed good affinity with H1N1 and H3N2 subtype influenza viruses of swine and human origin. Furthermore, the recognized antigenic epitope of mAb 5F10 was identified to consist of the conserved amino acid motif 81EHPSA85 in the second flexible loop region of NP protein through screening the phage display peptide library. Collectively, the mAb 5F10 which recognizes the novel universal NP linear B-cell epitope of IAV with diverse origins and subtypes will be a powerful tool for NP protein-based structural, functional, and mechanistic studies, as well as the development of detection methods and universal vaccines for IAV. KEY POINTS: • A broad-spectrum mAb against various subtypes and sources of IAV was developed • The mAb possessed good reactivity in IFA, western blot, IP, and IHC assays • The mAb targeted a novel conserved linear B-cell epitope involving 81EHPSA85 on NP protein.

Diploid species phylogeny and evolutionary reticulation indicate early radiation of Ephedra in the Tethys coast.

Reconstructing a robust species phylogeny and disentangling the evolutionary and biogeographic history of the gymnosperm genus Ephedra, which has a large genome and rich polyploids, remain a big challenge. Here we reconstructed a transcriptome-based phylogeny of 19 diploid Ephedra species, and explored evolutionary reticulations in this genus represented by 50 diploid and polyploid species, using four low-copy nuclear and nine plastid genes. The diploid species phylogeny indicates that the Mediterranean species diverged first, and the remaining species split into three clades, including the American species (Clade A), E. rhytidosperma, and all other Asian species (Clade B). The single-gene trees placed E. rhytidosperma sister to Clade A, Clade B, or Clades A + B in similar proportions, suggesting that radiation and gene flow likely occurred in the early evolution of Ephedra. In addition, reticulate evolution occurred not only among the deep nodes, but also in the recently evolved South American species, which further caused difficulty in phylogenetic reconstruction. Moreover, we found that allopolyploid speciation was pervasive in Ephedra. Our study also suggests that Ephedra very likely originated in the Tethys coast during the late Cretaceous, and the South American Ephedra species have a single origin by dispersal from Mexico or North America.

Zoonotic Threat of G4 Genotype Eurasian Avian-Like Swine Influenza A(H1N1) Viruses, China, 2020.

We investigated genetic and biologic characteristics of 2 Eurasian avian-like H1N1 swine influenza viruses from pigs in China that belong to the predominant G4 genotype. One swine isolate exhibited strikingly great homology to contemporaneous human Eurasian avian-like H1N1 isolates, preferential binding to the human-type receptor, and vigorous replication in mice without adaptation.

Identification and Molecular Mechanisms Study of Genes Associated with Osteoarthritis: A Comprehensive Bioinformatic Study of Cartilage and Synovium.

Osteoarthritis (OA) is a severe disease and has brought a massive burden to people's daily life. This study was performed to explore the hub genes associated with OA and predict potential biomarkers for OA. The cartilage (GSE114007) and synovial (GSE55235 and GSE55457) datasets downloaded from the Gene Expression Omnibus (GEO) database were used to screen the differentially expressed genes in OA compared with normal tissues. Then, based on weighted gene co-expression network analysis, the intersection genes between cartilage and synovium data were screened. The protein-protein interaction network and receiver operating characteristic (ROC) curve analysis were utilized to identify the OA-related hub genes. The gene ontology (GO), kyoto encyclopedia of genes and genomes (KEGG), and gene set variation analysis (GSVA) databases were used to explore the potential molecular mechanism of genes. Single-sample gene set enrichment analysis was performed to analyze the immune infiltration levels in OA synovium and cartilage tissues, respectively. A total of 131 intersection genes were screened. These genes were mainly enriched in the osteoclast differentiation and PI3K-Akt signaling pathways. Then, eight OA-related hub genes were further identified, including JUN, MYC, VEGFA, ATF3, NR4A1, BTG2, DUSP1, and JUNB. ROC curve analysis showed the area under the curve of these eight genes was > 0.6 in another OA dataset, suggesting their feasible capacity for predicting OA. Finally, we found that NR4A1 and BTG2 might be involved in multiple inflammatory responses of OA tissues. Our study identified some OA-related hub genes and revealed novel insights into the biological mechanism of OA, which provided a theoretical foundation for further experimental study.

Phylotranscriptomics reveals the evolutionary history of subtropical East Asian white pines: further insights into gymnosperm diversification.

Floristic composition within a geographic area is driven by a wide array of factors from local biotic interactions to biogeographical processes. Subtropical East Asia is a key biodiversity hotspot of the world, and harbors the most families of extant gymnosperms and a large number of endemic genera with ancient origins, but rare phylogenetic studies explored whether it served as a diversification center for gymnosperms. Here, we investigated the evolutionary and biogeographical history of subtropical East Asian white pines using an integrative approach that combines phylotranscriptomic and ecological analyses. Using 2,606 orthologous nuclear genes, we reconstructed a fully resolved and dated phylogeny of these species. Two main clades first diverged in the early Miocene, and by the late Miocene, all species appeared. Two white pines endemic to Taiwan Island experienced independent colonization events and regional extinction, which resulted in the present disjunctive distribution from mainland China. Ecological and biogeographical analyses indicate that the monsoon-driven assembly of evergreen broadleaved forests (EBLFs) might have significantly affected the diversification of subtropical East Asian white pines. Our study highlights the interactions of biotic and abiotic forces in the diversification and speciation of subtropical East Asian white pines. These findings indicate that subtropical East Asia is not only a floristic museum, but also a diversification center for gymnosperms. Our study also demonstrates the importance of phylotranscriptomics on species delimitation and biodiversity conservation, particularly for closely related species.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization.

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

New insights into the phylogeny and evolution of Podocarpaceae inferred from transcriptomic data.

Phylogenies of an increasing number of taxa have been resolved with the development of phylogenomics. However, the intergeneric relationships of Podocarpaceae, the second largest family of conifers comprising 19 genera and approximately 187 species mainly distributed in the Southern Hemisphere, have not been well disentangled in previous studies, even when genome-scale data sets were used. Here we used 993 nuclear orthologous groups (OGs) and 54 chloroplast OGs (genes), which were generated from 47 transcriptomes of Podocarpaceae and its sister group Araucariaceae, to reconstruct the phylogeny of Podocarpaceae. Our study completely resolved the intergeneric relationships of Podocarpaceae represented by all extant genera and revealed that topological conflicts among phylogenetic trees could be attributed to synonymous substitutions. Moreover, we found that two morphological traits, fleshy seed cones and flattened leaves, might be important for Podocarpaceae to adapt to angiosperm-dominated forests and thus could have promoted its species diversification. In addition, our results indicate that Podocarpaceae originated in Gondwana in the late Triassic and both vicariance and dispersal have contributed to its current biogeographic patterns. Our study provides the first robust transcriptome-based phylogeny of Podocarpaceae, an evolutionary framework important for future studies of this family.