Utilizing systems genetics to enhance understanding into molecular targets of skin cancer.

Despite progress made with immune checkpoint inhibitors and targeted therapies, skin cancer remains a significant public health concern in the United States. The intricacies of the disease, encompassing genetics, immune responses, and external factors, call for a comprehensive approach. Techniques in systems genetics, including transcriptional correlation analysis, functional pathway enrichment analysis, and protein-protein interaction network analysis, prove valuable in deciphering intricate molecular mechanisms and identifying potential diagnostic and therapeutic targets for skin cancer. Recent studies demonstrate the efficacy of these techniques in uncovering molecular processes and pinpointing diagnostic markers for various skin cancer types, highlighting the potential of systems genetics in advancing innovative therapies. While certain limitations exist, such as generalizability and contextualization of external factors, the ongoing progress in AI technologies provides hope in overcoming these challenges. By providing protocols and a practical example involving Braf, we aim to inspire early-career experimental dermatologists to adopt these tools and seamlessly integrate these techniques into their skin cancer research, positioning them at the forefront of innovative approaches in combating this devastating disease.

Utilizing a human TLR selective ligand in a humanized immune system mouse model to investigate human TLR4 signaling.

Mouse models with humanized immune systems are becoming increasingly prevalent in pharmaceutical research as a platform for preclinical testing with potential for greater translatability to clinical applications. However, the presence of both mouse and human cells that respond to TLR ligands poses a challenge for investigating therapeutic modalities targeting TLR signaling. AZ617 is a human TLR4 agonist, which has been shown in vitro to preferentially induce human cytokines via the TLR4 signaling pathway. We sought to examine the ability of AZ617 to preferentially induce human cytokines in CD34+ stem cell-engrafted NOG-EXL mice (huNOG-EXL), to determine its suitability as an in vivo human functional readout. AZ617 elicited a strong human TNFα and IL-6 response in vivo that demonstrated a 10- and 5-fold preference, respectively, over the mouse TNFα and IL-6. To assess efficacy of inhibiting a key protein in the TLR4 signaling pathway, PF-06650833, a small molecule inhibitor of IRAK4, was used as a tool molecule. PF-0660833 was found to effectively inhibit AZ617-induced human TNFα release in vitro. Likewise, PF-06650833 reduced AZ617-induced human TNFα in the huNOG-EXL mouse model, with a weaker effect on human IL-6. A longitudinal study tracking functionality of monocytes revealed that the ability of monocytes to respond to ex vivo stimuli was increased by 21 weeks after engraftment. Taken together, our data suggests that human selective TLR ligands could preferentially drive cytokine production from human cells in huNOG-EXL mice. This model will allow for investigation of pharmacological inhibition of human TLR signaling pathways in an in vivo model system.

Comprehensive phylogenomic time tree of bryophytes reveals deep relationships and uncovers gene incongruences in the last 500 million years of diversification.

PREMISE: Bryophytes form a major component of terrestrial plant biomass, structuring ecological communities in all biomes. Our understanding of the evolutionary history of hornworts, liverworts, and mosses has been significantly reshaped by inferences from molecular data, which have highlighted extensive homoplasy in various traits and repeated bursts of diversification. However, the timing of key events in the phylogeny, patterns, and processes of diversification across bryophytes remain unclear. METHODS: Using the GoFlag probe set, we sequenced 405 exons representing 228 nuclear genes for 531 species from 52 of the 54 orders of bryophytes. We inferred the species phylogeny from gene tree analyses using concatenated and coalescence approaches, assessed gene conflict, and estimated the timing of divergences based on 29 fossil calibrations. RESULTS: The phylogeny resolves many relationships across the bryophytes, enabling us to resurrect five liverwort orders and recognize three more and propose 10 new orders of mosses. Most orders originated in the Jurassic and diversified in the Cretaceous or later. The phylogenomic data also highlight topological conflict in parts of the tree, suggesting complex processes of diversification that cannot be adequately captured in a single gene-tree topology. CONCLUSIONS: We sampled hundreds of loci across a broad phylogenetic spectrum spanning at least 450 Ma of evolution; these data resolved many of the critical nodes of the diversification of bryophytes. The data also highlight the need to explore the mechanisms underlying the phylogenetic ambiguity at specific nodes. The phylogenomic data provide an expandable framework toward reconstructing a comprehensive phylogeny of this important group of plants.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions.

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

