RNA-Seq transcriptomics and pathway analyses reveal potential regulatory genes and molecular mechanisms in high- and low-residual feed intake in Nordic dairy cattle

BACKGROUND: The selective breeding of cattle with high-feed efficiencies (FE) is an important goal of beef and dairy cattle producers. Global gene expression patterns in relevant tissues can be used to study the functions of genes that are potentially involved in regulating FE. In the present study, high-throughput RNA sequencing data of liver biopsies from 19 dairy cows were used to identify differentially expressed genes (DEGs) between high- and low-FE groups of cows (based on Residual Feed Intake or RFI). Subsequently, a profile of the pathways connecting the DEGs to FE was generated, and a list of candidate genes and biomarkers was derived for their potential inclusion in breeding programmes to improve FE. RESULTS: The bovine RNA-Seq gene expression data from the liver was analysed to identify DEGs and, subsequently, identify the molecular mechanisms, pathways and possible candidate biomarkers of feed efficiency. On average, 57 million reads (short reads or short mRNA sequences < ~200 bases) were sequenced, 52 million reads were mapped, and 24,616 known transcripts were quantified according to the bovine reference genome. A comparison of the high- and low-RFI groups revealed 70 and 19 significantly DEGs in Holstein and Jersey cows, respectively. The interaction analysis (high vs. low RFI x control vs. high concentrate diet) showed no interaction effects in the Holstein cows, while two genes showed interaction effects in the Jersey cows. The analyses showed that DEGs act through certain pathways to affect or regulate FE, including steroid hormone biosynthesis, retinol metabolism, starch and sucrose metabolism, ether lipid metabolism, arachidonic acid metabolism and drug metabolism cytochrome P450. CONCLUSION: We used RNA-Seq-based liver transcriptomic profiling of high- and low-RFI dairy cows in two breeds and identified significantly DEGs, their molecular mechanisms, their interactions with other genes and functional enrichments of different molecular pathways. The DEGs that were identified were the CYP’s and GIMAP genes for the Holstein and Jersey cows, respectively, which are related to the primary immunodeficiency pathway and play a major role in feed utilization and the metabolism of lipids, sugars and proteins. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12864-017-3622-9) contains supplementary material, which is available to authorized users

Representativeness of microsatellite distributions in genomes, as revealed by 454 GS-FLX Titanium pyrosequencing

<p>Abstract</p> <p>Background</p> <p>Microsatellites are markers of choice in population genetics and genomics, as they provide useful insight into patterns and processes as diverse as genome evolutionary dynamics and demographic processes. The acquisition of microsatellites through multiplex-enriched libraries and 454 GS-FLX Titanium pyrosequencing is a promising new tool for the isolation of new markers in unknown genomes. This approach can also be used to evaluate the extent to which microsatellite-enriched libraries are representative of the genome from which they were isolated. In this study, we deciphered potential discrepancies in microsatellite content recovery for two reference genomes (<it>Apis mellifera </it>and <it>Danio rerio</it>), selected on the basis of their extreme heterogeneity in terms of the proportions and distributions of microsatellites on chromosomes.</p> <p>Results</p> <p>The <it>A. mellifera </it>genome, in particular, was found to be highly heterogeneous, due to extremely high rates of recombination, with hotspots, but the only bias consistently introduced into pyrosequenced multiplex-enriched libraries concerned sequence length, with the overrepresentation of sequences 160 to 320 bp in length. Other deviations from expected proportions or distributions of motifs on chromosomes were observed, but the significance and intensity of these deviations was mostly limited. Furthermore, no consistent adverse competition between multiplexed probes was observed during the motif enrichment phase.</p> <p>Conclusions</p> <p>This approach therefore appears to be a promising strategy for improving the development of microsatellites, as it introduces no major bias in terms of the proportions and distribution of microsatellites.</p

Genetic analysis of scattered populations of the Indian eri silkworm, Samia cynthia ricini Donovan: Differentiation of subpopulations

