Fat‐1 expression alleviates atherosclerosis in transgenic rabbits

Atherosclerosis is the main cause of cardiovascular diseases. The Fat‐1 gene can express the n‐3 fatty acid desaturase, which converts n‐6 polyunsaturated fatty acids (PUFA) to n‐3 PUFAs. The role of n‐3 PUFAs in atherosclerosis is widely debated. This study explored the effect of n‐3 PUFAs on atherosclerosis in rabbits. In this study, atherosclerosis was induced in Fat‐1 transgenic rabbits and their littermate (WT) rabbits by feeding a high‐cholesterol diet containing 0.3% cholesterol and 3% soybean oil for 16 weeks. Plasma lipid, fatty acid and pathological analyses of atherosclerotic lesions were conducted. Fatty acid composition in the liver and muscle showed that n‐3 PUFAs increased and n‐6 PUFAs decreased in the Fat‐1 group. Plasma high‐density lipoprotein cholesterol (HDL‐C) levels were significantly increased in the Fat‐1 group, and the atherosclerotic lesion area of the aortic arch in Fat‐1 transgenic rabbits was significantly reduced. Histological analysis showed that smooth muscle cells (SMCs) in atherosclerotic lesions decreased significantly. In conclusion, n‐3 PUFAs improve atherosclerosis in Fat‐1 transgenic rabbits, and this process may depend on the increase in plasma HDL‐C and the decrease in the amount of SMCs in atherosclerotic plaques.     

A genotyping‐in‐thousands by sequencing panel to inform invasive deer management using noninvasive fecal and hair samples

Studies in ecology, evolution, and conservation often rely on noninvasive samples, making it challenging to generate large amounts of high‐quality genetic data for many elusive and at‐risk species. We developed and optimized a Genotyping‐in‐Thousands by sequencing (GT‐seq) panel using noninvasive samples to inform the management of invasive Sitka black‐tailed deer (Odocoileus hemionus sitkensis) in Haida Gwaii (Canada). We validated our panel using paired high‐quality tissue and noninvasive fecal and hair samples to simultaneously distinguish individuals, identify sex, and reconstruct kinship among deer sampled across the archipelago, then provided a proof‐of‐concept application using field‐collected feces on SGang Gwaay, an island of high ecological and cultural value. Genotyping success across 244 loci was high (90.3%) and comparable to that of high‐quality tissue samples genotyped using restriction‐site associated DNA sequencing (92.4%), while genotyping discordance between paired high‐quality tissue and noninvasive samples was low (0.50%). The panel will be used to inform future invasive species operations in Haida Gwaii by providing individual and population information to inform management. More broadly, our GT‐seq workflow that includes quality control analyses for targeted SNP selection and a modified protocol may be of wider utility for other studies and systems where noninvasive genetic sampling is employed.     

High‐altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high‐altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high‐altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein‐coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high‐altitude vertebrates than in low‐altitude vertebrates. The seven rapidly evolving genes were shared by the high‐altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high‐altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high‐altitude vertebrates was higher than in low‐altitude vertebrates as their evolution requires more energy in a high‐altitude environment. Our study demonstrates the high‐altitude environment (low atmospheric O₂ levels) drives vertebrate evolution in mtDNA PCGs.     

Sex‐ and strain‐specific effects of mitochondrial uncoupling on age‐related metabolic diseases in high‐fat diet‐fed mice

Mild uncoupling of oxidative phosphorylation is an intrinsic property of all mitochondria and may have evolved to protect cells against the production of damaging reactive oxygen species. Therefore, compounds that enhance mitochondrial uncoupling are potentially attractive anti‐aging therapies; however, chronic ingestion is associated with a number of unwanted side effects. We have previously developed a controlled‐release mitochondrial protonophore (CRMP) that is functionally liver‐directed and promotes oxidation of hepatic triglycerides by causing a subtle sustained increase in hepatic mitochondrial inefficiency. Here, we sought to leverage the higher therapeutic index of CRMP to test whether mild mitochondrial uncoupling in a liver‐directed fashion could reduce oxidative damage and improve age‐related metabolic disease and lifespan in diet‐induced obese mice. Oral administration of CRMP (20 mg/[kg‐day] × 4 weeks) reduced hepatic lipid content, protein kinase C epsilon activation, and hepatic insulin resistance in aged (74‐week‐old) high‐fat diet (HFD)‐fed C57BL/6J male mice, independently of changes in body weight, whole‐body energy expenditure, food intake, or markers of hepatic mitochondrial biogenesis. CRMP treatment was also associated with a significant reduction in hepatic lipid peroxidation, protein carbonylation, and inflammation. Importantly, long‐term (49 weeks) hepatic mitochondrial uncoupling initiated late in life (94–104 weeks), in conjugation with HFD feeding, protected mice against neoplastic disorders, including hepatocellular carcinoma (HCC), in a strain and sex‐specific manner. Taken together, these studies illustrate the complex variation of aging and provide important proof‐of‐concept data to support further studies investigating the use of liver‐directed mitochondrial uncouplers to promote healthy aging in humans.     

