Modulation of matrix mineralization by Vwc2-like protein and its novel splicing isoforms

Genome-Wide Association Analysis of Reproductive Traits in Chinese Holstein Cattle

This study was to explore potential SNP loci for reproductive traits in Chinese Holstein cattle and identify candidate genes. Genome-wide Association Study based on mixed linear model was performed on 643 Holstein cattle using GeneSeek Bovine 50 K SNP chip. Our results detected forty significant SNP loci after Bonferroni correction. We identified five genes (VWC2L, STAT1, PPP3CA, LDB3, and CTNNA3) as being associated with pregnancy ratio of young cows, five genes (PAEP, ACOXL, EPAS1, GLRB, and MARVELD1) as being associated with pregnancy ratio of adult cows, and nine genes (PDE1B, SLCO1A2, ARHGAP26, ADAM10, APBB1, MON1B, COQ9, CDC42BPB, MARVELD1, and HPSE2) as being associated with daughter pregnancy rate. Our study may provide valuable insights into identifying genes related to reproductive traits and help promote the application of molecular breeding in dairy cows.

Estimation of direct and maternal genetic effects and annotation of potential candidate genes for weight and meat quality traits in a genotyped outdoor dual-purpose cattle breed

With regard to potential applications of genomic selection in small numbered breeds, we evaluated genomic models and focused on potential candidate gene annotations for weight and meat quality traits in the local Rotes Höhenvieh (RHV) breed. Traits included 6,003 birth weights (BWT), 5,719 200 d-weights (200dw), 4,594 365 d-weights (365dw), and 547 records for intramuscular fat content (IMF). A total of 581,304 SNP from 370 genotyped cattle with phenotypic records were included in genomic analyses. Model evaluations focused on single- and multiple-trait models with direct and with direct and maternal genetic effects. Genetic relationship matrices were based on pedigree (A-matrix), SNP markers (G-matrix), or both (H-matrix). Genome-wide association studies (GWASs) were carried out using linear mixed models to identify potential candidate genes for the traits of interest. De-regressed proofs (DRP) for direct and maternal genetic components were used as pseudo-phenotypes in the GWAS. Accuracies of direct breeding values were higher from models based on G or on H compared to A. Highest accuracies (> 0.89) were obtained for IMF with multiple-trait models using the G-matrix. Direct heritabilities with maternal genetic effects ranged from 0.62 to 0.66 for BWT, from 0.45 to 0.55 for 200dW, from 0.40 to 0.44 for 365dW, and from 0.48 to 0.75 for IMF. Maternal heritabilities for BWT, 200dW, and 365dW were in a narrow range from 0.21 to 0.24, 0.24 to 0.27, and 0.21 to 0.25, respectively, and from 0.25 to 0.65 for IMF. Direct genetic correlations among body weight traits were positive and favorable, and very similar from different models but showed a stronger variation with 0.31 (A), −0.13 (G), and 0.45 (H) between BWT and IMF. In gene annotations, we identified 6, 3, 1, and 6 potential candidate genes for direct genetic effect on BWT, 200dW, 365dW, and IMF traits, respectively. Regarding maternal genetic effects, four (SHROOM3, ZNF609, PECAM1, and TEX2) and two (TMEM182 and SEC11A) genes were detected as potential candidate genes for BWT and 365dW, respectively. Potential candidate genes for maternal effect on IMF were GRHL2, FGA, FGB, and CTNNA3. As the most important finding from a practical breeding perspective, a small number of genotyped RHV cattle enabled accurate breeding values for high heritability IMF.

From Mendel laws to whole genetic association study to decipher the swine mulefoot phenotype

The swine mulefoot (SM) is a rare condition characterized by a non-cloven hoof due to the partial or total fusion of the phalanges. No comprehensive study has been conducted to identify associated markers with this phenotype until now. We aimed to characterize the association between SNP and the mulefoot phenotype using a Genome-Wide Association Study (GWAS). An experimental population was produced using a half-sib mating where the male had the mulefoot phenotype and the females (n = 6) had cloven hoofs. The cross resulted in 27 (47%) animals with the mulefoot characteristic and 30 (53%) normal animals, indicating the possible dominant gene action. Animals were further genotyped using the Illumina PorcineSNP50k BeadChip, and SNPs were tested for associations. Twenty-nine SNPs located on the SSC15, SSC4, and SSCX were associated with the mulefoot phenotype (p-value <5 × 10^-5). Six markers were found in the intronic regions of VWC2L, CATIP, PDK3, PCYT1B, and POLA1 genes. The marker rs81277626, on SSC15:116,886,110 bp, is located in the Von Willebrand Factor C Domain (VWC2L), a possible functional candidate gene. The VWC2L is part of a biological process involved with the bone morphogenetic protein (BMP) signaling pathway, previously associated with syndactyly in other species. In conclusion, the identified markers suggest the involvement of the VWC2L gene in the SM phenotype in this population.

