The fossil record of appendicular muscle evolution in Synapsida on the line to mammals: Part II-Hindlimb

This paper is the second in a two-part series that charts the evolution of appendicular musculature along the mammalian stem lineage, drawing upon the exceptional fossil record of extinct synapsids. Here, attention is focused on muscles of the hindlimb. Although the hindlimb skeleton did not undergo as marked a transformation on the line to mammals as did the forelimb skeleton, the anatomy of extant tetrapods indicates that major changes to musculature have nonetheless occurred. To better understand these changes, this study surveyed the osteological evidence for muscular attachments in extinct mammalian and nonmammalian synapsids, two extinct amniote outgroups, and a large selection of extant mammals, saurians, and salamanders. Observations were integrated into an explicit phylogenetic framework, comprising 80 character-state complexes covering all muscles crossing the hip, knee, and ankle joints. These were coded for 33 operational taxonomic units spanning >330 Ma of tetrapod evolution, and ancestral state reconstruction was used to evaluate the sequence of muscular evolution along the stem lineage from Amniota to Theria. The evolutionary history of mammalian hindlimb musculature was complex, nonlinear, and protracted, with several instances of convergence and pulses of anatomical transformation that continued well into the crown group. Numerous traits typically regarded as characteristically "mammalian" have much greater antiquity than previously recognized, and for some traits, most synapsids are probably more reflective of the ancestral amniote condition than are extant saurians. More broadly, this study highlights the utility of the fossil record in interpreting the evolutionary appearance of distinctive anatomies.

Quantifying water-use efficiency in plant canopies with varying leaf angle and density distribution

BACKGROUND AND AIMS

Variation in architectural traits related to the spatial and angular distribution of leaf area can have considerable impacts on canopy-scale fluxes contributing to water-use efficiency (WUE). These architectural traits are frequent targets for crop improvement and for improving the understanding and predictions of net ecosystem carbon and water fluxes.

METHODS

A three-dimensional, leaf-resolving model along with a range of virtually generated hypothetical canopies were used to quantify interactions between canopy structure and WUE by examining its response to variation of leaf inclination independent of leaf azimuth, canopy heterogeneity, vegetation density and physiological parameters.

KEY RESULTS

Overall, increasing leaf area index (LAI), increasing the daily-averaged fraction of leaf area projected in the sun direction (Gavg) via the leaf inclination or azimuth distribution and increasing homogeneity had a similar effect on canopy-scale daily fluxes contributing to WUE. Increasing any of these parameters tended to increase daily light interception, increase daily net photosynthesis at low LAI and decrease it at high LAI, increase daily transpiration and decrease WUE. Isolated spherical crowns could decrease photosynthesis by ~60 % but increase daily WUE ≤130 % relative to a homogeneous canopy with equivalent leaf area density. There was no observed optimum in daily canopy WUE as LAI, leaf angle distribution or heterogeneity was varied. However, when the canopy was dense, a more vertical leaf angle distribution could increase both photosynthesis and WUE simultaneously.

CONCLUSIONS

Variation in leaf angle and density distributions can have a substantial impact on canopy-level carbon and water fluxes, with potential trade-offs between the two. These traits might therefore be viable target traits for increasing or maintaining crop productivity while using less water, and for improvement of simplified models. Increasing canopy density or decreasing canopy heterogeneity increases the impact of leaf angle on WUE and its dependent processes.

The scaling of metabolic traits differs among larvae and juvenile colonies of scleractinian corals

Body size profoundly affects organism fitness and ecosystem dynamics through the scaling of physiological traits. This study tests for variation in metabolic scaling and its potential drivers among corals differing in life history strategies and taxonomic identity. Data were compiled from published sources and augmented with empirical measurements of corals in Moorea, French Polynesia. The data compilation revealed metabolic isometry in broadcasted larvae, but size-independent metabolism in brooded larvae; empirical measures of Pocillopora acuta larvae also supported size-independent metabolism in brooded coral larvae. In contrast, for juvenile colonies (i.e., 1-4 cm diameter), metabolic scaling was isometric for Pocillopora spp., and negatively allometric for Porites spp. The scaling of biomass with surface area was isometric for Pocillopora spp., but positively allometric for Porites spp., suggesting the surface area:biomass ratio mediates metabolic scaling in these corals. The scaling of tissue biomass and metabolism were not affected by light treatment (i.e., either natural photoperiods or constant darkness) in both juvenile taxa. However, biomass was reduced by 9-15% in the juvenile corals from the light treatments and this coincided with higher metabolic scaling exponents, thus supporting the causal role of biomass in driving variation in scaling. This study shows that metabolic scaling is plastic in early life stages of corals, with intrinsic differences between life history strategy (i.e., brooded and broadcasted larvae) and taxa (i.e., Pocillopora spp. and Porites spp.), and acquired differences attributed to changes in area-normalized biomass.

Rhinoceros beetle (Trypoxylus dichotomus) cuticular hydrocarbons contain information about body size and sex

Japanese rhinoceros beetle (Trypoxylus dichotomus) males have exaggerated horns that are used to compete for territories. Larger males with larger horns tend to win these competitions, giving them access to females. Agonistic interactions include what appears to be assessment and often end without escalating to physical combat. However, it is unknown what information competitors use to assess each other. In many insect species chemical signals can carry a range of information, including social position, nutritional state, morphology, and sex. Specifically, cuticular hydrocarbons (CHCs), which are waxes excreted on the surface of insect exoskeletons, can communicate a variety of information. Here, we asked whether CHCs in rhinoceros beetles carry information about sex, body size, and condition that could be used by males during assessment behavior. Multivariate analysis of hydrocarbon composition revealed patterns associated with both sex and body size. We suggest that Rhinoceros beetles could be communicating information through CHCs that would explain behavioral decisions.

Triboelectric Energy Harvesting from Highly Conjugated Fused Aromatic Ladder Structure Under Extreme Environmental Conditions

Practical application of triboelectric nanogenerators (TENGs) has been challenging, particularly, under harsh environmental conditions. This work proposes a novel 3D-fused aromatic ladder (FAL) structure as a tribo-positive material for TENGs, to address these challenges. The 3D-FAL offers a unique materials engineering platform for tailored properties, such as high specific surface area and porosity, good thermal and mechanical stability, and tunable electronic properties. The fabricated 3D-FAL-based TENG reaches a maximum peak power density of 451.2 µW cm-2 at 5 Hz frequency. More importantly, the 3D-FAL-based TENG maintains stable output performance under harsh operating environments, over wide temperature (-45-100 °C) and humidity ranges (8.3-96.7% RH), representing the development of novel FAL for sustainable energy generation under challenging environmental conditions. Furthermore, the 3D-FAL-based TENG proves to be a promising device for a speed monitoring system engaging reconstruction in virtual reality in a snowy environment.

Granivore abundance shapes mutualism quality in plant-scatterhoarder interactions

Conditional mutualisms involve costs and benefits that vary with environmental factors, but mechanisms driving these dynamics remain poorly understood. Scatterhoarder-plant interactions are a prime example of this phenomenon, as scatterhoarders can either increase or reduce plant recruitment depending on the balance between seed dispersal and predation. We explored factors that drive the magnitude of net benefits for plants in this interaction using a mathematical model, with parameter values based on European beech (Fagus sylvatica) and yellow-necked mice (Apodemus flavicollis). We measured benefits as the percentage of germinating seeds, and examined how varying rodent survival (reflecting, e.g. changes in predation pressure), the rate of seed loss to other granivores, the abundance of alternative food resources, and changes in masting patterns affect the quality of mutualism. We found that increasing granivore abundance can degrade the quality of plant-scatterhoarder mutualism due to increased cache pilferage. Scatterhoarders are predicted to respond by increasing immediate consumption of gathered seeds, leading to higher costs and reduced benefits for plants. Thus, biotic changes that are detrimental to rodent populations can be beneficial for tree recruitment due to adaptive behavior of rodents. When scatterhoarder populations decline too drastically (< 5 individuals ha-1 ); however, tree recruitment may also suffer.

