A deletion containing a CTCF-element in intron 8 of the Bbs7 gene is partially responsible for juvenile obesity in the Berlin Fat Mouse.

The Berlin Fat Mouse Inbred (BFMI) line is a model for juvenile obesity. Previous studies on crosses between BFMI and C57Bl/6N (B6N) have identified a recessive defect causing juvenile obesity on chromosome 3 (jObes1). Bbs7 was identified as the most likely candidate gene for the observed effect. Comparative sequence analysis showed a 1578 bp deletion in intron 8 of Bbs7 in BFMI mice. A CTCF-element is located inside this deletion. To investigate the functional effect of this deletion, it was introduced into B6N mice using CRISPR/Cas9. Two mice containing the target deletion were obtained (B6N Bbs7emI8∆1 and Bbs7emI8∆2) and were subsequently mated to BFMI and B6N to generate two families suitable for complementation. Inherited alleles were determined and body composition was measured by quantitative magnetic resonance. Evidence for a partial complementation (13.1-15.1%) of the jObes1 allele by the CRISPR/Cas9 modified B6N Bbs7emI8∆1 and Bbs7emI8∆2 alleles was found. Mice carrying the complementation alleles had a 23-27% higher fat-to-lean ratio compared to animals which have a B6N allele (P(Bbs7emI8∆1) = 4.25 × 10-7; P(Bbs7emI8∆2) = 3.17 × 10-5). Consistent with previous findings, the recessive effect of the BFMI allele was also seen for the B6N Bbs7emI8∆1 and Bbs7emI8∆2 alleles. However, the effect size of the B6N Bbs7emI8∆1 and Bbs7emI8∆2 alleles was smaller than the BFMI allele, and thus showed only a partial complementation. Findings suggest additional variants near Bbs7 in addition to or interacting with the deletion in intron 8.

Identification of four novel QTL linked to the metabolic syndrome in the Berlin Fat Mouse.

BACKGROUND: The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with impaired glucose metabolism using the obese lines BFMI861-S1 and BFMI861-S2, which are genetically closely related, but differ in several traits. BFMI861-S1 is insulin resistant and stores ectopic fat in the liver, whereas BFMI861-S2 is insulin sensitive. METHODS: In generation 10, 397 males of an advanced intercross line (AIL) BFMI861-S1 × BFMI861-S2 were challenged with a high-fat, high-carbohydrate diet and phenotyped over 25 weeks. QTL-analysis was performed after selective genotyping of 200 mice using the GigaMUGA Genotyping Array. Additional 197 males were genotyped for 7 top SNPs in QTL regions. For the prioritization of positional candidate genes whole genome sequencing and gene expression data of the parental lines were used. RESULTS: Overlapping QTL for gonadal adipose tissue weight and blood glucose concentration were detected on chromosome (Chr) 3 (95.8-100.1 Mb), and for gonadal adipose tissue weight, liver weight, and blood glucose concentration on Chr 17 (9.5-26.1 Mb). Causal modeling suggested for Chr 3-QTL direct effects on adipose tissue weight, but indirect effects on blood glucose concentration. Direct effects on adipose tissue weight, liver weight, and blood glucose concentration were suggested for Chr 17-QTL. Prioritized positional candidate genes for the identified QTL were Notch2 and Fmo5 (Chr 3) and Plg and Acat2 (Chr 17). Two additional QTL were detected for gonadal adipose tissue weight on Chr 15 (67.9-74.6 Mb) and for body weight on Chr 16 (3.9-21.4 Mb). CONCLUSIONS: QTL mapping together with a detailed prioritization approach allowed us to identify candidate genes associated with traits of the metabolic syndrome. In addition, we provided evidence for direct and indirect genetic effects on blood glucose concentration in the insulin-resistant mouse line BFMI861-S1.

A 5' UTR Mutation Contributes to Down-Regulation of Bbs7 in the Berlin Fat Mouse.

