Investigation of the effects of phthalates on in vitro thyroid models with RNA-Seq and ATAC-Seq.

Introduction: Phthalates are a class of endocrine-disrupting chemicals that have been shown to negatively correlate with thyroid hormone serum levels in humans and to cause a state of hyperactivity in the thyroid. However, their mechanism of action is not well described at the molecular level. Methods: We analyzed the response of mouse thyroid organoids to the exposure to a biologically relevant dose range of the phthalates bis(2-ethylhexyl) phthalate (DEHP), di-iso-decylphthalate (DIDP), di-iso-nonylphthalate (DINP), and di-n-octylphthalate (DnOP) for 24 h and simultaneously analyzed mRNA and miRNA expression via RNA sequencing. Using the expression data, we performed differential expression and gene set enrichment analysis. We also exposed the human thyroid follicular epithelial cell line Nthy-ori 3-1 to 1 µM of DEHP or DINP for 5 days and analyzed changes in chromatin accessibility via ATAC-Seq. Results: Dose-series analysis showed how the expression of several genes increased or decreased at the highest dose tested. As expected with the low dosing scheme, the compounds induced a modest response on the transcriptome, as we identified changes in only mmu-miR-143-3p after DINP treatment and very few differentially expressed genes. No effect was observed on thyroid markers. Ing5, a component of histones H3 and H4 acetylation complexes, was consistently upregulated in three out of four conditions compared to control, and we observed a partial overlap among the genes differentially expressed by the treatments. Gene set enrichment analysis showed enrichment in the treatment samples of the fatty acid metabolism pathway and in the control of pathways related to the receptor signalling and extracellular matrix organization. ATAC-Seq analysis showed a general increase in accessibility compared to the control, however we did not identify significant changes in accessibility in the identified regions. Discussion: With this work, we showed that despite having only a few differentially expressed genes, diverse analysis methods could be applied to retrieve relevant information on phthalates, showing the potential of in vitro thyroid-relevant systems for the analysis of endocrine disruptors.

A novel orthotopic mouse model replicates human lung cancer cachexia.

INTRODUCTION: Cancer cachexia, highly prevalent in lung cancer, is a debilitating syndrome characterized by involuntary loss of skeletal muscle mass and is associated with poor clinical outcome, decreased survival and negative impact on tumour therapy. Various lung tumour-bearing animal models have been used to explore underlying mechanisms of cancer cachexia. However, these models do not simulate anatomical and immunological features key to lung cancer and associated muscle wasting. Overcoming these shortcomings is essential to translate experimental findings into the clinic. We therefore evaluated whether a syngeneic, orthotopic lung cancer mouse model replicates systemic and muscle-specific alterations associated with human lung cancer cachexia. METHODS: Immune competent, 11 weeks old male 129S2/Sv mice, were randomly allocated to either (1) sham control group or (2) tumour-bearing group. Syngeneic lung epithelium-derived adenocarcinoma cells (K-rasG12D ; p53R172HΔG ) were inoculated intrapulmonary into the left lung lobe of the mice. Body weight and food intake were measured daily. At baseline and weekly after surgery, grip strength was measured and tumour growth and muscle volume were assessed using micro cone beam CT imaging. After reaching predefined surrogate survival endpoint, animals were euthanized, and skeletal muscles of the lower hind limbs were collected for biochemical analysis. RESULTS: Two-third of the tumour-bearing mice developed cachexia based on predefined criteria. Final body weight (-13.7 ± 5.7%; P < 0.01), muscle mass (-13.8 ± 8.1%; P < 0.01) and muscle strength (-25.5 ± 10.5%; P < 0.001) were reduced in cachectic mice compared with sham controls and median survival time post-surgery was 33.5 days until humane endpoint. Markers for proteolysis, both ubiquitin proteasome system (Fbxo32 and Trim63) and autophagy-lysosomal pathway (Gabarapl1 and Bnip3), were significantly upregulated, whereas markers for protein synthesis (relative phosphorylation of Akt, S6 and 4E-BP1) were significantly decreased in the skeletal muscle of cachectic mice compared with control. The cachectic mice exhibited increased pentraxin-2 (P < 0.001) and CXCL1/KC (P < 0.01) expression levels in blood plasma and increased mRNA expression of IκBα (P < 0.05) in skeletal muscle, indicative for the presence of systemic inflammation. Strikingly, RNA sequencing, pathway enrichment and miRNA expression analyses of mouse skeletal muscle strongly mirrored alterations observed in muscle biopsies of patients with lung cancer cachexia. CONCLUSIONS: We developed an orthotopic model of lung cancer cachexia in immune competent mice. Because this model simulates key aspects specific to cachexia in lung cancer patients, it is highly suitable to further investigate the underlying mechanisms of lung cancer cachexia and to test the efficacy of novel intervention strategies.

