Genomic Evidence for Island Population Conversion Resolves Conflicting Theories of Polar Bear Evolution

Despite extensive genetic analysis, the evolutionary relationship between polar bears (Ursus maritimus) and brown bears (U. arctos) remains unclear. The two most recent comprehensive reports indicate a recent divergence with little subsequent admixture or a much more ancient divergence followed by extensive admixture. At the center of this controversy are the Alaskan ABC Islands brown bears that show evidence of shared ancestry with polar bears. We present an analysis of genome-wide sequence data for seven polar bears, one ABC Islands brown bear, one mainland Alaskan brown bear, and a black bear (U. americanus), plus recently published datasets from other bears. Surprisingly, we find clear evidence for gene flow from polar bears into ABC Islands brown bears but no evidence of gene flow from brown bears into polar bears. Importantly, while polar bears contributed <1% of the autosomal genome of the ABC Islands brown bear, they contributed 6.5% of the X chromosome. The magnitude of sex-biased polar bear ancestry and the clear direction of gene flow suggest a model wherein the enigmatic ABC Island brown bears are the descendants of a polar bear population that was gradually converted into brown bears via male-dominated brown bear admixture. We present a model that reconciles heretofore conflicting genetic observations. We posit that the enigmatic ABC Islands brown bears derive from a population of polar bears likely stranded by the receding ice at the end of the last glacial period. Since then, male brown bear migration onto the island has gradually converted these bears into an admixed population whose phenotype and genotype are principally brown bear, except at mtDNA and X-linked loci. This process of genome erosion and conversion may be a common outcome when climate change or other forces cause a population to become isolated and then overrun by species with which it can hybridize.     

p62/SQSTM1: A Missing Link between Protein Aggregates and the Autophagy Machinery

In eukaryotic cells short-lived proteins are degraded in a specific process by the ubiquitin-proteasome system (UPS), whereas long-lived proteins and damaged organelles are degraded by macroautophagy (hereafter referred to as autophagy). A growing body of evidence now suggests that autophagy is important for clearance of protein aggregates that form in cells as a consequence of ageing, oxidative stress, alterations that elevate the amounts of certain aggregation-prone proteins or expression of aggregating mutant variants of specific proteins. Autophagy is generally considered to be a non-specific, bulk degradation process. However, a recent study suggests that p62/SQSTM1 may link the recognition of polyubiquitinated protein aggregates to the autophagy machinery.1 This protein is able to polymerize via its N-terminal PB1 domain and to recognize polyubiquitin via its C-terminal UBA domain. It can also recruit the autophagosomal protein LC3 and co-localizes with many types of polyubiquitinated protein aggregates.1 Here we discuss possible implications of these findings and raise some questions for further investigation.     

Rational Design of a Fibroblast Growth Factor 21-Based Clinical Candidate,LY2405319

Fibroblast growth factor 21 is a novel hormonal regulator with the potential to treat a broad variety of metabolic abnormalities, such as type 2 diabetes, obesity, hepatic steatosis, and cardiovascular disease. Human recombinant wild type FGF21 (FGF21) has been shown to ameliorate metabolic disorders in rodents and non-human primates. However, development of FGF21 as a drug is challenging and requires re-engineering of its amino acid sequence to improve protein expression and formulation stability. Here we report the design and characterization of a novel FGF21 variant, LY2405319. To enable the development of a potential drug product with a once-daily dosing profile, in a preserved, multi-use formulation, an additional disulfide bond was introduced in FGF21 through Leu118Cys and Ala134Cys mutations. FGF21 was further optimized by deleting the four N-terminal amino acids, His-Pro-Ile-Pro (HPIP), which was subject to proteolytic cleavage. In addition, to eliminate an O-linked glycosylation site in yeast a Ser167Ala mutation was introduced, thus allowing large-scale, homogenous protein production in Pichia pastoris. Altogether re-engineering of FGF21 led to significant improvements in its biopharmaceutical properties. The impact of these changes was assessed in a panel of in vitro and in vivo assays, which confirmed that biological properties of LY2405319 were essentially identical to FGF21. Specifically, subcutaneous administration of LY2405319 in ob/ob and diet-induced obese (DIO) mice over 7–14 days resulted in a 25–50% lowering of plasma glucose coupled with a 10–30% reduction in body weight. Thus, LY2405319 exhibited all the biopharmaceutical and biological properties required for initiation of a clinical program designed to test the hypothesis that administration of exogenous FGF21 would result in effects on disease-related metabolic parameters in humans.     

