Next-generation sequencing identifies equine cartilage and subchondral bone miRNAs and suggests their involvement in osteochondrosis physiopathology.

BACKGROUND: MicroRNAs (miRNAs) are an abundant class of small single-stranded non-coding RNA molecules ranging from 18 to 24 nucleotides. They negatively regulate gene expression at the post-transcriptional level and play key roles in many biological processes, including skeletal development and cartilage maturation. In addition, miRNAs involvement in osteoarticular diseases has been proved and some of them were identified as suitable biomarkers for pathological conditions. Equine osteochondrosis (OC) is one of the most prevalent juvenile osteoarticular disorders in horses and represents a major concern for animal welfare and economic reasons. Its etiology and pathology remain controversial and biological pathways as well as molecular mechanisms involved in the physiopathology are still unclear. This study aims to investigate the potential role of miRNAs in equine osteochondrosis (OC) physiopathology.Short-read NGS technology (SOLID™, Life Technologies) was used to establish a comprehensive repertoire of miRNA expressed in either equine cartilage or subchondral bone. Undamaged cartilage and subchondral bone samples from healthy (healthy samples) and OC-affected (predisposed samples) 10-month Anglo-Arabian foals were analysed. Samples were also subjected or not to an experimental mechanical loading to evaluate the role of miRNAs in the regulation of mechano-transduction pathways. Predicted targets of annotated miRNAs were identified using miRmap. RESULTS: Epiphyseal cartilage and subchondral bone miRNome were defined, including about 300 new miRNAs. Differentially expressed miRNAs were identified between bone and cartilage from healthy and OC foals, as well as after an experimental mechanical loading. In cartilage, functional annotation of their predicted targets suggests a role in the maintenance of cartilage integrity through the control of cell cycle and differentiation, energy production and metabolism as well as extracellular matrix structure and dynamics. In bone, miRNA predicited targets were associated with osteoblasts and osteoclasts differentiation, though the regulation of energy production, vesicle transport and some growth factor signaling pathways. CONCLUSION: Taken together, our results suggest a role of miRNAs in equine OC physiopathology and in the cellular response to biomechanical stress in cartilage and bone. In silico target prediction and functional enrichment analysis provides new insight into OC molecular physiopathology.

Omics technologies provide new insights into the molecular physiopathology of equine osteochondrosis.

BACKGROUND: Osteochondrosis (OC(D)) is a juvenile osteo-articular disorder affecting several mammalian species. In horses, OC(D) is considered as a multifactorial disease and has been described as a focal disruption of endochondral ossification leading to the development of osteoarticular lesions. Nevertheless, OC(D) physiopathology is poorly understood. Affected horses may present joint swelling, stiffness and lameness. Thus, OC(D) is a major concern for the equine industry. Our study was designed as an integrative approach using omics technologies for the identification of constitutive defects in epiphyseal cartilage and/or subchondral bone associated with the development of primary lesions to further understand OC(D) pathology. This study compared samples from non-affected joints (hence lesion-free) from OC(D)-affected foals (n = 5, considered predisposed samples) with samples from OC-free foals (n = 5) considered as control samples. Consequently, results are not confounded by changes associated with the evolution of the lesion, but focus on altered constitutive molecular mechanisms. Comparative proteomics and micro computed tomography analyses were performed on predisposed and OC-free bone and cartilage samples. Metabolomics was also performed on synovial fluid from OC-free, OC(D)-affected and predisposed joints. RESULTS: Two lesion subtypes were identified: OCD (lesion with fragment) and OC (osteochondral defects). Modulated proteins were identified using omics technologies (2-DE proteomics) in cartilage and bone from affected foals compare to OC-free foals. These were associated with cellular processes including cell cycle, energy production, cell signaling and adhesion as well as tissue-specific processes such as chondrocyte maturation, extracellular matrix and mineral metabolism. Of these, five had already been identified in synovial fluid of OC-affected foals: ACTG1 (actin, gamma 1), albumin, haptoglobin, FBG (fibrinogen beta chain) and C4BPA (complement component 4 binding protein, alpha). CONCLUSION: This study suggests that OCD lesions may result from a cartilage defect whereas OC lesions may be triggered by both bone and cartilage defects, suggesting that different molecular mechanisms responsible for the equine osteochondrosis lesion subtypes and predisposition could be due to a defect in both bone and cartilage. This study will contribute to refining the definition of OC(D) lesions and may improve diagnosis and development of therapies for horses and other species, including humans.