Biparatopic nanobodies protect mice from lethal challenge with SARS‐CoV‐2 variants of concern

The ongoing COVID‐19 pandemic and the emergence of new SARS‐CoV‐2 variants of concern (VOCs) requires continued development of effective therapeutics. Recently, we identified high‐affinity neutralizing nanobodies (Nbs) specific for the receptor‐binding domain (RBD) of SARS‐CoV‐2. Taking advantage of detailed epitope mapping, we generate two biparatopic Nbs (bipNbs) targeting a conserved epitope outside and two different epitopes inside the RBD:ACE2 interface. Both bipNbs bind all currently circulating VOCs with high affinities and are capable to neutralize cellular infection with VOC B.1.351 (Beta) and B.1.617.2 (Delta) in vitro. To assess if the bipNbs NM1267 and NM1268 confer protection against SARS‐CoV‐2 infection in vivo, human ACE2 transgenic mice are treated intranasally before infection with a lethal dose of SARS‐CoV‐2 B.1, B.1.351 (Beta) or B.1.617.2 (Delta). Nb‐treated mice show significantly reduced disease progression and increased survival rates. Histopathological analyses further reveal a drastically reduced viral load and inflammatory response in lungs. These data suggest that both bipNbs are broadly active against a variety of emerging SARS‐CoV‐2 VOCs and represent easily applicable drug candidates.     

Engineering defensin α‐helix to produce high‐affinity SARS‐CoV‐2 spike protein binding ligands

The binding of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) spike protein to the angiotensin‐converting enzyme 2 (ACE2) receptor expressed on the host cells is a critical initial step for viral infection. This interaction is blocked through competitive inhibition by soluble ACE2 protein. Therefore, developing high‐affinity and cost‐effective ACE2 mimetic ligands that disrupt this protein–protein interaction is a promising strategy for viral diagnostics and therapy. We employed human and plant defensins, a class of small (2–5 kDa) and highly stable proteins containing solvent‐exposed alpha‐helix, conformationally constrained by two disulfide bonds. Therefore, we engineered the amino acid residues on the constrained alpha‐helix of defensins to mimic the critical residues on the ACE2 helix 1 that interact with the SARS‐CoV‐2 spike protein. The engineered proteins (h‐deface2, p‐deface2, and p‐deface2‐MUT) were soluble and purified to homogeneity with a high yield from a bacterial expression system. The proteins demonstrated exceptional thermostability (Tm 70.7°C), high‐affinity binding to the spike protein with apparent K _d values of 54.4 ± 11.3, 33.5 ± 8.2, and 14.4 ± 3.5 nM for h‐deface2, p‐deface2, and p‐deface2‐MUT, respectively, and were used in a diagnostic assay that detected SARS‐CoV‐2 neutralizing antibodies. This work addresses the challenge of developing helical ACE2 mimetics by demonstrating that defensins provide promising scaffolds to engineer alpha‐helices in a constrained form for designing of high‐affinity ligands.     

Characterization of an aminotransferase from Acanthamoeba polyphaga Mimivirus