Moving Domain Computational Fluid Dynamics to Interface with an Embryonic Model of Cardiac Morphogenesis

Peristaltic contraction of the embryonic heart tube produces time- and spatial-varying wall shear stress (WSS) and pressure gradients (∇P) across the atrioventricular (AV) canal. Zebrafish (Danio rerio) are a genetically tractable system to investigate cardiac morphogenesis. The use of Tg(fli1a:EGFP)^y1 transgenic embryos allowed for delineation and two-dimensional reconstruction of the endocardium. This time-varying wall motion was then prescribed in a two-dimensional moving domain computational fluid dynamics (CFD) model, providing new insights into spatial and temporal variations in WSS and ∇P during cardiac development. The CFD simulations were validated with particle image velocimetry (PIV) across the atrioventricular (AV) canal, revealing an increase in both velocities and heart rates, but a decrease in the duration of atrial systole from early to later stages. At 20-30 hours post fertilization (hpf), simulation results revealed bidirectional WSS across the AV canal in the heart tube in response to peristaltic motion of the wall. At 40-50 hpf, the tube structure undergoes cardiac looping, accompanied by a nearly 3-fold increase in WSS magnitude. At 110-120 hpf, distinct AV valve, atrium, ventricle, and bulbus arteriosus form, accompanied by incremental increases in both WSS magnitude and ∇P, but a decrease in bi-directional flow. Laminar flow develops across the AV canal at 20-30 hpf, and persists at 110-120 hpf. Reynolds numbers at the AV canal increase from 0.07±0.03 at 20-30 hpf to 0.23±0.07 at 110-120 hpf (p< 0.05, n=6), whereas Womersley numbers remain relatively unchanged from 0.11 to 0.13. Our moving domain simulations highlights hemodynamic changes in relation to cardiac morphogenesis; thereby, providing a 2-D quantitative approach to complement imaging analysis.     

Genetic variants associated with breast cancer risk for Ashkenazi Jewish women with strong family histories but no identifiable BRCA1/2 mutation

The ability to establish genetic risk models is critical for early identification and optimal treatment of breast cancer. For such a model to gain clinical utility, more variants must be identified beyond those discovered in previous genome-wide association studies (GWAS). This is especially true for women at high risk because of family history, but without BRCA1/2 mutations. This study incorporates three datasets in a GWAS analysis of women with Ashkenazi Jewish (AJ) homogeneous ancestry. Two independent discovery cohorts comprised 239 and 238 AJ women with invasive breast cancer or preinvasive ductal carcinoma in situ and strong family histories of breast cancer, but lacking the three BRCA1/2 founder mutations, along with 294 and 230 AJ controls, respectively. An independent, third cohort of 203 AJ cases with familial breast cancer history and 263 healthy controls of AJ women was used for validation. A total of 19 SNPs were identified as associated with familial breast cancer risk in AJ women. Among these SNPs, 13 were identified from a panel of 109 discovery SNPs, including an FGFR2 haplotype. In addition, six previously identified breast cancer GWAS SNPs were confirmed in this population. Seven of the 19 markers were significant in a multivariate predictive model of familial breast cancer in AJ women, three novel SNPs [rs17663555(5q13.2), rs566164(6q21), and rs11075884(16q22.2)], the FGFR2 haplotype, and three previously published SNPs [rs13387042(2q35), rs2046210(ESR1), and rs3112612(TOX3)], yielding moderate predictive power with an area under the curve (AUC) of the ROC (receiver-operator characteristic curve) of 0.74. Population-specific genetic variants in addition to variants shared with populations of European ancestry may improve breast cancer risk prediction among AJ women from high-risk families without founder BRCA1/2 mutations.     