High-throughput sequencing analyses of XX genital ridges lacking FOXL2 reveal DMRT1 up-regulation before SOX9 expression during the sex-reversal process in goats.

FOXL2 loss of function in goats leads to the early transdifferentiation of ovaries into testes, then to the full sex reversal of XX homozygous mutants. By contrast, Foxl2 loss of function in mice induces an arrest of follicle formation after birth, followed by complete female sterility. In order to understand the molecular role of FOXL2 during ovarian differentiation in the goat species, putative FOXL2 target genes were determined at the earliest stage of gonadal sex-specific differentiation by comparing the mRNA profiles of XX gonads expressing the FOXL2 protein or not. Of these 163 deregulated genes, around two-thirds corresponded to testicular genes that were up-regulated when FOXL2 was absent, and only 19 represented female-associated genes, down-regulated in the absence of FOXL2. FOXL2 should therefore be viewed as an antitestis gene rather than as a female-promoting gene. In particular, the key testis-determining gene DMRT1 was found to be up-regulated ahead of SOX9, thus suggesting in goats that SOX9 primary up-regulation may require DMRT1. Overall, our results equated to FOXL2 being an antitestis gene, allowing us to propose an alternative model for the sex-determination process in goats that differs slightly from that demonstrated in mice.

Involvement of mitochondrial dysfunction and ER-stress in the physiopathology of equine osteochondritis dissecans (OCD).

Osteochondrosis (OC) is a developmental bone disorder affecting several mammalian species including the horse. Equine OC is described as a focal disruption of endochondral ossification, leading to osteochondral lesions (osteochondritis dissecans, OCD) that may release free bodies within the joint. OCD lesions trigger joint swelling, stiffness and lameness and affects about 30% of the equine population. OCD is considered as multifactorial but its physiopathology is still poorly understood and genes involved in genetic predisposition are still unknown. Our study compared two healthy and two OC-affected 18-month-old French Trotters diagnosed with OCD lesions at the intermediate ridge of the distal tibia. A comparative shot-gun proteomic analysis of non-wounded cartilage and sub-chondral bone from healthy (healthy samples) and OC-affected foals (predisposed samples) identified 83 and 53 modulated proteins, respectively. These proteins are involved in various biological pathways including matrix structure and maintenance, protein biosynthesis, folding and transport, mitochondrial activity, energy and calcium metabolism. Transmission electron microscopy revealed typical features of mitochondrial swelling and ER-stress, such as large, empty mitochondria, and hyper-dilated rough endoplasmic reticulum, in the deep zone of both OC lesions and predisposed cartilage. Abnormal fibril organization surrounding chondrocytes and abnormal features at the ossification front were also observed. Combining these findings with quantitative trait loci and whole genome sequencing results identified about 140 functional candidate genes carrying putative damaging mutations in 30 QTL regions. In summary, our study suggests that OCD lesions may result from defective hypertrophic terminal differentiation associated with mitochondrial dysfunction and ER-stress, leading to impaired cartilage and bone biomechanical properties, making them prone to fractures. In addition, 11 modulated proteins and several candidate mutations located in QTL regions were identified, bringing new insight into the molecular physiopathology and genetic basis of OCD.

A method for proteomic analysis of equine subchondral bone and epiphyseal cartilage.

Proteomic analyses of cartilage and, to a lesser extent, of bone have long been impaired because of technical challenges related to their structure and biochemical properties. We have developed a unified method based on phenol extraction, 2DE, silver staining, and subsequent LC-MS/MS. This method proved to be efficient to characterize the proteome of equine cartilage and bone samples collected in vivo. Since proteins from several cellular compartments could be recovered, our procedure is mainly suitable for in situ molecular physiology studies focused on the cellular content of chondrocytes, osteoblasts, and osteoclasts as well as that of the extracellular matrix, with the exception of proteoglycans. Our method alleviates some drawbacks of cell culture that can mask physiological differences, as well as reduced reproducibility due to fractionation. Proteomic comparative studies between cartilage and bone samples from healthy and affected animals were thus achieved successfully. This achievement will contribute to increasing knowledge on the molecular mechanisms underlying the physiopathology of numerous osteoarticular diseases in horses and in humans.