CODA: a combo-Seq data analysis workflow.

The analysis of the combined mRNA and miRNA content of a biological sample can be of interest for answering several research questions, like biomarkers discovery, or mRNA-miRNA interactions. However, the process is costly and time-consuming, separate libraries need to be prepared and sequenced on different flowcells. Combo-Seq is a library prep kit that allows us to prepare combined mRNA-miRNA libraries starting from very low total RNA. To date, no dedicated bioinformatics method exists for the processing of Combo-Seq data. In this paper, we describe CODA (Combo-seq Data Analysis), a workflow specifically developed for the processing of Combo-Seq data that employs existing free-to-use tools. We compare CODA with exceRpt, the pipeline suggested by the kit manufacturer for this purpose. We also evaluate how Combo-Seq libraries analysed with CODA perform compared with conventional poly(A) and small RNA libraries prepared from the same samples. We show that using CODA more successfully trimmed reads are recovered compared with exceRpt, and the difference is more dramatic with short sequencing reads. We demonstrate how Combo-Seq identifies as many genes and fewer miRNAs compared to the standard libraries, and how miRNA validation favours conventional small RNA libraries over Combo-Seq. The CODA code is available at https://github.com/marta-nazzari/CODA.

Lead-exposure associated miRNAs in humans and Alzheimer's disease as potential biomarkers of the disease and disease processes.

Alzheimer's disease (AD) is a neurodegenerative disease that eventually affects memory and behavior. The identification of biomarkers based on risk factors for AD provides insight into the disease since the exact cause of AD remains unknown. Several studies have proposed microRNAs (miRNAs) in blood as potential biomarkers for AD. Exposure to heavy metals is a potential risk factor for onset and development of AD. Blood cells of subjects that are exposed to lead detected in the circulatory system, potentially reflect molecular responses to this exposure that are similar to the response of neurons. In this study we analyzed blood cell-derived miRNAs derived from a general population as proxies of potentially AD-related mechanisms triggered by lead exposure. Subsequently, we analyzed these mechanisms in the brain tissue of AD subjects and controls. A total of four miRNAs were identified as lead exposure-associated with hsa-miR-3651, hsa-miR-150-5p and hsa-miR-664b-3p being negatively and hsa-miR-627 positively associated. In human brain derived from AD and AD control subjects all four miRNAs were detected. Moreover, two miRNAs (miR-3651, miR-664b-3p) showed significant differential expression in AD brains versus controls, in accordance with the change direction of lead exposure. The miRNAs' gene targets were validated for expression in the human brain and were found enriched in AD-relevant pathways such as axon guidance. Moreover, we identified several AD relevant transcription factors such as CREB1 associated with the identified miRNAs. These findings suggest that the identified miRNAs are involved in the development of AD and might be useful in the development of new, less invasive biomarkers for monitoring of novel therapies or of processes involved in AD development.

Dynamic plastid and mitochondrial genomes in Chaetopeltidales (Chlorophyceae) and characterization of a new chlorophyte taxon.

PREMISE: Chaetopeltidales is a poorly characterized order in the Chlorophyceae, with only two plastid and no mitochondrial genomes published. Here we describe a new taxon in Chaetopeltidales, Gormaniella terricola gen. et sp. nov. and characterize both of its organellar genomes. METHODS: Gormaniella terricola was inadvertently isolated from a surface-sterilized hornwort thallus. Light microscopy was used to characterize its vegetative morphology. Organellar genomes were assembled, annotated, and analyzed using a variety of software packages. RESULTS: The mitochondrial genome (66,927 bp) represents the first complete mitochondrial genome published for Chaetopeltidales. The chloroplast genome, measuring 428,981 bp, is one of the largest plastid genomes published to date and shares this large size and an incredible number of short, dispersed repeats with the other sequenced chloroplast genomes in Chaetopeltidales. Despite these shared features, the chloroplast genomes of Chaetopeltidales appear to be highly rearranged when compared to one another, with numerous inversions, translocations, and duplications, suggesting a particularly dynamic chloroplast genome. Both the chloroplast and mitochondrial genomes of G. terricola contain a number of mobile group I and group II introns, which appear to have invaded separately. Three of the introns within the mitochondrial genome encode homing endonucleases that are phylogenetically nested within those found in fungi, rather than algae, suggesting a possible case of horizontal gene transfer. CONCLUSIONS: These results help to shed light on a poorly understood group of algae and their unusual organellar genomes, raising additional questions about the unique patterns of genome evolution within Chaetopeltidales.