The spatial network of skeletal proteins in a stony coral

Coral skeletons are materials composed of inorganic aragonitic fibres and organic molecules including proteins, sugars and lipids that are highly organized to form a solid biomaterial upon which the animals live. The skeleton contains tens of proteins, all of which are encoded in the animal genome and secreted during the biomineralization process. While recent advances are revealing the functions and evolutionary history of some of these proteins, how they are spatially arranged in the skeleton is unknown. Using a combination of chemical cross-linking and high-resolution tandem mass spectrometry, we identify, for the first time, the spatial interactions of the proteins embedded within the skeleton of the stony coral Stylophora pistillata. Our subsequent network analysis revealed that several coral acid-rich proteins are invariably associated with carbonic anhydrase(s), alpha-collagen, cadherins and other calcium-binding proteins. These spatial arrangements clearly show that protein-protein interactions in coral skeletons are highly coordinated and are key to understanding the formation and persistence of coral skeletons through time.

FAM20C directly binds to and phosphorylates Periostin

It is widely accepted that FAM20C functions as a Golgi casein kinase and has large numbers of kinase substrates within the secretory pathway. It has been previously reported that FAM20C is required for maintenance of healthy periodontal tissues. However, there has been no report that any extracellular matrix molecules expressed in periodontal tissues are indeed substrates of FAM20C. In this study, we sought to identify the binding partner(s) of FAM20C. FAM20C wild-type (WT) and its kinase inactive form D478A proteins were generated. These proteins were electrophoresed and the Coomassie Brilliant Blue (CBB)-positive bands were analyzed to identify FAM20C-binding protein(s) by Mass Spectrometry (MS) analysis. Periostin was found by the analysis and the binding between FAM20C and Periostin was investigated in cell cultures and in vitro. We further determined the binding region(s) within Periostin responsible for FAM20C-binding. Immunolocalization of FAM20C and Periostin was examined using mouse periodontium tissues by immunohistochemical analysis. In vitro kinase assay was performed using Periostin and FAM20C proteins to see whether FAM20C phosphorylates Periostin in vitro. We identified Periostin as one of FAM20C-binding proteins by MS analysis. Periostin interacted with FAM20C in a kinase-activity independent manner and the binding was direct in vitro. We further identified the binding domain of FAM20C in Periostin, which was mapped within Fasciclin (Fas) I domain 1-4 of Periostin. Immunolocalization of FAM20C was observed in periodontal ligament (PDL) extracellular matrix where that of Periostin was also immunostained in murine periodontal tissues. FAM20C WT, but not D478A, phosphorylated Periostin in vitro. Consistent with the overlapped expression pattern of FAM20C and Periostin, our data demonstrate for the first time that Periostin is a direct FAM20C-binding partner and that FAM20C phosphorylates Periostin in vitro.

Genomic inbreeding trends, influential sire lines and selection in the global Thoroughbred horse population

The Thoroughbred horse is a highly valued domestic animal population under strong selection for athletic phenotypes. Here we present a high resolution genomics-based analysis of inbreeding in the population that may form the basis for evidence-based discussion amid concerns in the breeding industry over the increasing use of small numbers of popular sire lines, which may accelerate a loss of genetic diversity. In the most comprehensive globally representative sample of Thoroughbreds to-date (n = 10,118), including prominent stallions (n = 305) from the major bloodstock regions of the world, we show using pan-genomic SNP genotypes that there has been a highly significant decline in global genetic diversity during the last five decades (FIS R2 = 0.942, P = 2.19 × 10-13; FROH R2 = 0.88, P = 1.81 × 10-10) that has likely been influenced by the use of popular sire lines. Estimates of effective population size in the global and regional populations indicate that there is some level of regional variation that may be exploited to improve global genetic diversity. Inbreeding is often a consequence of selection, which in managed animal populations tends to be driven by preferences for cultural, aesthetic or economically advantageous phenotypes. Using a composite selection signals approach, we show that centuries of selection for favourable athletic traits among Thoroughbreds acts on genes with functions in behaviour, musculoskeletal conformation and metabolism. As well as classical selective sweeps at core loci, polygenic adaptation for functional modalities in cardiovascular signalling, organismal growth and development, cellular stress and injury, metabolic pathways and neurotransmitters and other nervous system signalling has shaped the Thoroughbred athletic phenotype. Our results demonstrate that genomics-based approaches to identify genetic outcrosses will add valuable objectivity to augment traditional methods of stallion selection and that genomics-based methods will be beneficial to actively monitor the population to address the marked inbreeding trend.

