Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

Biomechanical Phenotyping of Chronic Low Back Pain: Protocol for BACPAC

OBJECTIVE

Biomechanics represents the common final output through which all biopsychosocial constructs of back pain must pass, making it a rich target for phenotyping. To exploit this feature, several sites within the NIH Back Pain Consortium (BACPAC) have developed biomechanics measurement and phenotyping tools. The overall aims of this article were to: 1) provide a narrative review of biomechanics as a phenotyping tool; 2) describe the diverse array of tools and outcome measures that exist within BACPAC; and 3) highlight how leveraging these technologies with the other data collected within BACPAC could elucidate the relationship between biomechanics and other metrics used to characterize low back pain (LBP).

METHODS

The narrative review highlights how biomechanical outcomes can discriminate between those with and without LBP, as well as among levels of severity of LBP. It also addresses how biomechanical outcomes track with functional improvements in LBP. Additionally, we present the clinical use case for biomechanical outcome measures that can be met via emerging technologies.

RESULTS

To answer the need for measuring biomechanical performance, our "Results" section describes the spectrum of technologies that have been developed and are being used within BACPAC.

CONCLUSION AND FUTURE DIRECTIONS

The outcome measures collected by these technologies will be an integral part of longitudinal and cross-sectional studies conducted in BACPAC. Linking these measures with other biopsychosocial data collected within BACPAC increases our potential to use biomechanics as a tool for understanding the mechanisms of LBP, phenotyping unique LBP subgroups, and matching these individuals with an appropriate treatment paradigm.

Deformation Rate-Adaptive Conducting Polymers and Composites

Materials are more easily damaged during accidents that involve rapid deformation. Here, a design strategy is described for electronic materials comprised of conducting polymers that defies this orthodox property, making their extensibility and toughness dynamically adaptive to deformation rates. This counterintuitive property is achieved through a morphology of interconnected nanoscopic core-shell micelles, where the chemical interactions are stronger within the shells than the cores. As a result, the interlinked shells retain material integrity under strain, while the rate of dissociation of the cores controls the extent of micelle elongation, which is a process that adapts to deformation rates. A prototype based on polyaniline shows a 7.5-fold increase in ultimate elongation and a 163-fold increase in toughness when deformed at increasing rates from 2.5 to 10 000% min-1 . This concept can be generalized to other conducting polymers and highly conductive composites to create "self-protective" soft electronic materials with enhanced durability under dynamic movement or deformation.

Dual chemical labeling enables nucleotide-resolution mapping of DNA abasic sites and common alkylation damage in human mitochondrial DNA

Mitochondrial DNA (mtDNA) modifications play an emerging role in innate immunity and inflammatory diseases. Nonetheless, relatively little is known regarding the locations of mtDNA modifications. Such information is critically important for deciphering their roles in mtDNA instability, mtDNA-mediated immune and inflammatory responses, and mitochondrial disorders. The affinity probe-based enrichment of lesion-containing DNA represents a key strategy for sequencing DNA modifications. Existing methods are limited in the enrichment specificity of abasic (AP) sites, a prevalent DNA modification and repair intermediate. Herein, we devise a novel approach, termed dual chemical labeling-assisted sequencing (DCL-seq), for mapping AP sites. DCL-seq features two designer compounds for enriching and mapping AP sites specifically at single-nucleotide resolution. For proof of principle, we mapped AP sites in mtDNA from HeLa cells under different biological conditions. The resulting AP site maps coincide with mtDNA regions with low TFAM (mitochondrial transcription factor A) coverage and with potential G-quadruplex-forming sequences. In addition, we demonstrated the broader applicability of the method in sequencing other DNA modifications in mtDNA, such as N7-methyl-2'-deoxyguanosine and N3-methyl-2'-deoxyadenosine, when coupled with a lesion-specific repair enzyme. Together, DCL-seq holds the promise to sequence multiple DNA modifications in various biological samples.

Optogenetic control of Wnt signaling models cell-intrinsic embryogenic patterning using 2D human pluripotent stem cell culture

In embryonic stem cell (ESC) models for early development, spatially and temporally varying patterns of signaling and cell types emerge spontaneously. However, mechanistic insight into this dynamic self-organization is limited by a lack of methods for spatiotemporal control of signaling, and the relevance of signal dynamics and cell-to-cell variability to pattern emergence remains unknown. Here, we combine optogenetic stimulation, imaging and transcriptomic approaches to study self-organization of human ESCs (hESC) in two-dimensional (2D) culture. Morphogen dynamics were controlled via optogenetic activation of canonical Wnt/β-catenin signaling (optoWnt), which drove broad transcriptional changes and mesendoderm differentiation at high efficiency (>99% cells). When activated within cell subpopulations, optoWnt induced cell self-organization into distinct epithelial and mesenchymal domains, mediated by changes in cell migration, an epithelial to mesenchymal-like transition and TGFβ signaling. Furthermore, we demonstrate that such optogenetic control of cell subpopulations can be used to uncover signaling feedback mechanisms between neighboring cell types. These findings reveal that cell-to-cell variability in Wnt signaling is sufficient to generate tissue-scale patterning and establish a hESC model system for investigating feedback mechanisms relevant to early human embryogenesis.

Association of nativity with survival among adults with hepatocellular carcinoma

INTRODUCTION

Immigrants comprise a considerable proportion of those diagnosed with hepatocellular carcinoma (HCC) in the United States. Nativity or birthplace affects incidence and risk factors for HCC, but little is known about its influence on survival after diagnosis.

METHODS

We identified 51 533 adults with HCC with available birthplace in the California Cancer Registry between 1988 and 2017. HCC cases were categorized as foreign born or US born and stratified by mutually exclusive race and ethnicity groups. Primary outcome was all-cause mortality. Race and ethnicity-specific Cox regression propensity score-weighted models evaluated the relationship between nativity and death as well as region of birth among foreign-born patients.

RESULTS

A total of 40% of all HCC cases were foreign born, and 92.2%, 45.2%, 9.1%, and 5.8% of Asian/Pacific Islander (API), Hispanic, White, and Black patients were foreign born, respectively. Five-year survival rates were higher in foreign-born patients compared with US-born patients: 12.9% vs 9.6% for White patients, 11.7% vs 9.8% for Hispanic patients, 12.8% vs 8.1% for Black patients, and 16.4% vs 12.4% for API patients. Nativity was associated with survival, with better survival in foreign-born patients: White patients: hazard ratio (HR) = 0.86 (95% confidence interval [CI] = 0.81 to 0.90), Hispanic patients: HR = 0.90 (95% CI = 0.86 to 0.93), Black patients: HR = 0.89 (95% CI = 0.76 to 1.05), and API patients: HR = 0.94 (95% CI = 0.88 to 1.00). Among foreign-born patients, lower mortality was observed in those from Central and South America compared with Mexico for Hispanic patients, East Asia compared with Southeast Asia for API patients, and East Europe and Greater Middle East compared with West/South/North Europe for White patients.

CONCLUSION

Foreign-born patients with HCC have better survival than US-born patients. Further investigation into the mechanisms of this survival disparity by nativity is needed.

Saccharomyces cerevisiae DNA polymerase IV overcomes Rad51 inhibition of DNA polymerase δ in Rad52-mediated direct-repeat recombination