Genome sequence for the blue-flowered Andean shrub Iochroma cyaneum reveals extensive discordance across the berry clade of Solanaceae.

The tomato (Solanum lycopersicum L.) family, Solanaceae, is a model clade for a wide range of applied and basic research questions. Currently, reference-quality genomes are available for over 30 species from seven genera, and these include numerous crops as well as wild species [e.g., Jaltomata sinuosa (Miers) Mione and Nicotiana attenuata Torr. ex S. Watson]. Here we present the genome of the showy-flowered Andean shrub Iochroma cyaneum (Lindl.) M. L. Green, a woody lineage from the tomatillo (Physalis philadelphica Lam.) subfamily Physalideae. The assembled size of the genome (2.7 Gb) is more similar in size to pepper (Capsicum annuum L.) (2.6 Gb) than to other sequenced diploid members of the berry clade of Solanaceae [e.g., potato (Solanum tuberosum L.), tomato, and Jaltomata]. Our assembly recovers 92% of the conserved orthologous set, suggesting a nearly complete genome for this species. Most of the genomic content is repetitive (69%), with Gypsy elements alone accounting for 52% of the genome. Despite the large amount of repetitive content, most of the 12 I. cyaneum chromosomes are highly syntenic with tomato. Bayesian concordance analysis provides strong support for the berry clade, including I. cyaneum, but reveals extensive discordance along the backbone, with placement of chili pepper and Jaltomata being highly variable across gene trees. The I. cyaneum genome contributes to a growing wealth of genomic resources in Solanaceae and underscores the need for expanded sampling of diverse berry genomes to dissect major morphological transitions.

Transcriptomics analysis of human iPSC-derived dopaminergic neurons reveals a novel model for sporadic Parkinson's disease.

Parkinson's disease (PD) is a progressive, neurodegenerative disease affecting over 1% of the population beyond 65 years of age. Although some PD cases are inheritable, the majority of PD cases occur in a sporadic manner. Risk factors comprise next to heredity, age, and gender also exposure to neurotoxins from for instance pesticides and herbicides. As PD is characterized by a loss of dopaminergic neurons in the substantia nigra, it is nearly impossible to access and extract these cells from patients for investigating disease mechanisms. The emergence of induced pluripotent stem (iPSC) technology allows differentiating and growing human dopaminergic neurons, which can be used for in vitro disease modeling. Here, we differentiated human iPSCs into dopaminergic neurons, and subsequently exposed the cells to increasing concentrations of the neurotoxin MPP+. Temporal transcriptomics analysis revealed a strong time- and dose-dependent response with genes over-represented across pathways involved in PD etiology such as "Parkinson's Disease", "Dopaminergic signaling" and "calcium signaling". Moreover, we validated this disease model by showing robust overlap with a meta-analysis of transcriptomics data from substantia nigra from post-mortem PD patients. The overlap included genes linked to e.g. mitochondrial dysfunction, neuron differentiation, apoptosis and inflammation. Our data shows, that MPP+-induced, human iPSC-derived dopaminergic neurons present molecular perturbations as observed in the etiology of PD. Therefore we propose iPSC-derived dopaminergic neurons as a foundation for a novel sporadic PD model to study the pathomolecular mechanisms of PD, but also to screen for novel anti-PD drugs and to develop and test new treatment strategies.

Re-evaluating the Systematics of Dendrolycopodium Using Restriction-Site Associated DNA-Sequencing.

The genus Dendrolycopodium (Lycopodiaceae) includes four to five species across North America and East Asia. Species identification in Dendrolycopodium is difficult due to limited or inconsistent characters. In addition, plants with intermediate morphologies regularly occur, potentially indicative of interspecific hybridization. To determine the species relationships in Dendrolycopodium and investigate the existence of hybrids, we generated a draft genome assembly for D. obscurum and carried out double-digest restriction-site associated DNA sequencing (RADSeq) on 86 Dendrolycopodium specimens. Our sampling includes all the described species and 11 individuals with intermediate morphology. We find that the genus can be divided into four clades that largely correspond to the described taxa, as well as evidence of interspecific hybridization. Within these clades, our STRUCTURE analysis suggests that there are multiple finer subgroups, with evidence of hybridization and introgression between these subgroups. Given the limited availability of specimens collected from Asia, the status of the various Asian species remains uncertain and will require further study. In summary, our study confirms several hybrid relationships in Dendrolycopodium and provides a clear phylogenetic framework for future taxonomic revision.