Deforestation and exploitation has led to the fragmentation of habitats and scattering of populations of the economically important eri silkworm, Samia cynthia ricini, in north-east India. Genetic analysis of 15 eri populations, using ISSR markers, showed 98% inter-population, and 23% to 58% intra-population polymorphism. Nei’s genetic distance between populations increased significantly with altitude (R2 = 0.71) and geographic distance (R2 = 0.78). On the dendrogram, the lower and upper Assam populations were clustered separately, with intermediate grouping of those from Barpathar and Chuchuyimlang, consistent with geographical distribution. The Nei’s gene diversity index was 0.350 in total populations and 0.121 in subpopulations. The genetic differentiation estimate (Gst) was 0.276 among scattered populations. Neutrality tests showed deviation of 118 loci from Hardy-Weinberg equilibrium. The number of loci that deviated from neutrality increased with altitude (R2 = 0.63). Test of linkage disequilibrium showed greater contribution of variance among eri subpopulations to total variance. D’2IS exceeded D’2ST, showed significant contribution of random genetic drift to the increase in variance of disequilibrium in subpopulations. In the Lakhimpur population, the peripheral part was separated from the core by a genetic distance of 0.260. Patchy habitats promoted low genetic variability, high linkage disequilibrium and colonization by new subpopulations. Increased gene flow and habitat-area expansion are required to maintain higher genetic variability and conservation of the original S. c. ricini gene pool

A Novel Approach for Mining Polymorphic Microsatellite Markers In Silico

An important emerging application of high-throughput 454 sequencing is the isolation of molecular markers such as microsatellites from genomic DNA. However, few studies have developed microsatellites from cDNA despite the added potential for targeting candidate genes. Moreover, to develop microsatellites usually requires the evaluation of numerous primer pairs for polymorphism in the focal species. This can be time-consuming and wasteful, particularly for taxa with low genetic diversity where the majority of primers often yield monomorphic polymerase chain reaction (PCR) products. Transcriptome assemblies provide a convenient solution, functional annotation of transcripts allowing markers to be targeted towards candidate genes, while high sequence coverage in principle permits the assessment of variability in silico. Consequently, we evaluated fifty primer pairs designed to amplify microsatellites, primarily residing within transcripts related to immunity and growth, identified from an Antarctic fur seal (Arctocephalus gazella) transcriptome assembly. In silico visualization was used to classify each microsatellite as being either polymorphic or monomorphic and to quantify the number of distinct length variants, each taken to represent a different allele. The majority of loci (n = 36, 76.0%) yielded interpretable PCR products, 23 of which were polymorphic in a sample of 24 fur seal individuals. Loci that appeared variable in silico were significantly more likely to yield polymorphic PCR products, even after controlling for microsatellite length measured in silico. We also found a significant positive relationship between inferred and observed allele number. This study not only demonstrates the feasibility of generating modest panels of microsatellites targeted towards specific classes of gene, but also suggests that in silico microsatellite variability may provide a useful proxy for PCR product polymorphism

An analysis of automatic gender classification

<p>[EN]Different researches suggest that inner facial features are not the only discriminative features for tasks such as person identification or gender classification. Indeed, they have shown an influence of features which are part of the local face context, such as hair, on these tasks. However, object-centered approaches which ignore local context dominate the research in computational vision based facial analysis. In this paper, we performed an analysis to study which areas and which resolutions are diagnostic for the gender classification problem. We first demonstrate the importance of contextual features in human observers for gender classification using a psychophysical ”bubbles” technique.</p

Polymorphic microsatellite loci for the endangered pine hoverfly Blera fallax (Diptera: Syrphidae)

We describe eleven polymorphic microsatellite loci developed for Blera fallax using &lsquo;next generation&rsquo; 454 whole genome shotgun sequencing, along with conditions for three multiplex PCR reactions. We tested allelic variation on forty B. fallax individuals from Scotland and Sweden. Allelic richness and expected heterozygosity were 3.03 (&plusmn;0.274) and 0.391 (&plusmn;0.057) respectively. The number of alleles per locus ranged from 2 to 6. We anticipate that these loci will assist conservation management by allowing the monitoring of translocated populations, estimating effective population size, and assessing population structure and dispersal in Scotland and across Europe