Saccharomyces cerevisiae DNA polymerase IV (Pol4) like its homolog, human DNA polymerase lambda (Polλ), is involved in Non-Homologous End-Joining and Microhomology-Mediated Repair. Using genetic analysis, we identified an additional role of Pol4 also in homology-directed DNA repair, specifically in Rad52-dependent/Rad51-independent direct-repeat recombination. Our results reveal that the requirement for Pol4 in repeat recombination was suppressed by the absence of Rad51, suggesting that Pol4 counteracts the Rad51 inhibition of Rad52-mediated repeat recombination events. Using purified proteins and model substrates, we reconstituted in vitro reactions emulating DNA synthesis during direct-repeat recombination and show that Rad51 directly inhibits Polδ DNA synthesis. Interestingly, although Pol4 was not capable of performing extensive DNA synthesis by itself, it aided Polδ in overcoming the DNA synthesis inhibition by Rad51. In addition, Pol4 dependency and stimulation of Polδ DNA synthesis in the presence of Rad51 occurred in reactions containing Rad52 and RPA where DNA strand-annealing was necessary. Mechanistically, yeast Pol4 displaces Rad51 from ssDNA independent of DNA synthesis. Together our in vitro and in vivo data suggest that Rad51 suppresses Rad52-dependent/Rad51-independent direct-repeat recombination by binding to the primer-template and that Rad51 removal by Pol4 is critical for strand-annealing dependent DNA synthesis.

Reference genome of an iconic lizard in western North America, Blainville's horned lizard Phrynosoma blainvillii

Genome assemblies are increasingly being used to identify adaptive genetic variation that can help prioritize the population management of protected species. This approach may be particularly relevant to species like Blainville's horned lizard, Phrynosoma blainvillii, due to its specialized diet on noxious harvester ants, numerous adaptative traits for avoiding predation (e.g. cranial horns, dorsoventrally compressed body, cryptic coloration, and blood squirting from the orbital sinuses), and status as Species of Special Concern in California. Rangewide decline since the early 20th century, the basis of its conservation status, has been driven mainly by habitat conversion, over-collecting, and invasion of a non-native ant that displaces its native ant prey base. Here, we report on a scaffold-level genome assembly for P. blainvillii as part of the California Conservation Genomics Project (CCGP), produced using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology. The de novo assembly has 78 scaffolds, a total length of ~2.21 Gb, a scaffold N50 length of ~352 Mb, and BUSCO score of 97.4%. This is the second species of Phrynosoma for which a reference genome has been assembled and represents a considerable improvement in terms of contiguity and completeness. Combined with the landscape genomics data being compiled by the CCGP, this assembly will help strategize efforts to maintain and/or restore local genetic diversity, where interventions like genetic rescue, translocation, and strategic land preservation may be the only means by which P. blainvillii and other low-vagility species can survive in the fragmented habitats of California.

Reference genome of the Virginia rail, Rallus limicola

The Virginia rail, Rallus limicola, is a member of the family Rallidae, which also includes many other species of secretive and poorly studied wetland birds. It is recognized as a single species throughout its broad distribution in North America where it is exploited as a game bird, often with generous harvest limits, despite a lack of systematic population surveys and evidence of declines in many areas due to wetland loss and degradation. To help advance understanding of the phylogeography, biology, and ecology of this elusive species, we report the first reference genome assembly for the Virginia rail, produced as part of the California Conservation Genomics Project (CCGP). We produced a de novo genome assembly using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology with an estimated sequencing error rate of 0.191%. The assembly consists of 1,102 scaffolds spanning 1.39 Gb, with a contig N50 of 11.0 Mb, scaffold N50 of 25.3 Mb, largest contig of 45 Mb, and largest scaffold of 128.4 Mb. It has a high BUSCO completeness score of 96.9% and represents the first genome assembly available for the genus Rallus. This genome assembly will help resolve questions about the complex evolutionary history of rails and evaluate the potential of rails for adaptive evolution in the face of growing threats from climate change and habitat loss and fragmentation. It will also provide a valuable resource for rail conservation efforts by quantifying Virginia rail vagility, population connectivity, and effective population sizes.

Reference genome of the black rail, Laterallus jamaicensis

The black rail, Laterallus jamaicensis, is one of the most secretive and poorly understood birds in the Americas. Two of its five subspecies breed in North America: the Eastern black rail (L. j. jamaicensis), found primarily in the southern and mid-Atlantic states, and the California black rail (L. j. coturniculus), inhabiting California and Arizona, are recognized across the highly disjunct distribution. Population declines, due primarily to wetland loss and degradation, have resulted in conservation status listings for both subspecies. To help advance understanding of the phylogeography, biology, and ecology of this elusive species, we report the first reference genome assembly for the black rail, produced as part of the California Conservation Genomics Project (CCGP). We produced a de novo genome assembly using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology with an estimated sequencing error rate of 0.182%. The assembly consists of 964 scaffolds spanning 1.39 Gb, with a contig N50 of 7.4 Mb, scaffold N50 of 21.4 Mb, largest contig of 44.8 Mb, and largest scaffold of 101.2 Mb. The assembly has a high BUSCO completeness score of 96.8% and represents the first genome assembly available for the genus Laterallus. This genome assembly can help resolve questions about the complex evolutionary history of rails, assess black rail vagility and population connectivity, estimate effective population sizes, and evaluate the potential of rails for adaptive evolution in the face of growing threats from climate change, habitat loss and fragmentation, and disease.

A highly contiguous genome assembly for the California quail (Callipepla californica)

The California quail (Callipepla californica) is an iconic native bird of scrub and oak woodlands in California and the Baja Peninsula of Mexico. Here, we report a draft reference assembly for the species generated from PacBio HiFi long read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 321 scaffolds totaling 1.08 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 5.5 Mb, scaffold N50 of 19.4 Mb, and BUSCO completeness score of 96.5%. Transposable elements (TEs) occupy 16.5% of the genome, more than previous Odontophoridae quail assemblies but in line with estimates of TE content for recent long-read assemblies of chicken and Peking duck. Together these metrics indicate that the present assembly is more complete than prior reference assemblies generated for Odontophoridae quail. This reference will serve as an essential resource for studies on local adaptation, phylogeography, and conservation genetics in this species of significant biological and recreational interest.

Reference genome for the American rubyspot damselfly, Hetaerina americana

Damselflies and dragonflies (Order: Odonata) play important roles in both aquatic and terrestrial food webs and can serve as sentinels of ecosystem health and predictors of population trends in other taxa. The habitat requirements and limited dispersal of lotic damselflies make them especially sensitive to habitat loss and fragmentation. As such, landscape genomic studies of these taxa can help focus conservation efforts on watersheds with high levels of genetic diversity, local adaptation, and even cryptic endemism. Here, as part of the California Conservation Genomics Project (CCGP), we report the first reference genome for the American rubyspot damselfly, Hetaerina americana, a species associated with springs, streams and rivers throughout California. Following the CCGP assembly pipeline, we produced two de novo genome assemblies. The primary assembly includes 1,630,044,487 base pairs, with a contig N50 of 5.4 Mb, a scaffold N50 of 86.2 Mb, and a BUSCO completeness score of 97.6%. This is the seventh Odonata genome to be made publicly available and the first for the subfamily Hetaerininae. This reference genome fills an important phylogenetic gap in our understanding of Odonata genome evolution, and provides a genomic resource for a host of interesting ecological, evolutionary, and conservation questions for which the rubyspot damselfly genus Hetaerina is an important model system.

Reference genome of the long-jawed orb-weaver, Tetragnatha versicolor (Araneae: Tetragnathidae)

Climate-driven changes in hydrological regimes are of global importance and are particularly significant in riparian ecosystems. Riparian ecosystems in California provide refuge to many native and vulnerable species within a xeric landscape. California Tetragnatha spiders play a key role in riparian ecosystems, serving as a link between terrestrial and aquatic elements. Their tight reliance on water paired with the widespread distributions of many species make them ideal candidates to better understand the relative role of waterways versus geographic distance in shaping the population structure of riparian species. To assist in better understanding population structure, we constructed a reference genome assembly for Tetragnatha versicolor using long-read sequencing, scaffolded with proximity ligation Omni-C data. The near-chromosome-level assembly is comprised of 174 scaffolds spanning 1.06 Gb pairs, with a scaffold N50 of 64.1 Mb pairs and BUSCO completeness of 97.6%. This reference genome will facilitate future study of T. versicolor population structure associated with the rapidly changing environment of California.