Monodopsis and Vischeria Genomes Shed New Light on the Biology of Eustigmatophyte Algae.

Members of eustigmatophyte algae, especially Nannochloropsis and Microchloropsis, have been tapped for biofuel production owing to their exceptionally high lipid content. Although extensive genomic, transcriptomic, and synthetic biology toolkits have been made available for Nannochloropsis and Microchloropsis, very little is known about other eustigmatophytes. Here we present three near-chromosomal and gapless genome assemblies of Monodopsis strains C73 and C141 (60 Mb) and Vischeria strain C74 (106 Mb), which are the sister groups to Nannochloropsis and Microchloropsis in the order Eustigmatales. These genomes contain unusually high percentages of simple repeats, ranging from 12% to 21% of the total assembly size. Unlike Nannochloropsis and Microchloropsis, long interspersed nuclear element repeats are abundant in Monodopsis and Vischeria and might constitute the centromeric regions. We found that both mevalonate and nonmevalonate pathways for terpenoid biosynthesis are present in Monodopsis and Vischeria, which is different from Nannochloropsis and Microchloropsis that have only the latter. Our analysis further revealed extensive spliced leader trans-splicing in Monodopsis and Vischeria at 36-61% of genes. Altogether, the high-quality genomes of Monodopsis and Vischeria not only serve as the much-needed outgroups to advance Nannochloropsis and Microchloropsis research, but also shed new light on the biology and evolution of eustigmatophyte algae.

The diversity and community structure of symbiotic cyanobacteria in hornworts inferred from long-read amplicon sequencing.

PREMISE: Nitrogen-fixing endosymbioses with cyanobacteria have evolved independently in five very different plant lineages. Expanding knowledge of these symbioses promises to improve the understanding of symbiosis evolution and broaden the toolkit for agricultural engineering to reduce artificial fertilizer use. Here we focused on hornworts, a bryophyte lineage in which all members host cyanobacteria, and investigated factors shaping the diversity of their cyanobiont communities. METHODS: We sampled hornworts and adjacent soils in upstate New York throughout the hornwort growing season. We included all three sympatric hornwort species in the area, allowing us to directly compare partner selectivity. To profile cyanobacteria communities, we established a metabarcoding protocol targeting rbcL-X with PacBio long reads. RESULTS: The hornwort cyanobionts detected were phylogenetically diverse, including clades that do not contain other known plant symbionts. We found significant overlap between hornwort cyanobionts and soil cyanobacteria, a pattern not previously reported in other plant-cyanobacteria symbioses. Cyanobiont communities differed between host plants only centimeters apart, but we did not detect an effect of sampling time or host species on the cyanobacterial community structure. CONCLUSIONS: This study expands the phylogenetic diversity of known symbiotic cyanobacteria. Our analyses suggest that hornwort cyanobionts have a tight connection to the soil background, and we found no evidence that time within growing season, host species, or distance at the scale of meters strongly govern cyanobacteria community assembly. This study provides a critical foundation for further study of the ecology, evolution, and interaction dynamics of plant-cyanobacteria symbiosis.

A novel thylakoid-less isolate fills a billion-year gap in the evolution of Cyanobacteria.

Cyanobacteria have played pivotal roles in Earth's geological history, especially during the rise of atmospheric oxygen. However, our ability to infer the early transitions in Cyanobacteria evolution has been limited by their extremely lopsided tree of life-the vast majority of extant diversity belongs to Phycobacteria (or "crown Cyanobacteria"), while its sister lineage, Gloeobacteria, is depauperate and contains only two closely related species of Gloeobacter and a metagenome-assembled genome. Here, we describe a new cultured member of Gloeobacteria, Anthocerotibacter panamensis, isolated from a tropical hornwort. Anthocerotibacter diverged from Gloeobacter over 1.4 Ga ago and has low 16S rDNA identities with environmental samples. Our ultrastructural, physiological, and genomic analyses revealed that this species possesses a unique combination of traits that are exclusively shared with either Gloeobacteria or Phycobacteria. For example, similar to Gloeobacter, it lacks thylakoids and circadian clock genes, but the carotenoid biosynthesis pathway is typical of Phycobacteria. Furthermore, Anthocerotibacter has one of the most reduced gene sets for photosystems and phycobilisomes among Cyanobacteria. Despite this, Anthocerotibacter is capable of oxygenic photosynthesis under a wide range of light intensities, albeit with much less efficiency. Given its key phylogenetic position, distinct trait combination, and availability as a culture, Anthocerotibacter opens a new window to further illuminate the dawn of oxygenic photosynthesis.