Mapping of quantitative trait loci affecting classical scrapie incubation time in a population comprising several generations of scrapie-infected sheep.

Although susceptibility to scrapie is largely controlled by the PrP gene, the role of other genes that affect scrapie resistance in sheep is now confirmed. Following the detection of quantitative trait loci (QTL) on chromosomes 6 and 18 in a half-sib family with an ARQ/VRQ susceptible PrP genotype, the whole pedigree of a naturally infected flock was investigated to confirm these QTL regions in different PrP genotypes. The present study has allowed us to confirm the QTL on chromosome 18, and to demonstrate the QTL effects in several PrP genotypes.

High-resolution comparative mapping among man, cattle and mouse suggests a role for repeat sequences in mammalian genome evolution.

BACKGROUND: Comparative mapping provides new insights into the evolutionary history of genomes. In particular, recent studies in mammals have suggested a role for segmental duplication in genome evolution. In some species such as Drosophila or maize, transposable elements (TEs) have been shown to be involved in chromosomal rearrangements. In this work, we have explored the presence of interspersed repeats in regions of chromosomal rearrangements, using an updated high-resolution integrated comparative map among cattle, man and mouse. RESULTS: The bovine, human and mouse comparative autosomal map has been constructed using data from bovine genetic and physical maps and from FISH-mapping studies. We confirm most previous results but also reveal some discrepancies. A total of 211 conserved segments have been identified between cattle and man, of which 33 are new segments and 72 correspond to extended, previously known segments. The resulting map covers 91% and 90% of the human and bovine genomes, respectively. Analysis of breakpoint regions revealed a high density of species-specific interspersed repeats in the human and mouse genomes. CONCLUSION: Analysis of the breakpoint regions has revealed specific repeat density patterns, suggesting that TEs may have played a significant role in chromosome evolution and genome plasticity. However, we cannot rule out that repeats and breakpoints accumulate independently in the few same regions where modifications are better tolerated. Likewise, we cannot ascertain whether increased TE density is the cause or the consequence of chromosome rearrangements. Nevertheless, the identification of high density repeat clusters combined with a well-documented repeat phylogeny should highlight probable breakpoints, and permit their precise dating. Combining new statistical models taking the present information into account should help reconstruct ancestral karyotypes.

LIPH expression in skin and hair follicles of normal coat and Rex rabbits.

Natural mutations in the LIPH gene were shown to be responsible for hair growth defects in humans and for the rex short hair phenotype in rabbits. In this species, we identified a single nucleotide deletion in LIPH (1362delA) introducing a stop codon in the C-terminal region of the protein. We investigated the expression of LIPH between normal coat and rex rabbits during critical fetal stages of hair follicle genesis, in adults and during hair follicle cycles. Transcripts were three times less expressed in both fetal and adult stages of the rex rabbits than in normal rabbits. In addition, the hair growth cycle phases affected the regulation of the transcription level in the normal and mutant phenotypes differently. LIPH mRNA and protein levels were higher in the outer root sheath (ORS) than in the inner root sheath (IRS), with a very weak signal in the IRS of rex rabbits. In vitro transfection shows that the mutant protein has a reduced lipase activity compared to the wild type form. Our results contribute to the characterization of the LIPH mode of action and confirm the crucial role of LIPH in hair production.

Overexpression of miR-30b in the developing mouse mammary gland causes a lactation defect and delays involution.

BACKGROUND: MicroRNA (miRNA) are negative regulators of gene expression, capable of exerting pronounced influences upon the translation and stability of mRNA. They are potential regulators of normal mammary gland development and of the maintenance of mammary epithelial progenitor cells. This study was undertaken to determine the role of miR-30b on the establishment of a functional mouse mammary gland. miR-30b is a member of the miR-30 family, composed of 6 miRNA that are highly conserved in vertebrates. It has been suggested to play a role in the differentiation of several cell types. METHODOLOGY/PRINCIPAL FINDINGS: The expression of miR-30b was found to be regulated during mammary gland development. Transgenic mice overexpressing miR-30b in mammary epithelial cells were used to investigate its role. During lactation, mammary histological analysis of the transgenic mice showed a reduction in the size of alveolar lumen, a defect of the lipid droplets and a growth defect of pups fed by transgenic females. Moreover some mammary epithelial differentiated structures persisted during involution, suggesting a delay in the process. The genes whose expression was affected by the overexpression of miR-30b were characterized by microarray analysis. CONCLUSION/SIGNIFICANCE: Our data suggests that miR-30b is important for the biology of the mammary gland and demonstrates that the deregulation of only one miRNA could affect lactation and involution.

