Bayesian spatial modeling of moose count data: increasing estimator efficiency and exploring ecological hypotheses

Moose management throughout much of Alaska and Canada relies on aerial count data, which are commonly collected using the geospatial population estimator (GSPE) protocol. The GSPE uses a model-based analytical approach and finite-population block kriging to estimate abundance from a collection of sampled survey units. Widespread implementation and well-documented analytical software have resulted in reliable estimates of moose abundance, density, and composition across a large proportion of their range. Analysis is conducted almost exclusively using the GSPE software, which fits a fixed model structure to data collected within a single year. The downside of this approach to analysis is that the fixed model structure is inefficient for estimation, leading to more field effort than would otherwise be necessary to achieve a desired level of estimator precision. We developed a more easily modified and flexible Bayesian spatial general additive model approach (BSG) that accommodates spatial and temporal covariates (e.g., habitat characteristics, trend), multiple survey events, prior information, and incomplete detection. Using a series of 6 GSPE surveys conducted in Yukon-Charley Rivers National Preserve, Alaska, USA, from 2003–2019, we established the equivalence of the 2 approaches under similar model structures. We then extended the BSG to demonstrate how a more comprehensive approach to analysis can affect estimator precision and be used to assess ecological relationships. The precision of annual abundance estimators from the BSG were improved by an average of 43% over those based on the standard GSPE analysis, highlighting the very real costs of assuming a fixed (i.e., suboptimal) model structure. The population increased at a rate of 2.3%/year (95% CrI = 0.8–3.8%), and the increase was largely explained by a parallel increase in wildfire extent (i.e., high quality moose habitat). These results suggest that our approach could be used to increase estimator efficiency or decrease future survey costs without any modifications to the basic protocol. While modification of the GSPE software is possible, practitioners may find the BSG approach more convenient for quickly developing model structures for a particular application, thereby allowing them to extract more information from existing and future datasets.     

Progesterone Receptor Membrane Component 1 Is a Functional Part of the Glucagon-like Peptide-1 (GLP-1) Receptor Complex in Pancreatic β Cells

