Dramatic genome-wide reprogramming of mRNA in hypometabolic muscle

In response to seasonal droughts, the green striped burrowing frog Cyclorana alboguttata enters a reversible hypometabolic state called aestivation where heart rate and oxygen consumption can be reduced despite warm (>25C°) ambient temperatures. With a view to understanding molecular mechanisms we profiled aestivating versus control gastrocnemius muscle using mRNA sequencing. This indicated an extensive metabolic reprogramming, with nearly a quarter of the entire transcriptome (3996 of 16,960 mRNA) exhibiting a nominal >2-fold change. Consistent with a physiological adaptation to spare carbohydrate reserves, carbohydrate catabolism was systemically downregulated. A 630-fold downregulation of ENO3 encoding the enolase enzyme was most striking. The 590 frog orthologs of mRNA encoding the mitoproteome were, viewed as a population, significantly downregulated during aestivation, although not to the same extent as mRNA encoding carbohydrate catabolism. Prominent examples include members of the TCA cycle (IDH2), electron transport chain (NDUFA6), the ATP synthase complex (ATP5F1B) and ADP/ATP intracellular transport (SLC25A4). Moreover, mRNA derived from the mt genome itself (e.g. mt-ND1) were also downregulated. Most prominent among the upregulated mRNA are those encoding aspects of regulated proteolysis including the proteosome (e.g. PSME4L), peptidases (USP25), atrogins (FBXO32) and ubiquitination (VCP). Finally, we note the ∼5-fold upregulation of the mRNA EIFG3 that encodes part of the EIF4F complex. This possesses global control of protein synthesis. Given protein synthesis is repressed in aestivating frogs this indicates the skeletal musculature is poised for accelerated translation of mRNA upon emergence, supporting a strategy to rapidly restore function when the summer rains come.

Differential voluntary feed intake and whole transcriptome profiling in the hypothalamus of young sheep offered CP and phosphorus-deficient diets

A reduction in voluntary feed intake is observed in ruminants consuming nutrient-deficient diets, such as those with a low CP or P content, and has been attributed to active metabolic regulation, rather than a physical constraint. The hypothalamus is the key integrator of feed intake regulation in mammals. The objectives of this experiment were to (1) establish a model of metabolic feed intake regulation in ruminants consuming diets of variable CP and P content, and (2) determine key biochemical pathways and influential points of regulation within the hypothalamus. Merino wethers [n = 40; 23.7 ± 1.4 kg liveweight (mean ± SD)] were fed one of five dietary treatments (n = 8/treatment) for 63 days in individual pens. The treatments included targeted combinations of high (H) and low (L) CP (110 and 55 g/kg DM) and high and low P (2.5 and 0.7 g/kg DM) with 9 MJ metabolisable energy (ME) per kg DM which were fed ad libitum (UMEI; unrestricted ME intake) resulting in four experimental diets (HCP-HP-UMEI, LCP-HP-UMEI, HCP-LP-UMEI and LCP-LP-UMEI). An additional nutritional treatment (HCP-HP-RMEI) restricted intake of the HCP-HP diet to an equivalent ME intake of wethers consuming the LCP-LP-UMEI treatment. Wethers offered the LCP-HP-UMEI, HCP-LP-UMEI and LCP-LP-UMEI treatments consumed 42, 32 and 49% less total DM (P ≤ 0.05), respectively than the HCP-HP-UMEI treatment, and this was not attributable to any physical limitation of the rumen. Plasma concentrations of urea nitrogen and inorganic phosphate indicated that these nutrient deficiencies were successfully established. To assess potential mechanisms, RNA-seq was conducted on samples from the arcuate nucleus (ARC), ventromedial hypothalamus and lateral hypothalamus of the wethers, yielding a total of 301, 8 and 148 differentially expressed genes across all pairwise comparisons, respectively. The expression of NPY, AGRP and CARTPT, known for their regulatory role in mammalian feed intake regulation, had a similar transcriptional response in the ARC of wethers consuming nutrient-deficient treatments and those consuming a ME-restricted treatment, despite these wethers expressing behaviours indicative of satiated and hungry states, respectively. In addition, genes involved with glycolysis (TPI1), the citric acid cycle (CS, OGDH, GLUD1, GOT1) and oxidative phosphorylation (COX5A, ATP5MC1, ATP5F1B, ATP5MC3) were downregulated in the ARC of wethers fed a nutrient deficient (LCP-LP-UMEI) relative to the non-deficient (HCP-HP-UMEI) treatment. In summary, a model for voluntary feed intake restriction was established to determine genome-wide molecular changes in the hypothalamus of young ruminants.

'Spikeopathy': COVID-19 Spike Protein Is Pathogenic, from Both Virus and Vaccine mRNA

The COVID-19 pandemic caused much illness, many deaths, and profound disruption to society. The production of 'safe and effective' vaccines was a key public health target. Sadly, unprecedented high rates of adverse events have overshadowed the benefits. This two-part narrative review presents evidence for the widespread harms of novel product COVID-19 mRNA and adenovectorDNA vaccines and is novel in attempting to provide a thorough overview of harms arising from the new technology in vaccines that relied on human cells producing a foreign antigen that has evidence of pathogenicity. This first paper explores peer-reviewed data counter to the 'safe and effective' narrative attached to these new technologies. Spike protein pathogenicity, termed 'spikeopathy', whether from the SARS-CoV-2 virus or produced by vaccine gene codes, akin to a 'synthetic virus', is increasingly understood in terms of molecular biology and pathophysiology. Pharmacokinetic transfection through body tissues distant from the injection site by lipid-nanoparticles or viral-vector carriers means that 'spikeopathy' can affect many organs. The inflammatory properties of the nanoparticles used to ferry mRNA; N1-methylpseudouridine employed to prolong synthetic mRNA function; the widespread biodistribution of the mRNA and DNA codes and translated spike proteins, and autoimmunity via human production of foreign proteins, contribute to harmful effects. This paper reviews autoimmune, cardiovascular, neurological, potential oncological effects, and autopsy evidence for spikeopathy. With many gene-based therapeutic technologies planned, a re-evaluation is necessary and timely.

Multi-breed genomic predictions and functional variants for fertility of tropical bulls

Worldwide, most beef breeding herds are naturally mated. As such, the ability to identify and select fertile bulls is critically important for both productivity and genetic improvement. Here, we collected ten fertility-related phenotypes for 6,063 bulls from six tropically adapted breeds. Phenotypes were comprised of four bull conformation traits and six traits directly related to the quality of the bull's semen. We also generated high-density DNA genotypes for all the animals. In total, 680,758 single nucleotide polymorphism (SNP) genotypes were analyzed. The genomic correlation of the same trait observed in different breeds was positive for scrotal circumference and sheath score on most breed comparisons, but close to zero for the percentage of normal sperm, suggesting a divergent genetic background for this trait. We confirmed the importance of a breed being present in the reference population to the generation of accurate genomic estimated breeding values (GEBV) in an across-breed validation scenario. Average GEBV accuracies varied from 0.19 to 0.44 when the breed was not included in the reference population. The range improved to 0.28 to 0.59 when the breed was in the reference population. Variants associated with the gene HDAC4, six genes from the spermatogenesis-associated (SPATA) family of proteins, and 29 transcription factors were identified as candidate genes. Collectively these results enable very early in-life selection for bull fertility traits, supporting genetic improvement strategies currently taking place within tropical beef production systems. This study also improves our understanding of the molecular basis of male fertility in mammals.

Seasonal resource selection of an arboreal habitat specialist in a human-dominated landscape: a case study using red panda

Human dominated landscapes provide heterogeneous wildlife habitat. Conservation of habitat specialists, like red pandas Ailurus fulgens, inhabiting such landscapes is challenging. Therefore information on resource use across spatial and temporal scales could enable informed-decision making with better conservation outcomes. We aimed to examine the effect of geo-physical, vegetation, and disturbance variables on fine-scale habitat selection of red pandas in one such landscape. We equipped 10 red pandas with GPS collars in eastern Nepal in 2019 and monitored them for one year. Our analysis was based on a generalized-linear-mixed model. We found the combined effect of geo-physical, vegetation and disturbance variables resulted in differences in resource selection of red pandas and that the degree of response to these variables varied across seasons. Human disturbances, especially road and cattle herding activities, affected habitat utilization throughout the year whereas other variables were important only during restricted periods. For instance, geo-physical variables were influential in the premating and cub-rearing seasons while vegetation variables were important in all seasons other than premating. Red pandas selected steeper slopes with high solar insolation in the premating season while they occupied elevated areas and preferred specific aspects in the cub-rearing season. Furthermore, the utilized areas had tall bamboo in the birthing and cub-rearing seasons while they also preferred diverse tree species and high shrub cover in the latter. Our study demonstrates the significance of season-specific management, suggests the importance of specific types of vegetation during biologically crucial periods, and emphasizes the necessity to minimize disturbances throughout the year.

Space use, interaction and recursion in a solitary specialized herbivore: a red panda case study

Better understanding of ecology is crucial for the success of endangered species conservation programmes. Little information is available on space use, conspecific interactions and recursions by one such species, the red panda Ailurus fulgens. To address this deficiency, we used GPS telemetry to examine their home range, core area, home-range overlap, dynamic interactions and recursive movement, and investigated the effect of sex, age and body mass on these behaviours across seasons. The median annual home range size was 1.41 km2, with nearly a quarter of this range being used as the core area. Sex and reproductive status were the key determinants of space use patterns on a seasonal scale, while body mass and age remained significant correlates for the core area. The home range of males was nearly double that of females, likely because of the polygynous mating system in red pandas. Females avoided overlapping home ranges, while males overlapped home range with up to 4 females, and neighbouring males overlapped nearly half of their ranges. We found rare interactions between males and females outside the mating season. Red pandas showed site fidelity within their territory, with seasonal variation across sex classes. We also observed high individual variation in patterns of both space use and recursive movement. Taken together, these results suggest that differences in biological requirements across seasons determine red panda space use patterns, conspecific interactions and recursion. However, forage availability and quality, climatic factors, disturbances and habitat fragmentation are also likely to influence these behaviours, and these factors need to be investigated.