Acanthamoeba polyphaga Mimivirus, a complex virus that infects amoeba, was first reported in 2003. It is now known that its DNA genome encodes for nearly 1,000 proteins including enzymes that are required for the biosynthesis of the unusual sugar 4‐amino‐4,6‐dideoxy‐d‐glucose, also known as d‐viosamine. As observed in some bacteria, the pathway for the production of this sugar initiates with a nucleotide‐linked sugar, which in the Mimivirus is thought to be UDP‐d‐glucose. The enzyme required for the installment of the amino group at the C‐4′ position of the pyranosyl moiety is encoded in the Mimivirus by the L136 gene. Here, we describe a structural and functional analysis of this pyridoxal 5′‐phosphate‐dependent enzyme, referred to as L136. For this analysis, three high‐resolution X‐ray structures were determined: the wildtype enzyme/pyridoxamine 5′‐phosphate/dTDP complex and the site‐directed mutant variant K185A in the presence of either UDP‐4‐amino‐4,6‐dideoxy‐d‐glucose or dTDP‐4‐amino‐4,6‐dideoxy‐d‐glucose. Additionally, the kinetic parameters of the enzyme utilizing either UDP‐d‐glucose or dTDP‐d‐glucose were measured and demonstrated that L136 is efficient with both substrates. This is in sharp contrast to the structurally related DesI from Streptomyces venezuelae, whose three‐dimensional architecture was previously reported by this laboratory. As determined in this investigation,DesI shows a profound preference in its catalytic efficiency for the dTDP‐linked sugar substrate. This difference can be explained in part by a hydrophobic patch in DesI that is missing in L136. Notably, the structure of L136 reported here represents the first three‐dimensional model for a virally encoded PLP‐dependent enzyme and thus provides new information on sugar aminotransferases in general.     

To minimize foraging time,use high‐efficiency,energy‐expensive search and capture methods when food is abundant but low‐efficiency,low‐cost methods during food shortages

Based on a mathematical model, I show that the amount of food in the habitat determines which among alternative methods for search of prey, respectively, for pursuit‐and‐capture give the shortest daily foraging time. The higher the locomotor activity, the higher the rate of energy expenditure and the larger the habitat space a predator can search for prey per time unit. Therefore, I assume that the more efficient a foraging method is, the higher its rate of energy expenditure. Survival selection favors individuals that use foraging methods that cover their energy needs in the shortest possible time. Therefore, I take the optimization criterion to be minimization of the daily foraging time or, equivalently, maximization of the rate of net energy gain. When time is limiting and food is in short supply, as during food bottleneck periods, low‐efficiency, low‐cost foraging methods give shorter daily foraging times than high‐efficiency, energy‐expensive foraging methods. When time is limiting, food is abundant and energy needs are large, as during reproduction, high‐efficiency high‐cost foraging methods give shorter daily foraging times than low‐efficiency low‐cost foraging methods. When time is not limiting, food is abundant, and energy needs are small, the choice of foraging method is not critical. Small animals have lower rates of energy expenditure for locomotion than large animals. At a given food density and with similar diet, small animals are therefore more likely than large ones to minimize foraging time by using high‐efficiency energy‐expansive foraging methods and to exploit patches and sites that require energy‐demanding locomotion modes. Survival selection takes place at food shortages, while low‐efficiency low‐cost foraging methods are used, whereas reproduction selection occurs when food is abundant and high‐efficiency energy‐expensive foraging methods do better. In seasonal environments, selection therefore acts on different foraging methods at different times. Morphological adaptation to one method may oppose adaptation to another. Such conflicts select against foraging and morphological specialization and tend to give species‐poor communities of year‐round resident generalists. But a stable year‐round food supply favors specialization, niche narrowing, and dense species packing.     

High‐resolution analysis of the conformational transition of pro‐apoptotic Bak at the lipid membrane

Permeabilization of the outer mitochondrial membrane by pore‐forming Bcl2 proteins is a crucial step for the induction of apoptosis. Despite a large set of data suggesting global conformational changes within pro‐apoptotic Bak during pore formation, high‐resolution structural details in a membrane environment remain sparse. Here, we used NMR and HDX‐MS (Hydrogen deuterium exchange mass spectrometry) in lipid nanodiscs to gain important high‐resolution structural insights into the conformational changes of Bak at the membrane that are dependent on a direct activation by BH3‐only proteins. Furthermore, we determined the first high‐resolution structure of the Bak transmembrane helix. Upon activation, α‐helix 1 in the soluble domain of Bak dissociates from the protein and adopts an unfolded and dynamic potentially membrane‐bound state. In line with this finding, comparative protein folding experiments with Bak and anti‐apoptotic BclxL suggest that α‐helix 1 in Bak is a metastable structural element contributing to its pro‐apoptotic features. Consequently, mutagenesis experiments aimed at stabilizing α‐helix 1 yielded Bak variants with delayed pore‐forming activity. These insights will contribute to a better mechanistic understanding of Bak‐mediated membrane permeabilization.     

