Novel equine tissue miRNAs and breed-related miRNA expressed in serum.

BACKGROUND: MiRNAs regulate multiple genes at the post-transcriptional level and therefore play an important role in many biological processes. It has been suggested that miRNA exported outside the cells contribute to inter-cellular communication. Consequently, circulating miRNAs are of particular interest and are promising biomarkers for many diseases. The number of miRNAs annotated in the horse genome is much lower compared to model organisms like human and mouse. We therefore aimed to identify novel equine miRNAs for tissue types and breed in serum. RESULTS: We analysed 71 small RNA-seq libraries derived from nine tissues (gluteus medius, platysma, masseter muscle, heart, liver, cartilage, bone, total blood and serum) using miRDeep2 and miRdentify tools. Known miRNAs represented between 2.3 and 62.9 % of the reads in 71 libraries. A total of 683 novel miRNAs were identified. Breed and tissue type affected the number of miRNAs detected and interestingly, affected its average intensity. A total of 50 miRNAs in serum proved to be potential biomarkers to differentiate specific breed types, of which miR-122, miR-200, miR-483 were over-expressed and miR-328 was under-expressed in ponies compared to Warmbloods. The different miRNAs profiles, as well as the differences in their expression levels provide a foundation for more hypotheses based on the novel miRNAs discovered. CONCLUSIONS: We identified 683 novel equine miRNAs expressed in seven solid tissues, blood and serum. Additionally, our approach evidenced that such data supported identification of specific miRNAs as markers of functions related to breeds or disease tissues.

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.

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.

Advancing responsible genomic analyses of ancient mollusc shells.

The analysis of the DNA entrapped in ancient shells of molluscs has the potential to shed light on the evolution and ecology of this very diverse phylum. Ancient genomics could help reconstruct the responses of molluscs to past climate change, pollution, and human subsistence practices at unprecedented temporal resolutions. Applications are however still in their infancy, partly due to our limited knowledge of DNA preservation in calcium carbonate shells and the need for optimized methods for responsible genomic data generation. To improve ancient shell genomic analyses, we applied high-throughput DNA sequencing to 27 Mytilus mussel shells dated to ~111-6500 years Before Present, and investigated the impact, on DNA recovery, of shell imaging, DNA extraction protocols and shell sub-sampling strategies. First, we detected no quantitative or qualitative deleterious effect of micro-computed tomography for recording shell 3D morphological information prior to sub-sampling. Then, we showed that double-digestion and bleach treatment of shell powder prior to silica-based DNA extraction improves shell DNA recovery, also suggesting that DNA is protected in preservation niches within ancient shells. Finally, all layers that compose Mytilus shells, i.e., the nacreous (aragonite) and prismatic (calcite) carbonate layers, with or without the outer organic layer (periostracum) proved to be valuable DNA reservoirs, with aragonite appearing as the best substrate for genomic analyses. Our work contributes to the understanding of long-term molecular preservation in biominerals and we anticipate that resulting recommendations will be helpful for future efficient and responsible genomic analyses of ancient mollusc shells.

Gene expression profiling in equine polysaccharide storage myopathy revealed inflammation, glycogenesis inhibition, hypoxia and mitochondrial dysfunctions.

BACKGROUND: Several cases of myopathies have been observed in the horse Norman Cob breed. Muscle histology examinations revealed that some families suffer from a polysaccharide storage myopathy (PSSM). It is assumed that a gene expression signature related to PSSM should be observed at the transcriptional level because the glycogen storage disease could also be linked to other dysfunctions in gene regulation. Thus, the functional genomic approach could be conducted in order to provide new knowledge about the metabolic disorders related to PSSM. We propose exploring the PSSM muscle fiber metabolic disorders by measuring gene expression in relationship with the histological phenotype. RESULTS: Genotypying analysis of GYS1 mutation revealed 2 homozygous (AA) and 5 heterozygous (GA) PSSM horses. In the PSSM muscles, histological data revealed PAS positive amylase resistant abnormal polysaccharides, inflammation, necrosis, and lipomatosis and active regeneration of fibers. Ultrastructural evaluation revealed a decrease of mitochondrial number and structural disorders. Extensive accumulation of an abnormal polysaccharide displaced and partially replaced mitochondria and myofibrils. The severity of the disease was higher in the two homozygous PSSM horses.Gene expression analysis revealed 129 genes significantly modulated (p < 0.05). The following genes were up-regulated over 2 fold: IL18, CTSS, LUM, CD44, FN1, GST01. The most down-regulated genes were the following: mitochondrial tRNA, SLC2A2, PRKCalpha, VEGFalpha. Data mining analysis showed that protein synthesis, apoptosis, cellular movement, growth and proliferation were the main cellular functions significantly associated with the modulated genes (p < 0.05). Several up-regulated genes, especially IL18, revealed a severe muscular inflammation in PSSM muscles. The up-regulation of glycogen synthase kinase-3 (GSK3beta) under its active form could be responsible for glycogen synthase (GYS1) inhibition and hypoxia-inducible factor (HIF1alpha) destabilization. CONCLUSION: The main disorders observed in PSSM muscles could be related to mitochondrial dysfunctions, glycogenesis inhibition and the chronic hypoxia of the PSSM muscles.