Movement and dispersal of a habitat specialist in human-dominated landscapes: a case study of the red panda

BACKGROUND

Habitat specialists living in human-dominated landscapes are likely to be affected by habitat fragmentation and human disturbances more than generalists. But there is a paucity of information on their response to such factors. We examined the effect of these factors on movement patterns of red pandas Ailurus fulgens, a habitat and diet specialist that inhabits the eastern Himalaya.

METHODS

We equipped 10 red pandas (six females, four males) with GPS collars and monitored them from September 2019 to March 2020 in Ilam, eastern Nepal. We collected habitat and disturbance data over four seasons. We considered geophysical covariates, anthropogenic factors and habitat fragmentation metrics, and employed linear -mixed models and logistic regression to evaluate the effect of those variables on movement patterns.

RESULTS

The median daily distance travelled by red pandas was 756 m. Males travelled nearly 1.5 times further than females (605 m). Males and sub-adults travelled more in the mating season while females showed no seasonal variation for their daily distance coverage. Red pandas were relatively more active during dawn and morning than the rest of the day, and they exhibited seasonal variation in distance coverage on the diel cycle. Both males and females appeared to be more active in the cub-rearing season, yet males were more active in the dawn in the birthing season. Two sub-adult females dispersed an average of 21 km starting their dispersal with the onset of the new moon following the winter solstice. The single subadult male did not disperse. Red pandas avoided roads, small-habitat patches and large unsuitable areas between habitat patches. Where connected habitat with high forest cover was scarce the animals moved more directly than when habitat was abundant.

CONCLUSIONS

Our study indicates that this habitat specialist is vulnerable to human disturbances and habitat fragmentation. Habitat restoration through improving functional connectivity may be necessary to secure the long-term conservation of specialist species in a human-dominated landscape. Regulation of human activities should go in parallel to minimize disturbances during biologically crucial life phases. We recommend habitat zonation to limit human activities and avoid disturbances, especially livestock herding and road construction in core areas.

Effect of disturbances and habitat fragmentation on an arboreal habitat specialist mammal using GPS telemetry: a case of the red panda

CONTEXT

Habitat specialists residing in human-modified landscapes are likely to be more vulnerable to disturbance because of a functional reliance on very particular habitat features. However, there have been few studies designed to specifically address that issue.

OBJECTIVES

This study aimed to explore how the red panda, an iconic endangered habitat specialist, behaves when faced with disturbances and habitat fragmentation. In particular, we attempted to examine the effect of anthropogenic disturbances and fragmentation on home-range size, activity patterns, and recursion.

METHODS

Using GPS telemetry we monitored 10 red pandas and documented disturbances using camera trapping for one year in eastern Nepal. We performed spatial analysis, analysed activity patterns and evaluated the effect of habitat fragmentation and disturbances on home-range size and residence time using Linear Mixed Models.

RESULTS

Home-range size increased in areas with low availability of forest cover whilst home ranges were smaller in areas with a high road density. Red pandas spent more time in large habitat patches away from roads and cattle stations. Crossing rates suggested that roads acted as a barrier for movement across their habitat. Red pandas also partitioned their activity to minimize interactions with disturbances.

CONCLUSIONS

Red pandas seem to make a trade-off to co-exist in human-dominated landscapes which may have adverse long-term effects on their survival. This indicates that current patterns of habitat fragmentation and forest exploitation may be adversely affecting red panda conservation efforts and that landscape-scale effects should be considered when planning conservation actions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10980-021-01357-w.

First record of Bengal Tiger, Panthera tigris tigris Linnaeus, 1758 (Felidae), in eastern Nepal

We report the first record of a Bengal Tiger, Panthera tigris tigris, in eastern Nepal in 2020 based on photographic evidence. We documented this evidence at 3,165 m a.s.l., which makes it the highest elevation record of a tiger in Nepal. We recorded this evidence in one of 46 trail cameras deployed for monitoring Red Pandas in the Panchthar-Ilam-Taplejung (PIT) area. The PIT area, which has non-protected status, borders India in the east. Our finding supports the importance of transboundary conservation, which will benefit local and flagship wildlife in the PIT area.

Improved Trapping and Handling of an Arboreal, Montane Mammal: Red Panda Ailurus fulgens

Simple Summary

Capture and handling is essential to study some biological and ecological properties of free-ranging animals. However, capturing an arboreal and cryptic species such as the red panda is challenging due to the difficult terrain, their elusive nature, and potential risks to human and animal safety. We developed and successfully tested a protocol for tracking, capture, immobilization, and handling of red pandas. This method could also be used, with some modifications, for other arboreal species. This study extends the known range of body weight and length of free-ranging red pandas. We also report some new morphometric data that could serve as a guide for field identification.

Abstract

It is sometimes essential to have an animal in the hand to study some of their ecological and biological characteristics. However, capturing a solitary, cryptic, elusive arboreal species such as the red panda in the wild is challenging. We developed and successfully tested a protocol for tracking, trapping, immobilization, and handling of red pandas in the wild in eastern Nepal. We established a red panda sighting rate of 0.89 panda/day with a capture success rate of 0.6. We trapped and collared one animal in 3.7 days. On average, we took nearly 136 (range 50–317) min to capture an animal after spotting it. Further processing was completed in 38.5 (21–70) min. Before capture, we found it difficult to recognize the sex of the red panda and to differentiate sub-adults above six months from adults. However, body weight, body length, tail length, shoulder height, and chest girth can be used for diagnosis, as these attributes are smaller in sub-adults. Our method is a welfare-friendly way of trapping and handling wild red pandas. We report new morphometric data that could serve as a guide for field identification.

A conserved haplotype in Wagyu cattle contains RAB4A whose encoded protein regulates glucose trafficking in muscle and fat cells

The Wagyu breed of taurine cattle possess favourable genetics for intramuscular fat (IMF) but genomic loci associated with the trait remain under characterised. Here, we report the identification of a previously unidentified genomic region possessing a particular haplotype structure in Wagyu. Through deployment of a genome-wide haplotype detection analysis that captures regions conserved in a target population but not other populations we screened 100 individual Wagyu and contrasted them with 100 individuals from two independent comparison breeds, Charolais and Angus, using high-density SNPs. An extreme level of Wagyu conservation was assigned to a single genomic window (spanning genomic coordinates BTA28:41 088-300 265 bp). In fact, a five-SNP region spanning 27 096 bp is almost perfectly conserved among the 100 Wagyu individuals assayed and partially overlaps RAB4A. Focussing in, two consecutive SNPs (genomic coordinates 236 949 and 239 950) are apparently fixed within the Wagyu (BB and AA respectively), but at mixed frequencies in the other two breeds. These SNPs are located in the two introns straddling exon 7. In a separate analysis using the 1000 Bulls database, we found that, coincident with exon 7 of RAB4A first allele frequencies were highest in the high IMF Japanese Native (Wagyu) breeds (0.78) and lowest in the low IMF indicine breeds (Nelore and Brahman), with intermediate marbling breeds (Angus and Charolais) assigned intermediate rankings (0.42). RAB4A is known to encode a protein that regulates intracellular trafficking of the insulin-regulated glucose transporter GLUT4. RAB4A can be considered an attractive new positional candidate for IMF development.

A low-density SNP genotyping panel for the accurate prediction of cattle breeds

