Advantage of ostarine over raloxifene and their combined treatments for muscle of estrogen-deficient rats

PURPOSE

Selective androgen (ostarine, OST) and estrogen (raloxifene, RAL) receptor modulators with improved tissue selectivity have been developed as alternatives to hormone replacement therapy. We investigated the combined effects of OST and RAL on muscle tissue in an estrogen-deficient rat model of postmenopausal conditions.

METHODS

Three-month-old Sprague Dawley rats were divided into groups: (1) untreated non-ovariectomized rats (Non-OVX), (2) untreated ovariectomized rats (OVX), (3) OVX rats treated with OST, (4) OVX rats treated with RAL, (5) OVX rats treated with OST and RAL. Both compounds were administered in the diet. The average dose received was 0.6 ± 0.1 mg for OST and 11.1 ± 1.2 mg for RAL per kg body weight/day. After thirteen weeks, rat activity, muscle weight, structure, gene expression, and serum markers were analyzed.

RESULTS

OST increased muscle weight, capillary ratio, insulin-like growth factor 1 (Igf-1) expression, serum phosphorus, uterine weight. RAL decreased muscle weight, capillary ratio, food intake, serum calcium and increased Igf-1 and Myostatin expression, serum follicle stimulating hormone (FSH). OST + RAL increased muscle nucleus ratio, uterine weight, serum phosphorus, FSH and luteinizing hormone and decreased body and muscle weight, serum calcium. Neither treatment changed muscle fiber size. OVX increased body and muscle weight, decreased uterine weight, serum calcium and magnesium.

CONCLUSION

OST had beneficial effects on muscle in OVX rats. Side effects of OST on the uterus and serum electrolytes should be considered before using it for therapeutic purposes. RAL and RAL + OST had less effect on muscle and showed endocrinological side effects on pituitary-gonadal axis.

Therapeutic Consequences of Targeting the IGF-1/PI3K/AKT/FOXO3 Axis in Sarcopenia: A Narrative Review

The high prevalence of sarcopenia in an aging population has an underestimated impact on quality of life by increasing the risk of falls and subsequent hospitalization. Unfortunately, the application of the major established key therapeutic-physical activity-is challenging in the immobile and injured sarcopenic patient. Consequently, novel therapeutic directions are needed. The transcription factor Forkhead-Box-Protein O3 (FOXO3) may be an option, as it and its targets have been observed to be more highly expressed in sarcopenic muscle. In such catabolic situations, Foxo3 induces the expression of two muscle specific ubiquitin ligases (Atrogin-1 and Murf-1) via the PI3K/AKT pathway. In this review, we particularly evaluate the potential of Foxo3-targeted gene therapy. Foxo3 knockdown has been shown to lead to increased muscle cross sectional area, through both the AKT-dependent and -independent pathways and the reduced impact on the two major downstream targets Atrogin-1 and Murf-1. Moreover, a Foxo3 reduction suppresses apoptosis, activates satellite cells, and initiates their differentiation into muscle cells. While this indicates a critical role in muscle regeneration, this mechanism might exhaust the stem cell pool, limiting its clinical applicability. As systemic Foxo3 knockdown has also been associated with risks of inflammation and cancer progression, a muscle-specific approach would be necessary. In this review, we summarize the current knowledge on Foxo3 and conceptualize a specific and targeted therapy that may circumvent the drawbacks of systemic Foxo3 knockdown. This approach presumably would limit the side effects and enable an activity-independent positive impact on skeletal muscle.

Effects of ligandrol as a selective androgen receptor modulator in a rat model for osteoporosis

INTRODUCTION

The selective androgen receptor modulator ligandrol (LGD-4033 or VK5211) has been shown to improve muscle tissue. In the present study, the effect of ligandrol on bone tissue was investigated in ovariectomized rat model.

MATERIALS AND METHODS

Three-month-old Sprague Dawley rats were either ovariectomized (OVX, n = 60) or left intact (NON-OVX, n = 15). After 9 weeks, OVX rats were divided into four groups: untreated OVX (n = 15) group and three OVX groups (each of 15 rats) treated with ligandrol orally at doses of 0.03, 0.3, or 3 mg/kg body weight. After five weeks, lumbar vertebral bodies (L), tibiae, and femora were examined using micro-computed tomographical, biomechanical, ashing, and gene expression analyses.