Glucagon-like peptide-1 (GLP-1) is an incretin hormone that regulates glucose homeostasis. Because of their direct stimulation of insulin secretion from pancreatic β cells, GLP-1 receptor (GLP-1R) agonists are now important therapeutic options for the treatment of type 2 diabetes. To better understand the mechanisms that control the insulinotropic actions of GLP-1, affinity purification and mass spectrometry (AP-MS) were employed to uncover potential proteins that functionally interact with the GLP-1R. AP-MS performed on Chinese hamster ovary cells or MIN6 β cells, both expressing the human GLP-1R, revealed 99 proteins potentially associated with the GLP-1R. Three novel GLP-1R interactors (PGRMC1, Rab5b, and Rab5c) were further validated through co-immunoprecipitation/immunoblotting, fluorescence resonance energy transfer, and immunofluorescence. Functional studies revealed that overexpression of PGRMC1, a novel cell surface receptor that associated with liganded GLP-1R, enhanced GLP-1-induced insulin secretion (GIIS) with the most robust effect. Knockdown of PGRMC1 in β cells decreased GIIS, indicative of positive interaction with GLP-1R. To gain insight mechanistically, we demonstrated that the cell surface PGRMC1 ligand P4-BSA increased GIIS, whereas its antagonist AG-205 decreased GIIS. It was then found that PGRMC1 increased GLP-1-induced cAMP accumulation. PGRMC1 activation and GIIS induced by P4-BSA could be blocked by inhibition of adenylyl cyclase/EPAC signaling or the EGF receptor–PI3K signal transduction pathway. These data reveal a dual mechanism for PGRMC1-increased GIIS mediated through cAMP and EGF receptor signaling. In conclusion, we identified several novel GLP-1R interacting proteins. PGRMC1 expressed on the cell surface of β cells was shown to interact with the activated GLP-1R to enhance the insulinotropic actions of GLP-1.Glucagon-like peptide-1 (GLP-1)¹ is a gastrointestinal hormone secreted by intestinal L cells upon food intake that is best known for its role in controlling glucose homeostasis. Acting through its cognate glucagon-like peptide-1 receptor (GLP-1R), GLP-1 has several important physiological and pharmacological functions. GLP-1 is best known for enhancing glucose-stimulated insulin secretion (GSIS) from the pancreatic β cells. Importantly, the insulinotropic properties of GLP-1 are maintained in patients with type 2 diabetes (1), which is characterized by insufficient insulin secretion from pancreatic β cells and an inability to maintain glucose homeostasis. Therefore, therapeutic strategies targeting GLP-1R have been developed to treat type 2 diabetes (2, 3). In addition to augmenting insulin secretion, GLP-1 has been known to improve glucose sensing, proinsulin biosynthesis, survival, and proliferation of β cells (3, 4) in a variety of experimental models. GLP-1 also has several extrapancreatic effects, including actions on the central nervous system to inhibit food intake (5), the stomach to decrease gastric emptying and gastric acid secretion (6), and the lungs to stimulate secretion of macromolecules from airways (7). Additionally, GLP-1 has an effect on the heart and possibly the kidney to modulate blood pressure and heart rate (8, 9).The GLP-1R is a member of the B1 family of G protein–coupled receptors (secretin receptor family). In mammals, GLP-1R is expressed in multiple tissues, including pancreatic β cells and δ cells (10), hypothalamus, lung, stomach, heart, kidney (11), and thyroid (12), which in part explains its diverse actions. Upon ligand binding, the GLP-1R is capable of coupling to diverse cell signal transduction pathways, but it is best known for its actions on G protein Gs α and adenylate cyclase activity to increase intracellular cAMP. It is known that other proteins can affect GLP-1R activity in addition to G proteins, including β-arrestin and caveolin, which affect receptor internalization and trafficking. β-Arrestin 1 is also required for proper GLP-1-stimulated cAMP production (13–15). More recently, it was shown that another B1 family member, gastric inhibitory polypeptide receptor heterodimerizes with GLP-1R, decreasing GLP-1-induced β-arrestin recruitment and mobilization (16). Very recently, our group identified several novel potential GLP-1R interactors using a membrane-based split-ubiquitin yeast two-hybrid (MYTH) assay (17). Three β cell–expressing membrane-bound interactors, solute carrier family 15 member 4 (SLC15A4), amyloid β A4 precursor-like protein 1 (APLP1), and adaptor-related protein complex 2 subunit mu (AP2M1), were further selected for individual knockdown in mouse insulinoma (MIN6) β cells using small interfering RNAs (siRNAs). GLP-1-induced insulin secretion was significantly enhanced when these genes were silenced, suggesting that these interactor proteins attenuate GLP-1R activity. These findings demonstrated that GLP-1R protein interactions are complex and the interactors can have measurable effects on receptor trafficking and downstream signaling. Such interactions may in part explain the diverse tissue-specific effects of GLP-1 and offer avenues for controlling GLP-1 actions in a tissue-selective manner.Although the MYTH system is well established (18) and has been applied to study G protein–coupled receptor interactomes (17), it is limited on two fronts. Firstly, it must be performed in yeast which is not an ideal representation of the mammalian system. Secondly, it is technically difficult to activate the receptor in MYTH, thus, effects of ligand stimulation on the receptor interactome cannot be assessed. Recently, affinity purification–mass spectrometry (AP-MS) has become a powerful tool for discovering and examining novel protein–protein interactions, including those between membrane-bound proteins in mammalian cells (19–21). In the current study, we applied AP-MS to discover novel GLP-1R interactors and employed a human GLP-1R harboring a FLAG^® epitope. GLP-1R-Flag was expressed in either Chinese hamster ovary (CHO) cells or MIN6 β cells, and interactors were studied in the presence or absence of GLP-1.     

Epstein-Barr virus nuclear antigen 3C (EBNA3C) interacts with the metabolism sensing C-terminal binding protein (CtBP) repressor to upregulate host genes

Epstein-Barr virus (EBV) infection is associated with the development of specific types of lymphoma and some epithelial cancers. EBV infection of resting B-lymphocytes in vitro drives them to proliferate as lymphoblastoid cell lines (LCLs) and serves as a model for studying EBV lymphomagenesis. EBV nuclear antigen 3C (EBNA3C) is one of the genes required for LCL growth and previous work has suggested that suppression of the CDKN2A encoded tumor suppressor p16^INK4A and possibly p14^ARF is central to EBNA3C’s role in this growth transformation. To directly assess whether loss of p16 and/or p14 was sufficient to explain EBNA3C growth effects, we used CRISPR/Cas9 to disrupt specific CDKN2A exons in EBV transformed LCLs. Disruption of p16 specific exon 1α and the p16/p14 shared exon 2 were each sufficient to restore growth in the absence of EBNA3C. Using EBNA3C conditional LCLs knocked out for either exon 1α or 2, we identified EBNA3C induced and repressed genes. By trans-complementing with EBNA3C mutants, we determined specific genes that require EBNA3C interaction with RBPJ or CtBP for their regulation. Unexpectedly, interaction with the CtBP repressor was required not only for repression, but also for EBNA3C induction of many host genes. Contrary to previously proposed models, we found that EBNA3C does not recruit CtBP to the promoters of these genes. Instead, our results suggest that CtBP is bound to these promoters in the absence of EBNA3C and that EBNA3C interaction with CtBP interferes with the repressive function of CtBP, leading to EBNA3C mediated upregulation.     