Maternal upbringing and selective breeding for voluntary exercise behavior modify patterns of DNA methylation and expression of genes in the mouse brain

Selective breeding has been utilized to study the genetic basis of exercise behavior, but research suggests that epigenetic mechanisms, such as DNA methylation, also contribute to this behavior. In a previous study, we demonstrated that the brains of mice from a genetically selected high runner (HR) line have sex-specific changes in DNA methylation patterns in genes known to be genomically imprinted compared to those from a non-selected control (C) line. Through cross-fostering, we also found that maternal upbringing can modify the DNA methylation patterns of additional genes. Here, we identify an additional set of genes in which DNA methylation patterns and gene expression may be altered by selection for increased wheel-running activity and maternal upbringing. We performed bisulfite sequencing and gene expression assays of 14 genes in the brain and found alterations in DNA methylation and gene expression for Bdnf, Pde4d and Grin2b. Decreases in Bdnf methylation correlated with significant increases in Bdnf gene expression in the hippocampus of HR compared to C mice. Cross-fostering also influenced the DNA methylation patterns for Pde4d in the cortex and Grin2b in the hippocampus, with associated changes in gene expression. We also found that the DNA methylation patterns for Atrx and Oxtr in the cortex and Atrx and Bdnf in the hippocampus were further modified by sex. Together with our previous study, these results suggest that DNA methylation and the resulting change in gene expression may interact with early-life influences to shape adult exercise behavior.

The name of the game: palaeoproteomics and radiocarbon dates further refine the presence and dispersal of caprines in eastern and southern Africa

We report the first large-scale palaeoproteomics research on eastern and southern African zooarchaeological samples, thereby refining our understanding of early caprine (sheep and goat) pastoralism in Africa. Assessing caprine introductions is a complicated task because of their skeletal similarity to endemic wild bovid species and the sparse and fragmentary state of relevant archaeological remains. Palaeoproteomics has previously proved effective in clarifying species attributions in African zooarchaeological materials, but few comparative protein sequences of wild bovid species have been available. Using newly generated type I collagen sequences for wild species, as well as previously published sequences, we assess species attributions for elements originally identified as caprine or 'unidentifiable bovid' from 17 eastern and southern African sites that span seven millennia. We identified over 70% of the archaeological remains and the direct radiocarbon dating of domesticate specimens allows refinement of the chronology of caprine presence in both African regions. These results thus confirm earlier occurrences in eastern Africa and the systematic association of domesticated caprines with wild bovids at all archaeological sites. The combined biomolecular approach highlights repeatability and accuracy of the methods for conclusive contribution in species attribution of archaeological remains in dry African environments.

Pathogen-induced m6A dynamics affect plant immunity

Posttranscriptional regulation of mRNA mediated by methylation at the N6 position of adenine (N6-methyladenosine [m6A]) has profound effects on transcriptome regulation in plants. Focused studies across eukaryotes offer glimpses into the processes governed by m6A throughout developmental and disease states. However, we lack an understanding of the dynamics and the regulatory potential of m6A during biotic stress in plants. Here, we provide a comprehensive look into the effects of m6A on both the short-term and long-term responses to pathogen signaling in Arabidopsis (Arabidopsis thaliana). We demonstrate that m6A-deficient plants are more resistant to bacterial and fungal pathogen infections and have altered immune responses. Furthermore, m6A deposition is specifically coordinated on transcripts involved in defense and immunity prior to and proceeding the pathogen signal flagellin. Consequently, the dynamic modulation of m6A on specific stress-responsive transcripts is correlated with changes in abundance and cleavage of these transcripts. Overall, we show that the m6A methylome is regulated prior to and during simulated and active pathogen stress and functions in the coordination and balancing of normal growth and pathogen responses.

A Comprehensive Study of Techniques to Optimize the Extraction of Lipids from the Autotrophic Strain of the Microalgae Chlorella vulgaris

Microalgae represent a promising source of triglycerides and free fatty acids, synthesized in the form of lipids, for use in renewable fuels and chemicals. One challenge is the ability to efficiently recover the lipids from within the microalgae cell. Although various techniques have been studied individually, a comprehensive study of extraction techniques using consistent experimental and analytical methodology is missing. This study aims to provide this unifying comparison using the common microalgae strain Chlorella vulgaris. The factors that were surveyed and then optimized to achieve maximum extraction efficiency included the solvent type; mechanical pre-treatment using a ball mill at a variety of grinding speeds; microalgae-to-solvent ratio; extraction facilitated by microwave; extraction facilitated by sonication; extraction facilitated using increased temperature; and extraction facilitated by in situ transesterification to convert the lipids into esters prior to extraction. The optimum conditions determined during these studies were utilizing methanol as the solvent, with ball mill pretreatment at a grinding speed of 500 rpm, and a 1:9 microalgae to solvent ratio. When used in combination with microwave-assisted extraction at a temperature of 140 °C, approximately 24 wt% of the initial lipids were recovered. Recoveries of over 70 wt% were obtained without a microwave at extraction temperatures of over 200 °C.

Historical biogeography supports Point Conception as the site of turnover between temperate East Pacific ichthyofaunas

The cold temperate and subtropical marine faunas of the Northeastern Pacific meet within California as part of one of the few eastern boundary upwelling ecosystems in the world. Traditionally, it is believed that Point Conception is the precise site of turnover between these two faunas due to sharp changes in oceanographic conditions. However, evidence from intraspecific phylogeography and species range terminals do not support this view, finding stronger biogeographic breaks elsewhere along the coast. Here I develop a new application of historical biogeographic approaches to uncover sites of transition between faunas without needing an a priori hypothesis of where these occur. I used this approach to determine whether the point of transition between northern and southern temperate faunas occurs at Point Conception or elsewhere within California. I also examined expert-vetted latitudinal range data of California fish species from the 1970s and the 2020s to assess how biogeography could change with the backdrop of climate change. The site of turnover was found to occur near Point Conception, in concordance with the traditional view. I suggest that recent species- and population-level processes could be expected to give signals of different events from historical biogeography, possibly explaining the discrepancy across studies. Species richness of California has increased since the 1970s, mostly due to species's ranges expanding northward from Baja California (Mexico). Range shifts under warming conditions seem to be increasing the disparity between northern and southern faunas of California, creating a more divergent biogeography.

Toward certifiable optimal motion planning for medical steerable needles

Medical steerable needles can follow 3D curvilinear trajectories to avoid anatomical obstacles and reach clinically significant targets inside the human body. Automating steerable needle procedures can enable physicians and patients to harness the full potential of steerable needles by maximally leveraging their steerability to safely and accurately reach targets for medical procedures such as biopsies. For the automation of medical procedures to be clinically accepted, it is critical from a patient care, safety, and regulatory perspective to certify the correctness and effectiveness of the planning algorithms involved in procedure automation. In this paper, we take an important step toward creating a certifiable optimal planner for steerable needles. We present an efficient, resolution-complete motion planner for steerable needles based on a novel adaptation of multi-resolution planning. This is the first motion planner for steerable needles that guarantees to compute in finite time an obstacle-avoiding plan (or notify the user that no such plan exists), under clinically appropriate assumptions. Based on this planner, we then develop the first resolution-optimal motion planner for steerable needles that further provides theoretical guarantees on the quality of the computed motion plan, that is, global optimality, in finite time. Compared to state-of-the-art steerable needle motion planners, we demonstrate with clinically realistic simulations that our planners not only provide theoretical guarantees but also have higher success rates, have lower computation times, and result in higher quality plans.