Genomic tools to better define breed composition in agriculturally important species have sparked scientific and commercial industry interest. Knowledge of breed composition can inform multiple scientifically important decisions of industry application including DNA marker-assisted selection, identification of signatures of selection, and inference of product provenance to improve supply chain integrity. Genomic tools are expensive but can be economized by deploying a relatively small number of highly informative single-nucleotide polymorphisms (SNP) scattered evenly across the genome. Using resources from the 1000 Bull Genomes Project we established calibration (more stringent quality criteria; N = 1,243 cattle) and validation (less stringent; N = 864) data sets representing 17 breeds derived from both taurine and indicine bovine subspecies. Fifteen successively smaller panels (from 500,000 to 50 SNP) were built from those SNP in the calibration data that increasingly satisfied 2 criteria, high differential allele frequencies across the breeds as measured by average Euclidean distance (AED) and high uniformity (even spacing) across the physical genome. Those SNP awarded the highest AED were in or near genes previously identified as important signatures of selection in cattle such as LCORL, NCAPG, KITLG, and PLAG1. For each panel, the genomic breed composition (GBC) of each animal in the validation dataset was estimated using a linear regression model. A systematic exploration of the predictive accuracy of the various sized panels was then undertaken on the validation population using 3 benchmarking approaches: (1) % error (expressed relative to the estimated GBC made from over 1 million SNP), (2) % breed misassignment (expressed relative to each individual's breed recorded), and (3) Shannon's entropy of estimated GBC across the 17 target breeds. Our analyses suggest that a panel of just 250 SNP represents an adequate balance between accuracy and cost-only modest gains in accuracy are made as one increases panel density beyond this point.

Visual Opsin Diversity in Sharks and Rays

The diversity of color vision systems found in extant vertebrates suggests that different evolutionary selection pressures have driven specializations in photoreceptor complement and visual pigment spectral tuning appropriate for an animal's behavior, habitat, and life history. Aquatic vertebrates in particular show high variability in chromatic vision and have become important models for understanding the role of color vision in prey detection, predator avoidance, and social interactions. In this study, we examined the capacity for chromatic vision in elasmobranch fishes, a group that have received relatively little attention to date. We used microspectrophotometry to measure the spectral absorbance of the visual pigments in the outer segments of individual photoreceptors from several ray and shark species, and we sequenced the opsin mRNAs obtained from the retinas of the same species, as well as from additional elasmobranch species. We reveal the phylogenetically widespread occurrence of dichromatic color vision in rays based on two cone opsins, RH2 and LWS. We also confirm that all shark species studied to date appear to be cone monochromats but report that in different species the single cone opsin may be of either the LWS or the RH2 class. From this, we infer that cone monochromacy in sharks has evolved independently on multiple occasions. Together with earlier discoveries in secondarily aquatic marine mammals, this suggests that cone-based color vision may be of little use for large marine predators, such as sharks, pinnipeds, and cetaceans.

Genome-Wide Co-Expression Distributions as a Metric to Prioritize Genes of Functional Importance

Genome-wide gene expression analysis are routinely used to gain a systems-level understanding of complex processes, including network connectivity. Network connectivity tends to be built on a small subset of extremely high co-expression signals that are deemed significant, but this overlooks the vast majority of pairwise signals. Here, we developed a computational pipeline to assign to every gene its pair-wise genome-wide co-expression distribution to one of 8 template distributions shapes varying between unimodal, bimodal, skewed, or symmetrical, representing different proportions of positive and negative correlations. We then used a hypergeometric test to determine if specific genes (regulators versus non-regulators) and properties (differentially expressed or not) are associated with a particular distribution shape. We applied our methodology to five publicly available RNA sequencing (RNA-seq) datasets from four organisms in different physiological conditions and tissues. Our results suggest that genes can be assigned consistently to pre-defined distribution shapes, regarding the enrichment of differential expression and regulatory genes, in situations involving contrasting phenotypes, time-series, or physiological baseline data. There is indeed a striking additional biological signal present in the genome-wide distribution of co-expression values which would be overlooked by currently adopted approaches. Our method can be applied to extract further information from transcriptomic data and help uncover the molecular mechanisms involved in the regulation of complex biological process and phenotypes.

Association analysis of loci implied in "buffering" epistasis

The existence of buffering mechanisms is an emerging property of biological networks, and this results in the buildup of robustness through evolution. So far, there are no explicit methods to find loci implied in buffering mechanisms. However, buffering can be seen as interaction with genetic background. Here we develop this idea into a tractable model for quantitative genetics, in which the buffering effect of one locus with many other loci is condensed into a single statistical effect, multiplicative on the total additive genetic effect. This allows easier interpretation of the results and simplifies the problem of detecting epistasis from quadratic to linear in the number of loci. Using this formulation, we construct a linear model for genome-wide association studies that estimates and declares the significance of multiplicative epistatic effects at single loci. The model has the form of a variance components, norm reaction model and likelihood ratio tests are used for significance. This model is a generalization and explanation of previous ones. We test our model using bovine data: Brahman and Tropical Composite animals, phenotyped for body weight at yearling and genotyped at high density. After association analysis, we find a number of loci with buffering action in one, the other, or both breeds; these loci do not have a significant statistical additive effect. Most of these loci have been reported in previous studies, either with an additive effect or as footprints of selection. We identify buffering epistatic SNPs present in or near genes reported in the context of signatures of selection in multi-breed cattle population studies. Prominent among these genes are those associated with fertility (INHBA, TSHR, ESRRG, PRLR, and PPARG), growth (MSTN, GHR), coat characteristics (KIT, MITF, PRLR), and heat resistance (HSPA6 and HSPA1A). In these populations, we found loci that have a nonsignificant statistical additive effect but a significant epistatic effect. We argue that the discovery and study of loci associated with buffering effects allow attacking the difficult problems, among others, of the release of maintenance variance in artificial and natural selection, of quick adaptation to the environment, and of opposite signs of marker effects in different backgrounds. We conclude that our method and our results generate promising new perspectives for research in evolutionary and quantitative genetics based on the study of loci that buffer effect of other loci.

Dynamics of Gene Co-expression Networks in Time-Series Data: A Case Study in Drosophila melanogaster Embryogenesis

Co-expression networks tightly coordinate the spatiotemporal patterns of gene expression unfolding during development. Due to the dynamic nature of developmental processes simply overlaying gene expression patterns onto static representations of co-expression networks may be misleading. Here, we aim to formally quantitate topological changes of co-expression networks during embryonic development using a publicly available Drosophila melanogaster transcriptome data set comprising 14 time points. We deployed a network approach which inferred 10 discrete co-expression networks by smoothly sliding along from early to late development using 5 consecutive time points per window. Such an approach allows changing network structure, including the presence of hubs, modules and other topological parameters to be quantitated. To explore the dynamic aspects of gene expression captured by our approach, we focused on regulator genes with apparent influence over particular aspects of development. Those key regulators were selected using a differential network algorithm to contrast the first 7 (early) with the last 7 (late) developmental time points. This assigns high scores to genes whose connectivity to abundant differentially expressed target genes has changed dramatically between states. We have produced a list of key regulators – some increasing (e.g., Tusp, slbo, Sidpn, DCAF12, and chinmo) and some decreasing (Rfx, bap, Hmx, Awh, and mld) connectivity during development – which reflects their role in different stages of embryogenesis. The networks we have constructed can be explored and interpreted within Cytoscape software and provide a new systems biology approach for the Drosophila research community to better visualize and interpret developmental regulation of gene expression.

The Efficacy of Sunitinib Treatment of Renal Cancer Cells Is Associated with the Protein PHAX In Vitro

Anti-angiogenic agents, such as the multi-tyrosine kinase inhibitor sunitinib, are key first line therapies for metastatic clear cell renal cell carcinoma (ccRCC), but their mechanism of action is not fully understood. Here, we take steps towards validating a computational prediction based on differential transcriptome network analysis that phosphorylated adapter RNA export protein (PHAX) is associated with sunitinib drug treatment. The regulatory impact factor differential network algorithm run on patient tissue samples suggests PHAX is likely an important regulator through changes in genome-wide network connectivity. Immunofluorescence staining of patient tumours showed strong localisation of PHAX to the microvasculature consistent with the anti-angiogenic effect of sunitinib. In normal kidney tissue, PHAX protein abundance was low but increased with tumour grade (G1 vs. G3/4; p < 0.01), consistent with a possible role in cancer progression. In organ culture, ccRCC cells had higher levels of PHAX protein expression than normal kidney cells, and sunitinib increased PHAX protein expression in a dose dependent manner (untreated vs. 100 µM; p < 0.05). PHAX knockdown in a ccRCC organ culture model impacted the ability of sunitinib to cause cancer cell death (p < 0.0001 untreated vs. treated), suggesting a role for PHAX in mediating the efficacy of sunitinib.

Gene expression identifies metabolic and functional differences between intramuscular and subcutaneous adipocytes in cattle

Background

This study used a genome-wide screen of gene expression to better understand the metabolic and functional differences between commercially valuable intramuscular fat (IMF) and commercially wasteful subcutaneous (SC) fat depots in Bos taurus beef cattle.

Results

We confirmed many findings previously made at the biochemical level and made new discoveries. The fundamental lipogenic machinery, such as ACACA and FASN encoding the rate limiting Acetyl CoA carboxylase and Fatty Acid synthase were expressed at 1.6–1.8 fold lower levels in IMF, consistent with previous findings. The FA elongation pathway including the rate limiting ELOVL6 was also coordinately downregulated in IMF compared to SC as expected. A 2-fold lower expression in IMF of ACSS2 encoding Acetyl Coenzyme A synthetase is consistent with utilisation of less acetate for lipogenesis in IMF compared to SC as previously determined using radioisotope incorporation. Reduced saturation of fat in the SC depot is reflected by 2.4 fold higher expression of the SCD gene encoding the Δ9 desaturase enzyme. Surprisingly, CH25H encoding the cholesterol 25 hydroxylase enzyme was ~ 36 fold upregulated in IMF compared to SC. Moreover, its expression in whole muscle tissue appears representative of the proportional representation of bovine marbling adipocytes. This suite of observations prompted quantification of a set of oxysterols (oxidised forms of cholesterol) in the plasma of 8 cattle exhibiting varying IMF. Using Liquid Chromatography-Mass Spectrometry (LC-MS) we found the levels of several oxysterols were significantly associated with multiple marbling measurements across the musculature, but (with just one exception) no other carcass phenotypes.

Conclusions

These data build on our molecular understanding of ruminant fat depot biology and suggest oxysterols represent a promising circulating biomarker for cattle marbling.

RAPID COMMUNICATION: A haplotype information theory method reveals genes of evolutionary interest in European vs. Asian pigs

Asian and European wild boars were independently domesticated ca. 10,000 yr ago. Since the 17th century, Chinese breeds have been imported to Europe to improve the genetics of European animals by introgression of favorable alleles, resulting in a complex mosaic of haplotypes. To interrogate the structure of these haplotypes further, we have run a new haplotype segregation analysis based on information theory, namely compression efficiency (CE). We applied the approach to sequence data from individuals from each phylogeographic region (n = 23 from Asia and Europe) including a number of major pig breeds. Our genome-wide CE is able to discriminate the breeds in a manner reflecting phylogeography. Furthermore, 24,956 nonoverlapping sliding windows (each comprising 1,000 consecutive SNP) were quantified for extent of haplotype sharing within and between Asia and Europe. The genome-wide distribution of extent of haplotype sharing was quite different between groups. Unlike European pigs, Asian pigs haplotype sharing approximates a normal distribution. In line with this, we found the European breeds possessed a number of genomic windows of dramatically higher haplotype sharing than the Asian breeds. Our CE analysis of sliding windows captures some of the genomic regions reported to contain signatures of selection in domestic pigs. Prominent among these regions, we highlight the role of a gene encoding the mitochondrial enzyme LACTB which has been associated with obesity, and the gene encoding MYOG a fundamental transcriptional regulator of myogenesis. The origin of these regions likely reflects either a population bottleneck in European animals, or selective targets on commercial phenotypes reducing allelic diversity in particular genes and/or regulatory regions.

MicroRNA profiling associated with muscle growth in modern broilers compared to an unselected chicken breed

Background

Genetically selected modern broiler chickens have acquired outstanding production efficiency through rapid growth and improved feed efficiency compared to unselected chicken breeds. Recently, we analyzed the transcriptome of breast muscle tissues obtained from modern pedigree male (PeM) broilers (rapid growth and higher efficiency) and foundational Barred Plymouth Rock (BPR) chickens (slow growth and poorer efficiency). This study was designed to investigate microRNAs that play role in rapid growth of the breast muscles in modern broiler chickens.

Results

In this study, differential abundance of microRNA (miRNA) was analyzed in breast muscle of PeM and BPR chickens and the results were integrated with differentially expressed (DE) mRNA in the same tissues. A total of 994 miRNA were identified in PeM and BPR chicken lines from the initial analysis of small RNA sequencing data. After filtering and statistical analyses, the results showed miR-2131-5p, miR-221-5p, miR-126-3p, miR-146b-5p, miR-10a-5p, let-7b, miR-125b-5p, and miR-146c-5p up-regulated whereas miR-206 down-regulated in PeM compared to BPR breast muscle. Based on inhibitory regulations of miRNAs on the mRNA abundance, our computational analysis using miRDB, an online software, predicated that 118 down-regulated mRNAs may be targeted by the up-regulated miRNAs, while 35 up-regulated mRNAs appear to be due to a down-regulated miRNA (i.e., miR-206). Functional network analyses of target genes of DE miRNAs showed their involvement in calcium signaling, axonal guidance signaling, and NRF2-mediated oxidative stress response pathways suggesting their involvement in breast muscle growth in chickens.

Conclusion

From the integrated analyses of differentially expressed miRNA-mRNA data, we were able to identify breast muscle specific miRNAs and their target genes whose concerted actions can contribute to rapid growth and higher feed efficiency in modern broiler chickens. This study provides foundation data for elucidating molecular mechanisms that govern muscle growth in chickens.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5061-7) contains supplementary material, which is available to authorized users.