Glucagon signalling in the dorsal vagal complex is sufficient and necessary for high‐protein feeding to regulate glucose homeostasis in vivo

High‐protein feeding acutely lowers postprandial glucose concentration compared to low‐protein feeding, despite a dichotomous rise of circulating glucagon levels. The physiological role of this glucagon rise has been largely overlooked. We here first report that glucagon signalling in the dorsal vagal complex (DVC) of the brain is sufficient to lower glucose production by activating a Gcgr–PKA–ERK–K_ATP channel signalling cascade in the DVC of rats in vivo. We further demonstrate that direct blockade of DVC Gcgr signalling negates the acute ability of high‐ vs. low‐protein feeding to reduce plasma glucose concentration, indicating that the elevated circulating glucagon during high‐protein feeding acts in the brain to lower plasma glucose levels. These data revise the physiological role of glucagon and argue that brain glucagon signalling contributes to glucose homeostasis during dietary protein intake.     

Modeling population size independent tissue epigenomes by ChIL‐seq with single thin sections

Recent advances in genome‐wide technologies have enabled analyses using small cell numbers of even single cells. However, obtaining tissue epigenomes with cell‐type resolution from large organs and tissues still remains challenging, especially when the available material is limited. Here, we present a ChIL‐based approach for analyzing the diverse cellular dynamics at the tissue level using high‐depth epigenomic data. “ChIL for tissues” allows the analysis of a single tissue section and can reproducibly generate epigenomic profiles from several tissue types, based on the distribution of target epigenomic states, tissue morphology, and number of cells. The proposed method enabled the independent evaluation of changes in cell populations and gene activation in cells from regenerating skeletal muscle tissues, using a statistical model of RNA polymerase II distribution on gene loci. Thus, the integrative analyses performed using ChIL can elucidate in vivo cell‐type dynamics of tissues.     

Tacrolimus alleviates LPS‐induced AKI by inhibiting TLR4/MyD88/NF‐κB signalling in mice

Lipopolysaccharide (LPS)‐induced sepsis‐associated acute kidney injury (SA‐AKI) is a model of clinical serious care syndrome, with high morbidity and mortality. Tacrolimus (TAC), a novel immunosuppressant that inhibits inflammatory response, plays a pivotal role in kidney diseases. In this study, LPS treated mice and cultured podocytes were used as the models of SA‐AKI in vivo and in vitro, respectively. Medium‐ and high‐dose TAC administration significantly attenuated renal function and renal pathological manifestations at 12, 24 and 48 h after LPS treatment in mice. Moreover, the Toll‐like receptor 4 (TLR4)/myeloid differential protein‐88 (MyD88)/nuclear factor‐kappa (NF‐κB) signalling pathway was also dramatically inhibited by medium‐ and high‐dose TAC administration at 12, 24 and 48 h of LPS treatment mice. In addition, TAC reversed LPS‐induced podocyte cytoskeletal injury and podocyte migratory capability. Our findings indicate that TAC has protective effects against LPS‐induced AKI by inhibiting TLR4/MyD88/NF‐κB signalling pathway and podocyte dysfunction, providing another potential therapeutic effects for the LPS‐induced SA‐AKI.     

Genome‐wide analysis reveals demographic and life‐history patterns associated with habitat modification in landlocked,deep‐spawning sockeye salmon (Oncorhynchus nerka)