The repertoire of family A-peptide GPCRs in archaic hominins.

Given the importance of G-protein coupled receptors in the regulation of many physiological functions, deciphering the relationships between genotype and phenotype in past and present hominin GPCRs is of main interest to understand the evolutionary process that contributed to the present-day variability in human traits and health. Here, we carefully examined the publicly available genomic and protein sequence databases of the archaic hominins (Neanderthal and Denisova) to draw up the catalog of coding variations in GPCRs for peptide ligands, in comparison with living humans. We then searched in the literature the functional changes, phenotypes and risk of disease possibly associated with the detected variants. Our survey suggests that Neanderthal and Denisovan hominins were likely prone to lower risk of obesity, to enhanced platelet aggregation in response to thrombin, to better response to infection, to less anxiety and aggressiveness and to favorable sociability. While some archaic variants were likely advantageous in the past, they might be responsible for maladaptive disorders today in the context of modern life and/or specific regional distribution. For example, an archaic haplotype in the neuromedin receptor 2 is susceptible to confer risk of diabetic nephropathy in type 1 diabetes in present-day Europeans. Paying attention to the pharmacological properties of some of the archaic variants described in this study may be helpful to understand the variability of therapeutic efficacy between individuals or ethnic groups.

Molecular characterization and chromosomal assignment of equine cartilage derived retinoic acid sensitive protein (CD-RAP)/melanoma inhibitory activity (MIA).

Cartilage-derived retinoic acid sensitive protein (CD-RAP) also known as melanoma inhibitory activity (MIA) has already been established as a marker for chondrocyte differentiation and a number of cancerous conditions in humans. Studies have also shown that CD-RAP/MIA is a potential marker of joint disease. The objective of this study was to characterize the equine CD-RAP/MIA gene and thus make it available as a marker in cartilage research and clinical studies. Gene analysis revealed that the equine gene (GenBank accession no. EF679787) consists of four exons and three introns, and the homology to the human gene is 90% for the translated region. The upstream sequence includes regulatory elements and putative transcription factor binding sites previously described in the human and murine promoter regions. The deduced amino acid sequence consists of 130 aa including a signal peptide of 23 aa, and has a 91% identity to the human protein. Using radiation hybrid mapping, the CD-RAP/MIA gene was localized to the p arm of equine chromosome 10 (ECA10p), which is in accordance with prediction based on the current human-equine comparative map. Gene expression studies showed expression of CD-RAP/MIA mRNA in articular cartilage and chondrocytes from horses with no signs of joint disease. The expression decreased as the cells dedifferentiated in monolayer culture. We also identified an equine CD-RAP/MIA splice variant similar to that reported in humans. The CD-RAP/MIA protein was detected in equine synovial fluid, serum and culture medium from chondrocyte cultures. In conclusion, CD-RAP/MIA is expressed in equine cartilage and chondrocytes, and the protein can be detected in equine serum, synovial fluid and in culture medium from chondrocyte cultures. The equine gene and resulting protein share great homology with the human gene, making future studies on CD-RAP/MIAs potential as a marker in joint disease possible using the equine joint as a model.

Oral health status in historic population: Macroscopic and metagenomic evidence.

Recent developments in High-Throughput DNA sequencing (HTS) technologies and ancient DNA (aDNA) research have opened access to the characterization of the microbial communities within past populations. Most studies have, however, relied on the analysis of dental calculus as one particular material type particularly prone to the molecular preservation of ancient microbial biofilms and potential of entire teeth for microbial characterization, both of healthy communities and pathogens in ancient individuals, remains overlooked. In this study, we used shotgun sequencing to characterize the bacterial composition from historical subjects showing macroscopic evidence of oral pathologies. We first carried out a macroscopic analysis aimed at identifying carious or periodontal diseases in subjects belonging to a French rural population of the 18th century AD. We next examined radiographically six subjects showing specific, characteristic dental pathologies and applied HTS shotgun sequencing to characterize the microbial communities present in and on the dental material. The presence of Streptococcus mutans and also Rothia dentocariosa, Actinomyces viscosus, Porphyromonas gingivalis, Tannerella forsythia, Pseudoramibacter alactolyticus, Olsenella uli and Parvimonas micra was confirmed through the presence of typical signatures of post-mortem DNA damage at an average depth-of-coverage ranging from 0.5 to 7X, with a minimum of 35% (from 35 to 93%) of the positions in the genome covered at least once. Each sampled tooth showed a specific bacterial signature associated with carious or periodontal pathologies. This work demonstrates that from a healthy independent tooth, without visible macroscopic pathology, we can identify a signature of specific pathogens and deduce the oral health status of an individual.