Data compression can discriminate broilers by selection line, detect haplotypes, and estimate genetic potential for complex phenotypes

Accurately establishing the relationships among individuals lays the foundation for genetic analyses such as genome-wide association studies and identification of selection signatures. Of particular interest to the poultry industry are estimates of genetic merit based on molecular data. These estimates can be commercially exploited in marker-assisted breeding programs to accelerate genetic improvement. Here, we test the utility of a new method we have recently developed to estimate animal relatedness and applied it to genetic parameter estimation in commercial broilers. Our approach is based on the concept of data compression from information theory. Using the real-world compressor gzip to estimate normalized compression distance (NCD) we have built compression-based relationship matrices (CRM) for 988 chickens from 4 commercial broiler lines-2 male and 2 female lines. For all pairs of individuals, we found a strong negative relationship between the commonly used genomic relationship matrix (GRM) and NCD. This reflects the fact that "similarity" is the inverse of "distance." The CRM explained more genetic variation than the corresponding GRM in 2 of 3 phenotypes, with corresponding improvements in accuracy of genomic-enabled predictions of breeding value. A sliding-window version of the analysis highlighted haplotype regions of the genome apparently under selection in a line-specific manner. In the male lines, we retrieved high population-specific scores for IGF-1 and a cognate receptor, INSR. For the female lines, we detected an extreme score for a region containing a reproductive hormone receptor (GNRHR). We conclude that our compression-based method is a valid approach to established relationships and identify regions under selective pressure in commercial lines of broiler chickens.

Molecular regulation of high muscle mass in developing Blonde d'Aquitaine cattle foetuses

The Blonde d'Aquitaine (BA) is a French cattle breed with enhanced muscularity, partly attributable to a MSTN mutation. The BA m. Semitendinosus has a faster muscle fibre isoform phenotype comprising a higher proportion of fast type IIX fibres compared to age-matched Charolais (CH). To better understand the molecular network of modifications in BA compared to CH muscle, we assayed the transcriptomes of the m. Semitendinosus at 110, 180, 210 and 260 days postconception (dpc). We used a combination of differential expression (DE) and regulatory impact factors (RIF) to compare and contrast muscle gene expression between the breeds. Prominently developmentally regulated genes in both breeds reflected the replacement of embryonic myosin isoforms (MYL4, MYH3) with adult isoforms (MYH1) and the upregulation of mitochondrial metabolism (CKMT2, AGXT2L1) in preparation for birth. However, the transition to a fast, glycolytic muscle phenotype in the MSTN mutant BA is detectable through downregulation of various slow twitch subunits (TNNC1, MYH7, TPM3, CSRP3) beyond 210 dpc, and a small but consistent genome-wide reduction in mRNA encoding the mitoproteome. Across the breeds, NRIP2 is the regulatory gene possessing a network change most similar to that of MSTN.

Summary: We explored the molecular basis of high muscle mass in Blonde d’Aquitaine cattle, and found that mRNA encoding slow fibres and mitochondrial proteins were downregulated in this breed.

Proteogenomics Reveals Enriched Ribosome Assembly and Protein Translation in Pectoralis major of High Feed Efficiency Pedigree Broiler Males

Background: In production animal agriculture, the cost of feed represents 60–70% of the total cost of raising an animal to market weight. Thus, development of viable biomarkers for feed efficiency (FE, g gain/g feed) to assist in genetic selection of breeding stock remains an important goal in commercial breeding programs.

Methods: Global gene (cDNA microarray, RNAseq) and protein expression (shotgun proteomics) analyses have been conducted on breast muscle samples obtained from pedigree broiler males (PedM) exhibiting high and low FE phenotypes. Using the entire datasets (i.e., no cutoffs for significance or fold difference in expression) the number of genes or proteins that were expressed numerically higher or lower in the high FE compared to the low FE phenotype for key terms or functions, e.g., ribosomal, mitochondrial ribosomal, tRNA, RNA binding motif, RNA polymerase, small nuclear ribonucleoprotein, and protein tyrosine phosphatase, were determined. Bionomial distribution analysis (exact) was then conducted on these datasets to determine significance between numerically up or down expression.

Results: Processes associated with mitochondrial proteome expression (e.g., mitochondrial ribosomal proteins, mitochondrial transcription, mitochondrial tRNA, and translation) were enriched in breast muscle from the high FE compared to the low FE pedigree male broiler phenotype. Furthermore, the high FE phenotype exhibited enrichment of ribosome assembly (e.g., RNA polymerase, mitochondrial and cytosolic ribosomes, small, and heterogeneous nuclear ribonucleoproteins), as well as nuclear transport and protein translation processes compared to the low FE phenotype. Quality control processes (proteosomes and autophagy) were also enriched in the high FE phenotype. In contrast, the low FE phenotype exhibited enrichment of cytoskeletal proteins, protein tyrosine phosphatases, and tyrosine kinases compared to the high FE phenotype. These results suggest that processes of mitochondrial and cytosolic ribosomal construction, activity, and protein translation would be enhanced in high FE breast muscle, and that phosphorylation of tyrosine moieties of proteins could be prolonged in the high compared to low FE phenotype. The results indicate the presence of a proteogenomic architecture that could enhance ribosome construction, protein translation, and quality control processes and contribute to the phenotypic expression of feed efficiency in this PedM broiler model.

Progesterone signalling in broiler skeletal muscle is associated with divergent feed efficiency

Background

We contrast the pectoralis muscle transcriptomes of broilers selected from within a single genetic line expressing divergent feed efficiency (FE) in an effort to improve our understanding of the mechanistic basis of FE.

Results

Application of a virtual muscle model to gene expression data pointed to a coordinated reduction in slow twitch muscle isoforms of the contractile apparatus (MYH15, TPM3, MYOZ2, TNNI1, MYL2, MYOM3, CSRP3, TNNT2), consistent with diminishment in associated slow machinery (myoglobin and phospholamban) in the high FE animals. These data are in line with the repeated transition from red slow to white fast muscle fibres observed in agricultural species selected on mass and FE. Surprisingly, we found that the expression of 699 genes encoding the broiler mitoproteome is modestly–but significantly–biased towards the high FE group, suggesting a slightly elevated mitochondrial content. This is contrary to expectation based on the slow muscle isoform data and theoretical physiological capacity arguments. Reassuringly, the extreme 40 most DE genes can successfully cluster the 12 individuals into the appropriate FE treatment group. Functional groups contained in this DE gene list include metabolic proteins (including opposing patterns of CA3 and CA4), mitochondrial proteins (CKMT1A), oxidative status (SEPP1, HIG2A) and cholesterol homeostasis (APOA1, INSIG1). We applied a differential network method (Regulatory Impact Factors) whose aim is to use patterns of differential co-expression to detect regulatory molecules transcriptionally rewired between the groups. This analysis clearly points to alterations in progesterone signalling (via the receptor PGR) as the major driver. We show the progesterone receptor localises to the mitochondria in a quail muscle cell line.

Conclusions

Progesterone is sometimes used in the cattle industry in exogenous hormone mixes that lead to a ~20% increase in FE. Because the progesterone receptor can localise to avian mitochondria, our data continue to point to muscle mitochondrial metabolism as an important component of the phenotypic expression of variation in broiler FE.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0396-2) contains supplementary material, which is available to authorized users.

Chicken muscle mitochondrial content appears co-ordinately regulated and is associated with performance phenotypes

Mitochondrial content is a fundamental cellular bioenergetic phenotype. Previous work has hypothesised possible links between variation in muscle mitochondrial content and animal performance. However, no population screens have been performed in any production species. Here, we have designed a high throughput molecular approach to estimate mitochondrial content in commercial broilers. Technical validity was established using several approaches, including its performance in monoclonal DF-1 cells, cross-tissue comparisons in tissues with differing metabolic demands (white fat<breast muscle<drumstick muscle<heart muscle) and, as a negative control, a near absence of mtDNA amplification from whole blood. We screened breast muscle and thigh muscle in 80 birds individually phenotyped for 11 growth and development traits. Substantial individual variation (fivefold) was discovered in both breast and thigh muscle mitochondrial content. Interestingly, across birds we detected a very strong positive relationship between breast and thigh content (correlation coefficient 0.61; P<0.0001), consistent with coordinate regulatory control across the musculature. Further, breast muscle mitochondrial content is negatively correlated with breast muscle yield (−0.27; P=0.037), abdominal fat content (−0.31; P=0.017) and carcass yield (−0.26; P=0.045). Therefore, low breast muscle mitochondrial content is associated with more muscular birds possessing higher abdominal fat, the latter being in line with biomedical models of obesity. Finally, thigh mitochondrial content is negatively correlated with the bow out leg defect (−0.30; P=0.011). Overall, our data point to mitochondrial content as a promising consideration in predictive modelling of production traits.