The Bardet-Biedl Syndrome 7 (Bbs7) gene was identified as the most likely candidate gene causing juvenile obesity in the Berlin Fat Mouse Inbred (BFMI) line. Bbs7 expression is significantly lower in the brain, adipose tissue, and liver of BFMI mice compared to lean C57BL/6NCrl (B6N) mice. A DNA sequence comparison between BFMI and B6N revealed 16 sequence variants in the Bbs7 promoter region. Here, we tested if these mutations contribute to the observed differential expression of Bbs7. In a cell-based dual-luciferase assay, we compared the effects of the BFMI and the B6N haplotypes of different regions of the Bbs7 promotor on the reporter gene expression. A single-nucleotide polymorphism (SNP) was identified causing a significant reduction in the reporter gene expression. This SNP (rs29947545) is located in the 5' UTR of Bbs7 at Chr3:36.613.350. The SNP is not unique to BFMI mice but also occurs in several other mouse strains, where the BFMI allele is not associated with lower Bbs7 transcript amounts. Thus, we suggest a compensatory mutation in the other mouse strains that keeps Bbs7 expression at the normal level. This compensatory mechanism is missing in BFMI mice and the cell lines tested.

Effects of DGAT1 on milk performance in Sudanese Butana × Holstein crossbred cattle.

The improvement of milk production of indigenous Sudanese cattle such as Bos indicus Butana and its cross with Holstein is a major goal of the Sudanese government to ensure sufficient healthy nutrition in the country. In this study, we investigated the K232A polymorphism of diacylglycerol acyltransferase (DGAT1), a well-known modulator of milk production in other breeds. We determined allele frequencies and the allele effects on milk production. Therefore, 93 purebred Butana and 203 Butana × Holstein crossbred cattle were genotyped using competitive allele-specific PCR assays. Association analysis was performed using a linear mixed model in R. In purebred Butana cattle, the lysine DGAT1 protein variant K232, which is found to be associated with higher fat and protein contents, as well as higher fat yield was highly frequent at 0.929, while its frequency in Butana × Holstein crossbred cattle was 0.394. Significant effects were found on milk yield (P = 7.6 × 10-20), fat yield (P = 2.2 × 10-17), protein yield (P = 2.0 × 10-19) and lactose yield (P = 4.0 × 10-18) in crossbred cattle. As expected, the protein variant K232 was disadvantageous since it was decreasing milk, protein, and lactose yields by 1.741 kg, 0.063 kg and 0.084 kg, respectively. No significant effects were found for milk fat, protein, and lactose contents. The high frequency of the lysine DGAT1 protein variant K232 in Butana cattle could contribute to their high milk fat content in combination with low milk yield. In Butana × Holstein crossbred cattle, the DGAT1 marker can be used for effective selection and thus genetic improvement of milk production.

Validation of somatic cell score-associated SNPs from Holstein cattle in Sudanese Butana and Butana × Holstein crossbred cattle.

The Bos indicus zebu cattle Butana is the most commonly used indigenous dairy cattle breed in Sudan. In the last years, high-yielding Holstein dairy cattle were introgressed into Butana cattle to improve their milk yield and simultaneously keep their good adaption to extreme environmental conditions. With the focus on the improvement of milk production, other problems arose such as an increased susceptibility to mastitis. Thus, genetic selection for mastitis resistance should be considered to maintain healthy and productive cows. In this study, we tested 10 single nucleotide polymorphisms (SNPs) which had been associated with somatic cell score (SCS) in Holstein cattle for association with SCS in 37 purebred Butana and 203 Butana × Holstein crossbred cattle from Sudan. Animals were genotyped by competitive allele-specific PCR assays and association analysis was performed using a linear mixed model. All 10 SNPs were segregating in the crossbred Butana × Holstein populations, but only 8 SNPs in Sudanese purebred Butana cattle. The SNP on chromosome 13 was suggestively associated with SCS in the Butana × Holstein crossbred population (rs109441194, 13:79,365,467, PBF = 0.054) and the SNP on chromosome 19 was significantly associated with SCS in both populations (rs41257403, 19:50,027,458, Butana: PBF = 0.003, Butana × Holstein: PBF = 6.2 × 10-16). The minor allele of both SNPs showed an increase in SCS. Therefore, selection against the disadvantageous minor allele could be used for genetic improvement of mastitis resistance in the studied populations. However, investigations in a bigger population and across the whole genome are needed to identify additional genomic loci.

Design and performance of a bovine 200 k SNP chip developed for endangered German Black Pied cattle (DSN).