Characterization of the equine glycogen debranching enzyme gene (AGL): Genomic and cDNA structure, localization, polymorphism and expression.

Glycogen debranching enzyme (AGL) is a multifunctional enzyme acting in the glycogen degradation pathway. In humans, the AGL activity deficiency causes a type III glycogen storage disease (Cori-Forbes disease). One particularity of AGL gene expression lies in the multiple alternative splicing in its 5' region. The AGL gene was localized on ECA5q14-q15. The sequence of the equine cDNA was determined to be 7.5 kb in length with an open reading frame of 4602 bp. The gene is 69 kb long and contains 35 exons. The equine AGL gene has an ubiquitous expression and presents five tissue-dependent cDNA variants arising from alternative splicing of the first exons. The equine skeletal muscle and heart contain four out of six variants previously described in humans and the equine liver express three of these four human variants. We identified a new alternative splicing variant expressed in equine skeletal and heart muscles. All these mRNA variants most probably encode only two different protein isoforms of 1533 and 1377 amino-acids. Four SNPs were detected in the mRNA. The equine in silico promoter sequence reveals a structure similar to those of other mammalian species. The disposition of the transcription factor biding sites does not correlate to the transcription start sites of tissue-specific variants.

Endurance Exercise Ability in the Horse: A Trait with Complex Polygenic Determinism.

Endurance horses are able to run at more than 20 km/h for 160 km (in bouts of 30-40 km). This level of performance is based on intense aerobic metabolism, effective body heat dissipation and the ability to endure painful exercise. The known heritabilities of endurance performance and exercise-related physiological traits in Arabian horses suggest that adaptation to extreme endurance exercise is influenced by genetic factors. The objective of the present genome-wide association study (GWAS) was to identify single nucleotide polymorphisms (SNPs) related to endurance racing performance in 597 Arabian horses. The performance traits studied were the total race distance, average race speed and finishing status (qualified, eliminated or retired). We used three mixed models that included a fixed allele or genotype effect and a random, polygenic effect. Quantile-quantile plots were acceptable, and the regression coefficients for actual vs. expected log10p-values ranged from 0.865 to 1.055. The GWAS revealed five significant quantitative trait loci (QTL) corresponding to 6 SNPs on chromosomes 6, 1, 7, 16, and 29 (two SNPs) with corrected p-values from 1.7 × 10-6 to 1.8 × 10-5. Annotation of these 5 QTL revealed two genes: sortilin-related VPS10-domain-containing receptor 3 (SORCS3) on chromosome 1 is involved in protein trafficking, and solute carrier family 39 member 12 (SLC39A12) on chromosome 29 is active in zinc transport and cell homeostasis. These two coding genes could be involved in neuronal tissues (CNS). The other QTL on chromosomes 6, 7, and 16 may be involved in the regulation of the gene expression through non-coding RNAs, CpG islands and transcription factor binding sites. On chromosome 6, a new candidate equine long non-coding RNA (KCNQ1OT1 ortholog: opposite antisense transcript 1 of potassium voltage-gated channel subfamily Q member 1 gene) was predicted in silico and validated by RT-qPCR in primary cultures of equine myoblasts and fibroblasts. This lncRNA could be one element of the cardiac rhythm regulation. Our GWAS revealed that equine performance during endurance races is a complex polygenic trait, and is partially governed by at least 5 QTL: two coding genes involved in neuronal tissues and three other loci with many regulatory functions such as slowing down heart rate.

Shotgun sequencing of the mitochondrial genome of the Aldabra giant tortoise (Aldabrachelys gigantea).

The complete mitochondrial genome of the Aldabra giant tortoise [Aldabrachelys gigantea (Schweigger, 1812): Reptilia, Testudines, Testudinidae] was sequenced using a shotgun approach on an Illumina HiSeq 2500 platform (Illumina Inc., San Diego, CA). This genome was 16 467 bp long and presents the typical organization found in vertebrates. The mean coverage of sequencing was 116×. A phylogenetic analysis of the Testudinidae confirms the placement of Aldabrachelys in an Indian Ocean group (including Madagascar). This mitogenome constitutes a reference for ancient DNA analyses of the extinct Madagascan lineages of Aldabrachelys.

Sequence of goat cyclin T1 cDNA, gene organisation and expression analysis.