Comparing in vitro human liver models to in vivo human liver using RNA-Seq.

The liver plays an important role in xenobiotic metabolism and represents a primary target for toxic substances. Many different in vitro cell models have been developed in the past decades. In this study, we used RNA-sequencing (RNA-Seq) to analyze the following human in vitro liver cell models in comparison to human liver tissue: cancer-derived cell lines (HepG2, HepaRG 3D), induced pluripotent stem cell-derived hepatocyte-like cells (iPSC-HLCs), cancerous human liver-derived assays (hPCLiS, human precision cut liver slices), non-cancerous human liver-derived assays (PHH, primary human hepatocytes) and 3D liver microtissues. First, using CellNet, we analyzed whether these liver in vitro cell models were indeed classified as liver, based on their baseline expression profile and gene regulatory networks (GRN). More comprehensive analyses using non-differentially expressed genes (non-DEGs) and differential transcript usage (DTU) were applied to assess the coverage for important liver pathways. Through different analyses, we noticed that 3D liver microtissues exhibited a high similarity with in vivo liver, in terms of CellNet (C/T score: 0.98), non-DEGs (10,363) and pathway coverage (highest for 19 out of 20 liver specific pathways shown) at the beginning of the incubation period (0 h) followed by a decrease during long-term incubation for 168 and 336 h. PHH also showed a high degree of similarity with human liver tissue and allowed stable conditions for a short-term cultivation period of 24 h. Using the same metrics, HepG2 cells illustrated the lowest similarity (C/T: 0.51, non-DEGs: 5623, and pathways coverage: least for 7 out of 20) with human liver tissue. The HepG2 are widely used in hepatotoxicity studies, however, due to their lower similarity, they should be used with caution. HepaRG models, iPSC-HLCs, and hPCLiS ranged clearly behind microtissues and PHH but showed higher similarity to human liver tissue than HepG2 cells. In conclusion, this study offers a resource of RNA-Seq data of several biological replicates of human liver cell models in vitro compared to human liver tissue.

Anthoceros genomes illuminate the origin of land plants and the unique biology of hornworts.

Hornworts comprise a bryophyte lineage that diverged from other extant land plants >400 million years ago and bears unique biological features, including a distinct sporophyte architecture, cyanobacterial symbiosis and a pyrenoid-based carbon-concentrating mechanism (CCM). Here, we provide three high-quality genomes of Anthoceros hornworts. Phylogenomic analyses place hornworts as a sister clade to liverworts plus mosses with high support. The Anthoceros genomes lack repeat-dense centromeres as well as whole-genome duplication, and contain a limited transcription factor repertoire. Several genes involved in angiosperm meristem and stomatal function are conserved in Anthoceros and upregulated during sporophyte development, suggesting possible homologies at the genetic level. We identified candidate genes involved in cyanobacterial symbiosis and found that LCIB, a Chlamydomonas CCM gene, is present in hornworts but absent in other plant lineages, implying a possible conserved role in CCM function. We anticipate that these hornwort genomes will serve as essential references for future hornwort research and comparative studies across land plants.

Complete Genomes of Symbiotic Cyanobacteria Clarify the Evolution of Vanadium-Nitrogenase.

Plant endosymbiosis with nitrogen-fixing cyanobacteria has independently evolved in diverse plant lineages, offering a unique window to study the evolution and genetics of plant-microbe interaction. However, very few complete genomes exist for plant cyanobionts, and therefore little is known about their genomic and functional diversity. Here, we present four complete genomes of cyanobacteria isolated from bryophytes. Nanopore long-read sequencing allowed us to obtain circular contigs for all the main chromosomes and most of the plasmids. We found that despite having a low 16S rRNA sequence divergence, the four isolates exhibit considerable genome reorganizations and variation in gene content. Furthermore, three of the four isolates possess genes encoding vanadium (V)-nitrogenase (vnf), which is uncommon among diazotrophs and has not been previously reported in plant cyanobionts. In two cases, the vnf genes were found on plasmids, implying possible plasmid-mediated horizontal gene transfers. Comparative genomic analysis of vnf-containing cyanobacteria further identified a conserved gene cluster. Many genes in this cluster have not been functionally characterized and would be promising candidates for future studies to elucidate V-nitrogenase function and regulation.