Root dynamics of Carex stricta-dominated tussock meadows

Background and aims

The roots of tussock-forming plants contribute to the formation of microtopographic features in many ecosystems, but the dynamics of such roots are poorly understood. We examined the spatial heterogeneity of tussock fine root dynamics to investigate allocation patterns and the role of root productivity in the persistence of tussock structures.

Methods

We compared the spatial variability of fine root (<1 mm, 1–2 mm) density, biomass, % live, allocation, turnover rate (using bomb 14C), and productivity of four Carex stricta Lam.-dominated tussock meadows in the upper Midwest, USA (3 reference, 1 restored site).

Results

Relative to underlying microsites, tussocks were warm, dry, and high in root density, productivity, % live biomass, and turnover. Root productivity averaged 649 g?m^?2 yr^?1 (±208) in reference sites, comprised 57 % (±10) of total net production, and was concentrated in tussocks (70 %?±?4). Root turnover rate averaged 0.63 yr^?1 (±0.08), but tussocks had ~50 % faster root turnover than the underlying soil, and <1 mm roots turned over ~40 % faster than 1–2 mm roots.

Conclusions

Our detailed analysis of the spatial heterogeneity of tussock root dynamics suggests that high allocation and elevated turnover of tussock roots facilitates organic matter accumulation and tussock persistence over time.     

Effects of small interfering RNA-mediated hepatic glucagon receptor inhibition on lipid metabolism in db/db mice

Hepatic glucose overproduction is a major characteristic of type 2 diabetes. Because glucagon is a key regulator for glucose homeostasis, antagonizing the glucagon receptor (GCGR) is a possible therapeutic strategy for the treatment of diabetes mellitus. To study the effect of hepatic GCGR inhibition on the regulation of lipid metabolism, we generated siRNA-mediated GCGR knockdown (si-GCGR) in the db/db mouse. The hepatic knockdown of GCGR markedly reduced plasma glucose levels; however, total plasma cholesterol was increased. The detailed lipid analysis showed an increase in the LDL fraction, and no change in VLDL HDL fractions. Further studies showed that the increase in LDL was the result of over-expression of hepatic lipogenic genes and elevated de novo lipid synthesis. Inhibition of hepatic glucagon signaling via siRNA-mediated GCGR knockdown had an effect on both glucose and lipid metabolism in db/db mice.     

A TRPC1 Protein-dependent Pathway Regulates Osteoclast Formation and Function

Ca²⁺ signaling is essential for bone homeostasis and skeletal development. Here, we show that the transient receptor potential canonical 1 (TRPC1) channel and the inhibitor of MyoD family, I-mfa, function antagonistically in the regulation of osteoclastogenesis. I-mfa null mice have an osteopenic phenotype characterized by increased osteoclast numbers and surface, which are normalized in mice lacking both Trpc1 and I-mfa. In vitro differentiation of pre-osteoclasts derived from I-mfa-deficient mice leads to an increased number of mature osteoclasts and higher bone resorption per osteoclast. These parameters return to normal levels in osteoclasts derived from double mutant mice. Consistently, whole cell currents activated in response to the depletion of intracellular Ca²⁺ stores are larger in pre-osteoclasts derived from I-mfa knock-out mice compared with currents in wild type mice and normalized in cells derived from double mutant mice, suggesting a cell-autonomous effect of I-mfa on TRPC1 in these cells. A new splice variant of TRPC1 (TRPC1ϵ) was identified in early pre-osteoclasts. Heterologous expression of TRPC1ϵ in HEK293 cells revealed that it is unique among all known TRPC1 isoforms in its ability to amplify the activity of the Ca²⁺ release-activated Ca²⁺ (CRAC) channel, mediating store-operated currents. TRPC1ϵ physically interacts with Orai1, the pore-forming subunit of the CRAC channel, and I-mfa is recruited to the TRPC1ϵ-Orai1 complex through TRPC1ϵ suppressing CRAC channel activity. We propose that the positive and negative modulation of the CRAC channel by TRPC1ϵ and I-mfa, respectively, fine-tunes the dynamic range of the CRAC channel regulating osteoclastogenesis.