Plastid phylogenomics and molecular evolution of Thismiaceae (Dioscoreales).

PREMISE: Species in Thismiaceae can no longer photosynthesize and instead obtain carbon from soil fungi. Here we infer Thismiaceae phylogeny using plastid genome data and characterize the molecular evolution of this genome. METHODS: We assembled five Thismiaceae plastid genomes from genome skimming data, adding to previously published data for phylogenomic inference. We investigated plastid-genome structural changes, considering locally colinear blocks (LCBs). We also characterized possible shifts in selection pressure in retained genes by considering changes in the ratio of nonsynonymous to synonymous changes (ω). RESULTS: Thismiaceae experienced two major pulses of gene loss around the early diversification of the family, with subsequent scattered gene losses across descendent lineages. In addition to massive size reduction, Thismiaceae plastid genomes experienced occasional inversions, and there were likely two independent losses of the plastid inverted repeat (IR) region. Retained plastid genes remain under generally strong purifying selection (ω << 1), with significant and sporadic weakening or strengthening in several instances. The bifunctional trnE-UUC gene of Thismia huangii may retain a secondary role in heme biosynthesis, despite a probable loss of functionality in protein translation. Several cis-spliced group IIA introns have been retained, despite the loss of the plastid intron maturase, matK. CONCLUSIONS: We infer that most gene losses in Thismiaceae occurred early and rapidly, following the initial loss of photosynthesis in its stem lineage. As a species-rich, fully mycoheterotrophic lineage, Thismiaceae provide a model system for uncovering the unique and divergent ways in which plastid genomes evolve in heterotrophic plants.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages.

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes.

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

Epithelial dysfunction is prevented by IL-22 treatment in a Citrobacter rodentium-induced colitis model that shares similarities with inflammatory bowel disease.

Inflammatory bowel disease (IBD) is characterized by a dysregulated intestinal epithelial barrier leading to breach of barrier immunity. Here we identified similar protein expression changes between IBD and Citrobacter rodentium-infected FVB mice with respect to dysregulation of solute transporters as well as components critical for intestinal barrier integrity. We attribute the disease associated changes in the model to the emergence of undifferentiated intermediate intestinal epithelial cells. Prophylactic treatment with IL-22.Fc in C. rodentium-infected FVB mice reduced disease severity and rescued the mice from lethality. Multi-omics and solute analyses revealed that IL-22.Fc treatment prevented disease-associated changes including disruption of the solute transporter machinery and restored proper physiological functions of the intestine, respectively. Taken together, we established the disease relevance of the C. rodentium-induced colitis model to IBD, demonstrated the protective role of IL-22 in amelioration of epithelial dysfunction and elucidated the molecular mechanisms with IL-22's effect on intestinal epithelial cells.

Phylogenomics and plastome evolution of a Brazilian mycoheterotrophic orchid, Pogoniopsis schenckii.

PREMISE: Pogoniopsis likely represents an independent photosynthesis loss in orchids. We use phylogenomic data to better identify the phylogenetic placement of this fully mycoheterotrophic taxon, and investigate its molecular evolution. METHODS: We performed likelihood analysis of plastid and mitochondrial phylogenomic data to localize the position of Pogoniopsis schenckii in orchid phylogeny, and investigated the evolution of its plastid genome. RESULTS: All analyses place Pogoniopsis in subfamily Epidendroideae, with strongest support from mitochondrial data, which also place it near tribe Sobralieae with moderately strong support. Extreme rate elevation in Pogoniopsis plastid genes broadly depresses branch support; in contrast, mitochondrial genes are only mildly rate elevated and display very modest and localized reductions in bootstrap support. Despite considerable genome reduction, including loss of photosynthesis genes and multiple translation apparatus genes, gene order in Pogoniopsis plastomes is identical to related autotrophs, apart from moderately shifted inverted repeat (IR) boundaries. All cis-spliced introns have been lost in retained genes. Two plastid genes (accD, rpl2) show significant strengthening of purifying selection. A retained plastid tRNA gene (trnE-UUC) of Pogoniopsis lacks an anticodon; we predict that it no longer functions in translation but retains a secondary role in heme biosynthesis. CONCLUSIONS: Slowly evolving mitochondrial genes clarify the placement of Pogoniopsis in orchid phylogeny, a strong contrast with analysis of rate-elevated plastome data. We documented the effects of the novel loss of photosynthesis: for example, despite massive gene loss, its plastome is fully colinear with other orchids, and it displays only moderate shifts in selective pressure in retained genes.