Predicting the shape, size, and placement of adult human pubic symphyses

OBJECTIVES

When reconstructing fossil pelves, the articulation of the pelvic bones largely relies on subjective decisions by researchers. Different positionings at the pubic symphysis can affect the overall morphology of the pelvis and the subsequent biological interpretation associated with that individual or species. This study aims to reduce this subjectivity using quantitative models to predict pubic symphysis morphology.

METHODS

We collected 3D landmarks and semilandmarks on the pubic symphysis and adjacent aspects on the CT scans of 103 adults. Using geometric morphometrics we, (1) quantified pubic symphysis morphology, (2) trained simple and two-stage least-squares linear regression models to predict pubic symphysis shape, and (3) assessed the shape variation in the sample. The model with the lowest prediction error was identified as the best model. Principal components analysis was used to explore the effects of each variable on shape and hypothetical shapes were generated from the model to illustrate these effects.

RESULTS

The best model is a two-stage least-squares model that predicts pubic symphysis size at the first stage using additive effects of sex and age, then subsequently interacts pubic symphysis size with sex and age at the second stage to predict pubic symphysis shape. Other models with low prediction errors included variables reflecting pelvic size and breadth.

CONCLUSION

Linear regression modeling can be used to systematically predict pubic symphysis morphology. This method can be used in addition to other techniques to improve fossil reconstructions by more accurately estimating the morphology of this region of the pelvis.

Exposure to quasi-ultrafine particulate matter accelerates memory impairment and Alzheimer's disease-like neuropathology in the AppNL-G-F knock-in mouse model

Exposure to traffic-related air pollution consisting of particulate matter (PM) is associated with cognitive decline leading to Alzheimer's disease (AD). In this study, we sought to examine the neurotoxic effects of exposure to ultrafine PM and how it exacerbates neuronal loss and AD-like neuropathology in wildtype (WT) mice and a knock-in mouse model of AD (AppNL-G-F/+-KI) when the exposure occurs at a prepathologic stage or at a later age with the presence of neuropathology. AppNL-G-F/+-KI and WT mice were exposed to concentrated ultrafine PM from local ambient air in Irvine, California, for 12 weeks, starting at 3 or 9 months of age. Particulate matter-exposed animals received concentrated ultrafine PM up to 8 times above the ambient levels, whereas control animals were exposed to purified air. Particulate matter exposure resulted in a marked impairment of memory tasks in prepathologic AppNL-G-F/+-KI mice without measurable changes in amyloid-β pathology, synaptic degeneration, and neuroinflammation. At aged, both WT and AppNL-G-F/+-KI mice exposed to PM showed a significant memory impairment along with neuronal loss. In AppNL-G-F/+-KI mice, we also detected an increased amyloid-β buildup and potentially harmful glial activation including ferritin-positive microglia and C3-positive astrocytes. Such glial activation could promote the cascade of degenerative consequences in the brain. Our results suggest that exposure to PM impairs cognitive function at both ages while exacerbation of AD-related pathology and neuronal loss may depend on the stage of pathology, aging, and/or state of glial activation. Further studies will be required to unveil the neurotoxic role of glial activation activated by PM exposure.

Distinct neural activations correlate with maximization of reward magnitude versus frequency

Choice selection strategies and decision-making are typically investigated using multiple-choice gambling paradigms that require participants to maximize expected value of rewards. However, research shows that performance in such paradigms suffers from individual biases towards the frequency of gains such that users often choose smaller frequent gains over larger rarely occurring gains, also referred to as melioration. To understand the basis of this subjective tradeoff, we used a simple 2-choice reward task paradigm in 186 healthy human adult subjects sampled across the adult lifespan. Cortical source reconstruction of simultaneously recorded electroencephalography suggested distinct neural correlates for maximizing reward magnitude versus frequency. We found that activations in the parahippocampal and entorhinal areas, which are typically linked to memory function, specifically correlated with maximization of reward magnitude. In contrast, maximization of reward frequency was correlated with activations in the lateral orbitofrontal cortices and operculum, typical areas involved in reward processing. These findings reveal distinct neural processes serving reward frequency versus magnitude maximization that can have clinical translational utility to optimize decision-making.

The non-canonical Wnt receptor Ror2 is required for cartilage cell polarity and morphogenesis of the craniofacial skeleton in zebrafish

Non-canonical/β-catenin-independent Wnt signaling plays crucial roles in tissue/cell polarity in epithelia, but its functions have been less well studied in mesenchymal tissues, such as the skeleton. Mutations in non-canonical Wnt signaling pathway genes cause human skeletal diseases such as Robinow syndrome and Brachydactyly Type B1, which disrupt bone growth throughout the endochondral skeleton. Ror2 is one of several non-canonical Wnt receptor/co-receptors. Here, we show that ror2-/- mutant zebrafish have craniofacial skeletal defects, including disruptions of chondrocyte polarity. ror1-/- mutants appear to be phenotypically wild type, but loss of both ror1 and ror2 leads to more severe cartilage defects, indicating partial redundancy. Skeletal defects in ror1/2 double mutants resemble those of wnt5b-/- mutants, suggesting that Wnt5b is the primary Ror ligand in zebrafish. Surprisingly, the proline-rich domain of Ror2, but not its kinase domain, is required to rescue its function in mosaic transgenic experiments in ror2-/- mutants. These results suggest that endochondral bone defects in ROR-related human syndromes reflect defects in cartilage polarity and morphogenesis.

Disparities in cancer incidence by rurality in California

BACKGROUND

Cancer rates in rural areas across the United States have different patterns than in urban areas. This study examines associations between rurality and incidence for the top 5 cancers in California and evaluates whether these associations vary jointly by sex, race, and ethnicity.

METHODS

We used 2015-2019 California Cancer Registry data to compare incidence rate ratios (IRRs) and trends for breast, prostate, lung, colorectal, and skin (melanoma) cancers. We leveraged census tract aggregation zones and 7 levels of percentage rural population (0%, >0% to <10%, 10% to <20%, 20% to <30%, 30% to <40%, 40% to <50%, and 50+%).

RESULTS

Zones with higher proportions of rural population were significantly associated with lower incidence of female breast cancer and prostate cancer, though the trends were not statistically significant overall. Zones with higher proportions of rural population were significantly associated with higher incidence of lung cancer and melanoma. There were no statistically significant trends for colorectal cancer overall. Comparing areas with 50% and over rural population with areas with 0% rural population, the IRR for lung cancer in Hispanic females was higher (IRR = 1.43, 95% confidence interval [CI] = 1.17 to 1.74) than in Hispanic males (IRR = 0.90, 95% CI = 0.72 to 1.11). Also, in areas with 50% or more rural population, the IRR for melanoma was higher in Hispanic females (IRR = 1.75, 95% CI = 1.23 to 2.45) than non-Hispanic White females (IRR = 0.87, 95% CI = 0.80 to 0.95).

CONCLUSIONS

Our findings show that rurality is associated with cancer incidence and underscore the importance of jointly examining rural disparities with sex, race, and ethnicity by cancer site.

Reference genome of the Woolly Sculpin, Clinocottus analis

Sculpins (Family Cottidae) are generally cold-temperate intertidal reef fishes most commonly found in the North Pacific. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Woolly Sculpin, Clinocottus analis, to establish a genomic model for understanding phylogeographic structure of inshore marine taxa along the California coast. These patterns, in turn, should further inform the design of marine protected areas using dispersal models based on genomic data. The small genome of C. analis is typical of marine fishes at less than 1 Gb (genome size = 538 Mb), and our assembly is near-chromosome level (contig N50 = 9.1 Mb, scaffold N50 = 21 Mb, BUSCO completeness = 97.9%). Within the context of the CCGP, the Woolly Sculpin genome will be used as a reference for future whole-genome resequencing projects aimed at enhancing our knowledge of the population structure of the species, and efficacy of marine protected areas across the state.