Exploring Conflicts in Whole Genome Phylogenetics: A Case Study Within Manakins (Aves: Pipridae)

Some phylogenetic problems remain unresolved even when large amounts of sequence data are analyzed and methods that accommodate processes such as incomplete lineage sorting are employed. In addition to investigating biological sources of phylogenetic incongruence, it is also important to reduce noise in the phylogenomic dataset by using appropriate filtering approach that addresses gene tree estimation errors. We present the results of a case study in manakins, focusing on the very difficult clade comprising the genera Antilophia and Chiroxiphia. Previous studies suggest that Antilophia is nested within Chiroxiphia, though relationships among Antilophia+Chiroxiphia species have been highly unstable. We extracted more than 11,000 loci (ultra-conserved elements and introns) from whole genomes and conducted analyses using concatenation and multispecies coalescent methods. Topologies resulting from analyses using all loci differed depending on the data type and analytical method, with 2 clades (Antilophia+Chiroxiphia and Manacus+Pipra+Machaeopterus) in the manakin tree showing incongruent results. We hypothesized that gene trees that conflicted with a long coalescent branch (e.g., the branch uniting Antilophia+Chiroxiphia) might be enriched for cases of gene tree estimation error, so we conducted analyses that either constrained those gene trees to include monophyly of Antilophia+Chiroxiphia or excluded these loci. While constraining trees reduced some incongruence, excluding the trees led to completely congruent species trees, regardless of the data type or model of sequence evolution used. We found that a suite of gene metrics (most importantly the number of informative sites and likelihood of intralocus recombination) collectively explained the loci that resulted in non-monophyly of Antilophia+Chiroxiphia. We also found evidence for introgression that may have contributed to the discordant topologies we observe in Antilophia+Chiroxiphia and led to deviations from expectations given the multispecies coalescent model. Our study highlights the importance of identifying factors that can obscure phylogenetic signal when dealing with recalcitrant phylogenetic problems, such as gene tree estimation error, incomplete lineage sorting, and reticulation events. [Birds; c-gene; data type; gene estimation error; model fit; multispecies coalescent; phylogenomics; reticulation].

Bioinspired stability enhancement in deuterium-substituted organic-inorganic hybrid perovskite solar cells

In hybrid perovskite solar cells (PSCs), the reaction of hydrogens (H) located in the amino group of the organic A-site cations with their neighboring halides plays a central role in degradation. Inspired by the retarded biological activities of cells in heavy water, we replaced the light H atom with its abundant, twice-as-heavy, nonradioactive isotope, deuterium (D) to hamper the motion of H. This D substitution retarded the formation kinetics of the detrimental H halides in Pb-based PSCs, as well as the H bond-mediated oxidation of Sn²⁺ in Sn-Pb-based narrow-bandgap PSCs, evidenced by accelerated stability studies. A computational study indicated that the zero point energy of D-based formamidinium (FA) is lower than that of pristine FA. In addition, the smaller increase in entropy in D-based FA than in pristine FA accounts for the increased formation free energy of the Sn²⁺ vacancies, which leads to the retarded oxidation kinetics of Sn²⁺. In this study, we show that substituting active H with D in organic cations is an effective way to enhance the stability of PSCs without sacrificing photovoltaic (PV) performance. This approach is also adaptable to other stabilizing methods.

Phylogenomics reveals deep relationships and diversification within phylactolaemate bryozoans

Bryozoans are mostly sessile colonial invertebrates that inhabit all kinds of aquatic ecosystems. Extant bryozoan species fall into two clades with one of them, Phylactolaemata, being the only exclusively freshwater clade. Phylogenetic relationships within the class Phylactolaemata have long been controversial owing to their limited distinguishable characteristics that reflect evolutionary relationships. Here, we present the first phylogenomic analysis of Phylactolaemata using transcriptomic data combined with dense taxon sampling of six families to better resolve the interrelationships and to estimate divergence time. Using maximum-likelihood and Bayesian inference approaches, we recovered a robust phylogeny for Phylactolaemata in which the interfamilial relationships are fully resolved. We show Stephanellidae is the sister taxon of all other phylactolaemates and confirm that Lophopodidae represents the second offshoot within the phylactolaemate tree. Plumatella fruticosa clearly falls outside Plumatellidae as previous investigations have suggested, and instead clusters with Pectinatellidae and Cristatellidae as the sister taxon of Fredericellidae. Our results demonstrate that cryptic speciation is very likely in F. sultana and in two species of Plumatella (P. repens and P. casmiana). Divergence time estimates show that Phylactolaemata appeared at the end of the Ediacaran and started to diverge in the Silurian, although confidence intervals were large for most nodes. The radiation of most extant phylactolaemate families occurred mainly in the Palaeogene and Neogene highlighting post-extinction diversification.

Disentangling the developmental origins of a novel phenotype: enhancement versus reversal of environmentally induced gene expression

Increasing evidence suggests that many novel traits might have originated via plasticity-led evolution (PLE). Yet, little is known of the developmental processes that underpin PLE, especially in its early stages. One such process is 'phenotypic accommodation', which occurs when, in response to a change in the environment, an organism experiences adjustments across variable parts of its phenotype that improve its fitness. Here, we asked if environmentally induced changes in gene expression are enhanced or reversed during phenotypic accommodation of a novel, complex phenotype in spadefoot toad tadpoles (Spea multiplicata). More genes than expected were affected by both the environment and phenotypic accommodation in the liver and brain. However, although phenotypic accommodation primarily reversed environmentally induced changes in gene expression in liver tissue, it enhanced these changes in brain tissue. Thus, depending on the tissue, phenotypic accommodation may either minimize functional disruption via reversal of gene expression patterns or promote novelty via enhancement of existing expression patterns. Our study thereby provides insights into the developmental origins of a novel phenotype and the incipient stages of PLE.

Molecular mechanisms underlying host-induced gene silencing

Host-induced gene silencing (HIGS) refers to the silencing of genes in pathogens and pests by expressing homologous double-stranded RNAs (dsRNA) or artificial microRNAs (amiRNAs) in the host plant. The discovery of such trans-kingdom RNA silencing has enabled the development of RNA interference-based approaches for controlling diverse crop pathogens and pests. Although HIGS is a promising strategy, the mechanisms by which these regulatory RNAs translocate from plants to pathogens, and how they induce gene silencing in pathogens, are poorly understood. This lack of understanding has led to large variability in the efficacy of various HIGS treatments. This variability is likely due to multiple factors, such as the ability of the target pathogen or pest to take up and/or process RNA from the host, the specific genes and target sequences selected in the pathogen or pest for silencing, and where, when, and how the dsRNAs or amiRNAs are produced and translocated. In this review, we summarize what is currently known about the molecular mechanisms underlying HIGS, identify key unanswered questions, and explore strategies for improving the efficacy and reproducibility of HIGS treatments in the control of crop diseases.

Substrate and Process Engineering for Biocatalytic Synthesis and Facile Purification of Human Milk Oligosaccharides