BACKGROUND: German Black Pied cattle (DSN) are an endangered dual-purpose breed which was largely replaced by Holstein cattle due to their lower milk yield. DSN cattle are kept as a genetic reserve with a current herd size of around 2500 animals. The ability to track sequence variants specific to DSN could help to support the conservation of DSN's genetic diversity and to provide avenues for genetic improvement. RESULTS: Whole-genome sequencing data of 304 DSN cattle were used to design a customized DSN200k SNP chip harboring 182,154 variants (173,569 SNPs and 8585 indels) based on ten selection categories. We included variants of interest to DSN such as DSN unique variants and variants from previous association studies in DSN, but also variants of general interest such as variants with predicted consequences of high, moderate, or low impact on the transcripts and SNPs from the Illumina BovineSNP50 BeadChip. Further, the selection of variants based on haplotype blocks ensured that the whole-genome was uniformly covered with an average variant distance of 14.4 kb on autosomes. Using 300 DSN and 162 animals from other cattle breeds including Holstein, endangered local cattle populations, and also a Bos indicus breed, performance of the SNP chip was evaluated. Altogether, 171,978 (94.31%) of the variants were successfully called in at least one of the analyzed breeds. In DSN, the number of successfully called variants was 166,563 (91.44%) while 156,684 (86.02%) were segregating at a minor allele frequency > 1%. The concordance rate between technical replicates was 99.83 ± 0.19%. CONCLUSION: The DSN200k SNP chip was proved useful for DSN and other Bos taurus as well as one Bos indicus breed. It is suitable for genetic diversity management and marker-assisted selection of DSN animals. Moreover, variants that were segregating in other breeds can be used for the design of breed-specific customized SNP chips. This will be of great value in the application of conservation programs for endangered local populations in the future.

Sex-specific genetic architecture in response to American and ketogenic diets.

BACKGROUND/OBJECTIVES: There is a growing appreciation for individual responses to diet. In a previous study, mouse strain-specific responses to American and ketogenic diets were observed. In this study, we searched for genetic variants underlying differences in the responses to American and ketogenic diets between C57BL/6J (B6) and FVB/NJ (FVB) mouse strains. RESULTS: Genetic mapping of fat and lean mass gain revealed QTLs on Chromosome (Chr) 1 at 191.6 Mb (Fmgq1) (P < 0.001, CI = 180.2-194.4 Mb), Chr5 at 73.7 Mb (Fmgq2, Lmgq1) (P < 0.001, CI = 66.1-76.6 Mb), and Chr7 at 40.5 Mb (Fmgq3) (P < 0.01, CI = 36.6-44.5 Mb). Analysis of serum HDL cholesterol concentration identified a significant (P < 0.001, CI = 160.6-176.1 Mb) QTL on Chr1 at 168.6 Mb (Hdlq1). Causal network inference suggests that HDL cholesterol and fat mass gain are both linked to Fmgq1. CONCLUSIONS: Strong sex effects were identified at both Fmgq2 and Lmgq1, which are also diet-dependent. Interestingly, Fmgq2 and Fmgq3 affect fat gain directly, while Fmgq1 influences fat gain directly and via an intermediate change in serum cholesterol. These results demonstrate how precision nutrition will be advanced through the integration of genetic variation and sex in physiological responses to diets varied in carbohydrate composition.

Immunity-related GTPase induces lipophagy to prevent excess hepatic lipid accumulation.

BACKGROUND & AIMS: Currently, only a few genetic variants explain the heritability of fatty liver disease. Quantitative trait loci (QTL) analysis of mouse strains has identified the susceptibility locus Ltg/NZO (liver triglycerides from New Zealand obese [NZO] alleles) on chromosome 18 as associating with increased hepatic triglycerides. Herein, we aimed to identify genomic variants responsible for this association. METHODS: Recombinant congenic mice carrying 5.3 Mbp of Ltg/NZO were fed a high-fat diet and characterized for liver fat. Bioinformatic analysis, mRNA profiles and electrophoretic mobility shift assays were performed to identify genes responsible for the Ltg/NZO phenotype. Candidate genes were manipulated in vivo by injecting specific microRNAs into C57BL/6 mice. Pulldown coupled with mass spectrometry-based proteomics and immunoprecipitation were performed to identify interaction partners of IFGGA2. RESULTS: Through positional cloning, we identified 2 immunity-related GTPases (Ifgga2, Ifgga4) that prevent hepatic lipid storage. Expression of both murine genes and the human orthologue IRGM was significantly lower in fatty livers. Accordingly, liver-specific suppression of either Ifgga2 or Ifgga4 led to a 3-4-fold greater increase in hepatic fat content. In the liver of low-fat diet-fed mice, IFGGA2 localized to endosomes/lysosomes, while on a high-fat diet it associated with lipid droplets. Pulldown experiments and proteomics identified the lipase ATGL as a binding partner of IFGGA2 which was confirmed by co-immunoprecipitation. Both proteins partially co-localized with the autophagic marker LC3B. Ifgga2 suppression in hepatocytes reduced the amount of LC3B-II, whereas overexpression of Ifgga2 increased the association of LC3B with lipid droplets and decreased triglyceride storage. CONCLUSION: IFGGA2 interacts with ATGL and protects against hepatic steatosis, most likely by enhancing the binding of LC3B to lipid droplets. LAY SUMMARY: The genetic basis of non-alcoholic fatty liver disease remains incompletely defined. Herein, we identified members of the immunity-related GTPase family in mice and humans that act as regulators of hepatic fat accumulation, with links to autophagy. Overexpression of the gene Ifgga2 was shown to reduce hepatic lipid storage and could be a therapeutic target for the treatment of fatty liver disease.