Mitochondrial genomic data are effective at placing mycoheterotrophic lineages in plant phylogeny.

Fully mycoheterotrophic plants can be difficult to place in plant phylogeny due to elevated substitution rates associated with photosynthesis loss. This potentially limits the effectiveness of downstream analyses of mycoheterotrophy that depend on accurate phylogenetic inference. Although mitochondrial genomic data sets are rarely used in plant phylogenetics, theory predicts that they should be resilient to long-branch artefacts, thanks to their generally slow evolution, coupled with limited rate elevation in heterotrophs. We examined the utility of mitochondrial genomes for resolving contentious higher-order placements of mycoheterotrophic lineages in two test cases: monocots (focusing on Dioscoreales) and Ericaceae. We find Thismiaceae to be distantly related to Burmanniaceae in the monocot order Dioscoreales, conflicting with current classification schemes based on few gene data sets. We confirm that the unusual Afrothismia is related to Taccaceae-Thismiaceae, with a corresponding independent loss of photosynthesis. In Ericaceae we recovered the first well supported relationships among its five major lineages: mycoheterotrophic Ericaceae are not monophyletic, as pyroloids are inferred to be sister to core Ericaceae, and monotropoids to arbutoids. Genes recovered from mitochondrial genomes collectively resolved previously ambiguous mycoheterotroph higher-order relationships. We propose that mitochondrial genomic data should be considered in standardised gene panels for inferring overall plant phylogeny.

Incomplete lineage sorting and local extinction shaped the complex evolutionary history of the Paleogene relict conifer genus, Chamaecyparis (Cupressaceae).

Inferring accurate biogeographic history of plant taxa with an East Asia (EA)-North America (NA) is usually hindered by conflicting phylogenies and a poor fossil record. The current distribution of Chamaecyparis (false cypress; Cupressaceae) with four species in EA, and one each in western and eastern NA, and its relatively rich fossil record, make it an excellent model for studying the EA-NA disjunction. Here we reconstruct phylogenomic relationships within Chamaecyparis using > 1400 homologous nuclear and 61 plastid genes. Our phylogenomic analyses using concatenated and coalescent approaches revealed strong cytonuclear discordance and conflicting topologies between nuclear gene trees. Incomplete lineage sorting (ILS) and hybridization are possible explanations of conflict; however, our coalescent analyses and simulations suggest that ILS is the major contributor to the observed phylogenetic discrepancies. Based on a well-resolved species tree and four fossil calibrations, the crown lineage of Chamaecyparis is estimated to have originated in the upper Cretaceous, followed by diversification events in the early and middle Paleogene. Ancestral area reconstructions suggest that Chamaecyparis had an ancestral range spanning both EA and NA. Fossil records further indicate that this genus is a relict of the "boreotropical" flora, and that local extinctions of European species were caused by global cooling. Overall, our results unravel a complex evolutionary history of a Paleogene relict conifer genus, which may have involved ILS, hybridization and the extinction of local species.

Macroevolutionary dynamics in the transition of angiosperms to aquatic environments.