RESULTS

In the 3-mg ligandrol group, bone structural properties were improved (trabecular number: 38 ± 8 vs. 35 ± 7 (femur), 26 ± 7 vs. 22 ± 6 (L), 12 ± 5 vs. 6 ± 3 (tibia) and serum phosphorus levels (1.81 ± 0.17 vs.1.41 ± 0.17 mmol/l), uterus (0.43 ± 0.04 vs. 0.11 ± 0.02 g), and heart (1.13 ± 0.11 vs. 1.01 ± 0.08 g) weights were increased compared to the OVX group. Biomechanical parameters were not changed. Low and medium doses did not affect bone tissue and had fewer side effects. Body weight and food intake were not affected by ligandrol; OVX led to an increase in these parameters and worsened all bone parameters.

CONCLUSION

Ligandrol at high dose showed a subtle anabolic effect on structural properties without any improvement in biomechanical properties of osteoporotic bones. Considering side effects of ligandrol at this dose, its further investigation for the therapy of postmenopausal osteoporosis should be reevaluated.

Treatment of osteoporosis using a selective androgen receptor modulator ostarine in an orchiectomized rat model

PURPOSE

The selective androgen receptor modulator ostarine has been shown to have advantageous effects on skeletal tissue properties, reducing muscle wasting and improving physical function in males. However, data on effects in male osteoporosis remain limited. In this study, the effects of ostarine on osteoporotic bone were evaluated in a rat model of male osteoporosis and compared with those of testosterone treatments.

METHODS

Eight-month-old male Sprague-Dawley rats were either non-orchiectomized to serve as a healthy control (Non-Orx, Group 1) or orchiectomized (Orx, Groups 2-6) and then grouped (n = 15/group): (1) Non-Orx, (2) Orx, (3) Ostarine Therapy, (4) Testosterone Therapy, (5) Ostarine Prophylaxis and (6) Testosterone Prophylaxis. Prophylaxis treatments started directly after orchiectomy and continued for 18 weeks, whereas Therapy treatments were initiated 12 weeks after Orx. Ostarine and Testosterone were applied orally at daily doses of 0.4 and 50 mg/kg body weight, respectively. The lumbar vertebral bodies and femora were analyzed using biomechanical, micro-CT, ashing, and gene expression analyses.

RESULTS

Ostarine Prophylaxis showed positive effects in preventing osteoporotic changes in cortical and trabecular bone (femoral trabecular density: 26.01 ± 9.1% vs. 20.75 ± 1.2% in Orx and in L4: 16.3 ± 7.3% vs 11.8 ± 2.9% in Orx); biomechanical parameters were not affected; prostate weight was increased (0.62 ± 0.13 g vs 0.18 ± 0.07 g in Orx). Ostarine Therapy increased solely the cortical density of the femur (1.25 ± 0.03 g/cm3 vs. 1.18 ± 0.04 g/cm3 in Orx); other bone parameters remained unaffected. Testosteron Prophylaxis positively influenced cortical density in femur (1.24 ± 0.05 g/cm3 vs. 1.18 ± 0.04 g/cm3 in Orx); Test. Therapy did not change any bony parameters.

CONCLUSION

Ostarine Prophylaxis could be further investigated as a preventative treatment for male osteoporosis, but an androgenic effect on the prostate should be taken into consideration, and combination therapies with other anti-osteoporosis agents could be considered.

Foxo3 Knockdown Mediates Decline of Myod1 and Myog Reducing Myoblast Conversion to Myotubes

Sarcopenia has a high prevalence among the aging population. Sarcopenia is of tremendous socioeconomic importance because it can lead to falls and hospitalization, subsequently increasing healthcare costs while limiting quality of life. In sarcopenic muscle fibers, the E3 ubiquitin ligase F-Box Protein 32 (Fbxo32) is expressed at substantially higher levels, driving ubiquitin-proteasomal muscle protein degradation. As one of the key regulators of muscular equilibrium, the transcription factor Forkhead Box O3 (FOXO3) can increase the expression of Fbxo32, making it a possible target for the regulation of this detrimental pathway. To test this hypothesis, murine C2C12 myoblasts were transduced with AAVs carrying a plasmid for four specific siRNAs against Foxo3. Successfully transduced myoblasts were selected via FACS cell sorting to establish single clone cell lines. Sorted myoblasts were further differentiated into myotubes and stained for myosin heavy chain (MHC) by immunofluorescence. The resulting area was calculated. Myotube contractions were induced by electrical stimulation and quantified. We found an increased Foxo3 expression in satellite cells in human skeletal muscle and an age-related increase in Foxo3 expression in older mice in silico. We established an in vitro AAV-mediated FOXO3 knockdown on protein level. Surprisingly, the myotubes with FOXO3 knockdown displayed a smaller myotube size and a lower number of nuclei per myotube compared to the control myotubes (AAV-transduced with a functionless control plasmid). During differentiation, a lower level of FOXO3 reduced the expression Fbxo32 within the first three days. Moreover, the expression of Myod1 and Myog via ATM and Tp53 was reduced. Functionally, the Foxo3 knockdown myotubes showed a higher contraction duration and time to peak. Early Foxo3 knockdown seems to terminate the initiation of differentiation due to lack of Myod1 expression, and mediates the inhibition of Myog. Subsequently, the myotube size is reduced and the excitability to electrical stimulation is altered.

Changes of the subchondral bone microchannel network in early osteoarthritis

OBJECTIVE

We have identified a 3D network of subchondral microchannels that connects the deep zone of cartilage to the bone marrow (i.e., cartilage-bone marrow microchannel connectors; CMMC). However, the pathological significance of CMMC is largely unknown. Here, we quantitatively evaluated how the CMMC microarchitecture is related to cartilage condition, as well as regional differences in early idiopathic osteoarthritis (OA).

METHODS

Two groups of cadaveric female human femoral heads (intact cartilage vs early cartilage lesions) were identified, and a biopsy-based high-resolution micro-CT imaging was employed. Subchondral bone (SB) thickness, CMMC number, maximum and minimum CMMC size, and the CMMC morphology were quantified and compared between the two groups. The effect of joint's region and cartilage condition was examined on each dependent variable.

RESULTS

The CMMC number and morphology were affected by region of the joint, but not by cartilage condition. On the other hand, the minimum and maximum CMMC size was changed by both the location on the joint, as well as the cartilage condition. The smallest CMMC were consistently detected at the load-bearing region (LBR) of the joint. Compared to non-pathological subjects, the size of the microchannels was enlarged in early OA, most noticeably at the non-load-bearing region (NLBR) and the peripheral rim (PR) of the femoral head. Furthermore, subchondral bone thinning was observed in early OA as a localized occurrence linked with areas of partial chondral defect.

CONCLUSION

Our data point to an enlargement of the SB microchannel network, and a collective structural deterioration of SB in early idiopathic OA.

A combined treatment with selective androgen and estrogen receptor modulators prevents bone loss in orchiectomized rats

PURPOSE

Enobosarm (EN), a selective androgen receptor modulator and raloxifene (RAL), a selective estrogen receptor modulator, have been shown to improve bone tissue in osteoporotic males. The present study evaluated the effects of a combination therapy of EN and RAL on bone properties in orchiectomized rats compared to the respective single treatments.

METHODS

Eight-month-old male Sprague-Dawley rats were either left intact (Non-Orx) or orchiectomized (Orx). The Orx rats were divided into four groups (n = 15 each): 1) Orx, 2) EN treatment (Orx + EN), 3) RAL treatment (Orx + RAL), 4) combined treatment (Orx + EN + RAL). EN and RAL (0.4 mg and 7 mg/kg body weight/day) were applied immediately after Orx with a soy-free pelleted diet for up to 18 weeks. The lumbar spine and femora were examined by micro-CT, biomechanical, histomorphological, ashing, and gene expression analyses.