Innovation in process development is essential for applying biocatalysis in industrial and laboratory production of organic compounds, including beneficial carbohydrates such as human milk oligosaccharides (HMOs). HMOs have attracted increasing attention for their potential application as key ingredients in products that can improve human health. To efficiently access HMOs through biocatalysis, a combined substrate and process engineering strategy is developed, namely multistep one-pot multienzyme (MSOPME) design. The strategy allows access to a pure tagged HMO in a single reactor with a single C18-cartridge purification process, despite the length of the target. Its efficiency is demonstrated in the high-yielding (71-91 %) one-pot synthesis of twenty tagged HMOs (83-155 mg), including long-chain oligosaccharides with or without fucosylation or sialylation up to nonaoses from a lactoside without the isolation of the intermediate oligosaccharides. Gram-scale synthesis of an important HMO derivative - tagged lacto-N-fucopentaose-I (LNFP-I) - proceeds in 84 % yield. Tag removal is carried out in high efficiency (94-97 %) without the need for column purification to produce the desired natural HMOs with a free reducing end. The method can be readily adapted for large-scale synthesis and automation to allow quick access to HMOs, other glycans, and glycoconjugates.

Electro-assembly of a dynamically adaptive molten fibril state for collagen

Collagen is a biological building block that is hierarchically assembled into diverse morphological structures that, in some cases, is dynamically adaptive in response to external cues and in other cases forms static terminal structures. Technically, there is limited capabilities to guide the emergence of collagen's hierarchical organization to recapitulate the richness of biological structure and function. Here, we report an electro-assembly pathway to create a dynamically adaptive intermediate molten fibril state for collagen. Structurally, this intermediate state is composed of partially aligned and reversibly associating fibrils with limited hierarchical structure. These molten fibrils can be reversibly reconfigured to offer dynamic properties such as stimuli-stiffening, stimuli-contracting, self-healing, and self-shaping. Also, molten fibrils can be guided to further assemble to recapitulate the characteristic hierarchical structural features of native collagen (e.g., aligned fibers with D-banding). We envision that the electro-assembly of collagen fibrils will provide previously unidentified opportunities for tailored collagen-based biomedical materials.

Phytolith evidence for the pastoral origins of multi-cropping in Mesopotamia (ancient Iraq)

Multi-cropping was vital for provisioning large population centers across ancient Eurasia. In Southwest Asia, multi-cropping, in which grain, fodder, or forage could be reliably cultivated during dry summer months, only became possible with the translocation of summer grains, like millet, from Africa and East Asia. Despite some textual sources suggesting millet cultivation as early as the third millennium BCE, the absence of robust archaeobotanical evidence for millet in semi-arid Mesopotamia (ancient Iraq) has led most archaeologists to conclude that millet was only grown in the region after the mid-first millennium BCE introduction of massive, state-sponsored irrigation systems. Here, we present the earliest micro-botanical evidence of the summer grain broomcorn millet (Panicum miliaceum) in Mesopotamia, identified using phytoliths in dung-rich sediments from Khani Masi, a mid-second millennium BCE site located in northern Iraq. Taphonomic factors associated with the region's agro-pastoral systems have likely made millet challenging to recognize using conventional macrobotanical analyses, and millet may therefore have been more widespread and cultivated much earlier in Mesopotamia than is currently recognized. The evidence for pastoral-related multi-cropping in Bronze Age Mesopotamia provides an antecedent to first millennium BCE agricultural intensification and ties Mesopotamia into our rapidly evolving understanding of early Eurasian food globalization.

Correlated mechanochemical maps of Arabidopsis thaliana primary cell walls using atomic force microscope infrared spectroscopy

Spatial heterogeneity in composition and organisation of the primary cell wall affects the mechanics of cellular morphogenesis. However, directly correlating cell wall composition, organisation and mechanics has been challenging. To overcome this barrier, we applied atomic force microscopy coupled with infrared (AFM-IR) spectroscopy to generate spatially correlated maps of chemical and mechanical properties for paraformaldehyde-fixed, intact Arabidopsis thaliana epidermal cell walls. AFM-IR spectra were deconvoluted by non-negative matrix factorisation (NMF) into a linear combination of IR spectral factors representing sets of chemical groups comprising different cell wall components. This approach enables quantification of chemical composition from IR spectral signatures and visualisation of chemical heterogeneity at nanometer resolution. Cross-correlation analysis of the spatial distribution of NMFs and mechanical properties suggests that the carbohydrate composition of cell wall junctions correlates with increased local stiffness. Together, our work establishes new methodology to use AFM-IR for the mechanochemical analysis of intact plant primary cell walls.

Causation without correlation: parasite-mediated frequency-dependent selection and infection prevalence

Parasite-mediated selection is thought to maintain host genetic diversity for resistance. We might thus expect to find a strong positive correlation between host genetic diversity and infection prevalence across natural populations. Here, we used computer simulations to examine host-parasite coevolution in 20 simi-isolated clonal populations across a broad range of values for both parasite virulence and parasite fecundity. We found that the correlation between host genetic diversity and infection prevalence can be significantly positive for intermediate values of parasite virulence and fecundity. But the correlation can also be weak and statistically non-significant, even when parasite-mediated frequency-dependent selection is the sole force maintaining host diversity. Hence correlational analyses of field populations, while useful, might underestimate the role of parasites in maintaining host diversity.

A chromosome-scale genome assembly and karyotype of the ctenophore Hormiphora californensis

Here, we present a karyotype, a chromosome-scale genome assembly, and a genome annotation from the ctenophore Hormiphora californensis (Ctenophora: Cydippida: Pleurobrachiidae). The assembly spans 110 Mb in 44 scaffolds and 99.47% of the bases are contained in 13 scaffolds. Chromosome micrographs and Hi-C heatmaps support a karyotype of 13 diploid chromosomes. Hi-C data reveal three large heterozygous inversions on chromosome 1, and one heterozygous inversion shares the same gene order found in the genome of the ctenophore Pleurobrachia bachei. We find evidence that H. californensis and P. bachei share thirteen homologous chromosomes, and the same karyotype of 1n = 13. The manually curated PacBio Iso-Seq-based genome annotation reveals complex gene structures, including nested genes and trans-spliced leader sequences. This chromosome-scale assembly is a useful resource for ctenophore biology and will aid future studies of metazoan evolution and phylogenetics.

Exploring multiuser virtual teaching simulation as an alternative learning environment for student instructors

In this mixed-method study, we investigated the impact and design of a multiuser, virtual reality (VR) supported teaching simulation, in comparison with live classroom teaching simulation, on the participatory training of teaching and the teaching knowledge development of student instructors. A total of 40 university teaching assistants participated in a 4-h training session in which they were randomly assigned to a VR simulation or a live classroom simulation condition. The study indicated that the VR simulation better promoted the lab-teaching knowledge development than the live simulation, whereas the latter better fostered class-teaching knowledge development. All participants reported higher teaching self-efficacy after the training. The qualitative data indicated that domain-specific challenges and authentic environmental prompting in the VR simulation fostered both experiential and vicarious learning of teaching. However, VR participants lacked mutual engagement in collaborative role-playing. The study findings suggest that VR-based simulation can supplement and work as an alternative to the live classroom simulation to host participatory teaching development.

PIF4 negatively modulates cold tolerance in tomato anthers via temperature-dependent regulation of tapetal cell death