Angiosperm lineages in aquatic environments are characterized by high structural and functional diversity, and wide distributions. A long-standing evolutionary riddle is what processes have caused the relatively low diversity of aquatic angiosperms compared to their terrestrial relatives. We use diversification and ancestral reconstruction models with a comprehensive > 10 000 genus angiosperm phylogeny to elucidate the macroevolutionary dynamics associated with transitions of terrestrial plants to water. Our study reveals that net diversification rates are significantly lower in aquatic than in terrestrial angiosperms due to lower speciation and higher extinction. Shifts from land to water started early in angiosperm evolution, but most events were concentrated during the last c. 25 million years. Reversals to a terrestrial habitat started only 40 million years ago, but occurred at much higher rates. Within aquatic angiosperms, the estimated pattern is one of gradual accumulation of lineages, and relatively low and constant diversification rates throughout the Cenozoic. Low diversification rates, together with infrequent water transitions, account for the low diversity of aquatic angiosperms today. The stressful conditions and small global surface of the aquatic habitat available for angiosperms are hypothesized to explain this pattern.

Meta-analysis of salt marsh vegetation impacts and recovery: a synthesis following the Deepwater Horizon oil spill.

Marine oil spills continue to be a global issue, heightened by spill events such as the 2010 Deepwater Horizon spill in the Gulf of Mexico, the largest marine oil spill in US waters and among the largest worldwide, affecting over 1,000 km of sensitive wetland shorelines, primarily salt marshes supporting numerous ecosystem functions. To synthesize the effects of the oil spill on foundational vegetation species in the salt marsh ecosystem, Spartina alterniflora and Juncus roemerianus, we performed a meta-analysis using data from 10 studies and 255 sampling sites over seven years post-spill. We examined the hypotheses that the oil spill reduced plant cover, stem density, vegetation height, aboveground biomass, and belowground biomass, and tracked the degree of effects temporally to estimate recovery time frames. All plant metrics indicated impacts from oiling, with 20-100% maximum reductions depending on oiling level and marsh zone. Peak reductions of ~70-90% in total plant cover, total aboveground biomass, and belowground biomass were observed for heavily oiled sites at the marsh edge. Both Spartina and Juncus were impacted, with Juncus affected to a greater degree. Most plant metrics had recovery time frames of three years or longer, including multiple metrics with incomplete recovery over the duration of our data, at least seven years post-spill. Belowground biomass was particularly concerning, because it declined over time in contrast with recovery trends in most aboveground metrics, serving as a strong indicator of ongoing impact, limited recovery, and impaired resilience. We conclude that the Deepwater Horizon spill had multiyear impacts on salt marsh vegetation, with full recovery likely to exceed 10 years, particularly in heavily oiled marshes, where erosion may preclude full recovery. Vegetation impacts and delayed recovery is likely to have exerted substantial influences on ecosystem processes and associated species, especially along heavily oiled shorelines. Our synthesis affords a greater understanding of ecosystem impacts and recovery following the Deepwater Horizon oil spill, and informs environmental impact analysis, contingency planning, emergency response, damage assessment, and restoration efforts related to oil spills.

A new carnivorous plant lineage (Triantha) with a unique sticky-inflorescence trap.

Carnivorous plants consume animals for mineral nutrients that enhance growth and reproduction in nutrient-poor environments. Here, we report that Triantha occidentalis (Tofieldiaceae) represents a previously overlooked carnivorous lineage that captures insects on sticky inflorescences. Field experiments, isotopic data, and mixing models demonstrate significant N transfer from prey to Triantha, with an estimated 64% of leaf N obtained from prey capture in previous years, comparable to levels inferred for the cooccurring round-leaved sundew, a recognized carnivore. N obtained via carnivory is exported from the inflorescence and developing fruits and may ultimately be transferred to next year's leaves. Glandular hairs on flowering stems secrete phosphatase, as seen in all carnivorous plants that directly digest prey. Triantha is unique among carnivorous plants in capturing prey solely with sticky traps adjacent to its flowers, contrary to theory. However, its glandular hairs capture only small insects, unlike the large bees and butterflies that act as pollinators, which may minimize the conflict between carnivory and pollination.

Investigating Target Gene Function in a CD40 Agonistic Antibody-induced Colitis Model using CRISPR/Cas9-based Technologies.