Summary: We have detected variation in muscle mitochondrial content across chickens. This variation relates to phenotypes and also indicates coordinate regulation of mitochondrial biogenesis across the musculature.

Disorder of endoplasmic reticulum calcium channel components is associated with the increased apoptotic potential in pale, soft, exudative pork

Eight pale, soft and exudative (PSE) and eight reddish-pink, firm and non-exudative (RFN) porcine longissimus muscle samples were selected based on pH and L* at 1h postmortem (PM), and drip loss at 24h PM, and used to evaluate the cellular calcium and apoptosis status. We found that SERCA1 was decreased, while IP3R was decreased in PSE meat (P<0.05), indicative of the overloaded sarcoplasmic calcium status. In PSE meat, the pro-apoptotic factor BAX was increased while the anti-apoptotic factor Bcl-2 was decreased (P<0.05). The significantly increased activity of caspase 3 and the expression of its cleavage fragment suggested higher apoptotic potential in PSE meat compared with RFN meat (P<0.05). Moreover, the significantly higher expression level of cytochrome C (P<0.05) suggests the important role of mitochondria during apoptosis appearance in PSE meat. Taken together, our data inferred that the calcium channel disorder present in PSE meat was associated with the increased apoptotic potential.

Compression distance can discriminate animals by genetic profile, build relationship matrices and estimate breeding values

BACKGROUND

Genetic relatedness is currently estimated by a combination of traditional pedigree-based approaches (i.e. numerator relationship matrices, NRM) and, given the recent availability of molecular information, using marker genotypes (via genomic relationship matrices, GRM). To date, GRM are computed by genome-wide pair-wise SNP (single nucleotide polymorphism) correlations.

RESULTS

We describe a new estimate of genetic relatedness using the concept of normalised compression distance (NCD) that is borrowed from Information Theory. Analogous to GRM, the resultant compression relationship matrix (CRM) exploits numerical patterns in genome-wide allele order and proportion, which are known to vary systematically with relatedness. We explored properties of the CRM in two industry cattle datasets by analysing the genetic basis of yearling weight, a phenotype of moderate heritability. In both Brahman (Bos indicus) and Tropical Composite (Bos taurus by Bos indicus) populations, the clustering inferred by NCD was comparable to that based on SNP correlations using standard principal component analysis approaches. One of the versions of the CRM modestly increased the amount of explained genetic variance, slightly reduced the 'missing heritability' and tended to improve the prediction accuracy of breeding values in both populations when compared to both NRM and GRM. Finally, a sliding window-based application of the compression approach on these populations identified genomic regions influenced by introgression of taurine haplotypes.

CONCLUSIONS

For these two bovine populations, CRM reduced the missing heritability and increased the amount of explained genetic variation for a moderately heritable complex trait. Given that NCD can sensitively discriminate closely related individuals, we foresee CRM having possible value for estimating breeding values in highly inbred populations.

Differences in muscle transcriptome among pigs phenotypically extreme for fatty acid composition

BACKGROUND

Besides having an impact on human health, the porcine muscle fatty acid profile determines meat quality and taste. The RNA-Seq technologies allowed us to explore the pig muscle transcriptome with an unprecedented detail. The aim of this study was to identify differentially-expressed genes between two groups of 6 sows belonging to an Iberian × Landrace backcross with extreme phenotypes according to FA profile.

RESULTS

We sequenced the muscle transcriptome acquiring 787.5 M of 75 bp paired-end reads. About 85.1% of reads were mapped to the reference genome. Of the total reads, 79.1% were located in exons, 6.0% in introns and 14.9% in intergenic regions, indicating expressed regions not annotated in the reference genome. We identified a 34.5% of the intergenic regions as interspersed repetitive regions. We predicted a total of 2,372 putative proteins. Pathway analysis with 131 differentially-expressed genes revealed that the most statistically-significant metabolic pathways were related with lipid metabolism. Moreover, 18 of the differentially-expressed genes were located in genomic regions associated with IMF composition in an independent GWAS study in the same genetic background. Thus, our results indicate that the lipid metabolism of FAs is differently modulated when the FA composition in muscle differs. For instance, a high content of PUFA may reduce FA and glucose uptake resulting in an inhibition of the lipogenesis. These results are consistent with previous studies of our group analysing the liver and the adipose tissue transcriptomes providing a view of each of the main organs involved in lipid metabolism.

CONCLUSIONS

The results obtained in the muscle transcriptome analysis increase the knowledge of the gene regulation of IMF deposition, FA profile and meat quality, in terms of taste and nutritional value. Besides, our results may be important in terms of human health.

A marker-derived gene network reveals the regulatory role of PPARGC1A, HNF4G, and FOXP3 in intramuscular fat deposition of beef cattle

High intramuscular fat (IMF) awards price premiums to beef producers and is associated with meat quality and flavor. Studying gene interactions and pathways that affect IMF might unveil causative physiological mechanisms and inform genomic selection, leading to increased accuracy of predictions of breeding value. To study gene interactions and pathways, a gene network was derived from genetic markers associated with direct measures of IMF, other fat phenotypes, feedlot performance, and a number of meat quality traits relating to body conformation, development, and metabolism that might be plausibly expected to interact with IMF biology. Marker associations were inferred from genomewide association studies (GWAS) based on high density genotypes and 29 traits measured on 10,181 beef cattle animals from 3 breed types. For the network inference, SNP pairs were assessed according to the strength of the correlation between their additive association effects across the 29 traits. The co-association inferred network was formed by 2,434 genes connected by 28,283 edges. Topological network parameters suggested a highly cohesive network, in which the genes are strongly functionally interconnected. Pathway and network analyses pointed towards a trio of transcription factors (TF) as key regulators of carcass IMF: PPARGC1A, HNF4G, and FOXP3. Importantly, none of these genes would have been deemed as significantly associated with IMF from the GWAS. Instead, a total of 313 network genes show significant co-association with the 3 TF. These genes belong to a wide variety of biological functions, canonical pathways, and genetic networks linked to IMF-related phenotypes. In summary, our GWAS and network predictions are supported by the current literature and suggest a cooperative role for the 3 TF and other interacting genes including CAPN6, STC2, MAP2K4, EYA1, COPS5, XKR4, NR2E1, TOX, ATF1, ASPH, TGS1, and TTPA as modulators of carcass and meat quality traits in beef cattle.

Information compression exploits patterns of genome composition to discriminate populations and highlight regions of evolutionary interest

Background

Genomic information allows population relatedness to be inferred and selected genes to be identified. Single nucleotide polymorphism microarray (SNP-chip) data, a proxy for genome composition, contains patterns in allele order and proportion. These patterns can be quantified by compression efficiency (CE). In principle, the composition of an entire genome can be represented by a CE number quantifying allele representation and order.

Results

We applied a compression algorithm (DEFLATE) to genome-wide high-density SNP data from 4,155 human, 1,800 cattle, 1,222 sheep, 81 dogs and 49 mice samples. All human ethnic groups can be clustered by CE and the clusters recover phylogeography based on traditional fixation index (F_ST) analyses. CE analysis of other mammals results in segregation by breed or species, and is sensitive to admixture and past effective population size. This clustering is a consequence of individual patterns such as runs of homozygosity. Intriguingly, a related approach can also be used to identify genomic loci that show population-specific CE segregation. A high resolution CE ‘sliding window’ scan across the human genome, organised at the population level, revealed genes known to be under evolutionary pressure. These include SLC24A5 (European and Gujarati Indian skin pigmentation), HERC2 (European eye color), LCT (European and Maasai milk digestion) and EDAR (Asian hair thickness). We also identified a set of previously unidentified loci with high population-specific CE scores including the chromatin remodeler SCMH1 in Africans and EDA2R in Asians. Closer inspection reveals that these prioritised genomic regions do not correspond to simple runs of homozygosity but rather compositionally complex regions that are shared by many individuals of a given population. Unlike F_ST, CE analyses do not require ab initio population comparisons and are amenable to the hemizygous X chromosome.

Conclusions

We conclude with a discussion of the implications of CE for a complex systems science view of genome evolution. CE allows one to clearly visualise the evolution of individual genomes and populations through a formal, mathematically-rigorous information space. Overall, CE makes a set of biological predictions, some of which are unique and await functional validation.

RNF14 is a regulator of mitochondrial and immune function in muscle

BACKGROUND

Muscle development and remodelling, mitochondrial physiology and inflammation are thought to be inter-related and to have implications for metabolism in both health and disease. However, our understanding of their molecular control is incomplete.

RESULTS

In this study we have confirmed that the ring finger 14 protein (RNF14), a poorly understood transcriptional regulator, influences the expression of both mitochondrial and immune-related genes. The prediction was based on a combination of network connectivity and differential connectivity in cattle (a non-model organism) and mice data sets, with a focus on skeletal muscle. They assigned similar probability to mammalian RNF14 playing a regulatory role in mitochondrial and immune gene expression. To try and resolve this apparent ambiguity we performed a genome-wide microarray expression analysis on mouse C2C12 myoblasts transiently transfected with two Rnf14 transcript variants that encode 2 naturally occurring but different RNF14 protein isoforms. The effect of both constructs was significantly different to the control samples (untransfected cells and cells transfected with an empty vector). Cluster analyses revealed that transfection with the two Rnf14 constructs yielded discrete expression signatures from each other, but in both cases a substantial set of genes annotated as encoding proteins related to immune function were perturbed. These included cytokines and interferon regulatory factors. Additionally, transfection of the longer transcript variant 1 coordinately increased the expression of 12 (of the total 13) mitochondrial proteins encoded by the mitochondrial genome, 3 of which were significant in isolated pair-wise comparisons (Mt-coxII, Mt-nd2 and mt-nd4l). This apparent additional mitochondrial function may be attributable to the RWD protein domain that is present only in the longer RNF14 isoform.