A reference genome assembly of the declining tricolored blackbird, Agelaius tricolor

The tricolored blackbird, Agelaius tricolor, is a gregarious species that forms enormous breeding and foraging colonies in wetland and agricultural habitats, primarily in California, USA. Once extremely abundant, species numbers have declined dramatically in the past century, largely due to losses of breeding and foraging habitats. Tricolored blackbirds are currently listed as Endangered by the IUCN, and Threatened under the California Endangered Species Act. Increased genetic information is needed to detail the evolutionary consequences of a species-wide bottleneck and inform conservation management. Here, we present a contiguous tricolored blackbird reference genome, assembled with PacBio HiFi long reads and Dovetail Omni-C data to generate a scaffold-level assembly containing multiple chromosome-length scaffolds. This genome adds a valuable resource for important evolutionary and conservation research on tricolored blackbirds and related species.

Reference genome of the Monkeyface Prickleback, Cebidichthys violaceus

Pricklebacks (Family Stichaeidae) are generally cold-temperate fishes most commonly found in the north Pacific. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Monkeyface Prickleback, Cebidichthys violaceus, to establish a genomic model for understanding phylogeographic patterns of marine organisms in California. These patterns, in turn, may inform the design of marine protected areas using dispersal models based on forthcoming population genomic data. The genome of C. violaceus is typical of many marine fishes at less than 1 Gb (genome size = 575.6 Mb), and our assembly is near-chromosome level (contig N50 = 1 Mb, scaffold N50 = 16.4 Mb, BUSCO completeness = 93.2%). Within the context of the CCGP, the genome will be used as a reference for future whole genome resequencing projects, enhancing our knowledge of the population structure of the species and more generally, the efficacy of marine protected areas as a primary conservation tool across California's marine ecosystems.

A chromosome-level reference genome for the Versatile Fairy Shrimp, Branchinecta lindahli

We present the novel reference genome of the Versatile Fairy Shrimp, Branchinecta lindahli. The Versatile Fairy Shrimp is a freshwater anostracan crustacean found across the western United States from Iowa to Oregon and from Alberta to Baja California. It is an ephemeral pool specialist, living in prairie potholes, irrigation ditches, tire treads, vernal pools, and other temporary freshwater wetlands. Anostracan fairy shrimp are facing global declines with 3 species in California on the Endangered Species list. This species was included in the California Conservation Genomics Project to provide an easily accessible reference genome, and to provide whole-genome resources for a generalist species, which may lead to new insights into Anostracan resiliency in the face of climate change. The final gapped genome comprises 15 chromosome-length scaffolds covering 98.63% of the 384.8 Mb sequence length, and an additional 55 unscaffolded contigs.

A draft reference genome of the Vernal Pool Fairy Shrimp, Branchinecta lynchi

We present the reference genome of the Vernal Pool Fairy Shrimp Branchinecta lynchi. This branchiopod crustacean is endemic to California's freshwater ephemeral ponds. It faces enormous habitat loss and fragmentation as urbanization and agriculture have fundamentally changed the vernal pool landscape over the past 3 centuries. The assembled genome consists of 22 chromosome-length scaffolds that account for 96.85% of the total sequence. One hundred and ninety-five unscaffolded contigs comprise the rest of the genome's 575.6 Mb length. The genome is substantially complete with a BUSCO score of 90.0%. There is no immediately identifiable sex chromosome, typical for this class of organism. This new resource will permit researchers to better understand the adaptive capacity of this imperiled species, as well as answer lingering questions about anostracan physiology, sex determination, and development.

Structural snapshots of base excision by the cancer-associated variant MutY N146S reveal a retaining mechanism

DNA glycosylase MutY plays a critical role in suppression of mutations resulted from oxidative damage, as highlighted by cancer-association of the human enzyme. MutY requires a highly conserved catalytic Asp residue for excision of adenines misinserted opposite 8-oxo-7,8-dihydroguanine (OG). A nearby Asn residue hydrogen bonds to the catalytic Asp in structures of MutY and its mutation to Ser is an inherited variant in human MUTYH associated with colorectal cancer. We captured structural snapshots of N146S Geobacillus stearothermophilus MutY bound to DNA containing a substrate, a transition state analog and enzyme-catalyzed abasic site products to provide insight into the base excision mechanism of MutY and the role of Asn. Surprisingly, despite the ability of N146S to excise adenine and purine (P) in vitro, albeit at slow rates, N146S-OG:P complex showed a calcium coordinated to the purine base altering its conformation to inhibit hydrolysis. We obtained crystal structures of N146S Gs MutY bound to its abasic site product by removing the calcium from crystals of N146S-OG:P complex to initiate catalysis in crystallo or by crystallization in the absence of calcium. The product structures of N146S feature enzyme-generated β-anomer abasic sites that support a retaining mechanism for MutY-catalyzed base excision.

Shared Decision-Making Concerning Anal Cancer Screening in Persons With Human Immunodeficiency Virus

BACKGROUND

Anal high-grade squamous intraepithelial lesion (aHSIL) is the immediate precursor of anal cancer. Anal cytology is a recommended screening test to identify aHSIL among people with human immunodeficiency virus (HIV; PWH). Heterogeneity of risk for invasive anal cancer among PWH suggests the value of a shared decision-making framework regarding screening.

METHODS

Using a longitudinal HIV cohort with a comprehensive anal cancer screening program, we estimated the adjusted probabilities of having aHSIL on the first anal cytology. We used logistic regression models with inverse probability weighting to account for differential screening in the cohort and to construct a predicted probability nomogram for aHSIL. Sensitivity analysis was performed to estimate aHSIL prevalence corrected for misclassification bias.

RESULTS

Of 8139 PWH under care between 2007 and 2020, 4105 (49.8%) underwent at least 1 anal cytology test. First-time cytology aHSIL was present in 502 (12.2%) PWH. The adjusted probability of having aHSIL varied from 5% to 18% depending on patient characteristics. Prespecified factors in the aHSIL prediction model included nadir CD4 cell count, ethnicity, race, age, sex, gender identity, and HIV risk factors. The ability of the model to discriminate cytological aHSIL was modest, with an area under the curve of 0.63 (95% confidence interval, .60-.65).

CONCLUSIONS

PWH are at increased risk for aHSIL and invasive anal cancer. Risk, however, varies by patient characteristics. Individual risk factor profiles predictive of aHSIL can be modeled and operationalized as nomograms to facilitate shared decision-making conversations concerning anal cancer screening.

Single cell evaluation of endocardial Hand2 gene regulatory networks reveals HAND2-dependent pathways that impact cardiac morphogenesis

The transcription factor HAND2 plays essential roles during cardiogenesis. Hand2 endocardial deletion (H2CKO) results in tricuspid atresia or double inlet left ventricle with accompanying intraventricular septum defects, hypo-trabeculated ventricles and an increased density of coronary lumens. To understand the regulatory mechanisms of these phenotypes, single cell transcriptome analysis of mouse E11.5 H2CKO hearts was performed revealing a number of disrupted endocardial regulatory pathways. Using HAND2 DNA occupancy data, we identify several HAND2-dependent enhancers, including two endothelial enhancers for the shear-stress master regulator KLF2. A 1.8 kb enhancer located 50 kb upstream of the Klf2 TSS imparts specific endothelial/endocardial expression within the vasculature and endocardium. This enhancer is HAND2-dependent for ventricular endocardium expression but HAND2-independent for Klf2 vascular and valve expression. Deletion of this Klf2 enhancer results in reduced Klf2 expression within ventricular endocardium. These data reveal that HAND2 functions within endocardial gene regulatory networks including shear-stress response.