The immune system functions to defend humans against foreign invaders such as bacteria and viruses. However, disorders of the immune system may lead to autoimmunity, inflammatory disease, and cancer. The inflammatory bowel diseases (IBD)-Crohn's disease (CD) and ulcerative colitis (UC)-are chronic diseases marked by relapsing intestinal inflammation. Although IBD is most prevalent in Western countries (1 in 1,000), incident rates are increasing around the world. Through association studies, researchers have linked hundreds of genes to the pathology of IBD. However, the elaborate pathology behind IBD and the high number of potential genes pose significant challenges in finding the best therapeutic targets. Additionally, the tools needed to functionally characterize each genetic association introduce many rate-limiting factors such as the generation of genetically modified mice for each gene. To investigate the therapeutic potential of target genes, a model system has been developed using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated endonuclease (Cas9)-based technologies and a cluster of differentiation 40 (CD40) agonistic antibody. The present study shows that CRISPR/Cas9-mediated editing in the immune system can be used to investigate the impact of genes in vivo. Limited to the hematopoietic compartment, this approach reliably edits the resulting reconstituted immune system. CRISPR/Cas9-edited mice are generated faster and are far less expensive than traditional genetically modified mice. Furthermore, CRISPR/Cas9 editing of mice has significant scientific advantages compared to generating and breeding genetically modified mice such as the ability to evaluate targets that are embryonic lethal. Using CD40 as a model target in the CD40 agonistic antibody-induced colitis model, this study demonstrates the feasibility of this approach.

Genomic Rearrangements and Sequence Evolution across Brown Algal Organelles.

Organellar genomes serve as useful models for genome evolution and contain some of the most widely used phylogenetic markers, but they are poorly characterized in many lineages. Here, we report 20 novel mitochondrial genomes and 16 novel plastid genomes from the brown algae. We focused our efforts on the orders Chordales and Laminariales but also provide the first plastid genomes (plastomes) from Desmarestiales and Sphacelariales, the first mitochondrial genome (mitome) from Ralfsiales and a nearly complete mitome from Sphacelariales. We then compared gene content, sequence evolution rates, shifts in genome structural arrangements, and intron distributions across lineages. We confirm that gene content is largely conserved in both organellar genomes across the brown algal tree of life, with few cases of gene gain or loss. We further show that substitution rates are generally lower in plastid than mitochondrial genes, but plastomes are more variable in gene arrangement, as mitomes tend to be colinear even among distantly related lineages (with exceptions). Patterns of intron distribution across organellar genomes are complex. In particular, the mitomes of several laminarialean species possess group II introns that have T7-like ORFs, found previously only in mitochondrial genomes of Pylaiella spp. (Ectocarpales). The distribution of these mitochondrial introns is inconsistent with vertical transmission and likely reflects invasion by horizontal gene transfer between lineages. In the most extreme case, the mitome of Hedophyllum nigripes is ∼40% larger than the mitomes of close relatives because of these introns. Our results provide substantial insight into organellar evolution across the brown algae.

Discordant Phylogenomic Placement of Hydnoraceae and Lactoridaceae Within Piperales Using Data From All Three Genomes.

Phylogenetic relationships within the magnoliid order Piperales have been studied extensively, yet the relationships of the monotypic family Lactoridaceae and the holoparasitic Hydnoraceae to the remainder of the order remain a matter of debate. Since the first confident molecular phylogenetic placement of Hydnoraceae among Piperales, different studies have recovered various contradictory topologies. Most phylogenetic hypotheses were inferred using only a few loci and have had incomplete taxon sampling at the genus level. Based on these results and an online survey of taxonomic opinion, the Angiosperm Phylogeny Group lumped both Hydnoraceae and Lactoridaceae in Aristolochiaceae; however, the latter family continues to have unclear relationships to the aforementioned taxa. Here we present extensive phylogenomic tree reconstructions based on up to 137 loci from all three subcellular genomes for all genera of Piperales. We infer relationships based on a variety of phylogenetic methods, explore instances of phylogenomic discordance between the subcellular genomes, and test alternative topologies. Consistent with these phylogenomic results and a consideration of the principles of phylogenetic classification, we propose to exclude Hydnoraceae and Lactoridaceae from the broad circumscription of Aristolochiaceae, and instead favor recognition of four monophyletic and morphologically well circumscribed families in the perianth-bearing Piperales: Aristolochiaceae, Asaraceae, Hydnoraceae, and Lactoridaceae, with a total of six families in the order.