CONCLUSIONS

RNF14 influences the expression of both mitochondrial and immune related genes in a skeletal muscle context, and has likely implications for the inter-relationship between bioenergetic status and inflammation.

Inferring the in vivo cellular program of developing bovine skeletal muscle from expression data

We outline an in vivo cellular program of bovine longissimus muscle development inferred from expression data from 60 days post conception to 3months postnatal. Analytic challenges included changes in cellular composition, ambiguous 'diagnostic' markers of cell type and contrasts between cattle human and mouse myogenesis. Nevertheless, the expression profiles of the myosin isoforms support slow and fast muscle fibres emanating from primary and secondary myogenesis respectively, while expression of the prenatal myosin subunits is down regulated prior to birth. Of the canonical pro-myogenic transcription factors (TF), MYF6 and MYF5 are negatively co-expressed, with MYF6 displaying higher expression in the post-natal samples and MYF5, MYOG, HES6 and PAX7 displaying higher expression in early development. A set of TFs (SIX1, EYA2 and DACH2) considered important in undifferentiated murine cells were equally abundant in differentiated bovine cells. An examination of mammalian regulators of fibre composition, muscle mass and muscle metabolism, underscored the roles of PPARGC1A, TGFβ signalling and the NHR4 Nuclear Hormone Receptors on bovine muscle development. Enriched among the most variably expressed genes from the entire data set were molecules regulating mitochondrial metabolism of carbohydrate (PDK4), fat (UCP3), protein (AGXT2L1) and high energy phosphate (CKMT2). The dramatic increase in the expression of these transcripts, which may enable the peri-natal transition to metabolic independence critical for new-born herbivores, provides surprising evidence for substantial developmental remodelling of muscle mitochondria and reflects changes in nutrient availability. Overall, despite differences in size, metabolism and physiology, the muscle structural subunit expression program appears very similar in ruminants, rodents and humans.

Gene expression phenotypes for lipid metabolism and intramuscular fat in skeletal muscle of cattle

Gene expression phenotypes were evaluated for intramuscular fat (IMF) in bovine skeletal muscle as an alternative to traditional estimates of IMF%. Gene expression data from a time course of LM development in high- and low-marbling Bos taurus cattle crosses were compared to identify genes involved in intramuscular adipocyte lipid metabolism with developmentally similar gene expression profiles. Three sets of genes were identified: triacylglyceride (TAG) synthesis and storage, fatty acid (FA) synthesis, and PPARγ-related genes. In an independent analysis in the LM of 48 Bos indicus cattle, TAG and FA gene sets were enriched in the top 100 genes of which expression was most correlated with IMF% (P = 1.2 × 10(-24) and 3.5 × 10(-9), respectively). In general, genes encoding enzymes involved in the synthesis of FA and TAG in the intramuscular adipocytes were present in the top 100 genes. In B. indicus, effects of a steroid hormone growth promotant (HGP), 2 experimental sites [New South Wales (NSW) and Western Australia (WA)], and 3 tenderness genotypes on the expression levels of genes in the TAG gene set and the correlation of gene expression with IMF% were investigated. Although correlation between expression of 12 individual TAG genes and IMF% was observed in HGP-treated animals in both experimental sites (mean r = 0.43), correlation was not observed for untreated animals at the NSW site (mean r = -0.07, P < 3 × 10(-6)). However, TAG genes showed an average 1.6-fold (P < 0.0004) reduction in expression in the LM of HGP-treated cattle relative to untreated cattle, an effect consistent across both experimental sites. Cattle possessing the favored tenderness calpain 1 and 3 and calpastatin alleles exhibited a greater (P = 0.008) reduction in expression in NSW (1.8-fold reduction, P = 0.0002) compared with WA (1.2-fold reduction, P = 0.03). Tenderness genotype had no impact (P > 0.05) on the correlation of TAG genes with IMF%. In general, the interactions among genotype, treatment and location, and TAG gene set gene expression were consistent with the interactions among the same factors and IMF% detected using 255 animals, of which the 48 in this study were a subset. Thus, the TAG gene set constitutes a gene expression phenotype able to predict effects of different genotypes and treatments on IMF% using much smaller groups than current approaches, even in animals with very low IMF%.

An Always Correlated gene expression landscape for ovine skeletal muscle, lessons learnt from comparison with an "equivalent" bovine landscape

Background

We have recently described a method for the construction of an informative gene expression correlation landscape for a single tissue, longissimus muscle (LM) of cattle, using a small number (less than a hundred) of diverse samples. Does this approach facilitate interspecies comparison of networks?

Findings

Using gene expression datasets from LM samples from a single postnatal time point for high and low muscling sheep, and from a developmental time course (prenatal to postnatal) for normal sheep and sheep exhibiting the Callipyge muscling phenotype gene expression correlations were calculated across subsets of the data comparable to the bovine analysis. An “Always Correlated” gene expression landscape was constructed by integrating the correlations from the subsets of data and was compared to the equivalent landscape for bovine LM muscle. Whilst at the high level apparently equivalent modules were identified in the two species, at the detailed level overlap between genes in the equivalent modules was limited and generally not significant. Indeed, only 395 genes and 18 edges were in common between the two landscapes.

Conclusions

Since it is unlikely that the equivalent muscles of two closely related species are as different as this analysis suggests, within tissue gene expression correlations appear to be very sensitive to the samples chosen for their construction, compounded by the different platforms used. Thus users need to be very cautious in interpretation of the differences. In future experiments, attention will be required to ensure equivalent experimental designs and use cross-species gene expression platform to enable the identification of true differences between different species.

Porcine tissue-specific regulatory networks derived from meta-analysis of the transcriptome

The processes that drive tissue identity and differentiation remain unclear for most tissue types. So are the gene networks and transcription factors (TF) responsible for the differential structure and function of each particular tissue, and this is particularly true for non model species with incomplete genomic resources. To better understand the regulation of genes responsible for tissue identity in pigs, we have inferred regulatory networks from a meta-analysis of 20 gene expression studies spanning 480 Porcine Affymetrix chips for 134 experimental conditions on 27 distinct tissues. We developed a mixed-model normalization approach with a covariance structure that accommodated the disparity in the origin of the individual studies, and obtained the normalized expression of 12,320 genes across the 27 tissues. Using this resource, we constructed a network, based on the co-expression patterns of 1,072 TF and 1,232 tissue specific genes. The resulting network is consistent with the known biology of tissue development. Within the network, genes clustered by tissue and tissues clustered by site of embryonic origin. These clusters were significantly enriched for genes annotated in key relevant biological processes and confirm gene functions and interactions from the literature. We implemented a Regulatory Impact Factor (RIF) metric to identify the key regulators in skeletal muscle and tissues from the central nervous systems. The normalization of the meta-analysis, the inference of the gene co-expression network and the RIF metric, operated synergistically towards a successful search for tissue-specific regulators. Novel among these findings are evidence suggesting a novel key role of ERCC3 as a muscle regulator. Together, our results recapitulate the known biology behind tissue specificity and provide new valuable insights in a less studied but valuable model species.

Beyond differential expression: the quest for causal mutations and effector molecules

High throughput gene expression technologies are a popular choice for researchers seeking molecular or systems-level explanations of biological phenomena. Nevertheless, there has been a groundswell of opinion that these approaches have not lived up to the hype because the interpretation of the data has lagged behind its generation. In our view a major problem has been an over-reliance on isolated lists of differentially expressed (DE) genes which – by simply comparing genes to themselves – have the pitfall of taking molecular information out of context. Numerous scientists have emphasised the need for better context. This can be achieved through holistic measurements of differential connectivity in addition to, or in replacement, of DE. However, many scientists continue to use isolated lists of DE genes as the major source of input data for common readily available analytical tools. Focussing this opinion article on our own research in skeletal muscle, we outline our resolutions to these problems – particularly a universally powerful way of quantifying differential connectivity. With a well designed experiment, it is now possible to use gene expression to identify causal mutations and the other major effector molecules with whom they cooperate, irrespective of whether they themselves are DE. We explain why, for various reasons, no other currently available experimental techniques or quantitative analyses are capable of reaching these conclusions.

Mitochondrial treason: a driver of pH decline rate in post-mortem muscle?

The rate and extent of post-mortem pH decline is an important determinant of meat quality, influencing among other things the incidence of pale, soft, exudative meat. In this perspective I argue two things. First, the accuracy of our modelling of pH decline rates might be improved by incorporation of a muscle mitochondrial content parameter. Second, that high mitochondrial content might help explain the surprisingly rapid pH decline that has sometimes been reported in red, oxidative muscles. The proposed physiological mechanism relates to the fact that mitochondria are known to act as net consumers, rather than producers, of ATP during periods of anoxia – a process called mitochondrial treason. To the best of my knowledge, the implications of the phenomenon of mitochondrial treason have never previously been applied or discussed in a meat science context.

Chronic exposure to anabolic steroids induces the muscle expression of oxytocin and a more than fiftyfold increase in circulating oxytocin in cattle

Molecular mechanisms in skeletal muscle associated with anabolic steroid treatment of cattle are unclear and we aimed to characterize transcriptional changes. Cattle were chronically exposed (68 ± 20 days) to a steroid hormone implant containing 200 mg trenbolone acetate and 20 mg estradiol (Revalor-H). Biopsy samples from 48 cattle (half treated) from longissimus dorsi (LD) muscle under local anesthesia were collected. Gene expression levels were profiled by microarray, covering 16,944 unique bovine genes: 121 genes were differentially expressed (DE) due to the implant (99.99% posterior probability of not being false positives). Among DE genes, a decrease in expression of a number of fat metabolism-associated genes, likely reflecting the lipid storage activity of intramuscular adipocytes, was observed. The expression of IGF1 and genes related to the extracellular matrix, slow twitch fibers, and cell cycle (including SOX8, a satellite cell marker) was increased in the treated muscle. Unexpectedly, a very large 21- (microarray) to 97 (real time quantitative PCR)-fold higher expression of the mRNA encoding the neuropeptide hormone oxytocin was observed in treated muscle. We also observed an ∼50-fold higher level of circulating oxytocin in the plasma of treated animals at the time of biopsy. Using a coexpression network strategy OXTR was identified as more likely than IGF1R to be a major mediator of the muscle response to Revalor-H. A re-investigation of in vivo cattle LD muscle samples during early to mid-fetal development identified a >128-fold increased expression of OXT, coincident with myofiber differentiation and fusion. We propose that oxytocin may be involved in mediating the anabolic effects of Revalor-H treatment.