A deletion in exon 9 of the LIPH gene is responsible for the rex hair coat phenotype in rabbits (Oryctolagus cuniculus).

The fur of common rabbits is constituted of 3 types of hair differing in length and diameter while that of rex animals is essentially made up of amazingly soft down-hair. Rex short hair coat phenotypes in rabbits were shown to be controlled by three distinct loci. We focused on the "r1" mutation which segregates at a simple autosomal-recessive locus in our rabbit strains. A positional candidate gene approach was used to identify the rex gene and the corresponding mutation. The gene was primo-localized within a 40 cM region on rabbit chromosome 14 by genome scanning families of 187 rabbits in an experimental mating scheme. Then, fine mapping refined the region to 0.5 cM (Z = 78) by genotyping an additional 359 offspring for 94 microsatellites present or newly generated within the first defined interval. Comparative mapping pointed out a candidate gene in this 700 kb region, namely LIPH (Lipase Member H). In humans, several mutations in this major gene cause alopecia, hair loss phenotypes. The rabbit gene structure was established and a deletion of a single nucleotide was found in LIPH exon 9 of rex rabbits (1362delA). This mutation results in a frameshift and introduces a premature stop codon potentially shortening the protein by 19 amino acids. The association between this deletion and the rex phenotype was complete, as determined by its presence in our rabbit families and among a panel of 60 rex and its absence in all 60 non-rex rabbits. This strongly suggests that this deletion, in a homozygous state, is responsible for the rex phenotype in rabbits.

Polymorphisms in the HSP90AA1 5' flanking region are associated with scrapie incubation period in sheep.

Susceptibility to scrapie is mainly controlled by point mutations at the PRNP locus. However, additional quantitative trait loci (QTL) have been identified across the genome including a region in OAR18. The gene which encodes the inducible form of the cytoplasmic Hsp90 chaperone (HSP90AA1) maps within this region and seems to be associated with the resistance/susceptibility to scrapie in sheep. Here, we have analyzed several polymorphisms which were previously described in the ovine HSP90AA1 5' flanking region and in intron 10 in two naturally scrapie infected Romanov sheep populations. First, we have studied 58 ARQ/VRQ animals pertaining to the sire family where the QTL influencing scrapie incubation period in OAR18 was detected. We have found a significant association between polymorphisms localized at -660 and -528 in the HSP90AA1 5' flanking region and the scrapie incubation period. These two polymorphisms have also been studied in a second sample constituted by 62 VRQ/VRQ sheep showing an extreme incubation period. Results are concordant with the first dataset. Finally, we have studied the HSP90AA1 expression in scrapie and control animals (N = 41) with different HSP90AA1 genotypes by real time PCR on blood samples. The HSP90AA1 expression rate was equivalent in CC(-600)AA(-528) and CG(-600)AG(-528) scrapie resistant animals (ARR/ARR) and was higher in their CC(-600)AA(-528) than in their CG(-600)AG(-528) scrapie susceptible counterparts (VRQ/VRQ). Our results support the hypothesis that the ovine HSP90AA1 gene acts as a modulator of scrapie susceptibility, contributing to the observed differences in the incubation period of scrapie infected animals with the same PRNP genotype.

The prion or the related Shadoo protein is required for early mouse embryogenesis.

The prion protein PrP has a key role in transmissible spongiform encephalopathies but its biological function remains largely unknown. Recently, a related protein, Shadoo, was discovered. Its biological properties and brain distribution partially overlap that of PrP. We report that the Shadoo-encoding gene knockdown in PrP-knockout mouse embryos results in a lethal phenotype, occurring between E8 and E11, not observed on the wild-type genetic background. It reveals that these two proteins play a shared, crucial role in mammalian embryogenesis, explaining the lack of severe phenotype in PrP-knockout mammals, an appreciable step towards deciphering the biological role of this protein family.