Polygenic risk of any, metastatic, and fatal prostate cancer in the Million Veteran Program

BACKGROUND

Genetic scores may provide an objective measure of prostate cancer risk and thus inform screening decisions. We evaluated whether a polygenic hazard score based on 290 genetic variants (PHS290) is associated with prostate cancer risk in a diverse population, including Black men, who have higher average risk of prostate cancer death but are often treated as a homogeneously high-risk group.

METHODS

This was a retrospective analysis of the Million Veteran Program, a national, population-based cohort study of US military veterans conducted 2011-2021. Cox proportional hazards analyses tested for association of genetic and other risk factors (including self-reported race and ethnicity and family history) with age at death from prostate cancer, age at diagnosis of metastatic (nodal or distant) prostate cancer, and age at diagnosis of any prostate cancer.

RESULTS

A total of 590 750 male participants were included. Median age at last follow-up was 69 years. PHS290 was associated with fatal prostate cancer in the full cohort and for each racial and ethnic group (P < .001). Comparing men in the highest 20% of PHS290 with those in the lowest 20% (based on percentiles from an independent training cohort), the hazard ratio for fatal prostate cancer was 4.42 (95% confidence interval = 3.91 to 5.02). When accounting for guideline-recommended risk factors (family history, race, and ethnicity), PHS290 remained a strong independent predictor of any, metastatic, and fatal prostate cancer.

CONCLUSIONS

PHS290 stratified US veterans of diverse ancestry for lifetime risk of prostate cancer, including metastatic and fatal cancer. Predicting genetic risk of lethal prostate cancer with PHS290 might inform individualized decisions about prostate cancer screening.

External Factors Modulating Vaping-Induced Thermal Degradation of Vitamin E Acetate

Despite previous studies indicating the thermal stability of vitamin E acetate (VEA) at low temperatures, VEA has been shown to readily decompose into various degradation products such as alkenes, long-chain alcohols, and carbonyls such as duroquinone (DQ) at vaping temperatures of <200 °C. While most models simulate the thermal decomposition of e-liquids under pyrolysis conditions, numerous factors, including vaping behavior, device construction, and the surrounding environment, may impact the thermal degradation process. In this study, we investigated the role of the presence of molecular oxygen (O2) and transition metals in promoting thermal oxidation of e-liquids, resulting in greater degradation than predicted by pure pyrolysis. Thermal degradation of VEA was performed in inert (N2) and oxidizing atmospheres (clean air) in the absence and presence of Ni-Cr and Cu-Ni alloy nanopowders, metals commonly found in the heating coil and body of e-cigarettes. VEA degradation was analyzed using thermogravimetric analysis (TGA) and gas chromatography/mass spectrometry (GC/MS). While the presence of O2 was found to significantly enhance the degradation of VEA at both high (356 °C) and low (176 °C) temperatures, the addition of Cu-Ni to oxidizing atmospheres was found to greatly enhance VEA degradation, resulting in the formation of numerous degradation products previously identified in VEA vaping emissions. O2 and Cu-Ni nanopowder together were also found to significantly increase the production of OH radicals, which has implications for e-liquid degradation pathways as well as the potential risk of oxidative damage to biological systems in real-world vaping scenarios. Ultimately, the results presented in this study highlight the importance of oxidation pathways in VEA thermal degradation and may aid in the prediction of thermal degradation products from e-liquids.

Neurog2 regulates Isl1 to modulate horizontal cell number

The population sizes of different retinal cell types vary between different strains of mice, and that variation can be mapped to genomic loci in order to identify its polygenic origin. In some cases, controlling genes act independently, whereas in other instances, they exhibit epistasis. Here, we identify an epistatic interaction revealed through the mapping of quantitative trait loci from a panel of recombinant inbred strains of mice. The population of retinal horizontal cells exhibits a twofold variation in number, mapping to quantitative trait loci on chromosomes 3 and 13, where these loci are shown to interact epistatically. We identify a prospective genetic interaction underlying this, mediated by the bHLH transcription factor Neurog2, at the chromosome 3 locus, functioning to repress the LIM homeodomain transcription factor Isl1, at the chromosome 13 locus. Using single and double conditional knockout mice, we confirm the countervailing actions of each gene, and validate in vitro a crucial role for two single nucleotide polymorphisms in the 5'UTR of Isl1, one of which yields a novel E-box, mediating the repressive action of Neurog2.

Activity-Based Sensing for Chemistry-Enabled Biology: Illuminating Principles, Probes, and Prospects for Boronate Reagents for Studying Hydrogen Peroxide

Activity-based sensing (ABS) offers a general approach that exploits chemical reactivity as a method for selective detection and manipulation of biological analytes. Here, we illustrate the value of this chemical platform to enable new biological discovery through a case study in the design and application of ABS reagents for studying hydrogen peroxide (H2O2), a major type of reactive oxygen species (ROS) that regulates a diverse array of vital cellular signaling processes to sustain life. Specifically, we summarize advances in the use of activity-based boronate probes for the detection of H2O2 featuring high molecular selectivity over other ROS, with an emphasis on tailoring designs in chemical structure to promote new biological principles of redox signaling.

Mental health help-seeking among Latina/o/x undocumented college students

OBJECTIVES

Informed by a social-ecological framework, this study nested undocumented students' individual mental health needs within micro-level campus factors and the macro-level immigration policy context to examine how these are associated with undocumented Latina/o/x college students' use of on-campus mental health services.

METHOD

A large-scale survey was administered to 1,277 undocumented college students attending 4-year public universities in California. Only Latina/o/x respondents were included in this study (N = 1,181). Fifty percent of students attended a UC system (n = 589). On average, students were 21.84 years old (SE = .15), and most were women (75.3%, n = 890).

RESULTS

Greater level of mental health symptoms and perceived mental health need, and greater use of campus-wide resources and undocumented student services predicted greater likelihood of using on-campus mental health services. Greater perceptions of social exclusion due to the immigration policy context predicted lower use of mental health services.

CONCLUSIONS

Results indicate that a greater use of resources and an inclusive campus environment, as well as efforts to minimize policy-related feelings of social exclusion, may facilitate undocumented students' professional mental health help-seeking. These findings emphasize the need to take multiple and multi-level ecological factors into account when considering mental health service use, particularly in the case of undocumented immigrants and likely other structurally marginalized groups. (PsycInfo Database Record (c) 2022 APA, all rights reserved).

Distinct Reproductive Risk Profiles for Intrinsic-Like Breast Cancer Subtypes: Pooled Analysis of Population-Based Studies

BACKGROUND

Reproductive factors have been shown to be differentially associated with risk of estrogen receptor (ER)-positive and ER-negative breast cancer. However, their associations with intrinsic-like subtypes are less clear.

METHODS

Analyses included up to 23 353 cases and 71 072 controls pooled from 31 population-based case-control or cohort studies in the Breast Cancer Association Consortium across 16 countries on 4 continents. Polytomous logistic regression was used to estimate the association between reproductive factors and risk of breast cancer by intrinsic-like subtypes (luminal A-like, luminal B-like, luminal B-HER2-like, HER2-enriched-like, and triple-negative breast cancer) and by invasiveness. All statistical tests were 2-sided.

RESULTS

