Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae).

BACKGROUND AND AIMS: Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico. METHODS: Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages. KEY RESULTS: The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages. CONCLUSIONS: The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

Genomic landscape of the global oak phylogeny.

The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny. We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 individuals representing nearly 250 Quercus species to infer a time-calibrated phylogeny of the world's oaks. We use a reversible-jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among-clade sampling biases. We then map the > 20 000 RAD-seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny. Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity. The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.

Sympatric parallel diversification of major oak clades in the Americas and the origins of Mexican species diversity.

Oaks (Quercus, Fagaceae) are the dominant tree genus of North America in species number and biomass, and Mexico is a global center of oak diversity. Understanding the origins of oak diversity is key to understanding biodiversity of northern temperate forests. A phylogenetic study of biogeography, niche evolution and diversification patterns in Quercus was performed using 300 samples, 146 species. Next-generation sequencing data were generated using the restriction-site associated DNA (RAD-seq) method. A time-calibrated maximum likelihood phylogeny was inferred and analyzed with bioclimatic, soils, and leaf habit data to reconstruct the biogeographic and evolutionary history of the American oaks. Our highly resolved phylogeny demonstrates sympatric parallel diversification in climatic niche, leaf habit, and diversification rates. The two major American oak clades arose in what is now the boreal zone and radiated, in parallel, from eastern North America into Mexico and Central America. Oaks adapted rapidly to niche transitions. The Mexican oaks are particularly numerous, not because Mexico is a center of origin, but because of high rates of lineage diversification associated with high rates of evolution along moisture gradients and between the evergreen and deciduous leaf habits. Sympatric parallel diversification in the oaks has shaped the diversity of North American forests.

Phylogeny and biogeography of East Asian evergreen oaks (Quercus section Cyclobalanopsis; Fagaceae): Insights into the Cenozoic history of evergreen broad-leaved forests in subtropical Asia.

The evolutionary history of Quercus section Cyclobalanopsis, a dominant lineage in East Asian evergreen broadleaved forests (EBLFs), has not been comprehensively studied using molecular tools. In this study, we reconstruct the first comprehensive phylogeny of this lineage using a genomic approach (restriction-site associated DNA sequencing, RAD-seq), sampling 35 of the ca. 90 species currently recognized, representing all main morphological groups of section Cyclobalanopsis. In addition, 10 other species of Quercus and two outgroups were also sampled. Divergence times were estimated using a relaxed clock model and two fossil calibrations. Ancestral areas and dispersal routes were inferred using statistical dispersal-vicariance analysis and the dispersal-extinction-cladogenesis (DEC) model. The phylogeny of Quercus section Cyclobalanopsis demonstrates the section to be monophyletic, comprising two main lineages and six subclades that are well supported by anatomical traits. Biogeographical reconstructions indicate that the wide northern hemisphere distribution of Quercus was disrupted in the Late Eocene, leading to the main extant groups at about 33 Ma. The earliest divergences in section Cyclobalanopsis correspond to the phased uplift of the Himalayas and lateral extrusion of Indochina at the transition of the Oligocene and Miocene, where the highest rate of diversification occurred in the late Miocene. Dispersal from Sino-Himalaya and the Palaeotropics to Sino-Japan in the Miocene was facilitated by the increased intensity of East Asian summer monsoons and by the Middle Miocene Climatic Optimum. Our results highlight the importance of climatic changes and Indo-Eurasian collision-induced tectonic activities from the Neogene onward to the spatial-temporal diversification patterns of Asian EBLF lineages.

RAD-seq linkage mapping and patterns of segregation distortion in sedges: meiosis as a driver of karyotypic evolution in organisms with holocentric chromosomes.

Meiotic drive, the class of meiotic mechanisms that drive unequal segregation of alleles among gametes, may be an important force in karyotype evolution. Its role in holocentric organisms, whose chromosomes lack localized centromeres, is poorly understood. We crossed two individuals of Carex scoparia (Cyperaceae) with different chromosome numbers (2n = 33II  = 66 × 2n = 32II  = 64) to obtain F1 individuals, which we then self-pollinated to obtain second-generation (F2) crosses. RAD-seq was performed for 191 individuals (including the parents, five F1 individuals and 184 F2 individuals). Our F2 linkage map based on stringent editing of the RAD-seq data set yielded 32 linkage groups. In the final map, 865 loci were located on a linkage map of 3966.99 cM (linkage groups ranged from 24.39 to 193.31 cM in length and contained 5-51 loci each). Three linkage groups exhibit more loci under segregation distortion than expected by chance; within linkage groups, loci exhibiting segregation distortion are clustered. This finding implicates meiotic drive in the segregation of chromosome variants, suggesting that selection of chromosome variants in meiosis may contribute to the establishment and fixation of chromosome variants in Carex, which is renowned for high chromosomal and species diversity. This is an important finding as previous studies demonstrate that chromosome divergence may play a key role in differentiation and speciation in Carex.

Chromosomal rearrangements in holocentric organisms lead to reproductive isolation by hybrid dysfunction: The correlation between karyotype rearrangements and germination rates in sedges.