Polygenic risk score for disability and insights into disability-related molecular mechanisms

Late life disability is a highly devastating condition affecting 20% or more of persons aged 65 years and older in the USA; it is an important determinant of acute medical and long-term care costs which represent a growing burden on national economies. Disability is a multifactorial trait that contributes substantially to decline of health/wellbeing. Accordingly, gaining insights into the genetics of disability could help in identifying molecular mechanisms of this devastating condition and age-related processes contributing to a large fraction of specific geriatric conditions, concordantly with geroscience. We performed a genome-wide association study of disability in a sample of 24,068 subjects from five studies with 12,550 disabled individuals. We identified 30 promising disability-associated polymorphisms in 19 loci at p < 10^-4; four of them attained suggestive significance, p < 10^-5. In contrast, polygenic risk scores aggregating effects of minor alleles of independent SNPs that were adversely or beneficially associated with disability showed highly significant associations in meta-analysis, p = 3.13 × 10^-45 and p = 5.60 × 10^-23, respectively, and were replicated in each study. The analysis of genetic pathways, related diseases, and biological functions supported the connections of genes for the identified SNPs with disabling and age-related conditions primarily through oxidative/nitrosative stress, inflammatory response, and ciliary signaling. We identified musculoskeletal system development, maintenance, and regeneration as important components of gene functions. The beneficial and adverse gene sets may be differently implicated in the development of musculoskeletal-related disability with the beneficial set characterized, e.g., by regulation of chondrocyte proliferation and bone formation, and the adverse set by inflammation and bone loss.

AMPA receptors and their minions: auxiliary proteins in AMPA receptor trafficking

To correctly transfer information, neuronal networks need to continuously adjust their synaptic strength to extrinsic stimuli. This ability, termed synaptic plasticity, is at the heart of their function and is, thus, tightly regulated. In glutamatergic neurons, synaptic strength is controlled by the number and function of AMPA receptors at the postsynapse, which mediate most of the fast excitatory transmission in the central nervous system. Their trafficking to, at, and from the synapse, is, therefore, a key mechanism underlying synaptic plasticity. Intensive research over the last 20 years has revealed the increasing importance of interacting proteins, which accompany AMPA receptors throughout their lifetime and help to refine the temporal and spatial modulation of their trafficking and function. In this review, we discuss the current knowledge about the roles of key partners in regulating AMPA receptor trafficking and focus especially on the movement between the intracellular, extrasynaptic, and synaptic pools. We examine their involvement not only in basal synaptic function, but also in Hebbian and homeostatic plasticity. Included in our review are well-established AMPA receptor interactants such as GRIP1 and PICK1, the classical auxiliary subunits TARP and CNIH, and the newest additions to AMPA receptor native complexes.

VWC2 Increases Bone Formation Through Inhibiting Activin Signaling

By a bioinformatics approach, we have identified a novel cysteine knot protein member, VWC2 (von Willebrand factor C domain containing 2) previously known as Brorin. Since Brorin has been proposed to function as a bone morphogenetic protein (BMP) antagonist, we investigated the binding of Brorin/VWC2 to several BMPs; however, none of the BMPs tested were bound to VWC2. Instead, the βA subunit of activin was found as a binding partner among transforming growth factor (TGF)-β superfamily members. Here, we show that Vwc2 gene expression is temporally upregulated early in osteoblast differentiation, VWC2 protein is present in bone matrix, and localized at osteoblasts/osteocytes. Activin A-induced Smad2 phosphorylation was inhibited in the presence of exogenous VWC2 in MC3T3-E1 osteoblast cell line and primary osteoblasts. The effect of VWC2 on ex vivo cranial bone organ cultures treated with activin A was investigated, and bone morphometric parameters decreased by activin A were restored with VWC2. When we further investigated the biological mechanism how VWC2 inhibited the effects of activin A on bone formation, we found that the effects of activin A on osteoblast cell growth, differentiation, and mineralization were reversed by VWC2. Taken together, a novel secretory protein, VWC2 promotes bone formation by inhibiting Activin-Smad2 signaling pathway.

Global differential gene expression in the pituitary gland and the ovaries of pre- and postpubertal Brahman heifers