A genome-wide association study for clinical mastitis in the dual-purpose German Black Pied cattle breed.

The dual-purpose German Black Pied Cattle (DSN) has become an endangered breed of approximately 2,550 registered cows in Germany. The breed is genetically related to Holstein-Friesian cattle because the old DSN breed contributed to the selection of the modern Holstein dairy cow. In dairy farms, breeders aim to improve animal health and well-being by reducing the number of mastitis cases, which would also reduce milk losses and treatment costs. On the genomic level, no markers associated with clinical mastitis have been reported in DSN. Therefore, we performed a genome-wide association study on 1,062 DSN cows using a univariate linear mixed model that included a relatedness matrix to correct for population stratification. Although the statistical power was limited by the small population size, 3 markers were significantly associated, and 2 additional markers showed a suggestive association with clinical mastitis. Those markers accounted for 1 to 3% of the variance of clinical mastitis in the examined DSN population. One marker was found in the intragenic region of NEURL1 on BTA26, and the other 4 markers in intergenic regions on BTA3, BTA6, and BTA9. Further analyses identified 23 positional candidate genes. Among them is BMPR1B, which has been previously associated with clinical mastitis in other dairy cattle breeds. The markers presented here can be used for selection for mastitis-resistant animals in the endangered DSN population, and can broadly contribute to a better understanding of mastitis determinants in dairy cattle breeds.

Whey protein polymorphisms in Sudanese goat breeds.

The aim of the present study was to assess genetic variation that is characteristic for Sudanese goat breeds in the milk whey protein genes (LALBA and BLG). Four Sudanese goat breeds were screened for variability in LALBA and BLG genes at the DNA level by comparative sequencing of five animals per breed. Sixteen SNPs were identified in LALBA: seven in the upstream region, six synonymous, and three in the 3´-UTR. Three novel synonymous SNPs in exon 2 (ss5197800003, ss5197800012, and ss5197800004) were found in Nubian, Desert, and Nilotic, but not in Taggar goats. One SNP in the promoter of LALBA (rs642745519) modifies a predicted transcription factor binding site for Tcfe2a. The SNPs in the 3'-UTR (rs657915405, rs641559728, and rs664225585) affect predicted miRNA target sites. With respect to haplotypes in the exonic region, haplotype LALBA-A is most frequent in Nubian, Desert, and Nilotic goats, while haplotype LALBA-D is prevalent in Taggar goats. In BLG, 30 SNPs were detected: eight in the upstream gene region, two synonymous, 17 intronic, and three in the 3'-UTR. Among the 30 identified SNPs, 15 were novel. Four of these novel SNPs were located in the upstream gene region, one was synonymous, and ten were intronic. The novel synonymous SNP (ss5197800017), located in exon 2, was only found in Nubian and Nilotic goats. The SNPs ss5197800010 and rs635615192 in the promoter are located in predicted binding sites of transcription factors (M6097, Elk3, Elf5, and GABPA). Among seven haplotypes detected in the coding region, haplotype BLG-A is most frequent in Nubian and Nilotic goats while haplotype BLG-B is most frequent in Desert and Taggar goats. The high variability in regulatory gene regions among Sudanese goats could potentially affect the quality and yield of whey proteins in goat milk and provide a wide resource for genetic improvement of milk production and milk technology characteristics.

High dosage of zinc modulates T-cells in a time-dependent manner within porcine gut-associated lymphatic tissue.