PREMISE OF THE STUDY: Understanding the drivers of speciation is a central task of evolutionary biology. Chromosomal rearrangements are known to play an important role in species diversification, but the role of rearrangements of holocentric chromosomes-chromosomes without localized centromeres-is poorly understood. METHODS: We made numerous artificial crosses between Carex scoparia individuals of different diploid chromosome numbers and, for comparison, between individuals of the same chromosome number. We studied chromosome pairing and chromosomal rearrangements in the F1 individuals using light microscopy. We then estimated germination rates as a function of geographic distance, genetic distance, chromosome number differences in parents, and pairing irregularities in F1 individuals, using generalized least squares to fit alternative regression models. KEY RESULTS: The most informative predictors of germination rates in the F1 generation are chromosome number differences and minimum number of chromosome pairing irregularities in the F1 individuals. Genetic and geographic distances between parents are not significant predictors. CONCLUSIONS: Holocentric chromosomal rearrangements play an important role in postzygotic reproductive isolation in Carex through F1 hybrid inviability and sterility. Hybrid dysfunction seems to be a suitable model for chromosomal speciation when there are several chromosomal rearrangements between parents. However, we have not tested the hypothesis that genome rearrangements may also play an important role in suppressing recombination between cytogenetically divergent populations.

Genotyping-by-sequencing as a tool to infer phylogeny and ancestral hybridization: a case study in Carex (Cyperaceae).

Determining phylogenetic relationships among very closely related species has remained a challenge for evolutionary biologists due to interlocus phylogenetic discordance and the difficulty of obtaining variable markers. Here, we used a Genotyping-by-Sequencing (GBS) approach to sample a reduced representation genomic data set and infer the phylogeny of seven closely related species in the genus Carex (Cyperaceae). Past attempts to reconstruct phylogenetic relationships among these species produced conflicting and poorly-supported results. We inferred a robust phylogeny based on >3000 GBS loci and >1300 SNPs (with a minimum sequence depth within individuals of 10) using maximum likelihood and Bayesian inference. We also tested for historical introgression using the D-statistic test. We compared these analyses with partitioned RAD analysis, which is designed to identify suboptimal trees reflecting secondary phylogenetic signal that may be obscured by the dominant signal in the data. Phylogenetic analyses yielded fully resolved trees with high support. We found two main clades, one grouping Carex scoparia populations and C. waponahkikensis, and a second clade grouping C. longii, C. vexans, C. suberecta and C. albolutescens. We detected marginally significant signals of introgression between C. scoparia and C. suberecta or C. albolutescens, and we rejected a hybrid origin hypothesis for C. waponahkikensis. Our results demonstrate the power of NGS data sets for resolving some of the most difficult phylogenetic challenges where traditional phylogenetic markers have failed.

COVID-19 and Perinatal Stress Experience - a Study Conducted as Part of the COVGEN Initiative.

Introduction During the COVID-19 pandemic, stress and anxiety in the population increased due to concerns about people's own health and that of their relatives, as well as changes in everyday life due to measures taken to reduce the infection rate. Pregnant women are particularly stressed. The present study examines how the COVID-19 pandemic affects the stress experience and mental health of pregnant women and mothers of newborns and how care could be optimized. Methods As part of the international COVGEN initiative ( https://www.covgen.org ) to investigate the effects of the COVID-19 pandemic on the peripartum period, pregnant and postpartum women were asked about their experience with stress using the COPE-IS (Coronavirus Perinatal Experiences - Impact Survey) questionnaire developed for this purpose and translated from the English. In addition, demographic data, pre-existing diseases, pregnancy complications and the care situation were recorded. The questionnaire was either administered as hardcopy to inpatients at the Department of Women's Health, University Hospital Tübingen, Germany, or online. All pregnant women and mothers who were pregnant or had given birth after the official start of the COVID-19 pandemic (11 March 2020) were eligible to participate. Results Complete data sets of n = 156 pregnant women and n = 221 postpartum women were available for evaluation. The general stress level assessed with the COPE-IS was significantly increased by the COVID-19 pandemic in both, pregnant and postpartum women, with pre-existing conditions such as respiratory diseases and pregnancy-related diseases like gestational diabetes adding to the stress. The subjectively perceived quality of care/support during pregnancy also influenced the stress level. Conclusions Fears of a COVID-19 infection and changes in preventive and aftercare services were a burden for the women surveyed. Intensified care during pregnancy and puerperium could help to stabilize the mental situation and reduce stress.

Tackling Rapid Radiations With Targeted Sequencing.