Impact of strong selection for the PrP major gene on genetic variability of four French sheep breeds(Open Access publication).

Effective selection on the PrP gene has been implemented since October 2001 in all French sheep breeds. After four years, the ARR "resistant" allele frequency increased by about 35% in young males. The aim of this study was to evaluate the impact of this strong selection on genetic variability. It is focussed on four French sheep breeds and based on the comparison of two groups of 94 animals within each breed: the first group of animals was born before the selection began, and the second, 3-4 years later. Genetic variability was assessed using genealogical and molecular data (29 microsatellite markers). The expected loss of genetic variability on the PrP gene was confirmed. Moreover, among the five markers located in the PrP region, only the three closest ones were affected. The evolution of the number of alleles, heterozygote deficiency within population, expected heterozygosity and the Reynolds distances agreed with the criteria from pedigree and pointed out that neutral genetic variability was not much affected. This trend depended on breed, i.e. on their initial states (population size, PrP frequencies) and on the selection strategies for improving scrapie resistance while carrying out selection for production traits.

Positional cloning of the PIS mutation in goats and its impact on understanding mammalian sex-differentiation.

In goats, the PIS (polled intersex syndrome) mutation is responsible for both the absence of horns in males and females and sex-reversal affecting exclusively XX individuals. The mode of inheritance is dominant for the polled trait and recessive for sex-reversal. In XX PIS-/- mutants, the expression of testis-specific genes is observed very precociously during gonad development. Nevertheless, a delay of 4-5 days is observed in comparison with normal testis differentiation in XY males. By positional cloning, we demonstrate that the PIS mutation is an 11.7-kb regulatory-deletion affecting the expression of two genes, PISRT1 and FOXL2 which could act synergistically to promote ovarian differentiation. The transcriptional extinction of these two genes leads, very early, to testis-formation in XX homozygous PIS-/- mutants. According to their expression profiles and bibliographic data, we propose that FOXL2 may be an ovary-differentiating gene, and the non-coding RNA PISRT1, an anti-testis factor repressing SOX9, a key regulator of testis differentiation. Under this hypothesis, SRY, the testis-determining factor would inhibit these two genes in the gonads of XY males, to ensure testis differentiation.

Foxl2 gene and the development of the ovary: a story about goat, mouse, fish and woman.

In this review, we describe recent results concerning the genetics of sex determination in mammals. Particularly, we developed the study of the FOXL2 gene and its implication in genetic anomalies in goats (PIS mutation) and humans (BPES). We present the expression of FOXL2 in the ovaries of different species.

Expression profiles and chromosomal localization of genes controlling meiosis and follicular development in the sheep ovary.

In female sheep fetuses, two of the most crucial stages of ovarian development are prophase of meiosis I and follicle formation. In the present study, sheep ovaries collected on Days 25, 38, 49, 56, 67, 75, 94, and 120 of gestation, at birth, and in adulthood were tested by reverse transcription-polymerase chain reaction (RT-PCR) for the expression of 14 genes known to be involved in the ovarian differentiation in diverse organisms. The aim of this study was to determine 1) the expression pattern of six genes involved in germ cell development or meiosis (DMC1, SPO11, MSH4, MSH5, DAZL, and Boule) and five ovary-derived factors (OVOL1, SIAH2, DIAPH2, FOXL2, and FGF9), 2) the onset of gene expression for several members of the bone morphogenetic protein (BMP) pathway involved in follicular development (GDF9, BMP15, BMPR-IB), and 3) the chromosomal localization of seven of these genes in the sheep genome. The RT-PCR analysis revealed that the two germline-specific genes, DAZL and Boule, were expressed between 49 and 94 days postcoitum (dpc) with a similar pattern to typical meiosis genes (DMC1, MSH4, and MSH5), suggesting their possible participation in prophase of meiosis I. GDF9 and OVOL1 gene transcription started at 56 dpc and extended until birth, while BMP15 presented a more restricted window of expression between 94 dpc and birth, corresponding to the formation of first growing follicles. The homologous ovine genes for SPO11, DMC1, MSH5, DAZL, FGF9, DIAPH2, and SIAH2 were located on OAR 13q21-22, 3q35, 20q22, 19q13, 10q15, Xq44, and 1q41-42, respectively. In sheep, quantitative trait loci affecting female reproductive capacities are currently being detected. The ontology and precise mapping of ovarian genes will be useful to identify potential candidate genes that might underlie these effects.