The cyclin T1 (Cyc T1) protein has been recently identified, associated with the cyclin-dependent kinase 9 (CDK 9), as to be involved in the transcriptional activation of the Human Immunodeficiency Virus type 1 (HIV-1) by the Tat protein. In this study, the sequence of the 7 kb goat Cyc T1 cDNA is reported as well as the exon/intron structure of the gene. Its observed ubiquitous expression is consistent with the promoter structure.

cDNA sequence of the horse (Equus caballus) LAMA3 gene and characterization of two intronic SNP markers.

Laminins are large heterotrimeric basement membrane glycoproteins composed of alpha, beta and gamma chains. The Laminin 5 isoform has an alpha3beta3gamma2 composition and is essential for the adhesion of basal keratinocytes to the underlying epithelial basement membrane where it is mainly located. Mutations in the genes coding for the 3 chains have been associated with a severe skin blistering disease, Herlitz's junctional epidermolysis bullosa (JEB), observed in different species as man, dog, cat and horse. In this study, we report the sequence of the 5.2 kb horse laminin alpha 3 cDNA (LAMA3) as well as the detection of two intronic SNPs. These data will be useful to further identify causal mutations for the disease in this gene.

Putative FLJ20436 gene characterisation in goat. Observed ubiquitous expression in goat and transgenic mice allowed to restrict the location of an hypothesised insulator element.

Eukaryotic genomes are organised into independent domains through the establishment of boundaries which allow to have distinct pattern of gene expression both during development and in differentiated cells. The previously reported site independent expression of the mammary-specific goat alpha-lactalbumin gene in transgenic mice suggested the existence of cis-regulatory elements located upstream of this gene. The nearby presence of a second ubiquitously expressed gene (the cyclin T1 gene) allowed to define two chromatin domains putatively separated by a boundary insulator-like element. In this study, the characterisation of a third putative gene (FLJ20436) present between the alpha-lactalbumin and the cyclin T1 loci is reported. It was found to be a functional gene, ubiquitously expressed both in goat and transgenic mice. A complex pattern of alternative splicing events was observed in several analysed tissues leading to various mRNAs and putative FLJ20436 proteins, as suggested in human and mouse species. This allowed us to assign this gene to one of the two hypothesised chromatin domains, refining the location of the potential insulator element.

Protein catabolism and high lipid metabolism associated with long-distance exercise are revealed by plasma NMR metabolomics in endurance horses.

During long distance endurance races, horses undergo high physiological and metabolic stresses. The adaptation processes involve the modulation of the energetic pathways in order to meet the energy demand. The aims were to evaluate the effects of long endurance exercise on the plasma metabolomic profiles and to investigate the relationships with the individual horse performances. The metabolomic profiles of the horses were analyzed using the non-dedicated methodology, NMR spectroscopy and statistical multivariate analysis. The advantage of this method is to investigate several metabolomic pathways at the same time in a single sample. The plasmas were obtained before exercise (BE) and post exercise (PE) from 69 horses competing in three endurance races at national level (130-160 km). Biochemical assays were also performed on the samples taken at PE. The proton NMR spectra were compared using the supervised orthogonal projection on latent structure method according to several factors. Among these factors, the race location was not significant whereas the effect of the race exercise (sample BE vs PE of same horse) was highly discriminating. This result was confirmed by the projection of unpaired samples (only BE or PE sample of different horses). The metabolomic profiles proved that protein, energetic and lipid metabolisms as well as glycoproteins content are highly affected by the long endurance exercise. The BE samples from finisher horses could be discriminated according to the racing speed based on their metabolomic lipid content. The PE samples could be discriminated according to the horse ranking position at the end of the race with lactate as unique correlated metabolite. As a conclusion, the metabolomic profiles of plasmas taken before and after the race provided a better understanding of the high energy demand and protein catabolism pathway that could expose the horses to metabolic disorders.

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.

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.

Ubiquitous expression of goat cyclin T1 in transgenic mice.

Ubiquitous gene expression has a variety of applications in transgenesis that include, for example cell lineage analyses in chimeras and gain-of-function related to xenotransplantations. Although several promoters have already been used to these aims, they often do not reliably or reproducibly target gene expression in mice. We have recently reported the site-independent expression of a bacterial artificial chromosome (BAC)-derived goat alpha-lactalbumin transgene in the mammary gland of mice and the subsequent localisation within the insert of this BAC of the cyclin T1 locus. This ubiquitously expressed gene encodes for a protein that acts as a co-factor for the HIV nuclear transcriptional activator. In the present paper, we report that the goat BAC transgene, which encompasses around 30 kb of the cyclin T1 promoter, also confers ubiquitous expression of this gene in the six transgenic mouse lines analysed. These results suggest that the cyclin T1 promoter could be a useful alternative to target ubiquitous gene expression in transgenics.