We present here the de novo genome assembly CerEla1.0 for the red deer, Cervus elaphus, an emblematic member of the natural megafauna of the Northern Hemisphere. Humans spread the species in the South. Today, the red deer is also a farm-bred animal and is becoming a model animal in biomedical and population studies. Stag DNA was sequenced at 74× coverage by Illumina technology. The ALLPATHS-LG assembly of the reads resulted in 34.7 × 10 scaffolds, 26.1 × 10 of which were utilized in Cer.Ela1.0. The assembly spans 3.4 Gbp. For building the red deer pseudochromosomes, a pre-established genetic map was used for main anchor points. A nearly complete co-linearity was found between the mapmarker sequences of the deer genetic map and the order and orientation of the orthologous sequences in the syntenic bovine regions. Syntenies were also conserved at the in-scaffold level. The cM distances corresponded to 1.34 Mbp uniformly along the deer genome. Chromosomal rearrangements between deer and cattle were demonstrated. 2.8 × 10 SNPs, 365 × 10 indels and 19368 protein-coding genes were identified in CerEla1.0, along with positions for centromerons. CerEla1.0 demonstrates the utilization of dual references, i.e., when a target genome (here C. elaphus) already has a pre-established genetic map, and is combined with the well-established whole genome sequence of a closely related species (here Bos taurus). Genome-wide association studies (GWAS) that CerEla1.0 (NCBI, MKHE00000000) could serve for are discussed.

BackgroundMangalicas are fatty type local/rare pig breeds with an increasing presence in the niche pork market in Hungary and in other countries. To explore their genetic resources, we have analysed data from next-generation sequencing of an individual male from each of three Mangalica breeds along with a local male Duroc pig. Structural variations, such as SNPs, INDELs and CNVs, were identified and particular genes with SNP variations were analysed with special emphasis on functions related to fat metabolism in pigs.ResultsMore than 60 Gb of sequence data were generated for each of the sequenced individuals, resulting in 11× to 19× autosomal median coverage. After stringent filtering, around six million SNPs, of which approximately 10% are novel compared to the dbSNP138 database, were identified in each animal. Several hundred thousands of INDELs and about 1,000 CNV gains were also identified. The functional annotation of genes with exonic, non-synonymous SNPs, which are common in all three Mangalicas but are absent in either the reference genome or the sequenced Duroc of this study, highlighted 52 genes in lipid metabolism processes. Further analysis revealed that 41 of these genes are associated with lipid metabolic or regulatory pathways, 49 are in fat-metabolism and fatness-phenotype QTLs and, with the exception of ACACA, ANKRD23, GM2A, KIT, MOGAT2, MTTP, FASN, SGMS1, SLC27A6 and RETSAT, have not previously been associated with fat-related phenotypes.ConclusionsGenome analysis of Mangalica breeds revealed that local/rare breeds could be a rich source of sequence variations not present in cosmopolitan/industrial breeds. The identified Mangalica variations may, therefore, be a very useful resource for future studies of agronomically important traits in pigs.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-761) contains supplementary material, which is available to authorized users.

Understanding the molecular mechanisms underlying bone development is a fundamental and fascinating problem in developmental biology, with significant medical implications. Here, we have identified the expression patterns for 36 genes that were characteristic or dominant in the consecutive cell differentiation zones (mesenchyme, precartilage, cartilage) of the tip section of the developing velvet antler of red deer Cervus elaphus. Two major functional groups of these genes clearly outlined: six genes linked to high metabolic demand and other five to tumor biology. Our study demonstrates the advantages of the antler as a source of mesenchymal markers, for distinguishing precartilage and cartilage by different gene expression patterns and for identifying genes involved in the robust bone development, a striking feature of the growing antler. Putative roles for "antler" genes that encode alpha-tropomyosine (tpm1), transgelin (tagln), annexin 2 (anxa2), phosphatidylethanolamine-binding protein (pebp) and apolipoprotein D (apoD) in intense but still controlled tissue proliferation are discussed.