Musical beauty and information compression: Complex to the ear but simple to the mind?

BACKGROUND

The biological origin of music, its universal appeal across human cultures and the cause of its beauty remain mysteries. For example, why is Ludwig Van Beethoven considered a musical genius but Kylie Minogue is not? Possible answers to these questions will be framed in the context of Information Theory.

PRESENTATION OF THE HYPOTHESIS

The entire life-long sensory data stream of a human is enormous. The adaptive solution to this problem of scale is information compression, thought to have evolved to better handle, interpret and store sensory data. In modern humans highly sophisticated information compression is clearly manifest in philosophical, mathematical and scientific insights. For example, the Laws of Physics explain apparently complex observations with simple rules. Deep cognitive insights are reported as intrinsically satisfying, implying that at some point in evolution, the practice of successful information compression became linked to the physiological reward system. I hypothesise that the establishment of this "compression and pleasure" connection paved the way for musical appreciation, which subsequently became free (perhaps even inevitable) to emerge once audio compression had become intrinsically pleasurable in its own right.

TESTING THE HYPOTHESIS

For a range of compositions, empirically determine the relationship between the listener's pleasure and "lossless" audio compression. I hypothesise that enduring musical masterpieces will possess an interesting objective property: despite apparent complexity, they will also exhibit high compressibility.

IMPLICATIONS OF THE HYPOTHESIS

Artistic masterpieces and deep Scientific insights share the common process of data compression. Musical appreciation is a parasite on a much deeper information processing capacity. The coalescence of mathematical and musical talent in exceptional individuals has a parsimonious explanation. Musical geniuses are skilled in composing music that appears highly complex to the ear yet transpires to be highly simple to the mind. The listener's pleasure is influenced by the extent to which the auditory data can be resolved in the simplest terms possible.

Genome-wide patterns of promoter sharing and co-expression in bovine skeletal muscle

BACKGROUND

Gene regulation by transcription factors (TF) is species, tissue and time specific. To better understand how the genetic code controls gene expression in bovine muscle we associated gene expression data from developing Longissimus thoracis et lumborum skeletal muscle with bovine promoter sequence information.

RESULTS

We created a highly conserved genome-wide promoter landscape comprising 87,408 interactions relating 333 TFs with their 9,242 predicted target genes (TGs). We discovered that the complete set of predicted TGs share an average of 2.75 predicted TF binding sites (TFBSs) and that the average co-expression between a TF and its predicted TGs is higher than the average co-expression between the same TF and all genes. Conversely, pairs of TFs sharing predicted TGs showed a co-expression correlation higher that pairs of TFs not sharing TGs. Finally, we exploited the co-occurrence of predicted TFBS in the context of muscle-derived functionally-coherent modules including cell cycle, mitochondria, immune system, fat metabolism, muscle/glycolysis, and ribosome. Our findings enabled us to reverse engineer a regulatory network of core processes, and correctly identified the involvement of E2F1, GATA2 and NFKB1 in the regulation of cell cycle, fat, and muscle/glycolysis, respectively.

CONCLUSION

The pivotal implication of our research is two-fold: (1) there exists a robust genome-wide expression signal between TFs and their predicted TGs in cattle muscle consistent with the extent of promoter sharing; and (2) this signal can be exploited to recover the cellular mechanisms underpinning transcription regulation of muscle structure and development in bovine. Our study represents the first genome-wide report linking tissue specific co-expression to co-regulation in a non-model vertebrate.

Activity, abundance, distribution and expression of Na+/K+-ATPase in the salt glands of Crocodylus porosus following chronic saltwater acclimation

Saltwater crocodiles, Crocodylus porosus, possess lingual salt glands which function to remove excess Na+ and Cl− accumulated as a consequence of living in salt water. Little is known about the nature of ion transport systems in C. porosus salt glands and how these systems respond to an osmotic challenge. In the present study, we examined the distribution and regulation of the Na+/K+-ATPase (NKA) pump, specifically the α-(catalytic) subunit in the salt glands of C. porosus chronically acclimated (6 months) to freshwater (FW) or 70% seawater (SW). We hypothesised that in the SW-acclimated C. porosus there would be an up-regulation of the abundance, activity and gene expression of the NKA transporter. NKA was immunolocalised to the lateral and basal membrane of secretory cells. As predicted, the NKA α-subunit was 2-fold more abundant in SW-acclimated C. porosus salt glands. NKA gene expression was also elevated in the salt glands of SW- vs FW-acclimated crocodiles. There was no increase in the specific activity of NKA in SW-acclimated animals and the in vitro rate of oxygen consumption by salt gland slices from SW-acclimated animals was not significantly different from that of FW-acclimated animals. The proportion of tissue oxygen consumption rate attributable to NKA activity was not different between SW- and FW-acclimated animals (approximately 50%). These data suggest that either chronic SW acclimation does not affect NKA in crocodile salt glands in the same manner as seen in other models or crocodiles possess the capacity to moderate NKA activity following prolonged exposure to SW.

Regulatory impact factors: unraveling the transcriptional regulation of complex traits from expression data

MOTIVATION

Although transcription factors (TF) play a central regulatory role, their detection from expression data is limited due to their low, and often sparse, expression. In order to fill this gap, we propose a regulatory impact factor (RIF) metric to identify critical TF from gene expression data.

RESULTS

To substantiate the generality of RIF, we explore a set of experiments spanning a wide range of scenarios including breast cancer survival, fat, gonads and sex differentiation. We show that the strength of RIF lies in its ability to simultaneously integrate three sources of information into a single measure: (i) the change in correlation existing between the TF and the differentially expressed (DE) genes; (ii) the amount of differential expression of DE genes; and (iii) the abundance of DE genes. As a result, RIF analysis assigns an extreme score to those TF that are consistently most differentially co-expressed with the highly abundant and highly DE genes (RIF1), and to those TF with the most altered ability to predict the abundance of DE genes (RIF2). We show that RIF analysis alone recovers well-known experimentally validated TF for the processes studied. The TF identified confirm the importance of PPAR signaling in adipose development and the importance of transduction of estrogen signals in breast cancer survival and sexual differentiation. We argue that RIF has universal applicability, and advocate its use as a promising hypotheses generating tool for the systematic identification of novel TF not yet documented as critical.

Skeletal muscle atrophy occurs slowly and selectively during prolonged aestivation in Cyclorana alboguttata (Günther 1867)

We investigated the effect of prolonged immobilisation of six and nine months duration on the morphology and antioxidant biochemistry of skeletal muscles in the amphibian aestivator Cyclorana alboguttata. We hypothesised that, in the event of atrophy occurring during aestivation,larger jumping muscles were more likely to be preserved over smaller non-jumping muscles. Whole muscle mass (g), muscle cross-sectional area (CSA)(μm2), water content (%) and myofibre number (per mm2) remained unchanged in the cruralis muscle after six to nine months of aestivation; however, myofibre area (μm2) was significantly reduced. Whole muscle mass, water content, myofibre number and myofibre CSA remained unchanged in the gastrocnemius muscle after six to nine months of aestivation. However, iliofibularis dry muscle mass, whole muscle CSA and myofibre CSA was significantly reduced during aestivation. Similarly,sartorius dry muscle mass, water content and whole muscle CSA was significantly reduced during aestivation. Endogenous antioxidants were maintained at control levels throughout aestivation in all four muscles. The results suggest changes to muscle morphology during aestivation may occur when lipid reserves have been depleted and protein becomes the primary fuel substrate for preserving basal metabolic processes. Muscle atrophy as a result of this protein catabolism may be correlated with locomotor function, with smaller non-jumping muscles preferentially used as a protein source during fasting over larger jumping muscles. Higher levels of endogenous antioxidants in the jumping muscles may confer a protective advantage against oxidative damage during aestivation; however, it is not clear whether they play a role during aestivation or upon resumption of normal metabolic activity.

Surviving the drought: burrowing frogs save energy by increasing mitochondrial coupling

During dormancy energy conservation is a key priority and as such dormant animals undergo a major metabolic depression to conserve their limited endogenous fuel supplies. Mitochondrial coupling efficiency, the efficiency with which mitochondria convert oxygen into ATP, significantly affects aerobic metabolism and thus to maximise energy savings during dormancy it has been hypothesised that coupling efficiency should increase. However, previous studies have shown coupling efficiency to be maintained or even to decrease. In this study we measured state 3 and state 4 mitochondrial respiration in the muscle of the burrowing frog, Cyclorana alboguttata and calculated the respiratory control ratio as a measure of coupling efficiency. After 7 months in aestivation, C. alboguttata significantly reduced oxygen consumption of isolated mitochondria by 83% and, unlike other dormant animals, the frogs appeared to decrease rates of proton leak to a greater extent than ATP synthesis, consistent with an increase in mitochondrial coupling efficiency. The significant energy savings observed at the mitochondrial level were reflected at higher levels of biological organisation, with tissue oxygen consumption depressed by as much as 81% and whole animal metabolic rate by 82%. Cyclorana alboguttata can survive in a dormant state for several years and we propose the hypothesis that energy efficiency is increased during aestivation.

A differential wiring analysis of expression data correctly identifies the gene containing the causal mutation