A novel experimental system using the liverwort Marchantia polymorpha and its fungal endophytes reveals diverse and context-dependent effects.

Fungal symbioses are ubiquitous in plants, but their effects have mostly been studied in seed plants. This study aimed to assess the diversity of fungal endophyte effects in a bryophyte and identify factors contributing to the variability of outcomes in these interactions. Fungal endophyte cultures and axenic liverwort clones were isolated from wild populations of the liverwort, Marchantia polymorpha. These collections were combined in a gnotobiotic system to test the effects of fungal isolates on the growth rates of hosts under laboratory conditions. Under the experimental conditions, fungi isolated from M. polymorpha ranged from aggressively pathogenic to strongly growth-promoting, but the majority of isolates caused no detectable change in host growth. Growth promotion by selected fungi depended on nutrient concentrations and was inhibited by coinoculation with multiple fungi. The M. polymorpha endophyte system expands the resources for this model liverwort. The experiments presented here demonstrate a wealth of diversity in fungal interactions even in a host reported to lack standard mycorrhizal symbiosis. In addition, they show that some known pathogens of vascular plants live in M. polymorpha and can confer benefits to this nonvascular host. This highlights the importance of studying endophyte effects across the plant tree of life.

Phylogenomics reveals a complex evolutionary history of lobed-leaf white oaks in western North America.

Species within the genus Quercus (oak) hybridize in complex patterns that have yet to be fully explored with phylogenomic data. Analyses to date have recovered reasonable divergent patterns, suggesting that the impact of introgression may not always be obvious in inferred oak phylogenies. We explore this phenomenon using RADseq data for 136 samples representing 54 oak species by conducting phylogenetic analyses designed to distinguish signals of lineage diversification and hybridization, focusing on the lobed-leaf species Quercus gambelii, Q. lobata, and Q. garryana in the context of a broad sampling of allied white oaks (Quercus section Quercus), and particularly the midwestern Q. macrocarpa. We demonstrate that historical introgressive hybridization between once sympatric species affects phylogeny estimation. Historical range expansion during periods of favorable climate likely explains our observations; analyses support genetic exchange between ancestral populations of Q. gambelii and Q. macrocarpa. We conclude that the genomic consequences of introgression caused the attraction of distant lineages in phylogenetic tree space, and that introgressive and divergent signals can be disentangled to produce a robust estimate of the phylogenetic history of the species.

The evolution and diversification of the red oaks of the California Floristic Province (Quercus section Lobatae, series Agrifoliae).

PREMISE OF THE STUDY: The California Floristic Province (CA-FP) is a unique and diverse region of floral endemism, yet the timing and nature of divergence and diversification of many lineages remain underexplored. We seek to elucidate the evolutionary history of the red oaks of the CA-FP, the Agrifoliae. METHODS: We collected PstI-associated RAD-seq data as well as morphometrics from individuals of the four species across their ranges, including varieties and hybrids. Phylogeny and divergence times were estimated. We analyzed morphological differentiation in over 70 plants using PCA and assessed species delimitation and admixture using genotype clustering analysis in over 40 plants. KEY RESULTS: We find that the Agrifoliae are monophyletic and sister to all other red oak species. Within the Agrifoliae, all species are supported, with Quercus kelloggii sister to a clade of subevergreen taxa: (Quercus agrifolia - (Q. parvula + Q. wislizeni)). Molecular and morphometric analyses are equivocal for named varieties. Notably, Q. parvula var. tamalpaisensis appears to be part of a hybrid swarm between Q. parvula and Q. wislizeni. Dating estimates were concordant with previous hypotheses and geological evidence, with diversification occurring between 10 and 20 million years ago. CONCLUSIONS: The Agrifoliae represent a geographically discrete, early-diverging red oak lineage that diversified during the period of drying and warming associated with Sierran uplift during the middle Miocene. Molecular differentiation within the clade supports the current taxonomy, including an east-west species level pattern (Q. parvula and Q. wislizeni) and north-south intraspecific patterns to some degree, although the latter require additional study.