Extreme temperature conditions seriously impair male reproductive development in plants; however, the molecular mechanisms underlying the response of anthers to extreme temperatures remain poorly described. The transcription factor phytochrome-interacting factor4 (PIF4) acts as a hub that integrates multiple signaling pathways to regulate thermosensory growth and architectural adaptation in plants. Here, we report that SlPIF4 in tomato (Solanum lycopersicum) plays a pivotal role in regulating cold tolerance in anthers. CRISPR (clustered regularly interspaced short palindromic repeats)-associated nuclease Cas9-generated SlPIF4 knockout mutants showed enhanced cold tolerance in pollen due to reduced temperature sensitivity of the tapetum, while overexpressing SlPIF4 conferred pollen abortion by delaying tapetal programmed cell death (PCD). SlPIF4 directly interacts with SlDYT1, a direct upstream regulator of SlTDF1, both of which (SlDYT1 and SlTDF1) play important roles in regulating tapetum development and tapetal PCD. Moderately low temperature (MLT) promotes the transcriptional activation of SlTDF1 by the SlPIF4-SlDYT1 complex, resulting in pollen abortion, while knocking out SlPIF4 blocked the MLT-induced activation of SlTDF1. Furthermore, SlPIF4 directly binds to the canonical E-box sequence in the SlDYT1 promoter. Collectively, these findings suggest that SlPIF4 negatively regulates cold tolerance in anthers by directly interacting with the tapetal regulatory module in a temperature-dependent manner. Our results shed light on the molecular mechanisms underlying the adaptation of anthers to low temperatures.

Measuring global multi-scale place connectivity using geotagged social media data

Shaped by human movement, place connectivity is quantified by the strength of spatial interactions among locations. For decades, spatial scientists have researched place connectivity, applications, and metrics. The growing popularity of social media provides a new data stream where spatial social interaction measures are largely devoid of privacy issues, easily assessable, and harmonized. In this study, we introduced a global multi-scale place connectivity index (PCI) based on spatial interactions among places revealed by geotagged tweets as a spatiotemporal-continuous and easy-to-implement measurement. The multi-scale PCI, demonstrated at the US county level, exhibits a strong positive association with SafeGraph population movement records (10% penetration in the US population) and Facebook's social connectedness index (SCI), a popular connectivity index based on social networks. We found that PCI has a strong boundary effect and that it generally follows the distance decay, although this force is weaker in more urbanized counties with a denser population. Our investigation further suggests that PCI has great potential in addressing real-world problems that require place connectivity knowledge, exemplified with two applications: (1) modeling the spatial spread of COVID-19 during the early stage of the pandemic and (2) modeling hurricane evacuation destination choice. The methodological and contextual knowledge of PCI, together with the open-sourced PCI datasets at various geographic levels, are expected to support research fields requiring knowledge in human spatial interactions.

High-Quality Genome Assembly and Comprehensive Transcriptome of the Painted Lady Butterfly Vanessa cardui

The painted lady butterfly, Vanessa cardui, has the longest migration routes, the widest hostplant diversity, and one of the most complex wing patterns of any insect. Due to minimal culturing requirements, easily characterized wing pattern elements, and technical feasibility of CRISPR/Cas9 genome editing, V. cardui is emerging as a functional genomics model for diverse research programs. Here, we report a high-quality, annotated genome assembly of the V. cardui genome, generated using 84× coverage of PacBio long-read data, which we assembled into 205 contigs with a total length of 425.4 Mb (N50 = 10.3 Mb). The genome was very complete (single-copy complete Benchmarking Universal Single-Copy Orthologs [BUSCO] 97%), with contigs assembled into presumptive chromosomes using synteny analyses. Our annotation used embryonic, larval, and pupal transcriptomes, and 20 transcriptomes across five different wing developmental stages. Gene annotations showed a high level of accuracy and completeness, with 14,437 predicted protein-coding genes. This annotated genome assembly constitutes an important resource for diverse functional genomic studies ranging from the developmental genetic basis of butterfly color pattern, to coevolution with diverse hostplants.

Domain-driven models yield better predictions at lower cost than reservoir computers in Lorenz systems