To understand genes, pathways, and networks related to puberty, we characterized the transcriptome of two tissues: the pituitary gland and ovaries. Samples were harvested from pre- and postpubertal Brahman heifers (same age group). Brahman heifers () are older at puberty compared with , a productivity issue. With RNA sequencing, we identified differentially expressed (DEx) genes and important transcription factors (TF) and predicted coexpression networks. The number of DEx genes detected in the pituitary gland was 284 ( < 0.05), and was the most DEx gene (fold change = 4.12, = 0.01). The gene promotes bone mineralization through transforming growth factor-β (TGFβ) signaling. Further studies of the link between bone mineralization and puberty could target . In ovaries, 3,871 genes were DEx ( < 0.05). Four highly DEx genes were noteworthy for their function: (a γ-aminobutyric acid [GABA] transporter), (), and () and its receptor . These genes had higher ovarian expression in postpubertal heifers. The GABA and its receptors and transporters were expressed in the ovaries of many mammals, suggesting a role for this pathway beyond the brain. The pathway has been known to influence the timing of puberty in rats, via modulation of GnRH. The effects of at the hypothalamus, pituitary gland, and ovaries have been documented. and its receptors are known factors in the release of GnRH, similar to and GABA, although their roles in ovarian tissue are less clear. Pathways previously related to puberty such as TGFβ signaling ( = 6.71 × 10), Wnt signaling ( = 4.1 × 10), and peroxisome proliferator-activated receptor (PPAR) signaling ( = 4.84 × 10) were enriched in our data set. Seven genes were identified as key TF in both tissues: , , , , , , and a novel gene. An ovarian subnetwork created with TF and significant ovarian DEx genes revealed five zinc fingers as regulators: , , , , and . Recent work of hypothalamic gene expression also pointed to zinc fingers as TF for bovine puberty. Although some zinc fingers may be ubiquitously expressed, the identification of DEx genes in common across tissues points to key regulators of puberty. The hypothalamus and pituitary gland had eight DEx genes in common. The hypothalamus and ovaries had 89 DEx genes in common. The pituitary gland and ovaries had 48 DEx genes in common. Our study confirmed the complexity of puberty and suggested further investigation on genes that code zinc fingers.

Molecular Therapy for Degenerative Disc Disease: Clues from Secretome Analysis of the Notochordal Cell-Rich Nucleus Pulposus

Degenerative disc disease (DDD) is associated with spinal pain often leading to long-term disability. However, the non-chondrodystrophic canine intervertebral disc is protected from the development of DDD, ostensibly due to its retention of notochordal cells (NC) in the nucleus pulposus (NP). In this study, we hypothesized that secretome analysis of the NC-rich NP will lead to the identification of key proteins that delay the onset of DDD. Using mass-spectrometry, we identified 303 proteins including components of TGFβ- and Wnt-signaling, anti-angiogeneic factors and proteins that inhibit axonal ingrowth in the bioactive fractions of serum free, notochordal cell derived conditioned medium (NCCM). Ingenuity Pathway Analysis revealed TGFβ1 and CTGF as major hubs in protein interaction networks. In vitro treatment with TGFβ1 and CTGF promoted the synthesis of healthy extra-cellular matrix proteins, increased cell proliferation and reduced cell death in human degenerative disc NP cells. A single intra-discal injection of recombinant TGFβ1 and CTGF proteins in a pre-clinical rat-tail disc injury model restored the NC and stem cell rich NP. In conclusion, we demonstrate the potential of TGFβ1 and CTGF to mitigate the progression of disc degeneration and the potential use of these molecules in a molecular therapy to treat the degenerative disc.

Genome Wide Association Analysis Reveals New Production Trait Genes in a Male Duroc Population

In this study, 796 male Duroc pigs were used to identify genomic regions controlling growth traits. Three production traits were studied: food conversion ratio, days to 100 KG, and average daily gain, using a panel of 39,436 single nucleotide polymorphisms. In total, we detected 11 genome-wide and 162 chromosome-wide single nucleotide polymorphism trait associations. The Gene ontology analysis identified 14 candidate genes close to significant single nucleotide polymorphisms, with growth-related functions: six for days to 100 KG (WT1, FBXO3, DOCK7, PPP3CA, AGPAT9, and NKX6-1), seven for food conversion ratio (MAP2, TBX15, IVL, ARL15, CPS1, VWC2L, and VAV3), and one for average daily gain (COL27A1). Gene ontology analysis indicated that most of the candidate genes are involved in muscle, fat, bone or nervous system development, nutrient absorption, and metabolism, which are all either directly or indirectly related to growth traits in pigs. Additionally, we found four haplotype blocks composed of suggestive single nucleotide polymorphisms located in the growth trait-related quantitative trait loci and further narrowed down the ranges, the largest of which decreased by ~60 Mb. Hence, our results could be used to improve pig production traits by increasing the frequency of favorable alleles via artificial selection.