Compared with nulliparous women, parous women had a lower risk of luminal A-like, luminal B-like, luminal B-HER2-like, and HER2-enriched-like disease. This association was apparent only after approximately 10 years since last birth and became stronger with increasing time (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.49 to 0.71; and OR = 0.36, 95% CI = 0.28 to 0.46 for multiparous women with luminal A-like tumors 20 to less than 25 years after last birth and 45 to less than 50 years after last birth, respectively). In contrast, parous women had a higher risk of triple-negative breast cancer right after their last birth (for multiparous women: OR = 3.12, 95% CI = 2.02 to 4.83) that was attenuated with time but persisted for decades (OR = 1.03, 95% CI = 0.79 to 1.34, for multiparous women 25 to less than 30 years after last birth). Older age at first birth (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) and breastfeeding (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) were associated with lower risk of triple-negative breast cancer but not with other disease subtypes. Younger age at menarche was associated with higher risk of all subtypes; older age at menopause was associated with higher risk of luminal A-like but not triple-negative breast cancer. Associations for in situ tumors were similar to luminal A-like.

CONCLUSIONS

This large and comprehensive study demonstrates a distinct reproductive risk factor profile for triple-negative breast cancer compared with other subtypes, with implications for the understanding of disease etiology and risk prediction.

Reference Genome of the Black Surfperch, Embiotoca jacksoni (Embiotocidae, Perciformes), a California Kelp Forest Fish That Lacks a Pelagic Larval Stage

Surfperches (Family Embiotocidae) are viviparous temperate reef fishes that brood their young. This life history trait translates into limited dispersal, strong population structure, and an unusually strong potential for local adaptation in a marine fish. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Black Surfperch, Embiotoca jacksoni, to establish a genomic model for understanding phylogeographic patterns of marine organisms in California. These patterns, in turn, may inform the design of marine protected areas using dispersal models based on genomic data. The genome of E. jacksoni is typical of marine fishes at less than 1Gb (genome size = 635 Mb), and our assembly is near-chromosome level (contig N50 = 6.5Mb, scaffold N50 = 15.5 Mb, BUSCO = 98.1%). Within the context of the CCGP, the genome will be used as a reference for future whole genome resequencing projects aimed at enhancing our knowledge of the population structure of the species, and efficacy of Marine Protected Areas across the state.

NPOmix: A machine learning classifier to connect mass spectrometry fragmentation data to biosynthetic gene clusters

Microbial specialized metabolites are an important source of and inspiration for many pharmaceuticals, biotechnological products and play key roles in ecological processes. Untargeted metabolomics using liquid chromatography coupled with tandem mass spectrometry is an efficient technique to access metabolites from fractions and even environmental crude extracts. Nevertheless, metabolomics is limited in predicting structures or bioactivities for cryptic metabolites. Efficiently linking the biosynthetic potential inferred from (meta)genomics to the specialized metabolome would accelerate drug discovery programs by allowing metabolomics to make use of genetic predictions. Here, we present a k-nearest neighbor classifier to systematically connect mass spectrometry fragmentation spectra to their corresponding biosynthetic gene clusters (independent of their chemical class). Our new pattern-based genome mining pipeline links biosynthetic genes to metabolites that they encode for, as detected via mass spectrometry from bacterial cultures or environmental microbiomes. Using paired datasets that include validated genes-mass spectral links from the Paired Omics Data Platform, we demonstrate this approach by automatically linking 18 previously known mass spectra (17 for which the biosynthesis gene clusters can be found at the MIBiG database plus palmyramide A) to their corresponding previously experimentally validated biosynthetic genes (e.g., via nuclear magnetic resonance or genetic engineering). We illustrated a computational example of how to use our Natural Products Mixed Omics (NPOmix) tool for siderophore mining that can be reproduced by the users. We conclude that NPOmix minimizes the need for culturing (it worked well on microbiomes) and facilitates specialized metabolite prioritization based on integrative omics mining.

Internal sensory neurons regulate stage-specific growth in Drosophila

Animals control their developmental schedule in accordance with internal states and external environments. In Drosophila larvae, it is well established that nutrient status is sensed by different internal organs, which in turn regulate production of insulin-like peptides and thereby control growth. In contrast, the impact of the chemosensory system on larval development remains largely unclear. Here, we performed a genetic screen to identify gustatory receptor (Gr) neurons regulating growth and development, and found that Gr28a-expressing neurons are required for proper progression of larval growth. Gr28a is expressed in a subset of peripheral internal sensory neurons, which directly extend their axons to insulin-producing cells (IPCs) in the central nervous system. Silencing of Gr28a-expressing neurons blocked insulin-like peptide release from IPCs and suppressed larval growth during the mid-larval period. These results indicate that Gr28a-expressing neurons promote larval development by directly regulating growth-promoting endocrine signaling in a stage-specific manner.

Functional metabolite reserves and lipid homeostasis revealed by the MA-10 Leydig cell metabolome

In Leydig cells, intrinsic factors that determine cellular steroidogenic efficiency is of functional interest to decipher and monitor pathophysiology in many contexts. Nevertheless, beyond basic regulation of cholesterol storage and mobilization, systems biology interpretation of the metabolite networks in steroidogenic function is deficient. To reconstruct and describe the different molecular systems regulating steroidogenesis, we profiled the metabolites in resting MA-10 Leydig cells. Our results identified 283-annotated components (82 neutral lipids, 154 membrane lipids, and 47 other metabolites). Neutral lipids were represented by an abundance of triacyglycerols (97.1%), and low levels of cholesterol esters (2.0%). Membrane lipids were represented by an abundance of glycerophospholipids (77.8%), followed by sphingolipids (22.2%). Acylcarnitines, nucleosides, amino acids and their derivatives were the other metabolite classes identified. Among nonlipid metabolites, we recognized substantial reserves of aspartic acid, choline, creatine, betaine, glutamine, homoserine, isoleucine, and pantothenic acid none of which have been previously considered as a requirement in steroidogenic function. Individually limiting use of betaine, choline, or pantothenic acid, during luteinizing hormone-induced steroidogenesis in MA-10 cells resulted in substantial decreases to acute steroidogenic capacity, explained by intermediary metabolite imbalances affecting homeostasis. As such, our dataset represents the current level of baseline characterization and unravels the functional resting state of steroidogenic MA-10 Leydig cells. In identifying metabolite stockpiles and causal mechanisms, these results serve to further comprehend the cellular setup and regulation of steroid biosynthesis.

Sleep loss leads to the withdrawal of human helping across individuals, groups, and large-scale societies

Humans help each other. This fundamental feature of homo sapiens has been one of the most powerful forces sculpting the advent of modern civilizations. But what determines whether humans choose to help one another? Across 3 replicating studies, here, we demonstrate that sleep loss represents one previously unrecognized factor dictating whether humans choose to help each other, observed at 3 different scales (within individuals, across individuals, and across societies). First, at an individual level, 1 night of sleep loss triggers the withdrawal of help from one individual to another. Moreover, fMRI findings revealed that the withdrawal of human helping is associated with deactivation of key nodes within the social cognition brain network that facilitates prosociality. Second, at a group level, ecological night-to-night reductions in sleep across several nights predict corresponding next-day reductions in the choice to help others during day-to-day interactions. Third, at a large-scale national level, we demonstrate that 1 h of lost sleep opportunity, inflicted by the transition to Daylight Saving Time, reduces real-world altruistic helping through the act of donation giving, established through the analysis of over 3 million charitable donations. Therefore, inadequate sleep represents a significant influential force determining whether humans choose to help one another, observable across micro- and macroscopic levels of civilized interaction. The implications of this effect may be non-trivial when considering the essentiality of human helping in the maintenance of cooperative, civil society, combined with the reported decline in sufficient sleep in many first-world nations.

Reconstructing the history of founder events using genome-wide patterns of allele sharing across individuals