Seasonal Shifts in Nocturnal Habitat Use by Coastal Bat Species

Sensitivity of bats to land use change depends on their foraging ecology, which varies among species based on ecomorphological traits. Additionally, because prey availability, vegetative clutter, and temperature change throughout the year, some species may display seasonal shifts in their nocturnal habitat use. In the Coastal Plain of South Carolina, USA, the northern long-eared bat (Myotis septentrionalis), southeastern myotis (Myotis austroriparius), tri-colored bat (Perimyotis subflavus), and northern yellow bat (Lasiurus intermedius) are species of conservation concern that are threatened by habitat loss. Our objective was to identify characteristics of habitat used by these species during their nightly active period and compare use between summer and winter. We conducted acoustic surveys at 125 sites during May–August and at 121 of the same 125 sites December–March 2018 and 2019 in upland forests, bottomland forests, fields, ponds, and salt marsh and used occupancy models to assess habitat use. The northern long-eared bat and southeastern myotis (i.e., myotis bats) used sites that were closer to hardwood stands, pine stands, and fresh water year-round. We did not identify any strong predictors of tri-colored bat habitat use in summer, but during winter they used bottomland forests, fields, and ponds more than salt marsh and upland forests. During summer and winter, northern yellow bats used sites close to fresh water and salt marsh. Additionally, during summer they used fields, ponds, and salt marsh more than upland and bottomland forests, but in winter they used bottomland forests, fields, and ponds more than upland forest and salt marsh. Our results highlight important land cover types for bats in this area (e.g., bottomland forests, ponds, and salt marsh), and that habitat use changes between seasons. Accounting for and understanding how habitat use changes throughout the year will inform managers about how critical habitat features may vary in their importance to bats throughout the year. © 2021 The Wildlife Society.     

Social structure and facultative mating systems of hoary marmots (Marmota caligata)

Mate-choice theory predicts different optimal mating systems depending on resource availability and habitat stability. Regions with limited resources are thought to promote monogamy. We tested predictions of monogamy in a social rodent, the hoary marmot (Marmota caligata), at the northern climatic extreme of its distribution. Mating systems, social structure and genetic relationships were investigated within and among neighbouring colonies of marmots within a 4 km(2) valley near Kluane National Park, Yukon, Canada, using 21 microsatellite loci. While both monogamous and polygynous populations of hoary marmots have been observed in the southern reaches of this species' range; northern populations of this species are thought to be predominantly monogamous. Contrary to previous studies, we did not find northern hoary marmot social groups to be predominantly monogamous; rather, the mating system seemed to be facultative, varying between monogamy and polygyny within, as well as among, social groups. These findings reveal that the mating systems within colonies of this species are more flexible than previously thought, potentially reflecting local variation in resource availability.     

Fluorogenic cephalosporin substrates for β-lactamase TEM-1

Cephalosporin was used to synthesize soluble and precipitating fluorogenic β-lactam substrates that demonstrated differential catalytic hydrolysis by three different subtypes of β-lactamase: TEM-1 (class A), p99 (class C), and a Bacillus cereus enzyme sold by Genzyme (class B). The most successful soluble substrate contained difluorofluorescein (Oregon Green 488) ligated to two cephalosporin moieties that, therefore, required two turnovers to produce the fluorescent Oregon Green 488 leaving group. The bis-cephalosporin modification was required so that the final reaction product was the Oregon Green 488 carboxylic acid rather than a less bright phenolic adduct of the dye. Hydrolysis in pH 5.5 Mes and pH 7.2 phosphate-buffered saline (PBS) buffers was similar, but in pH 8.0 Tris the hydrolysis rate nearly doubled. Activity of the β-lactamases on the various substrates was shown to depend highly on the linker between the cephalosporin and the fluorophore, with an allyl linker promoting faster turnover than a phenol ether linker. Measured K_m values for dichlorofluorescein and difluorofluorescein cephalosporin substrates were approximately the same as K_m values for penicillin G and ampicillin found in the literature (∼30–40 μM).     

4E10-Resistant HIV-1 Isolated from Four Subjects with Rare Membrane-Proximal External Region Polymorphisms