Zn serves as a powerful feed additive to reduce post-weaning diarrhoea in pigs. However, the mechanisms responsible for Zn-associated effects on the adaptive immune responses following feeding of a very high dosage of Zn remain elusive. In this study, we examined the T-cell response in gut-associated lymphatic tissues of seventy-two weaned piglets. Piglets received diets with 57 mg Zn/kg (low Zn concentration, LZn), 164 mg Zn/kg (medium Zn concentration, MZn) or 2425 mg Zn/kg (high Zn concentration, HZn) mg Zn/kg feed for 1, 2 or 4 weeks. We observed that feeding the HZn diet for 1 week increased the level of activated T-helper cells (CD4+ and CD8α dim) compared with feeding MZn and LZn (P<0·05). In addition, we observed higher transcript amounts of interferon γ and T-box 21 (TBET) in the HZn group compared with the MZn and LZn groups (P<0·05). A gene set enrichment analysis revealed an over-representation of genes associated with 'cytokine signalling in immune system'. Remarkably, feeding of a very high Zn dosage led to a switch in the immune response after 2 weeks. We detected higher relative cell counts of CD4+CD25high regulatory T-helper cells (P<0·05) and a higher expression of forkhead box P3 (FOXP3) transcripts (P<0·05). After 4 weeks of feeding a high-dosage Zn diet, the relative CD4+ T-cell count (P<0·05) and the relative CD8ß + T-cell count (P<0·1) were reduced compared with the MZn group. We hypothesise that after 1 week the cellular T-helper 1 response is switched on and after 2 weeks it is switched off, leading to decreased numbers of T-cells.

Cell Specific eQTL Analysis without Sorting Cells.

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn's disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus.

Whole genome population genetics analysis of Sudanese goats identifies regions harboring genes associated with major traits.

BACKGROUND: Sudan is endowed with a variety of indigenous goat breeds which are used for meat and milk production and which are well adapted to the local environment. The aim of the present study was to determine the genetic diversity and relationship within and between the four main Sudanese breeds of Nubian, Desert, Taggar and Nilotic goats. Using the 50 K SNP chip, 24 animals of each breed were genotyped. RESULTS: More than 96% of high quality SNPs were polymorphic with an average minor allele frequency of 0.3. In all breeds, no significant difference between observed (0.4) and expected (0.4) heterozygosity was found and the inbreeding coefficients (FIS) did not differ from zero. Fst coefficients for the genetic distance between breeds also did not significantly deviate from zero. In addition, the analysis of molecular variance revealed that 93% of the total variance in the examined population can be explained by differences among individuals, while only 7% result from differences between the breeds. These findings provide evidence for high genetic diversity and little inbreeding within breeds on one hand, and low diversity between breeds on the other hand. Further examinations using Nei's genetic distance and STRUCTURE analysis clustered Taggar goats distinct from the other breeds. In a principal component (PC) analysis, PC1 could separate Taggar, Nilotic and a mix of Nubian and Desert goats into three groups. The SNPs that contributed strongly to PC1 showed high Fst values in Taggar goat versus the other goat breeds. PCA allowed us to identify target genomic regions which contain genes known to influence growth, development, bone formation and the immune system. CONCLUSIONS: The information on the genetic variability and diversity in this study confirmed that Taggar goat is genetically different from the other goat breeds in Sudan. The SNPs identified by the first principal components show high Fst values in Taggar goat and allowed to identify candidate genes which can be used in the development of breed selection programs to improve local breeds and find genetic factors contributing to the adaptation to harsh environments.

Feeding of Enterococcus faecium NCIMB 10415 Leads to Intestinal miRNA-423-5p-Induced Regulation of Immune-Relevant Genes.

Probiotics are widely used in human and animal health, but little is known about the mode of action of probiotics. One possible mechanism at the molecular level could be an influence on microRNAs (miRNAs) and the related immune-relevant target genes. Here, we analyzed differential expression of miRNA and potential target genes of ileal and jejunal lymphatic tissues from Enterococcus faeciumNCIMB 10415-fed piglets versus untreated controls by using next-generation sequencing. We identified miR-423-5p as being greatly affected by the treatment group (2.32-fold;P= 0.014). Validation by reverse transcription-quantitative PCR (RT-qPCR) confirmed a significant upregulation of miR-423-5p (2.11-fold;P= 0.03) and, additionally, downregulation of the important immune-relevant immunoglobulin lambda light C region (IGLC) (0.61-fold;P= 0.03) and immunoglobulin kappa constant (IGKC) (0.69-fold;P= 0.04) target genes. Expression analysis of miR-423-5p and IGLC at different age points shows a clear anti correlated relationship. Luciferase reporter assays with a HeLa cell line verified IGLC as a target of miR-423-5p. The results provided evidence for an effect of feeding of E. faeciumon the expression of miR-423-5p and on the regulation of the IGLC gene through miR-423-5p. This might be a possible mode of action of E. faeciumon immune cell regulation in the small intestine.