Founder events play a critical role in shaping genetic diversity, fitness and disease risk in a population. Yet our understanding of the prevalence and distribution of founder events in humans and other species remains incomplete, as most existing methods require large sample sizes or phased genomes. Thus, we developed ASCEND that measures the correlation in allele sharing between pairs of individuals across the genome to infer the age and strength of founder events. We show that ASCEND can reliably estimate the parameters of founder events under a range of demographic scenarios. We then apply ASCEND to two species with contrasting evolutionary histories: ~460 worldwide human populations and ~40 modern dog breeds. In humans, we find that over half of the analyzed populations have evidence for recent founder events, associated with geographic isolation, modes of sustenance, or cultural practices such as endogamy. Notably, island populations have lower population sizes than continental groups and most hunter-gatherer, nomadic and indigenous groups have evidence of recent founder events. Many present-day groups––including Native Americans, Oceanians and South Asians––have experienced more extreme founder events than Ashkenazi Jews who have high rates of recessive diseases due their known history of founder events. Using ancient genomes, we show that the strength of founder events differs markedly across geographic regions and time––with three major founder events related to the peopling of Americas and a trend in decreasing strength of founder events in Europe following the Neolithic transition and steppe migrations. In dogs, we estimate extreme founder events in most breeds that occurred in the last 25 generations, concordant with the establishment of many dog breeds during the Victorian times. Our analysis highlights a widespread history of founder events in humans and dogs and elucidates some of the demographic and cultural practices related to these events.

A founder event occurs when small numbers of ancestral individuals give rise to a large fraction of the population. Founder events reduce genetic variation and increase the risk of recessive diseases. Despite their importance in evolutionary and disease studies, we still only have a limited comprehension of their prevalence and properties in humans and other species, as most existing methods require large sample sizes or phased genomes. Here, we present a flexible method, ASCEND, to infer the timing and the strength of founder events that is suitable for sparse datasets with few samples or limited coverage. ASCEND provides reliable estimates across a wide range of demographic scenarios. By applying it to data from two species (humans and dogs), we document a widespread history of recent founder events in both species and provide insights about the demographic processes related to these events. Our analysis helps to identify groups with strong founder events that should be prioritized for future studies as they offer a unique opportunity for biological discovery and reducing disease burden through mapping of recessive disease-causing genes and pathways, as previously shown in studies of Ashkenazi Jews and Finns.

Proteomic characterization of phagocytic primary human monocyte-derived macrophages

Macrophages play a vital role in the innate immune system, identifying and destroying unwanted cells. However, it has been difficult to attain a comprehensive understanding of macrophage protein abundance due to technical limitations. In addition, it remains unclear how changes in proteome composition are linked to phagocytic activity. In this study we developed methods to derive human macrophages and prepare them for mass spectrometry analysis in order to more-deeply understand the proteomic consequences of macrophage stimulation. Interferon gamma (IF-g), an immune stimulating cytokine, was used to induce macrophage activation, increasing phagocytosis of cancer cells by 2-fold. These conditions were used to perform comparative shotgun proteomics between resting macrophages and stimulated macrophages with increased phagocytic activity. Our analysis revealed that macrophages bias their protein production toward biological processes associated with phagocytosis and antigen processing in response to stimulation. We confirmed our findings by antibody-based western blotting experiments, validating both previously reported and novel proteins of interest. In addition to whole protein changes, we evaluated active protein synthesis by treating cells with the methionine surrogate probe homopropargylglycine (HPG). We saw increased rates of HPG incorporation during phagocytosis-inducing stimulation, suggesting protein synthesis rates are altered by stimulation. Together our findings provide the most comprehensive proteomic insight to date into primary human macrophages. We anticipate that this data can be used as a launchpoint to generate new hypotheses about innate immune function.

Perspectives on conformationally constrained peptide nucleic acid (PNA): insights into the structural design, properties and applications

Peptide nucleic acid or PNA is a synthetic DNA mimic that contains a sequence of nucleobases attached to a peptide-like backbone derived from N-2-aminoethylglycine. The semi-rigid PNA backbone acts as a scaffold that arranges the nucleobases in a proper orientation and spacing so that they can pair with their complementary bases on another DNA, RNA, or even PNA strand perfectly well through the standard Watson-Crick base-pairing. The electrostatically neutral backbone of PNA contributes to its many unique properties that make PNA an outstanding member of the xeno-nucleic acid family. Not only PNA can recognize its complementary nucleic acid strand with high affinity, but it does so with excellent specificity that surpasses the specificity of natural nucleic acids and their analogs. Nevertheless, there is still room for further improvements of the original PNA in terms of stability and specificity of base-pairing, direction of binding, and selectivity for different types of nucleic acids, among others. This review focuses on attempts towards the rational design of new generation PNAs with superior performance by introducing conformational constraints such as a ring or a chiral substituent in the PNA backbone. A large collection of conformationally rigid PNAs developed during the past three decades are analyzed and compared in terms of molecular design and properties in relation to structural data if available. Applications of selected modified PNA in various areas such as targeting of structured nucleic acid targets, supramolecular scaffold, biosensing and bioimaging, and gene regulation will be highlighted to demonstrate how the conformation constraint can improve the performance of the PNA. Challenges and future of the research in the area of constrained PNA will also be discussed.

Cbp1, a fungal virulence factor under positive selection, forms an effector complex that drives macrophage lysis

Intracellular pathogens secrete effectors to manipulate their host cells. Histoplasma capsulatum (Hc) is a fungal intracellular pathogen of humans that grows in a yeast form in the host. Hc yeasts are phagocytosed by macrophages, where fungal intracellular replication precedes macrophage lysis. The most abundant virulence factor secreted by Hc yeast cells is Calcium Binding Protein 1 (Cbp1), which is absolutely required for macrophage lysis. Here we take an evolutionary, structural, and cell biological approach to understand Cbp1 function. We find that Cbp1 is present only in the genomes of closely related dimorphic fungal species of the Ajellomycetaceae family that lead primarily intracellular lifestyles in their mammalian hosts (Histoplasma, Paracoccidioides, and Emergomyces), but not conserved in the extracellular fungal pathogen Blastomyces dermatitidis. We observe a high rate of fixation of non-synonymous substitutions in the Cbp1 coding sequences, indicating that Cbp1 is under positive selection. We determine the de novo structures of Hc H88 Cbp1 and the Paracoccidioides americana (Pb03) Cbp1, revealing a novel "binocular" fold consisting of a helical dimer arrangement wherein two helices from each monomer contribute to a four-helix bundle. In contrast to Pb03 Cbp1, we show that Emergomyces Cbp1 orthologs are unable to stimulate macrophage lysis when expressed in the Hc cbp1 mutant. Consistent with this result, we find that wild-type Emergomyces africanus yeast are able to grow within primary macrophages but are incapable of lysing them. Finally, we use subcellular fractionation of infected macrophages and indirect immunofluorescence to show that Cbp1 localizes to the macrophage cytosol during Hc infection, making this the first instance of a phagosomal human fungal pathogen directing an effector into the cytosol of the host cell. We additionally show that Cbp1 forms a complex with Yps-3, another known Hc virulence factor that accesses the cytosol. Taken together, these data imply that Cbp1 is a fungal virulence factor under positive selection that localizes to the cytosol to trigger host cell lysis.

Avian-associated Aspergillus fumigatus displays broad phylogenetic distribution, no evidence for host specificity, and multiple genotypes within epizootic events

Birds are highly susceptible to aspergillosis, which can manifest as a primary infection in both domestic and wild birds. Aspergillosis in wild birds causes mortalities ranging in scale from single animals to large-scale epizootic events. However, pathogenicity factors associated with aspergillosis in wild birds have not been examined. Specifically, it is unknown whether wild bird-infecting strains are host-adapted (i.e. phylogenetically related). Similarly, it is unknown whether epizootics are driven by contact with clonal strains that possess unique pathogenic or virulence properties, or by distinct and equally pathogenic strains. Here, we use a diverse collection of Aspergillus fumigatus isolates taken from aspergillosis-associated avian carcasses, representing 24 bird species from a wide geographic range, and representing individual bird mortalities as well as epizootic events. These isolates were sequenced and analyzed along with 130 phylogenetically diverse human clinical isolates to investigate the genetic diversity and phylogenetic placement of avian-associated A. fumigatus, the geographic and host distribution of avian isolates, evidence for clonal outbreaks among wild birds, and the frequency of azole resistance in avian isolates. We found that avian isolates were phylogenetically diverse, with no clear distinction from human clinical isolates, and no sign of host or geographic specificity. Avian isolates from the same epizootic events were diverse and phylogenetically distant, suggesting that avian aspergillosis is not contagious among wild birds and that outbreaks are likely driven by environmental spore loads or host comorbidities. Finally, all avian isolates were susceptible to Voriconazole and none contained the canonical azole resistance gene variants.