Human antibody 4E10 targets the highly conserved membrane-proximal external region (MPER) of the HIV-1 transmembrane glycoprotein, gp41, and has extraordinarily broad neutralizing activity. It is considered by many to be a prototype for vaccine development. In this study, we describe four subjects infected with viruses carrying rare MPER polymorphisms associated with resistance to 4E10 neutralization. In one case resistant virus carrying a W680G substitution was transmitted from mother to infant. We used site-directed mutagenesis to demonstrate that the W680G substitution is necessary for conferring the 4E10-resistant phenotype, but that it is not sufficient to transfer the phenotype to a 4E10-sensitive Env. Our third subject carried Envs with a W680R substitution causing variable resistance to 4E10, indicating that residues outside the MPER are required to confer the phenotype. A fourth subject possessed a F673L substitution previously associated with 4E10 resistance. For all three subjects with W680 polymorphisms, we observed additional residues in the MPER that co-varied with position 680 and preserved charged distributions across this region. Our data provide important caveats for vaccine development targeting the MPER. Naturally occurring Env variants described in our study also represent unique tools for probing the structure-function of HIV-1 envelope.     

Dioxane contributes to the altered conformation and oligomerization state of a designed engrailed homeodomain variant

Our goal was to compute a stable, full-sequence design of the Drosophila melanogaster engrailed homeodomain. Thermal and chemical denaturation data indicated the design was significantly more stable than was the wild-type protein. The data were also nearly identical to those for a similar, later full-sequence design, which was shown by NMR to adopt the homeodomain fold: a three-helix, globular monomer. However, a 1.65 A crystal structure of the design described here turned out to be of a completely different fold: a four-helix, rodlike tetramer. The crystallization conditions included approximately 25% dioxane, and subsequent experiments by circular dichroism and sedimentation velocity analytical ultracentrifugation indicated that dioxane increases the helicity and oligomerization state of the designed protein. We attribute at least part of the discrepancy between the target fold and the crystal structure to the presence of a high concentration of dioxane.     

Enumeration of Mycobacterium leprae using real-time PCR

Mycobacterium leprae is not cultivable in axenic media, and direct microscopic enumeration of the bacilli is complex, labor intensive, and suffers from limited sensitivity and specificity. We have developed a real-time PCR assay for quantifying M. leprae DNA in biological samples. Primers were identified to amplify a shared region of the multicopy repeat sequence (RLEP) specific to M. leprae and tested for sensitivity and specificity in the TaqMan format. The assay was specific for M. leprae and able to detect 10 fg of purified M. leprae DNA, or approximately 300 bacteria in infected tissues. We used the RLEP TaqMan PCR to assess the short and long-term growth results of M. leprae in foot pad tissues obtained from conventional mice, a gene knock-out mouse strain, athymic nude mice, as well as from reticuloendothelial tissues of M. leprae-infected nine-banded armadillos. We found excellent correlative results between estimates from RLEP TaqMan PCR and direct microscopic counting (combined r = 0.98). The RLEP TaqMan PCR permitted rapid analysis of batch samples with high reproducibility and is especially valuable for detection of low numbers of bacilli. Molecular enumeration is a rapid, objective and highly reproducible means to estimate the numbers of M. leprae in tissues, and application of the technique can facilitate work with this agent in many laboratories.     

Genome Analysis Traces Regional Dispersal of Rice in Taiwan and Southeast Asia

The dispersal of rice (Oryza sativa) following domestication influenced massive social and cultural changes across South, East, and Southeast (SE) Asia. The history of dispersal across islands of SE Asia, and the role of Taiwan and the Austronesian expansion in this process remain largely unresolved. Here, we reconstructed the routes of dispersal of O. sativa ssp. japonica rice to Taiwan and the northern Philippines using whole-genome resequencing of indigenous rice landraces coupled with archaeological and paleoclimate data. Our results indicate that japonica rice found in the northern Philippines diverged from Indonesian landraces as early as 3,500 years before present (BP). In contrast, rice cultivated by the indigenous peoples of the Taiwanese mountains has complex origins. It comprises two distinct populations, each best explained as a result of admixture between temperate japonica that presumably came from northeast Asia, and tropical japonica from the northern Philippines and mainland SE Asia, respectively. We find that the temperate japonica component of these indigenous Taiwan populations diverged from northeast Asia subpopulations at about 2,600 BP, whereas gene flow from the northern Philippines had begun before ∼1,300  BP. This coincides with a period of intensified trade established across the South China Sea. Finally, we find evidence for positive selection acting on distinct genomic regions in different rice subpopulations, indicating local adaptation associated with the spread of japonica rice.