Recent advances in computing algorithms and hardware have rekindled interest in developing high-accuracy, low-cost surrogate models for simulating physical systems. The idea is to replace expensive numerical integration of complex coupled partial differential equations at fine time scales performed on supercomputers, with machine-learned surrogates that efficiently and accurately forecast future system states using data sampled from the underlying system. One particularly popular technique being explored within the weather and climate modelling community is the echo state network (ESN), an attractive alternative to other well-known deep learning architectures. Using the classical Lorenz 63 system, and the three tier multi-scale Lorenz 96 system (Thornes T, Duben P, Palmer T. 2017 Q. J. R. Meteorol. Soc. 143, 897-908. (doi:10.1002/qj.2974)) as benchmarks, we realize that previously studied state-of-the-art ESNs operate in two distinct regimes, corresponding to low and high spectral radius (LSR/HSR) for the sparse, randomly generated, reservoir recurrence matrix. Using knowledge of the mathematical structure of the Lorenz systems along with systematic ablation and hyperparameter sensitivity analyses, we show that state-of-the-art LSR-ESNs reduce to a polynomial regression model which we call Domain-Driven Regularized Regression (D2R2). Interestingly, D2R2 is a generalization of the well-known SINDy algorithm (Brunton SL, Proctor JL, Kutz JN. 2016 Proc. Natl Acad. Sci. USA 113, 3932-3937. (doi:10.1073/pnas.1517384113)). We also show experimentally that LSR-ESNs (Chattopadhyay A, Hassanzadeh P, Subramanian D. 2019 (http://arxiv.org/abs/1906.08829)) outperform HSR ESNs (Pathak J, Hunt B, Girvan M, Lu Z, Ott E. 2018 Phys. Rev. Lett. 120, 024102. (doi:10.1103/PhysRevLett.120.024102)) while D2R2 dominates both approaches. A significant goal in constructing surrogates is to cope with barriers to scaling in weather prediction and simulation of dynamical systems that are imposed by time and energy consumption in supercomputers. Inexact computing has emerged as a novel approach to helping with scaling. In this paper, we evaluate the performance of three models (LSR-ESN, HSR-ESN and D2R2) by varying the precision or word size of the computation as our inexactness-controlling parameter. For precisions of 64, 32 and 16 bits, we show that, surprisingly, the least expensive D2R2 method yields the most robust results and the greatest savings compared to ESNs. Specifically, D2R2 achieves 68 × in computational savings, with an additional 2 × if precision reductions are also employed, outperforming ESN variants by a large margin. This article is part of the theme issue 'Machine learning for weather and climate modelling'.

Lessons from movement ecology for the return to work: Modeling contacts and the spread of COVID-19

Human behavior (movement, social contacts) plays a central role in the spread of pathogens like SARS-CoV-2. The rapid spread of SARS-CoV-2 was driven by global human movement, and initial lockdown measures aimed to localize movement and contact in order to slow spread. Thus, movement and contact patterns need to be explicitly considered when making reopening decisions, especially regarding return to work. Here, as a case study, we consider the initial stages of resuming research at a large research university, using approaches from movement ecology and contact network epidemiology. First, we develop a dynamical pathogen model describing movement between home and work; we show that limiting social contact, via reduced people or reduced time in the workplace are fairly equivalent strategies to slow pathogen spread. Second, we develop a model based on spatial contact patterns within a specific office and lab building on campus; we show that restricting on-campus activities to labs (rather than labs and offices) could dramatically alter (modularize) contact network structure and thus, potentially reduce pathogen spread by providing a workplace mechanism to reduce contact. Here we argue that explicitly accounting for human movement and contact behavior in the workplace can provide additional strategies to slow pathogen spread that can be used in conjunction with ongoing public health efforts.

The development of affiliative and coercive reproductive tactics in male chimpanzees

Like many animals, adult male chimpanzees often compete for a limited number of mates. They fight other males as they strive for status that confers reproductive benefits and use aggression to coerce females to mate with them. Nevertheless, small-bodied, socially immature adolescent male chimpanzees, who cannot compete with older males for status nor intimidate females, father offspring. We investigated how they do so through a study of adolescent and young adult males at Ngogo in Kibale National Park, Uganda. Adolescent males mated with nulliparous females and reproduced primarily with these first-time mothers, who are not preferred as mating partners by older males. Two other factors, affiliation and aggression, also influenced mating success. Specifically, the strength of affiliative bonds that males formed with females and the amount of aggression males directed toward females predicted male mating success. The effect of male aggression toward females on mating success increased as males aged, especially when they directed it toward females with whom they shared affiliative bonds. These results mirror sexual coercion in humans, which occurs most often between males and females involved in close, affiliative relationships.

Characterisation of charred organic matter in micromorphological thin sections by means of Raman spectroscopy

Burned or charred organic matter in anthropogenic combustion features may provide important clues about past human activities related to fire. To interpret archaeological hearths, a correct identification of the organic source material is key. In the present work, Raman spectroscopy is applied to characterise the structural properties of char produced in laboratory heating- and open-fire experiments. This reference data set is compared to analyses of three different archaeological sites with Middle Palaeolithic combustion contexts. The results show that it is possible to determine whether a charred fragment is the product of burning animal-derived matter (e.g. meat) or plant-derived matter (e.g. wood) by plotting a few Raman spectral parameters (i.e. position of G and D bands, and intensity ratios H _D/H _G and H _V/H _G) against one another. The most effective parameters for discriminating animal- from plant-derived matter are the position of the G band and the H _V/H _G intensity ratio. This method can be applied on raw sample material and on uncovered micromorphological thin sections. The latter greatly compliments micromorphology by providing information about char fragments without any clear morphological characteristics. This study is the first of its kind and may provide archaeologists with a robust new method to distinguish animal- from plant-derived char in thin sections.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12520-020-01263-3.

Matrix models and deformations of JT gravity

Recently, it has been found that Jackiw-Teitelboim (JT) gravity, which is a two-dimensional theory with bulk action - 1 / 2 ∫ d 2 x g ϕ ( R + 2 ) , is dual to a matrix model, that is, a random ensemble of quantum systems rather than a specific quantum mechanical system. In this article, we argue that a deformation of JT gravity with bulk action - 1 / 2 ∫ d 2 x g ( ϕ R + W ( ϕ ) ) is likewise dual to a matrix model. With a specific procedure for defining the path integral of the theory, we determine the density of eigenvalues of the dual matrix model. There is a simple answer if W(0) = 0, and otherwise a rather complicated answer.

The effect of sleep on public good contributions and punishment: Experimental evidence

We investigate the effect of a full week of sleep restriction (SR) vs. well-restedness (WR) on contributions in a common public good experiment, the voluntary contributions mechanism (VCM). We examine the effect of sleep manipulation on decisions regarding both contributions and punishment of non-contributors. Actigraphy devices are used to confirm that our random assignment to sleep condition generates significant differences in objective nightly sleep duration and sleepiness. We find that when punishment is unavailable public good contributions do not differ by SR/WR assignment. When punishment is available, we find evidence that SR subjects contribute more than WR subjects, respond more to the availability of punishment than do WR subjects, and that the availability of punishment significantly increases the contributions of SR but not WR subjects. Yet SR subjects do not punish others more or less than WR subjects. Our main findings are robust when considering compliance and sample selection. However, some findings are not robust to an alternative but less objective sleep control measure that is based partly on participants' self-identified optimal sleep levels.

River networks and funerary metal in the Bronze Age of the Carpathian Basin

Archaeologists use differences in metals from burial contexts to identify variation in social inequalities during the European Bronze Age. Many have argued that these social inequalities depended on access to, and control of, trade routes. In this paper, I model critical gateways in the Tisza river-a river system in the Carpathian Basin that might have enabled privileged access to metal in some areas but not others. I then evaluate the concentration of metal on different topological nodes of the river network in an attempt to understand what best explains the distribution of metals across this landscape. I do this by describing Bronze Age metal consumption and display in cemeteries from four micro-regions of the Tisza, and compare them with network 'betweenness centrality' values for locations along the river. I find support for the argument that favourably located river nodes had better access to metal in the earlier part of the Bronze Age.

Load Position and Weight Classification during Carrying Gait Using Wearable Inertial and Electromyographic Sensors

Lifting and carrying heavy objects is a major aspect of physically intensive jobs. Wearable sensors have previously been used to classify different ways of picking up an object, but have seen only limited use for automatic classification of load position and weight while a person is walking and carrying an object. In this proof-of-concept study, we thus used wearable inertial and electromyographic sensors for offline classification of different load positions (frontal vs. unilateral vs. bilateral side loads) and weights during gait. Ten participants performed 19 different carrying trials each while wearing the sensors, and data from these trials were used to train and evaluate classification algorithms based on supervised machine learning. The algorithms differentiated between frontal and other loads (side/none) with an accuracy of 100%, between frontal vs. unilateral side load vs. bilateral side load with an accuracy of 96.1%, and between different load asymmetry levels with accuracies of 75-79%. While the study is limited by a lack of electromyographic sensors on the arms and a limited number of load positions/weights, it shows that wearable sensors can differentiate between different load positions and weights during gait with high accuracy. In the future, such approaches could be used to control assistive devices or for long-term worker monitoring in physically demanding occupations.

Ontogenetic variability in crystallography and mosaicity of conodont apatite: implications for microstructure, palaeothermometry and geochemistry

X-ray diffraction data from Silurian conodonts belonging to various developmental stages of the species Dapsilodus obliquicostatus demonstrate changes in crystallography and degree of nanocrystallite ordering (mosaicity) in both lamellar crown tissue and white matter. The exclusive use of a single species in this study, combined with systematic testing of each element type at multiple locations, provided insight into microstructural and crystallographic differentiation between element type (Sa , Sb -c , M) as well as between juveniles and adults. A relative increase in the unit cell dimensions a/c ratio of nanocrystallites during growth was apparent in areas demonstrating single-crystal behaviour, but no such relationship was seen in dominantly polycrystalline areas. Systematic variations in mosaicity were identified, with mosaicity (as a proxy for disorder) increasing during growth, as well as along elements from tip to base. These results provide potential insight into the integrity of conodont apatite as a recorder of palaeoseawater chemistry, as well as demonstrate the need to consider the influence of ontogeny and element type on the use of conodonts in palaeothermometry and geochemical investigations.

Position-wise binding preference is important for miRNA target site prediction

MOTIVATION

It is a fundamental task to identify microRNAs (miRNAs) targets and accurately locate their target sites. Genome-scale experiments for miRNA target site detection are still costly. The prediction accuracies of existing computational algorithms and tools are often not up to the expectation due to a large number of false positives. One major obstacle to achieve a higher accuracy is the lack of knowledge of the target binding features of miRNAs. The published high-throughput experimental data provide an opportunity to analyze position-wise preference of miRNAs in terms of target binding, which can be an important feature in miRNA target prediction algorithms.

RESULTS

We developed a Markov model to characterize position-wise pairing patterns of miRNA-target interactions. We further integrated this model as a scoring method and developed a dynamic programming (DP) algorithm, MDPS (Markov model-scored Dynamic Programming algorithm for miRNA target site Selection) that can screen putative target sites of miRNA-target binding. The MDPS algorithm thus can take into account both the dependency of neighboring pairing positions and the global pairing information. Based on the trained Markov models from both miRNA-specific and general datasets, we discovered that the position-wise binding information specific to a given miRNA would benefit its target prediction. We also found that miRNAs maintain region-wise similarity in their target binding patterns. Combining MDPS with existing methods significantly improves their precision while only slightly reduces their recall. Therefore, position-wise pairing patterns have the promise to improve target prediction if incorporated into existing software tools.

AVAILABILITY AND IMPLEMENTATION

The source code and tool to calculate MDPS score is available at http://hulab.ucf.edu/research/projects/MDPS/index.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Duplication and functional divergence of a calcium sensor in the Brassicaceae

The presence of varied numbers of CALCINEURIN B-LIKE10 (CBL10) calcium sensor genes in species across the Brassicaceae and the demonstrated role of CBL10 in salt tolerance in Arabidopsis thaliana and Eutrema salsugineum provided a unique opportunity to determine if CBL10 function is modified in different species and linked to salt tolerance. Salinity effects on species growth and cross-species complementation were used to determine the extent of conservation and divergence of CBL10 function in four species representing major lineages within the core Brassicaceae (A. thaliana, E. salsugineum, Schrenkiella parvula, and Sisymbrium irio) as well as the first diverging lineage (Aethionema arabicum). Evolutionary and functional analyses indicate that CBL10 duplicated within expanded lineage II of the Brassicaceae and that, while portions of CBL10 function are conserved across the family, there are species-specific variations in CBL10 function. Paralogous CBL10 genes within a species diverged in expression and function probably contributing to the maintenance of the duplicated gene pairs. Orthologous CBL10 genes diverged in function in a species-specific manner, suggesting that functions arose post-speciation. Multiple CBL10 genes and their functional divergence may have expanded calcium-mediated signaling responses and contributed to the ability of certain members of the Brassicaceae to maintain growth in salt-affected soils.

A catastrophic tropical drought kills hydraulically vulnerable tree species

Drought-related tree mortality is now a widespread phenomenon predicted to increase in magnitude with climate change. However, the patterns of which species and trees are most vulnerable to drought, and the underlying mechanisms have remained elusive, in part due to the lack of relevant data and difficulty of predicting the location of catastrophic drought years in advance. We used long-term demographic records and extensive databases of functional traits and distribution patterns to understand the responses of 20-53 species to an extreme drought in a seasonally dry tropical forest in Costa Rica, which occurred during the 2015 El Niño Southern Oscillation event. Overall, species-specific mortality rates during the drought ranged from 0% to 34%, and varied little as a function of tree size. By contrast, hydraulic safety margins correlated well with probability of mortality among species, while morphological or leaf economics spectrum traits did not. This firmly suggests hydraulic traits as targets for future research.

Computational wrapping: A universal method to wrap 3D-curved surfaces with nonstretchable materials for conformal devices

This study starts from the counterintuitive question of how we can render conventional stiff, nonstretchable, and even brittle materials sufficiently conformable to fully wrap curved surfaces, such as spheres, without failure. Here, we extend the geometrical design method of computational origami to wrapping. Our computational wrapping approach provides a robust and reliable method for fabricating conformal devices for arbitrary curved surfaces with a computationally designed nonpolyhedral developable net. This computer-aided design transforms two-dimensional (2D)-based materials, such as Si wafers and steel sheets, into various targeted conformal structures that can fully wrap desired 3D structures without fracture or severe plastic deformation. We further demonstrate that our computational wrapping approach enables a design platform that can transform conventional nonstretchable 2D-based devices, such as electroluminescent lighting and flexible batteries, into conformal 3D curved devices.

Harvest models of small populations of a large carnivore using Bayesian forecasting

Harvesting large carnivores can be a management tool for meeting politically set goals for their desired abundance. However, harvesting carnivores creates its own set of conflicts in both society and among conservation professionals, where one consequence is a need to demonstrate that management is sustainable, evidence-based, and guided by science. Furthermore, because large carnivores often also have high degrees of legal protection, harvest quotas have to be carefully justified and constantly adjusted to avoid damaging their conservation status. We developed a Bayesian state-space model to support adaptive management of Eurasian lynx harvesting in Scandinavia. The model uses data from the annual monitoring of lynx abundance and results from long-term field research on lynx biology, which has provided detailed estimates of key demographic parameters. We used the model to predict the probability that the forecasted population size will be below or above the management objectives when subjected to different harvest quotas. The model presented here informs decision makers about the policy risks of alternative harvest levels. Earlier versions of the model have been available for wildlife managers in both Sweden and Norway to guide lynx harvest quotas and the model predictions showed good agreement with observations. We combined monitoring data with data on vital rates and were able to estimate unobserved additional mortality rates, which are most probably due to poaching. In both countries, the past quota setting strategy suggests that there has been a de facto threshold strategy with increasing proportion, which means that there is no harvest below a certain population size, but above this threshold there is an increasing proportion of the population harvested as the population size increases. The annual assessment of the monitoring results, the use of forecasting models, and a threshold harvest approach to quota setting will all reduce the risk of lynx population sizes moving outside the desired goals. The approach we illustrate could be adapted to other populations of mammals worldwide.

Non-Catalytic Subunits Facilitate Quaternary Organization of Plastidic Acetyl-CoA Carboxylase

Arabidopsis (Arabidopsis thaliana), like most dicotyledonous plants, expresses a multicomponent, heteromeric acetyl-CoA carboxylase (htACCase), which catalyzes the generation of the malonyl-CoA precursor of de novo fatty acid biosynthesis. This enzyme consists of four catalytic subunits: biotin carboxylase (BC), carboxyltransferase (CT)-α, CT-β, and biotin carboxyl carrier protein (BCCP1 or BCCP2). By coexpressing combinations of components in a bacterial expression system, we demonstrate noncatalytic BADCs facilitate the assembly and activation of BCCP-BADC-BC subcomplexes catalyzing the bicarbonate-dependent hydrolysis of ATP, which is the first half-reaction catalyzed by the htACCase enzyme. Although BADC proteins do not directly impact formation of the CT-αβ subcomplex, the BADC-facilitated BCCP-BADC-BC subcomplex can more readily interact with the CT-αβ subcomplex to facilitate the generation of malonyl-CoA. The Arabidopsis genome encodes three BADC isoforms (BADC1, BADC2, and BADC3), and BADC2 and BADC3 (rather than BADC1), in combination with BCCP1, best support this quaternary-structural organization and catalytic activation of the htACCase enzyme. Physiological genetic studies validate these attributes as Arabidopsis double mutants singularly expressing BADC2, BADC3, or BADC1 present increasingly greater deleterious impacts on morphological and biochemical phenotypes. Specifically, plants expressing only BADC2 develop normally, plants only expressing BADC3 suffer a stunted root-growth phenotype, and plants expressing only BADC1 are embryo-lethal. The latter phenotype may also be associated with the distinct suborganelle localization of BADC1 in plastids as compared to the localization of the other two BADC homologs. These finding can inspire novel strategies to improve the biological sources of fats and oils for dietary and industrial applications.

Conservation and Divergence in the Meiocyte sRNAomes of Arabidopsis, Soybean, and Cucumber

Meiosis is a critical process for sexual reproduction. During meiosis, genetic information on homologous chromosomes is shuffled through meiotic recombination to produce gametes with novel allelic combinations. Meiosis and recombination are orchestrated by several mechanisms including regulation by small RNAs (sRNAs). Our previous work in Arabidopsis (Arabidopsis thaliana) meiocytes showed that meiocyte-specific sRNAs (ms-sRNAs) have distinct characteristics, including positive association with the coding region of genes that are transcriptionally upregulated during meiosis. Here, we characterized the ms-sRNAs in two important crops, soybean (Glycine max) and cucumber (Cucumis sativus). Ms-sRNAs in soybean have the same features as those in Arabidopsis, suggesting that they may play a conserved role in eudicots. We also investigated the profiles of microRNAs (miRNAs) and phased secondary small interfering RNAs in the meiocytes of all three species. Two conserved miRNAs, miR390 and miR167, are highly abundant in the meiocytes of all three species. In addition, we identified three novel cucumber miRNAs. Intriguingly, our data show that the previously identified phased secondary small interfering RNA pathway involving soybean-specific miR4392 is more abundant in meiocytes. These results showcase the conservation and divergence of ms-sRNAs in flowering plants, and broaden our understanding of sRNA function in crop species.

A Systematic Literature Review of Animal-Assisted Interventions in Oncology (Part II): Theoretical Mechanisms and Frameworks

Animal-assisted interventions (AAIs) can improve patients' quality of life as complementary medical treatments. Part I of this 2-paper systematic review focused on the methods and results of cancer-related AAIs; Part II discusses the theories of the field's investigators. Researchers cite animal personality, physical touch, physical movement, distraction, and increased human interaction as sources of observed positive outcomes. These mechanisms then group under theoretical frameworks such as the social support hypothesis or the human-animal bond concept to fully explain AAI in oncology. The cognitive activation theory of stress, the science of unitary human beings, and the self-object hypothesis are additional frameworks mentioned by some researchers. We also discuss concepts of neurobiological transduction connecting mechanisms to AAI benefits. Future researchers should base study design on theories with testable hypotheses and use consistent terminology to report results. This review aids progress toward a unified theoretical framework and toward more holistic cancer treatments.

Controlling for baseline telomere length biases estimates of the rate of telomere attrition

Longitudinal studies have sought to establish whether environmental exposures such as smoking accelerate the attrition of individuals' telomeres over time. These studies typically control for baseline telomere length (TL) by including it as a covariate in statistical models. However, baseline TL also differs between smokers and non-smokers, and telomere attrition is spuriously linked to baseline TL via measurement error and regression to the mean. Using simulated datasets, we show that controlling for baseline TL overestimates the true effect of smoking on telomere attrition. This bias increases with increasing telomere measurement error and increasing difference in baseline TL between smokers and non-smokers. Using a meta-analysis of longitudinal datasets, we show that as predicted, the estimated difference in telomere attrition between smokers and non-smokers is greater when statistical models control for baseline TL than when they do not, and the size of the discrepancy is positively correlated with measurement error. The bias we describe is not specific to smoking and also applies to other exposures. We conclude that to avoid invalid inference, models of telomere attrition should not control for baseline TL by including it as a covariate. Many claims of accelerated telomere attrition in individuals exposed to adversity need to be re-assessed.

The history and impact of digitization and digital data mobilization on biodiversity research

The first two decades of the twenty-first century have seen a rapid rise in the mobilization of digital biodiversity data. This has thrust natural history museums into the forefront of biodiversity research, underscoring their central role in the modern scientific enterprise. The advent of mobilization initiatives such as the United States National Science Foundation's Advancing Digitization of Biodiversity Collections (ADBC), Australia's Atlas of Living Australia (ALA), Mexico's National Commission for the Knowledge and Use of Biodiversity (CONABIO), Brazil's Centro de Referência em Informação (CRIA) and China's National Specimen Information Infrastructure (NSII) has led to a rapid rise in data aggregators and an exponential increase in digital data for scientific research and arguably provide the best evidence of where species live. The international Global Biodiversity Information Facility (GBIF) now serves about 131 million museum specimen records, and Integrated Digitized Biocollections (iDigBio) in the USA has amassed more than 115 million. These resources expose collections to a wider audience of researchers, provide the best biodiversity data in the modern era outside of nature itself and ensure the primacy of specimen-based research. Here, we provide a brief history of worldwide data mobilization, their impact on biodiversity research, challenges for ensuring data quality, their contribution to scientific publications and evidence of the rising profiles of natural history collections.This article is part of the theme issue 'Biological collections for understanding biodiversity in the Anthropocene'.

Allosteric interactions in a birod model of DNA

Allosteric interactions between molecules bound to DNA at distant locations have been known for a long time. The phenomenon has been studied via experiments and numerical simulations, but a comprehensive understanding grounded in a theory of DNA elasticity remains a challenge. Here, we quantify allosteric interactions between two entities bound to DNA by using the theory of birods. We recognize that molecules bound to DNA cause local deformations that can be captured in a birod model which consists of two elastic strands interacting via an elastic web representing the basepairs. We show that the displacement field caused by bound entities decays exponentially with distance from the binding site. We compute the interaction energy between two proteins on DNA as a function of distance between them and find that it decays exponentially while oscillating with the periodicity of the double helix, in excellent agreement with experiments. The decay length of the interaction energy can be determined in terms of the mechanical properties of the strands and the webbing in our birod model, and it varies with the GC content of the DNA. Our model provides a framework for viewing allosteric interactions in DNA within the ambit of configurational forces of continuum elasticity.

An integrated RNAseq-1H NMR metabolomics approach to understand soybean primary metabolism regulation in response to Rhizoctonia foliar blight disease

BACKGROUND

Rhizoctonia solani AG1-IA is a devastating phytopathogen causing Rhizoctonia foliar blight (RFB) of soybean worldwide with yield losses reaching 60%. Plant defense mechanisms are complex and information from different metabolic pathways is required to thoroughly understand plant defense regulation and function. Combining information from different "omics" levels such as transcriptomics, metabolomics, and proteomics is required to gain insights into plant metabolism and its regulation. As such, we studied fluctuations in soybean metabolism in response to R. solani infection at early and late disease stages using an integrated transcriptomics-metabolomics approach, focusing on the regulation of soybean primary metabolism and oxidative stress tolerance.

RESULTS

Transcriptomics (RNAseq) and metabolomics (1H NMR) data were analyzed individually and by integration using bidirectional orthogonal projections to latent structures (O2PLS) to reveal possible links between the metabolome and transcriptome during early and late infection stages. O2PLS analysis detected 516 significant transcripts, double that reported in the univariate analysis, and more significant metabolites than detected in partial least squares discriminant analysis. Strong separation of treatments based on integration of the metabolomes and transcriptomes of the analyzed soybean leaves was revealed, similar trends as those seen in analyses done on individual datasets, validating the integration method being applied. Strong fluctuations of soybean primary metabolism occurred in glycolysis, the TCA cycle, photosynthesis and photosynthates in response to R. solani infection. Data were validated using quantitative real-time PCR on a set of specific markers as well as randomly selected genes. Significant increases in transcript and metabolite levels involved in redox reactions and ROS signaling, such as peroxidases, thiamine, tocopherol, proline, L-alanine and GABA were also recorded. Levels of ethanol increased 24 h post-infection in soybean leaves, and alcohol dehydrogenase (ADH) loss-of-function mutants of Arabidopsis thaliana had higher necrosis than wild type plants.

CONCLUSIONS

As a proof-of-concept, this study offers novel insights into the biological correlations and identification of candidate genes and metabolites that can be used in soybean breeding for resistance to R. solani AG1-IA infection. Additionally, these findings imply that alcohol and its associated gene product ADH may have important roles in plant resistance to R. solani AG1-IA causing foliar blight.