We applied an interdisciplinary approach to investigate kinship patterns and funerary practices during the middle Neolithic. Genetic studies, radiocarbon dating, and taphonomic analyses were used to examine two grave clusters from Krusza Zamkowa, Poland. To reconstruct kinship and determine biological sex, we extracted DNA from bones and teeth, analyzed mitochondrial genomes and nuclear SNPs using the HID-Ion AmpliSeq™ Identity panel generated on Illumina and Ion Torrent platforms, respectively. We further dated the material (AMS 14C) and to exclude aquatic radiocarbon reservoir effects, measures of carbon and nitrogen stable isotopes for diet reconstruction were used. We found distinct mitochondrial genomes belonging to haplogroups U5b2a1a, K1c and H3d in the first grave cluster, and excluded maternal kin patterns among the three analyzed individuals. In the second grave cluster one individual belonged to K1a4. However, we could not affiliate the second individual to a certain haplogroup due to the fragmented state of the mitochondrial genome. Although the individuals from the second grave cluster differ at position 6643, we believe that more data is needed to fully resolve this issue. We retrieved between 26 and 77 autosomal SNPs from three of the individuals. Based on kinship estimations, taking into account the allelic dropout distribution, we could not exclude first degree kin relation between the two individuals from the second grave cluster. We could, however, exclude a first degree kinship between these two individuals and an individual from the first grave cluster. Presumably, not only biological kinship, but also social relations played an important role in the funerary practice during this time period. We further conclude that the HID-Ion AmpliSeq™ Identity Panel may prove useful for first degree kin relation studies for samples with good DNA preservation, and that mitochondrial genome capture enrichment is a powerful tool for excluding direct maternal relationship in ancient individuals.

Hybrid enrichment is an increasingly popular approach for obtaining hundreds of loci for phylogenetic analysis across many taxa quickly and cheaply. The genes targeted for sequencing are typically single-copy loci, which facilitate a more straightforward sequence assembly and homology assignment process. However, this approach limits the inclusion of most genes of functional interest, which often belong to multi-gene families. Here, we demonstrate the feasibility of including large gene families in hybrid enrichment protocols for phylogeny reconstruction and subsequent analyses of molecular evolution, using a new set of bait sequences designed for the “portullugo” (Caryophyllales), a moderately sized lineage of flowering plants (~ 2200 species) that includes the cacti and harbors many evolutionary transitions to C$$_{mathrm{4}}$$ and CAM photosynthesis. Including multi-gene families allowed us to simultaneously infer a robust phylogeny and construct a dense sampling of sequences for a major enzyme of C$$_{mathrm{4}}$$ and CAM photosynthesis, which revealed the accumulation of adaptive amino acid substitutions associated with C$$_{mathrm{4}}$$ and CAM origins in particular paralogs. Our final set of matrices for phylogenetic analyses included 75–218 loci across 74 taxa, with ~ 50% matrix completeness across data sets. Phylogenetic resolution was greatly improved across the tree, at both shallow and deep levels. Concatenation and coalescent-based approaches both resolve the sister lineage of the cacti with strong support: Anacampserotaceae $$+$$ Portulacaceae, two lineages of mostly diminutive succulent herbs of warm, arid regions. In spite of this congruence, BUCKy concordance analyses demonstrated strong and conflicting signals across gene trees. Our results add to the growing number of examples illustrating the complexity of phylogenetic signals in genomic-scale data.

Smallpox holds a unique position in the history of medicine. It was the first disease for which a vaccine was developed and remains the only human disease eradicated by vaccination. Although there have been claims of smallpox in Egypt, India, and China dating back millennia [1–4], the timescale of emergence of the causative agent, variola virus (VARV), and how it evolved in the context of increasingly widespread immunization, have proven controversial [4–9]. In particular, some molecular-clock-based studies have suggested that key events in VARV evolution only occurred during the last two centuries [4–6] and hence in apparent conflict with anecdotal historical reports, although it is difficult to distinguish smallpox from other pustular rashes by description alone. To address these issues, we captured, sequenced, and reconstructed a draft genome of an ancient strain of VARV, sampled from a Lithuanian child mummy dating between 1643 and 1665 and close to the time of several documented European epidemics [1, 2, 10]. When compared to vaccinia virus, this archival strain contained the same pattern of gene degradation as 20th century VARVs, indicating that such loss of gene function had occurred before ca. 1650. Strikingly, the mummy sequence fell basal to all currently sequenced strains of VARV on phylogenetic trees. Molecular-clock analyses revealed a strong clock-like structure and that the timescale of smallpox evolution is more recent than often supposed, with the diversification of major viral lineages only occurring within the 18th and 19th centuries, concomitant with the development of modern vaccination.

Crassulacean acid metabolism (CAM) is a modified form of photosynthesis that has arisen independently at least 35 times in flowering plants. The occurrence of CAM is often correlated with shifts to arid, semiarid, or epiphytic habits, as well as transitions in leaf morphology (e.g. increased leaf thickness) and anatomy (e.g. increased cell size and packing). We assess shifts between C3 and CAM photosynthesis in the subfamily Agavoideae (Asparagaceae) through phylogenetic analysis of targeted loci captured from the nuclear and chloroplast genomes of over 60 species. Carbon isotope data was used as a proxy for mode of photosynthesis in extant species and ancestral states were estimated on the phylogeny. Ancestral character state mapping suggests three independent origins of CAM in the Agavoideae. CAM species differ from C3 species in climate space and are found to have thicker leaves with densely packed cells. C3 ancestors of CAM species show a predisposition toward CAM-like morphology. Leaf characteristics in the ancestral C3 species may have enabled the repeated evolution of CAM in the Agavoideae subfamily. Anatomical changes, including a tendency toward 3D venation, may have initially arisen in C3 ancestors in response to aridity as a way to increase leaf succulence for water storage.