WormQTLHD--a web database for linking human disease to natural variation data in C. elegans.

Interactions between proteins are highly conserved across species. As a result, the molecular basis of multiple diseases affecting humans can be studied in model organisms that offer many alternative experimental opportunities. One such organism-Caenorhabditis elegans-has been used to produce much molecular quantitative genetics and systems biology data over the past decade. We present WormQTL(HD) (Human Disease), a database that quantitatively and systematically links expression Quantitative Trait Loci (eQTL) findings in C. elegans to gene-disease associations in man. WormQTL(HD), available online at http://www.wormqtl-hd.org, is a user-friendly set of tools to reveal functionally coherent, evolutionary conserved gene networks. These can be used to predict novel gene-to-gene associations and the functions of genes underlying the disease of interest. We created a new database that links C. elegans eQTL data sets to human diseases (34 337 gene-disease associations from OMIM, DGA, GWAS Central and NHGRI GWAS Catalogue) based on overlapping sets of orthologous genes associated to phenotypes in these two species. We utilized QTL results, high-throughput molecular phenotypes, classical phenotypes and genotype data covering different developmental stages and environments from WormQTL database. All software is available as open source, built on MOLGENIS and xQTL workbench.

Pheno2Geno - High-throughput generation of genetic markers and maps from molecular phenotypes for crosses between inbred strains.

BACKGROUND: Genetic markers and maps are instrumental in quantitative trait locus (QTL) mapping in segregating populations. The resolution of QTL localization depends on the number of informative recombinations in the population and how well they are tagged by markers. Larger populations and denser marker maps are better for detecting and locating QTLs. Marker maps that are initially too sparse can be saturated or derived de novo from high-throughput omics data, (e.g. gene expression, protein or metabolite abundance). If these molecular phenotypes are affected by genetic variation due to a major QTL they will show a clear multimodal distribution. Using this information, phenotypes can be converted into genetic markers. RESULTS: The Pheno2Geno tool uses mixture modeling to select phenotypes and transform them into genetic markers suitable for construction and/or saturation of a genetic map. Pheno2Geno excludes candidate genetic markers that show evidence for multiple possibly epistatically interacting QTL and/or interaction with the environment, in order to provide a set of robust markers for follow-up QTL mapping. We demonstrate the use of Pheno2Geno on gene expression data of 370,000 probes in 148 A. thaliana recombinant inbred lines. Pheno2Geno is able to saturate the existing genetic map, decreasing the average distance between markers from 7.1 cM to 0.89 cM, close to the theoretical limit of 0.68 cM (with 148 individuals we expect a recombination every 100/148=0.68 cM); this pinpointed almost all of the informative recombinations in the population. CONCLUSION: The Pheno2Geno package makes use of genome-wide molecular profiling and provides a tool for high-throughput de novo map construction and saturation of existing genetic maps. Processing of the showcase dataset takes less than 30 minutes on an average desktop PC. Pheno2Geno improves QTL mapping results at no additional laboratory cost and with minimum computational effort. Its results are formatted for direct use in R/qtl, the leading R package for QTL studies. Pheno2Geno is freely available on CRAN under "GNU GPL v3". The Pheno2Geno package as well as the tutorial can also be found at: http://pheno2geno.nl .

The direction of cross affects [corrected] obesity after puberty in male but not female offspring.