RESULTS

EN exhibited an anabolic effect on bone, improving some of its parameters in Orx rats, but did not affect biomechanical properties. RAL exhibited antiresorptive activity, maintaining the biomechanical and trabecular parameters of Orx rats at the levels of Non-Orx rats. EN + RAL exerted a stronger effect than the single treatments, improving most of the bone parameters. Liver weight increased after all treatments; the kidney, prostate, and levator ani muscle weights increased after EN and EN + RAL treatments. BW was reduced due to a decreased food intake in the Orx + RAL group and due a reduced visceral fat weight in the Orx + EN + RAL group.

CONCLUSION

The EN + RAL treatment appeared to be promising in preventing male osteoporosis, but given the observed side effects on liver, kidney, and prostate weights, it requires further investigation.

Investigating the Microchannel Architectures Inside the Subchondral Bone in Relation to Estimated Hip Reaction Forces on the Human Femoral Head

The interplay between articular cartilage (AC) and subchondral bone (SB) plays a pivotal role in cartilage homeostasis and functionality. As direct connective pathways between the two are poorly understood, we examined the location-dependent characteristics of the 3D microchannel network within the SB that connects the basal cartilage layer to the bone marrow (i.e. cartilage-bone marrow microchannel connectors; CMMC). 43 measuring points were defined on five human cadaveric femoral heads with no signs of osteoarthritis (OA) (age ≤ 60), and cartilage-bone cylinders with diameters of 2.00 mm were extracted for high-resolution scanning (n = 215). The micro-CT data were categorized into three groups (load-bearing region: LBR, n = 60; non-load-bearing region: NLBR, n = 60; and the peripheral rim: PR, n = 95) based on a gait analysis estimation of the joint reaction force (young, healthy cohort with no signs of OA). At the AC-SB interface, the number of CMMC in the LBR was 1.8 times and 2.2 times higher compared to the NLBR, and the PR, respectively. On the other hand, the median Feret size of the CMMC were smallest in the LBR (55.2 µm) and increased in the NLBR (73.5 µm; p = 0.043) and the PR (89.1 µm; p = 0.043). AC thickness was positively associated with SB thickness (Pearson's r = 0.48; p < 1e-13), CMMC number. (r = 0.46; p < 1e-11), and circularity index (r = 0.61; p < 1e-38). In conclusion, our data suggest that regional differences in the microchannel architecture of SB might reflect regional differences in loading.

Predicting the disaster - The role of CRP in acetabular surgery

OBJECTIVES

Acetabular fractures represent a complex surgical challenge. Given the heterogenous fracture pattern, the patient characteristics and spectrum of complications demand individual solutions. Surgical site infections (SSI) threaten osteosynthesis, and early detection of them and treatment remain crucial. What is the value of postoperative C-reactive protein (CRP) in this group of patients as well as its normal course?

DESIGN & METHODS

115 patients with isolated fractures of the acetabulum were retrospectively evaluated. CRP, white blood cell count (WBC) and fracture patterns as well as patient characteristics were assessed for 20 days following operative fixation of the acetabular fracture (n = 71) and in fractures that were managed conservatively (n = 44).

RESULTS

Twelve patients suffered an infectious complication. With a one-phase decay, 70.55% of the variance of postoperative CRP kinetics was predicted. To anticipate maximum CRP as well as an infection, the preoperative CRP represented the best prognostic parameter. To predict an infection, the single variable "peak CRP value above 100 mg/l" resulted in a sensitivity and specificity of 91.67% and 36.21%, respectively. Combining a second peak of CRP with maximum CRP and day 5 CRP value for receiver-operating characteristic (ROC) analysis resulted in 83.3% and 88.1%, respectively.

CONCLUSIONS

Predicting surgical site infections after an acetabular fracture is most predictive when analyzing the maximum overall CRP, the second peak and the CRP after day 5. With a combination of these parameters, a sensitivity and specificity of 83.3% and 88.1% to detect an infection was achieved.

Current State of Bone Adhesives-Necessities and Hurdles

The vision of gluing two bone fragments with biodegradable and biocompatible adhesives remains highly fascinating and attractive to orthopedic surgeons. Possibly shorter operation times, better stabilization, lower infection rates, and unnecessary removal make this approach very appealing. After 30 years of research in this field, the first adhesive systems are now appearing in scientific reports that may fulfill the comprehensive requirements of bioadhesives for bone. For a successful introduction into clinical application, special requirements of the musculoskeletal system, challenges in the production of a bone adhesive, as well as regulatory hurdles still need to be overcome. In this article, we will give an overview of existing synthetic polymers, biomimetic, and bio-based adhesive approaches, review the regulatory hurdles they face, and discuss perspectives of how bone adhesives could be efficiently introduced into clinical application, including legal regulations.

Radiation-induced sensitivity of tissue-resident mesenchymal stem cells in the head and neck region

BACKGROUND

Tissue-resident mesenchymal stem cells (MSCs) possess the ability to migrate to areas of inflammation and promote the regeneration of damaged tissue. However, it remains unclear how radiation influences this capacity of MSC in the head and neck region.

METHODS

Two types of MSCs of the head and neck region (mucosa [mMSC] and parotid gland [pMSC]) were isolated, cultured and exposed to single radiation dosages of 2 Gy/day up to 10 days. Effects on morphology, colony forming ability, apoptosis, chemokine receptor expression, cytokine secretion, and cell migration were analyzed.

RESULTS

Although MSC preserved MSC-specific regenerative abilities and immunomodulatory properties following irradiation in our in vitro model, we found a deleterious impact on colony forming ability, especially in pMSC.

CONCLUSIONS

MSC exhibited robustness and activation upon radiation for the support of tissue regeneration, but lost their potential to replicate, thus possibly leading to depletion of the local MSC-pool after irradiation.

Effects of RANKL Knockdown by Virus-like Particle-Mediated RNAi in a Rat Model of Osteoporosis

Rebalancing of the RANKL/OPG system seems to be an effective treatment strategy in postmenopausal osteoporosis. Here, we evaluate the knockdown of RANKL by in-vivo-delivered siRNA in a rat model of osteoporosis. Virus-like-particles (VLPs) derived from polyoma JC virus were used for delivering RANKL siRNA in ovariectomized (OVX) rats. 48 rats were ovariectomized and treated with either 17β-estradiol (E2), VLPs containing RANKL siRNA (siRANKL), or VLPs containing non-cognate siRNA (siCtrl). All OVX groups were subdivided into the prophylaxis group (PG) and the therapy group (TG). The PG received treatment directly after being OVX for 10 weeks. The TG received treatment 5 weeks after being OVX for 5 weeks. Rats were sacrificed 10 weeks after being OVX. Bone and blood samples were analyzed. E2 and siRANKL showed a significant knockdown of RANKL mRNA. A protein knockdown was observed with E2 and siRANKL in the TG but not in the PG. No distinct improvements in biomechanical and morphological properties of the bones were observed after siRANKL treatment. In the PG, E2 protected the bone structure. We demonstrated successful mRNA and protein knockdown by VLP-mediated RNAi in vivo. Knockdown of membranous RANKL did not result in significant improvements of bone properties in this model of early-stage postmenopausal osteoporosis.

Effect of Strontium Ranelate on the Muscle and Vertebrae of Ovariectomized Rats

Osteoporosis is often accompanied by sarcopenia. The effect of strontium ranelate (SR) on muscle tissue has not been investigated sufficiently. In this study, the effect of different SR treatments on muscle was studied. Additionally, the lumbar vertebrae were analyzed. Three-month-old female rats were divided into five groups (n = 12): Group 1: untreated (NON-OVX); Group 2: ovariectomized and left untreated (OVX); Group 3: SR after OVX until the study ended (13 weeks, SR prophylaxis and therapy = pr+th); Group 4: OVX and SR for 8 weeks (SR prophylaxis = pr); Group 5: SR for 5 weeks from the 8 week after OVX (SR therapy = SR th). SR was applied in food (630 mg/kg body weight). The size of muscle fibers, capillary density, metabolic enzymes, and mRNA expression were assessed in soleus, gastrocnemius, and longissimus muscles. The vertebral bodies underwent micro-CT, biomechanical, and ashing analyses. In general, SR did not alter the muscle histological parameters. The changes in fiber size and capillary ratio were related to the body weight. Myostatin mRNA was decreased in Sr pr+th; protein expression was not changed. SR th led to increase in mRNA expression of vascular endothelial growth factor (Vegf-B). In lumbar spine, SR pr+th enhanced biomechanical properties, bone mineral density, trabecular area, density, and thickness and cortical density. The reduced calcium/phosphate ratio in the SR pr+th group indicates the replacement of calcium by strontium ions. SR has no adverse effects on muscle tissue and it shows a favorable time-dependent effect on vertebrae. A functional analysis of muscles could verify these findings.

The Impact of the CD9 Tetraspanin on Lentivirus Infectivity and Exosome Secretion

Efficient transduction tools are a hallmark for both research and therapy development. Here, we introduce new insights into the generation of lentiviral vectors with improved performance by utilizing producer cells with increased production rates of extracellular vesicles through CD9 overexpression. Most human cells secrete small vesicles from their surface (microvesicles) or intraluminal endosome-derived membranes (exosomes). In particular, enhanced levels of the tetraspanin CD9 result in significantly increased numbers of extracellular vesicles with exosome-like features that were secreted from four different human cell lines. Intriguingly, exosomes and their biogenesis route display similarities to lentivirus and we examined the impact of CD9 expression on release and infectivity of recombinant lentiviral vectors. Although the titers of released viral particles were not increased upon production in high CD9 cells, we observed improved performance in terms of both speed and efficiency of lentiviral gene delivery into numerous human cell lines, including HEK293, HeLa, SH-SY5Y, as well as B and T lymphocytes. Here, we demonstrate that enhanced CD9 enables lentiviral transduction in the absence of any pseudotyping viral glycoprotein or fusogenic molecule. Our findings indicate an important role of CD9 for lentiviral vector and exosome biogenesis and point out a remarkable function of this tetraspanin in membrane fusion, viral infectivity, and exosome-mediated horizontal information transfer.

In Vivo siRNA Delivery Using JC Virus-like Particles Decreases the Expression of RANKL in Rats

Bone remodeling requires a precise balance between formation and resorption. This complex process involves numerous factors that orchestrate a multitude of biochemical events. Among these factors are hormones, growth factors, vitamins, cytokines, and, most notably, osteoprotegerin (OPG) and the receptor activator for nuclear factor-kappaB ligand (RANKL). Inflammatory cytokines play a major role in shifting the RANKL/OPG balance toward excessive RANKL, resulting in osteoclastogenesis, which in turn initiates bone resorption, which is frequently associated with osteoporosis. Rebalancing RANKL/OPG levels may be achieved through either upregulation of OPG or through transient silencing of RANKL by means of RNA interference. Here, we describe the utilization of a viral capsid-based delivery system for in vivo and in vitro RNAi using synthetic small interfering RNA (siRNA) molecules in rat osteoblasts. Polyoma JC virus-derived virus-like particles are capable of delivering siRNAs to target RANKL in osteoblast cells both in vitro and in a rat in vivo system. Expression levels were monitored using quantitative real-time polymerase reaction and enzyme-linked immunosorbent assay after single and repeated injections over a 14-day period. Our data indicate that this is an efficient and safe route for in vivo delivery of gene modulatory tools to study important molecular factors in a rat osteoporosis model.

Mitogenomic phylogeny of the common long-tailed macaque (Macaca fascicularis fascicularis)

Background

Long-tailed macaques (Macaca fascicularis) are an important model species in biomedical research and reliable knowledge about their evolutionary history is essential for biomedical inferences. Ten subspecies have been recognized, of which most are restricted to small islands of Southeast Asia. In contrast, the common long-tailed macaque (M. f. fascicularis) is distributed over large parts of the Southeast Asian mainland and the Sundaland region. To shed more light on the phylogeny of M. f. fascicularis, we sequenced complete mitochondrial (mtDNA) genomes of 40 individuals from all over the taxon’s range, either by classical PCR-amplification and Sanger sequencing or by DNA-capture and high-throughput sequencing.

Results

Both laboratory approaches yielded complete mtDNA genomes from M. f. fascicularis with high accuracy and/or coverage. According to our phylogenetic reconstructions, M. f. fascicularis initially diverged into two clades 1.70 million years ago (Ma), with one including haplotypes from mainland Southeast Asia, the Malay Peninsula and North Sumatra (Clade A) and the other, haplotypes from the islands of Bangka, Java, Borneo, Timor, and the Philippines (Clade B). The three geographical populations of Clade A appear as paraphyletic groups, while local populations of Clade B form monophyletic clades with the exception of a Philippine individual which is nested within the Borneo clade. Further, in Clade B the branching pattern among main clades/lineages remains largely unresolved, most likely due to their relatively rapid diversification 0.93-0.84 Ma.

Conclusions

Both laboratory methods have proven to be powerful to generate complete mtDNA genome data with similarly high accuracy, with the DNA-capture and high-throughput sequencing approach as the most promising and only practical option to obtain such data from highly degraded DNA, in time and with relatively low costs. The application of complete mtDNA genomes yields new insights into the evolutionary history of M. f. fascicularis by providing a more robust phylogeny and more reliable divergence age estimations than earlier studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1437-0) contains supplementary material, which is available to authorized users.