Ontogenesis of female-to-male sex-reversal in XX polled goats.

The association of polledness and intersexuality in domestic goats (PIS mutation) made them a practical genetic model for studying mammalian female-to-male sex reversal. In this study, gonads from XX sex-reversed goats (PIS-/-) were thoroughly characterized at the molecular and histologic level from the first steps of gonadal differentiation (36 days post coitum [dpc]) to birth. The first histologic signs of gonadal sex reversal were detectable between 36 and 40 dpc (4-5 days later than the XY male) and were mainly characterized by the reduction of the ovarian cortex and the organization of seminiferous cords. As early as 36 dpc, aromatase (CYP19) gene expression was decreased in XX (PIS-/-) gonads, whereas genes normally up-regulated in males, such as SOX9 and AMH, showed an increased expression level from 40 dpc. Thereafter, steroidogenic cell precursors were affected, and at 56 dpc, WNT4 and 3beta-HSD were expressed in a male-specific manner in sex-reversed gonads. Another noticeable feature was a progressive disappearance of germ cells, clearly visible in testicular cords around 70 dpc where 50-75% of germ cells were absent in XX (PIS-/-) gonads. These observations indicated that the causal mutation of PIS acts very early in the sex-determining cascade and affects primarily the supporting cells of the gonad.

Cytogenetic localization of 136 genes in the horse: comparative mapping with the human genome.

The aim of this study was to increase the number of type I markers on the horse cytogenetic map and to improve comparison with maps of other species, thus facilitating positional candidate cloning studies. BAC clones from two different sources were FISH mapped: homologous horse BAC clones selected from our newly extended BAC library using consensus primer sequences and heterologous goat BAC clones. We report the localization of 136 genes on the horse cytogenetic map, almost doubling the number of cytogenetically mapped genes with 48 localizations from horse BAC clones and 88 from goat BAC clones. For the first time, genes were mapped to ECA13p, ECA29, and probably ECA30. A total of 284 genes are now FISH mapped on the horse chromosomes. Comparison with the human map defines 113 conserved segments that include new homologous segments not identified by Zoo-FISH on ECA7 and ECA13p.

A first generation bovine BAC-based physical map.

A first generation clone-based physical map for the bovine genome was constructed combining, fluorescent double digestion fingerprinting and sequence tagged site (STS) marker screening. The BAC clones were selected from an Inra BAC library (105,984 clones) and a part of the CHORI-240 BAC library (26,500 clones). The contigs were anchored using the screening information for a total of 1303 markers (451 microsatellites, 471 genes, 127 EST, and 254 BAC ends). The final map, which consists of 6615 contigs assembled from 100,923 clones, will be a valuable tool for genomic research in ruminants, including targeted marker production, positional cloning or targeted sequencing of regions of specific interest.

Sheep/human comparative map in a chromosome region involved in scrapie incubation time shows multiple breakpoints between human chromosomes 14 and 15 and sheep chromosomes 7 and 18.

A chromosome region involved in scrapie incubation time was identified on sheep chromosome 18 (OAR18). Since OAR18 (and OAR7) share conserved chromosome segments with human chromosomes HSA14 and HSA15, a dense map of type I markers was constructed by FISH mapping of bacterial artificial chromosomes containing genes located on these human chromosomes. In this study, we used the complete human sequence information (gene positions in megabases, Mb) to locate approximately one gene every 2 Mb on HSA15 (19 genes mapped between 19.51 and 66.02 Mb) and on HSA14 (11 genes between 73.24 and 102.62 Mb). Combined with previous work carried out in cattle and goats, our results made it possible to refine the comparative map between ruminants and humans for these two highly rearranged chromosomes (10 segments on HSA15 and 7 on HSA14). Furthermore, we identified relatively short intervals containing evolutionary breakpoints, which is a prerequisite to position them precisely. This work is also the first step in the cloning of the region involved in scrapie incubation period in sheep.