BACKGROUND: We investigated parent-of-origin and allele-specific expression effects on obesity and hepatic gene expression in reciprocal crosses between the Berlin Fat Mouse Inbred line (BFMI) and C57Bl/6NCrl (B6N). RESULTS: We found that F1-males with a BFMI mother developed 1.8 times more fat mass on a high fat diet at 10 weeks than F1-males of a BFMI father. The phenotype was detectable from six weeks on and was preserved after cross-fostering. RNA-seq data of liver provided evidence for higher biosynthesis and elongation of fatty acids (p = 0.00635) in obese male offspring of a BFMI mother versus lean offspring of a BFMI father. Furthermore, fatty acid degradation (p = 0.00198) and the peroxisome pathway were impaired (p = 0.00094). The circadian rhythm was affected as well (p = 0.00087). Among the highest up-regulated protein coding genes in obese males were Acot4 (1.82 fold, p = 0.022), Cyp4a10 (1.35 fold, p = 0.026) and Cyp4a14 (1.32 fold, p = 0.012), which hydroxylize fatty acids and which are known to be increased in liver steatosis. Obese males showed lower expression of the genetically imprinted and paternally expressed 3 (Peg3) gene (0.31 fold, p = 0.046) and higher expression of the androgen receptor (Ar) gene (2.38 fold, p = 0.068). Allelic imbalance was found for expression of ATP-binding cassette transporter gene Abca8b. Several of the differentially expressed genes contain estrogen response elements. CONCLUSIONS: Parent-of-origin effects during gametogenesis and/or fetal development in an obese mother epigenetically modify the transcription of genes that lead to enhanced fatty acid synthesis and impair ß-oxidation in the liver of male, but not female F1 offspring. Down-regulation of Peg3 could contribute to trigger this metabolic setting. At puberty, higher amounts of the androgen receptor and altered access to estrogen response elements in affected genes are likely responsible for male specific expression of genes that were epigenetically triggered. A suggestive lack of estrogen binding motifs was found for highly down-regulated genes in adult hepatocytes of obese F1 males (p = 0.074).

WormQTL--public archive and analysis web portal for natural variation data in Caenorhabditis spp.

Here, we present WormQTL (http://www.wormqtl.org), an easily accessible database enabling search, comparative analysis and meta-analysis of all data on variation in Caenorhabditis spp. Over the past decade, Caenorhabditis elegans has become instrumental for molecular quantitative genetics and the systems biology of natural variation. These efforts have resulted in a valuable amount of phenotypic, high-throughput molecular and genotypic data across different developmental worm stages and environments in hundreds of C. elegans strains. WormQTL provides a workbench of analysis tools for genotype-phenotype linkage and association mapping based on but not limited to R/qtl (http://www.rqtl.org). All data can be uploaded and downloaded using simple delimited text or Excel formats and are accessible via a public web user interface for biologists and R statistic and web service interfaces for bioinformaticians, based on open source MOLGENIS and xQTL workbench software. WormQTL welcomes data submissions from other worm researchers.

Bioinformatics tools and database resources for systems genetics analysis in mice--a short review and an evaluation of future needs.

During a meeting of the SYSGENET working group 'Bioinformatics', currently available software tools and databases for systems genetics in mice were reviewed and the needs for future developments discussed. The group evaluated interoperability and performed initial feasibility studies. To aid future compatibility of software and exchange of already developed software modules, a strong recommendation was made by the group to integrate HAPPY and R/qtl analysis toolboxes, GeneNetwork and XGAP database platforms, and TIQS and xQTL processing platforms. R should be used as the principal computer language for QTL data analysis in all platforms and a 'cloud' should be used for software dissemination to the community. Furthermore, the working group recommended that all data models and software source code should be made visible in public repositories to allow a coordinated effort on the use of common data structures and file formats.

Identifying genotype-by-environment interactions in the metabolism of germinating arabidopsis seeds using generalized genetical genomics.

A complex phenotype such as seed germination is the result of several genetic and environmental cues and requires the concerted action of many genes. The use of well-structured recombinant inbred lines in combination with "omics" analysis can help to disentangle the genetic basis of such quantitative traits. This so-called genetical genomics approach can effectively capture both genetic and epistatic interactions. However, to understand how the environment interacts with genomic-encoded information, a better understanding of the perception and processing of environmental signals is needed. In a classical genetical genomics setup, this requires replication of the whole experiment in different environmental conditions. A novel generalized setup overcomes this limitation and includes environmental perturbation within a single experimental design. We developed a dedicated quantitative trait loci mapping procedure to implement this approach and used existing phenotypical data to demonstrate its power. In addition, we studied the genetic regulation of primary metabolism in dry and imbibed Arabidopsis (Arabidopsis thaliana) seeds. In the metabolome, many changes were observed that were under both environmental and genetic controls and their interaction. This concept offers unique reduction of experimental load with minimal compromise of statistical power and is of great potential in the field of systems genetics, which requires a broad understanding of both plasticity and dynamic regulation.