Neuroanatomical correlations of visuospatial processing in primary progressive aphasia

Clinical phenotyping of primary progressive aphasia has largely focused on speech and language presentations, leaving other cognitive domains under-examined. This study investigated the diagnostic utility of visuospatial profiles and examined their neural basis among the three main primary progressive aphasia variants. We studied the neuropsychological performances of 118 primary progressive aphasia participants and 30 cognitively normal controls, across 11 measures of visuospatial cognition, and investigated their neural correlates via voxel-based morphometry analysis using visuospatial composite scores derived from principal component analysis. The principal component analysis identified three main factors: visuospatial-executive, visuospatial-memory and visuomotor components. Logopenic variant primary progressive aphasia performed significantly worst across all components; nonfluent/agrammatic variant primary progressive aphasia showed deficits in the visuospatial-executive and visuomotor components compared with controls; and the semantic variant primary progressive aphasia scored significantly lower than nonfluent/agrammatic variant primary progressive aphasia and control in the visuospatial-memory component. Grey matter volumes over the right parieto-occipital cortices correlated with visuospatial-executive performance; volumetric changes in the right anterior parahippocampal gyrus and amygdala were associated with visuospatial-memory function, and visuomotor composite scores correlated significantly with the grey matter volume at the right precentral gyrus. Discriminant function analysis identified three visuospatial measures: Visual Object and Space Perception and Benson figure copy and recall test, which classified 79.7% (94/118) of primary progressive aphasia into their specific variant. This study shows that each primary progressive aphasia variant also carries a distinctive visuospatial cognitive profile that corresponds with grey matter volumetric changes and in turn can be largely represented by their performance on the visuomotor, visuospatial-memory and executive functions.

Unique Transcriptomic Changes Underlie Hormonal Interactions During Mammary Histomorphogenesis in Female Pigs

Successful lactation and the risk for developing breast cancer depend on growth and differentiation of the mammary gland (MG) epithelium that is regulated by ovarian steroids (17β-estradiol [E] and progesterone [P]) and pituitary-derived prolactin (PRL). Given that the MG of pigs share histomorphogenic features present in the normal human breast, we sought to define the transcriptional responses within the MG of pigs following exposure to all combinations of these hormones. Hormone-ablated female pigs were administered combinations of E, medroxyprogesterone 17-acetate (source of P), and either haloperidol (to induce PRL) or 2-bromo-α-ergocryptine. We subsequently monitored phenotypic changes in the MG including mitosis, receptors for E and P (ESR1 and PGR), level of phosphorylated STAT5 (pSTAT5), and the frequency of terminal ductal lobular unit (TDLU) subtypes; these changes were then associated with all transcriptomic changes. Estrogen altered the expression of approximately 20% of all genes that were mostly associated with mitosis, whereas PRL stimulated elements of fatty acid metabolism and an inflammatory response. Several outcomes, including increased pSTAT5, highlighted the ability of E to enhance PRL action. Regression of transcriptomic changes against several MG phenotypes revealed 1669 genes correlated with proliferation, among which 29 were E inducible. Additional gene expression signatures were associated with TDLU formation and the frequency of ESR1 or PGR. These data provide a link between the hormone-regulated genome and phenome of the MG in a species having a complex histoarchitecture like that in the human breast, and highlight an underexplored synergy between the actions of E and PRL during MG development.

Instigators of COVID-19 in Immune Cells Are Increased in Tobacco Cigarette Smokers and Electronic Cigarette Vapers Compared With Nonsmokers

INTRODUCTION

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the virus responsible for the COVID-19 pandemic, gains entry into the host cell when its Spike protein is cleaved by host proteases TMPRSS2 and furin, thereby markedly increasing viral affinity for its receptor, angiotensin-converting enzyme-2 (ACE2). In rodent and diseased human lungs, tobacco cigarette (TCIG) smoke increases ACE2, but the effect of electronic cigarette vaping (ECIG) is unknown. It is unknown whether nicotine (in both TCIGs and ECIGs) or non-nicotine constituents unique to TCIG smoke increase expression of key proteins in COVID-19 pathogenesis.

METHODS

Immune (CD45+) cells collected before the pandemic in otherwise healthy young people, including TCIG smokers (n = 9), ECIG vapers (n = 12), or nonsmokers (NS) (n = 12), were studied. Using flow cytometry, expression of key proteins in COVID-19 pathogenesis were compared among these groups.

RESULTS

TCIG smokers and ECIG vapers had similar smoking or vaping burdens as indicated by similar plasma cotinine levels. TCIG smokers compared with NS had a significantly increased percentage of cells that were positive for ACE2 (10-fold, p < .001), TMPRSS2 (5-fold, p < .001), and ADAM17 (2.5-fold, p < .001). Additionally, the mean fluorescence intensity (MFI) consistently showed greater mean ACE2 (2.2-fold, p < .001), TMPRSS2 (1.5-fold, p < .001), furin (1.1-fold, p < .05), and ADAM17 (2-fold, p < .001) in TCIG smokers compared with NS. In ECIG vapers, furin MFI was increased (1.15-fold, p < .05) and TMPRSS2 MFI tended to be increased (1.1-fold, p = .077) compared with NS.

CONCLUSIONS

The finding that key instigators of COVID-19 infection are lower in ECIG vapers compared with TCIG smokers is intriguing and warrants additional investigation to determine if switching to ECIGs is an effective harm reduction strategy. However, the trend toward increased proteases in ECIG vapers remains concerning.

IMPLICATIONS

(1) This is the first human study to report a marked increase in proteins critical for COVID-19 infection, including ACE2, TMPRSS2, and ADAM17, in immune cells from healthy tobacco cigarette smokers without lung disease compared with e-cigarette vapers and nonsmokers. (2) These findings warrant additional investigation to determine whether switching to electronic cigarettes may be an effective harm reduction strategy in smokers addicted to nicotine who are unable or unwilling to quit. (3) The increase in proteases in electronic cigarette vapers remains concerning.

Evolutionary systems biology reveals patterns of rice adaptation to drought-prone agro-ecosystems

Rice (Oryza sativa) was domesticated around 10,000 years ago and has developed into a staple for half of humanity. The crop evolved and is currently grown in stably wet and intermittently dry agro-ecosystems, but patterns of adaptation to differences in water availability remain poorly understood. While previous field studies have evaluated plant developmental adaptations to water deficit, adaptive variation in functional and hydraulic components, particularly in relation to gene expression, has received less attention. Here, we take an evolutionary systems biology approach to characterize adaptive drought resistance traits across roots and shoots. We find that rice harbors heritable variation in molecular, physiological, and morphological traits that is linked to higher fitness under drought. We identify modules of co-expressed genes that are associated with adaptive drought avoidance and tolerance mechanisms. These expression modules showed evidence of polygenic adaptation in rice subgroups harboring accessions that evolved in drought-prone agro-ecosystems. Fitness-linked expression patterns allowed us to identify the drought-adaptive nature of optimizing photosynthesis and interactions with arbuscular mycorrhizal fungi. Taken together, our study provides an unprecedented, integrative view of rice adaptation to water-limited field conditions.