In phylogenetic studies across angiosperms, at various taxonomic levels, polytomies have persisted despite efforts to resolve them by increasing sampling of taxa and loci. The large amount of genomic data now available and statistical tools to analyze them provide unprecedented power for phylogenetic inference. Targeted sequencing has emerged as a strong tool for estimating species trees in the face of rapid radiations, lineage sorting, and introgression. Evolutionary relationships in Cyperaceae have been studied mostly using Sanger sequencing until recently. Despite ample taxon sampling, relationships in many genera remain poorly understood, hampered by diversification rates that outpace mutation rates in the loci used. The C4 Cyperus clade of the genus Cyperus has been particularly difficult to resolve. Previous studies based on a limited set of markers resolved relationships among Cyperus species using the C3 photosynthetic pathway, but not among C4 Cyperus clade taxa. We test the ability of two targeted sequencing kits to resolve relationships in the C4 Cyperus clade, the universal Angiosperms-353 kit and a Cyperaceae-specific kit. Sequences of the targeted loci were recovered from data generated with both kits and used to investigate overlap in data between kits and relative efficiency of the general and custom approaches. The power to resolve shallow-level relationships was tested using a summary species tree method and a concatenated maximum likelihood approach. High resolution and support are obtained using both approaches, but high levels of missing data disproportionately impact the latter. Targeted sequencing provides new insights into the evolution of morphology in the C4 Cyperus clade, demonstrating for example that the former segregate genus Alinula is polyphyletic despite its seeming morphological integrity. An unexpected result is that the Cyperus margaritaceus-Cyperus niveus complex comprises a clade separate from and sister to the core C4 Cyperus clade. Our results demonstrate that data generated with a family-specific kit do not necessarily have more power than those obtained with a universal kit, but that data generated with different targeted sequencing kits can often be merged for downstream analyses. Moreover, our study contributes to the growing consensus that targeted sequencing data are a powerful tool in resolving rapid radiations.

A nuclear DNA barcode for eastern North American oaks and application to a study of hybridization in an Arboretum setting.

DNA barcoding has proved difficult in a number of woody plant genera, including the ecologically important oak genus Quercus. In this study, we utilized restrictionsite-associated DNA sequencing (RAD-seq) to develop an economical single nucleotide polymorphism (SNP) DNA barcoding system that suffices to distinguish eight common, sympatric eastern North American white oak species. Two de novo clustering pipelines, PyRAD and Stacks, were used in combination with postclustering bioinformatic tools to generate a list of 291 potential SNPs, 80 of which were included in a barcoding toolkit that is easily implemented using MassARRAY mass spectrometry technology. As a proof-of-concept, we used the genotyping toolkit to infer potential hybridization between North American white oaks transplanted outside of their native range (Q. michauxii, Q. montana, Q muehlenbergii/Q. prinoides, and Q. stellata) into a horticultural collection surrounded by natural forests of locally native trees (Q. alba and Q. macrocarpa) in the living collection at The Morton Arboretum (Lisle, IL, USA). Phylogenetic and clustering analyses suggested low rates of hybridization between cultivated and native species, with the exception of one Q. michauxii mother tree, the acorns of which exhibited high admixture from either Q. alba or Q. stellata and Q. macrocarpa, and a hybrid between Q. stellata that appears to have backcrossed almost exclusively to Q. alba. Together, RAD-seq and MassARRAY technologies allow for efficient development and implementation of a multispecies barcode for one of the more challenging forest tree genera.

White-tailed deer are a biotic filter during community assembly, reducing species and phylogenetic diversity.

Community assembly entails a filtering process, where species found in a local community are those that can pass through environmental (abiotic) and biotic filters and successfully compete. Previous research has demonstrated the ability of white-tailed deer (Odocoileus virginianus) to reduce species diversity and favour browse-tolerant plant communities. In this study, we expand on our previous work by investigating deer as a possible biotic filter altering local plant community assembly. We used replicated 23-year-old deer exclosures to experimentally assess the effects of deer on species diversity (H'), richness (SR), phylogenetic community structure and phylogenetic diversity in paired browsed (control) and unbrowsed (exclosed) plots. Additionally, we developed a deer-browsing susceptibility index (DBSI) to assess the vulnerability of local species to deer. Deer browsing caused a 12 % reduction in H' and 17 % reduction in SR, consistent with previous studies. Furthermore, browsing reduced phylogenetic diversity by 63 %, causing significant phylogenetic clustering. Overall, graminoids were the least vulnerable to deer browsing based on DBSI calculations. These findings demonstrate that deer are a significant driver of plant community assembly due to their role as a selective browser, or more generally, as a biotic filter. This study highlights the importance of knowledge about the plant tree of life in assessing the effects of biotic filters on plant communities. Application of such knowledge has considerable potential to advance our understanding of plant community assembly.

The Bacillus subtilis spore coat protein interaction network.

Bacterial spores are surrounded by a morphologically complex, mechanically flexible protein coat, which protects the spore from toxic molecules. The interactions among the over 50 proteins that make up the coat remain poorly understood. We have used cell biological and protein biochemical approaches to identify novel coat proteins in Bacillus subtilis and describe the network of their interactions, in order to understand coat assembly and the molecular basis of its protective functions and mechanical properties. Our analysis characterizes the interactions between 32 coat proteins. This detailed view reveals a complex interaction network. A key feature of the network is the importance of a small subset of proteins that direct the assembly of most of the coat. From an analysis of the network topology, we propose a model in which low-affinity interactions are abundant in the coat and account, to a significant degree, for the coat's mechanical properties as well as structural variation between spores.