Transcription factor (TF) regulation is often post-translational. TF modifications such as reversible phosphorylation and missense mutations, which can act independent of TF expression level, are overlooked by differential expression analysis. Using bovine Piedmontese myostatin mutants as proof-of-concept, we propose a new algorithm that correctly identifies the gene containing the causal mutation from microarray data alone. The myostatin mutation releases the brakes on Piedmontese muscle growth by translating a dysfunctional protein. Compared to a less muscular non-mutant breed we find that myostatin is not differentially expressed at any of ten developmental time points. Despite this challenge, the algorithm identifies the myostatin ‘smoking gun’ through a coordinated, simultaneous, weighted integration of three sources of microarray information: transcript abundance, differential expression, and differential wiring. By asking the novel question “which regulator is cumulatively most differentially wired to the abundant most differentially expressed genes?” it yields the correct answer, “myostatin”. Our new approach identifies causal regulatory changes by globally contrasting co-expression network dynamics. The entirely data-driven ‘weighting’ procedure emphasises regulatory movement relative to the phenotypically relevant part of the network. In contrast to other published methods that compare co-expression networks, significance testing is not used to eliminate connections.

Evolution, development, and cancer are governed by regulatory circuits where the central nodes are transcription factors. Consequently, there is great interest in methods that can identify the causal mutation/perturbation responsible for any circuit rewiring. The most widely available high-throughput technology, the microarray, assays the transcriptome. However, many regulatory perturbations are post-transcriptional. This means that they are overlooked by traditional differential gene expression analysis. We hypothesised that by viewing biological systems as networks one could identify causal mutations and perturbations by examining those regulators whose position in the network changes the most. Using muscular myostatin mutant cattle as a proof-of-concept, we propose an analysis that succeeds based solely on microarray expression data from just 27 animals. Our analysis differs from competing network approaches in that we do not use significance testing to eliminate connections. All connections are contrasted, no matter how weak. Further, the identity of target genes is maintained throughout the analysis. Finally, the analysis is ‘weighted’ such that movement relative to the phenotypically most relevant part of the network is emphasised. By identifying the question to which myostatin is the answer, we present a comparison of network connectivity that is potentially generalisable.

Gene expression studies of developing bovine longissimus muscle from two different beef cattle breeds

BACKGROUND

The muscle fiber number and fiber composition of muscle is largely determined during prenatal development. In order to discover genes that are involved in determining adult muscle phenotypes, we studied the gene expression profile of developing fetal bovine longissimus muscle from animals with two different genetic backgrounds using a bovine cDNA microarray. Fetal longissimus muscle was sampled at 4 stages of myogenesis and muscle maturation: primary myogenesis (d 60), secondary myogenesis (d 135), as well as beginning (d 195) and final stages (birth) of functional differentiation of muscle fibers. All fetuses and newborns (total n = 24) were from Hereford dams and crossed with either Wagyu (high intramuscular fat) or Piedmontese (GDF8 mutant) sires, genotypes that vary markedly in muscle and compositional characteristics later in postnatal life.

RESULTS

We obtained expression profiles of three individuals for each time point and genotype to allow comparisons across time and between sire breeds. Quantitative reverse transcription-PCR analysis of RNA from developing longissimus muscle was able to validate the differential expression patterns observed for a selection of differentially expressed genes, with one exception. We detected large-scale changes in temporal gene expression between the four developmental stages in genes coding for extracellular matrix and for muscle fiber structural and metabolic proteins. FSTL1 and IGFBP5 were two genes implicated in growth and differentiation that showed developmentally regulated expression levels in fetal muscle. An abundantly expressed gene with no functional annotation was found to be developmentally regulated in the same manner as muscle structural proteins. We also observed differences in gene expression profiles between the two different sire breeds. Wagyu-sired calves showed higher expression of fatty acid binding protein 5 (FABP5) RNA at birth. The developing longissimus muscle of fetuses carrying the Piedmontese mutation shows an emphasis on glycolytic muscle biochemistry and a large-scale up-regulation of the translational machinery at birth. We also document evidence for timing differences in differentiation events between the two breeds.

CONCLUSION

Taken together, these findings provide a detailed description of molecular events accompanying skeletal muscle differentiation in the bovine, as well as gene expression differences that may underpin the phenotype differences between the two breeds. In addition, this study has highlighted a non-coding RNA, which is abundantly expressed and developmentally regulated in bovine fetal muscle.

Obese humans as economically designed feed converters: symmorphosis and low oxidative capacity skeletal muscle

Human obesity is considered a consequence of a thrifty or economic metabolism. In this hypothesis, we apply an established economic design theory, called symmorphosis, to help explain the known association between obesity and low oxidative capacity skeletal muscle. Symmorphosis reflects an engineering principle, and predicts that physiological systems are most economically designed when unnecessary spare capacity is eliminated. This is because the structural/functional adaptations accounting for spare capacity themselves bear energetic costs of construction, maintenance and load. As oxidation of feed energy occurs in mitochondria, and because skeletal muscle accounts for 30% of resting metabolism, we focus on skeletal muscle mitochondria. In the same way that the most economically designed elevator is supported by a cable that is strong enough, but not too strong, symmorphosis predicts that the most economically designed feed converters should have enough, but not too much mitochondrial oxidative (fuel burning) capacity. While ATP demand is clearly more efficiently met by oxidative (38 molecules of ATP) rather than glycolytic (2 molecules of ATP) metabolism, symmorphosis predicts that having excess oxidative capacity actually reduces feed efficiency. This inefficiency is manifest by having to maintain, ultimately using feed energy, the expensive inner mitochondrial proton gradient in the superfluous mitochondria. On this basis, we predict that established molecular controllers of mitochondrial biogenesis and oxidative capacity such as eNOS, SIN3 co-repressor, TFAM and PPARgamma may yield useful DNA markers and therapeutic targets for issues relating to frugal energetics, namely predisposition to obesity and starvation resilience.

The effects of saltwater acclimation on neurotransmitters in the lingual salt glands of the estuarine crocodile, Crocodylus porosus

INTRODUCTION

Most avian and reptilian salt glands display marked phenotypic plasticity when animals are exposed to hyperosmotic conditions. In addition, the activity of most salt glands is under considerable control by the nervous system and nerves containing cholinergic, adrenergic and peptidergic neurotransmitters have been identified in avian and reptilian salt gland tissues. The present study sought to determine whether the salt glands of the estuarine crocodile, Crocodylus porosus contain the peptidergic neurotransmitters SP, CGRP, VIP, and PACAP and the gaseous neurotransmitter, NO. In addition, we sought to determine whether there was any evidence for the adaptation of the C. porosus salt gland nervous system to hyperosmotic conditions.

METHODS

Salt glands from freshwater- and saltwater-acclimated C. porosus hatchlings were sectioned and examined immunohistochemically for neurotransmitters within the tissue.

RESULTS

Neurons containing SP, CGRP, VIP, PACAP and NO synthase were identified within C. porosus salt glands. There was no difference in the overall number (density) of neurons within SW-acclimated tissues when compared with FW-acclimated animals. However, there was a significant reduction in density of neurons containing SP and PACAP in SW-acclimated animals.

CONCLUSION

C. porosus salt glands display phenotypic plasticity following exposure to hyperosmotic conditions. In addition to cholinergic and adrenergic neurons, they contain a variety of peptidergic neurotransmitters and the gaseous neurotransmitter NO. Additionally, there appears to be some evidence of acclimation of the nervous system of C. porosus to hypersaline conditions, although the functional significance of these changes remains to be determined.

Skeletal muscle extracellular matrix remodelling after aestivation in the green striped burrowing frog, Cyclorana alboguttata

Connective tissue has recently been found to play a role in mediating mammalian skeletal muscle atrophy. We investigated connective tissue remodelling in the skeletal muscle of a species of the Australian burrowing frog, Cyclorana alboguttata. Despite being inactive whilst aestivating, the frog shows an inhibition of muscle atrophy. Connective tissue size and distribution was measured in histological sections of the cruralis muscle of control and aestivating C. alboguttata. Using a custom written software application we could detect no significant difference in any connective tissue morphological parameter between the two treatment groups. Biochemical measurements of gelatinase activity showed 2-fold higher activity in aestivating gastrocnemius muscle than in controls (p<0.001). We measured the messenger RNA transcript levels for C. alboguttata metalloproteinase 2 (MMP2) and tissue inhibitor of metalloproteinase 2 (TIMP2) in cruralis skeletal muscle using quantitative real-time PCR. The trend of reduced expression of the two genes in the aestivators did not meet statistical significance. This work indicates that aestivation in C. alboguttata leads to subtle and specific changes in some extracellular matrix remodelling factors. Their main impact is to maintain proportional representation of extracellular matrix components of skeletal muscle and therefore preserve the active frog phenotype.

Gene expression profiling of bovine skeletal muscle in response to and during recovery from chronic and severe undernutrition1

Gene expression profiles of LM from beef cattle that underwent significant postweaning undernutrition were studied using complementary DNA (cDNA) microarrays. After 114 d of undernutrition, the RNA from LM showed 2- to 6-fold less expression of many genes from the classes of muscle structural proteins, muscle metabolic enzymes, and extracellular matrix compared with animals on a rapid growth diet. The expression levels of these genes had mostly returned to pretreatment levels after 84 d of realimentation. The gene expression changes associated with undernutrition and BW loss showed an emphasis on downregulation of gene expression specific to fast-twitch fibers, typical of starving mammals, with a preferential atrophy of glycolytic fast-twitch fibers. We also identified a small group of genes that showed 2- to 5-fold elevated expression in LM after 114 d of undernutrition. Putative roles for these genes in atrophying skeletal muscle are regulation of myogenic differentiation (CSRP3), maintenance of mesenchymal stem cells (CYR61), modulation of membrane function (TM4SF2), prevention of oxidative damage (SESN1), and regulation of muscle protein degradation (SQSTM1). A significant increase in stearoyl-CoA desaturase (SCD) gene expression was observed in atrophying muscle, suggesting either that increased fatty acid synthesis is part of the tissue response to caloric restriction, or that SCD plays another role in energy metabolism in the mixed cellular environment of bovine skeletal muscle.