Human‐mediated habitat fragmentation in freshwater ecosystems can negatively impact genetic diversity, demography, and life history of native biota, while disrupting the behavior of species that are dependent on spatial connectivity to complete their life cycles. In the Alouette River system (British Columbia, Canada), dam construction in 1928 impacted passage of anadromous sockeye salmon (Oncorhynchus nerka), with the last records of migrants occurring in the 1930s. Since that time, O. nerka persisted as a resident population in Alouette Reservoir until experimental water releases beginning in 2005 created conditions for migration; two years later, returning migrants were observed for the first time in ~70 years, raising important basic and applied questions regarding life‐history variation and population structure in this system. Here, we investigated the genetic distinctiveness and population history of Alouette Reservoir O. nerka using genome‐wide SNP data (n = 7,709 loci) collected for resident and migrant individuals, as well as for neighboring anadromous sockeye salmon and resident kokanee populations within the Fraser River drainage (n = 312 individuals). Bayesian clustering and principal components analyses based on neutral loci revealed five distinct clusters, largely associated with geography, and clearly demonstrated that Alouette Reservoir resident and migrant individuals are genetically distinct from other O. nerka populations in the Fraser River drainage. At a finer level, there was no clear evidence for differentiation between Alouette Reservoir residents and migrants; although we detected eight high‐confidence outlier loci, they all mapped to sex chromosomes suggesting that differences were likely due to uneven sex ratios rather than life history. Taken together, these data suggest that contemporary Alouette Reservoir O. nerka represents a landlocked sockeye salmon population, constituting the first reported instance of deep‐water spawning behavior associated with this life‐history form. This finding punctuates the need for reassessment of conservation status and supports ongoing fisheries management activities in Alouette Reservoir.     

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.     

Classifying behavior from short‐interval biologging data: An example with GPS tracking of birds

1. Recent advances in digital data collection have spurred accumulation of immense quantities of data that have potential to lead to remarkable ecological insight, but that also present analytic challenges. In the case of biologging data from birds, common analytical approaches to classifying movement behaviors are largely inappropriate for these massive data sets.
2. We apply a framework for using K‐means clustering to classify bird behavior using points from short time interval GPS tracks. K‐means clustering is a well‐known and computationally efficient statistical tool that has been used in animal movement studies primarily for clustering segments of consecutive points. To illustrate the utility of our approach, we apply K‐means clustering to six focal variables derived from GPS data collected at 1–11 s intervals from free‐flying bald eagles (Haliaeetus leucocephalus) throughout the state of Iowa, USA. We illustrate how these data can be used to identify behaviors and life‐stage‐ and age‐related variation in behavior.
3. After filtering for data quality, the K‐means algorithm identified four clusters in >2 million GPS telemetry data points. These four clusters corresponded to three movement states: ascending, flapping, and gliding flight; and one non‐moving state: perching. Mapping these states illustrated how they corresponded tightly to expectations derived from natural history observations; for example, long periods of ascending flight were often followed by long gliding descents, birds alternated between flapping and gliding flight.
4. The K‐means clustering approach we applied is both an efficient and effective mechanism to classify and interpret short‐interval biologging data to understand movement behaviors. Furthermore, because it can apply to an abundance of very short, irregular, and high‐dimensional movement data, it provides insight into small‐scale variation in behavior that would not be possible with many other analytical approaches.

     

The alternative serotonin transporter promoter P2 impacts gene function in females with irritable bowel syndrome

Irritable bowel syndrome (IBS) is a gut‐brain disorder in which symptoms are shaped by serotonin acting centrally and peripherally. The serotonin transporter gene SLC6A4 has been implicated in IBS pathophysiology, but the underlying genetic mechanisms remain unclear. We sequenced the alternative P2 promoter driving intestinal SLC6A4 expression and identified single nucleotide polymorphisms (SNPs) that were associated with IBS in a discovery sample. Identified SNPs built different haplotypes, and the tagging SNP rs2020938 seems to associate with constipation‐predominant IBS (IBS‐C) in females. rs2020938 validation was performed in 1978 additional IBS patients and 6,038 controls from eight countries. Meta‐analysis on data from 2,175 IBS patients and 6,128 controls confirmed the association with female IBS‐C. Expression analyses revealed that the P2 promoter drives SLC6A4 expression primarily in the small intestine. Gene reporter assays showed a functional impact of SNPs in the P2 region. In silico analysis of the polymorphic promoter indicated differential expression regulation. Further follow‐up revealed that the major allele of the tagging SNP rs2020938 correlates with differential SLC6A4 expression in the jejunum and with stool consistency, indicating functional relevance. Our data consolidate rs2020938 as a functional SNP associated with IBS‐C risk in females, underlining the relevance of SLC6A4 in IBS pathogenesis.