BackgroundMangalica breeds are indigenous to Hungary and their breeding history dates back to about 200–250 years ago. They are fat-type pigs and have a rare curly hair phenotype. The aim of our study was to establish the relationships between these unique breeds and other European breeds.ResultsBased on a core sequence of 382 bp present in 2713 mitochondrial D-loop sequences from pigs belonging to 38 local breeds from nine countries, five cosmopolitan breeds and wild boars from 14 countries, we identified 164 haplotypes. More than half of the 2713 sequences belonged to either four haplotypes characteristic of continental European breeds or two haplotypes characteristic of British/cosmopolitan breeds; each haplotype is present in more than 100 individuals. Most Mangalica individuals belonged either to one of these common continental European haplotypes or to two Mangalica-specific haplotypes that were absent in all other breeds. In addition, we identified the ancestral mitochondrial D-loop signature present in these 2713 sequences and found that ~ 80% carried the European ancient signatures, ANC-Aside and ANC-Cside or their closely related signatures, while most of the remaining sequences carried a modern Asian signature, ANC-Easia. Mangalica individuals carried the ANC-Aside signature, but not the ANC-Cside or ANC-Easia signatures.ConclusionsIn all the Mangalica individuals, a unique ancient European signature was found in the mitochondrial DNA D-loop region, but they belonged almost exclusively to either certain very abundant European or two Mangalica-specific D-loop haplotypes. This indicates that the present-day Mangalica population in Hungary evolved either by introgression of other European breeds and wild boars or via total isolation after the divergence of European ancient porcine bloodlines.

Red deer is the most valuable game of the fauna in Hungary, and there is a strong need for genetic identification of individuals. For this purpose, 10 tetranucleotide STR markers were developed and amplified in two 5-plex systems. The study presented here includes the flanking region sequence analysis and the allele nomenclature of the 10 loci as well as the PCR optimization of the DeerPlex I and II. LD pairwise tests and cross-species similarity analyses showed the 10 loci to be independently inherited. Considerable levels of genetic differences between two subpopulations were recorded, and F(ST) was 0.034 using AMOVA. The average probability of identity (PI(ave)) was at the value of 2.6736 × 10(-15). This low value for PI(ave) nearly eliminates false identification. An illegal hunting case solved by DeerPlex is described herein. The calculated likelihood ratio (LR) illustrates the potential of the 10 red deer microsatellite markers for forensic investigations.

Osteoporosis attacks 10% of the population worldwide. Humans or even the model animals of the disease cannot recover from porous bone. Regeneration in skeletal elements is the unique feature of our newly investigated osteoporosis model, the red deer (Cervus elaphus) stag. Cyclic physiological osteoporosis is a consequence of the annual antler cycle. This phenomenon raises the possibility to identify genes involved in the regulation of bone mineral density on the basis of comparative genomics between deer and human. We compare gene expression activity of osteoporotic and regenerating rib bone samples versus autumn dwell control in red deer by microarray hybridization. Identified genes were tested on human femoral bone tissue from non-osteoporotic controls and patients affected with age-related osteoporosis. Expression data were evaluated by Principal Components Analysis and Canonical Variates Analysis. Separation of patients into a normal and an affected group based on ten formerly known osteoporosis reference genes was significantly improved by expanding the data with newly identified genes. These genes include IGSF4, FABP3, FABP4, FKBP2, TIMP2, TMSB4X, TRIB, and members of the Wnt signaling. This study supports that extensive comparative genomic analyses, here deer and human, provide a novel approach to identify new targets for human diagnostics and therapy.

Annual re-growth of deer antler represents a unique example of complete organ regeneration. Because antler mesenchymal cells retain their embryonic capacity to develop into cartilage or bone, studying antler development provides a natural system to follow gene expression changes during mesenchymal differentiation toward chondrogenic/osteogenic lineage. To identify novel genes involved either in early events of mesenchymal cell specialization or in robust bone development, we have introduced a 3 K heterologous microarray set-up (deer cDNA versus mouse template). Fifteen genes were differentially expressed; genes for housekeeping, regulatory functions (components of different signaling pathways, including FGF, TGFbeta, Wnt), and genes encoding members of the Polycomb group were represented. Expression dynamics for genes are visualized by an expression logo. The expression profile of the gene C21orf70 of unknown function is described along with the effects when over-expressed; furthermore the nuclear localization of the cognate protein is shown. In this report, we demonstrate the particular advantage of the velvet antler model in bone research for: (1) identification of mesenchymal and precartilaginous genes and (2) targeting genes upregulated in robust cartilage development.

The lactose operon of Escherichia coli is a paradigm system for quantitative understanding of gene regulation in prokaryotes. Yet, none of the many mathematical models built so far to study the dynamics of this system considered the fact that the Lac repressor regulates its own transcription by forming a transcriptional roadblock at the O3 operator site. Here we study the effect of autoregulation on intracellular LacI levels and also show that cAMP-CRP binding does not affect the efficiency of autoregulation. We built a mathematical model to study the role of LacI autoregulation in the lactose utilization system. Previously, it has been argued that negative autoregulation can significantly reduce noise as well as increase the speed of response. We show that the particular molecular mechanism, a transcriptional roadblock, used to achieve self-repression in the lac system does neither. Instead, LacI autoregulation balances two opposing states, one that allows quicker response to smaller pulses of external lactose, and the other that minimizes production costs in the absence of lactose.