QTL global meta-analysis: are trait determining genes clustered?

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when the effects of homozygosity for recessive mutations are studied in non-isogenic backgrounds, genes located proximal to the mutation on the same chromosome often appear over-represented among those genes identified as differentially expressed (DE). One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect functional responses to the mutation but, instead, result from an unequal distribution of expression quantitative trait loci (eQTLs) between sample groups of mutant or wild-type genotypes. This is problematic because eQTL expression differences are difficult to distinguish from genes that are DE due to functional responses to a mutation. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) are also observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between those sample groups subjected to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting-based approaches. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. Our method for DAR analysis requires only RNA-sequencing data, facilitating its application across new and existing datasets.

Meta-QTL analysis in wheat: progress, challenges and opportunities

Wheat, an important cereal crop globally, faces major challenges due to increasing global population and changing climates. The production and productivity are challenged by several biotic and abiotic stresses. There is also a pressing demand to enhance grain yield and quality/nutrition to ensure global food and nutritional security. To address these multifaceted concerns, researchers have conducted numerous meta-QTL (MQTL) studies in wheat, resulting in the identification of candidate genes that govern these complex quantitative traits. MQTL analysis has successfully unraveled the complex genetic architecture of polygenic quantitative traits in wheat. Candidate genes associated with stress adaptation have been pinpointed for abiotic and biotic traits, facilitating targeted breeding efforts to enhance stress tolerance. Furthermore, high-confidence candidate genes (CGs) and flanking markers to MQTLs will help in marker-assisted breeding programs aimed at enhancing stress tolerance, yield, quality and nutrition. Functional analysis of these CGs can enhance our understanding of intricate trait-related genetics. The discovery of orthologous MQTLs shared between wheat and other crops sheds light on common evolutionary pathways governing these traits. Breeders can leverage the most promising MQTLs and CGs associated with multiple traits to develop superior next-generation wheat cultivars with improved trait performance. This review provides a comprehensive overview of MQTL analysis in wheat, highlighting progress, challenges, validation methods and future opportunities in wheat genetics and breeding, contributing to global food security and sustainable agriculture.

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when homozygous mutations are studied in non-isogenic backgrounds, genes from the same chromosome as a mutation often appear over-represented among differentially expressed (DE) genes. One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect true biological responses to the mutation but, instead, result from differences in representation of expression quantitative trait loci (eQTLs) between sample groups selected on the basis of mutant or wild-type genotype. This is problematic when inclusion of spurious DE genes in a functional enrichment study results in incorrect inferences of mutation effect. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) can also be observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between groups of samples subject to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting of gene-level rankings. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. The DAR metric provides a solid foundation for addressing the eQTL issue in new and existing datasets because it relies solely on RNA-sequencing data.

Genetic architecture of divergence: the selfing syndrome in Ipomoea lacunosa

PREMISE

Highly selfing plant species frequently display a distinctive suite of traits termed the selfing syndrome. Here we tested the hypothesis that these traits are grouped into correlated evolutionary modules and determined the degree of independence between such modules.

METHODS

We evaluated phenotypic correlations and QTL overlaps in F2 offspring of a cross between the morning glories Ipomoea lacunosa and I. cordatotriloba and investigated how traits clustered into modules at both the phenotypic and genetic level. We then compared our findings to other QTL studies of the selfing syndrome.

RESULTS

In the I. lacunosa selfing syndrome, traits grouped into modules that displayed correlated evolution within but not between modules. QTL overlap predicted phenotypic correlations, and QTLs affecting the same trait module were significantly physically clustered in the genome. The genetic architecture of the selfing syndrome varied across systems, but the pattern of stronger within- than between-module correlation was widespread.

CONCLUSIONS

The genetic architecture we observe in the selfing syndrome is consistent with a growing understanding of floral morphological integration achieved via pleiotropy in clustered traits. This view of floral evolution is consistent with resource limitation or predation driving the evolution of the selfing syndrome, but invites further research into both the selective causes of the selfing syndrome and how genetic architecture itself evolves in response to changes in mating system.

Zebrafish Chromosome 14 Gene Differential Expression in the fmr1 h u2787 Model of Fragile X Syndrome

Zebrafish represent a valuable model for investigating the molecular and cellular basis of Fragile X syndrome (FXS). Reduced expression of the zebrafish FMR1 orthologous gene, fmr1, causes developmental and behavioural phenotypes related to FXS. Zebrafish homozygous for the hu2787 non-sense mutation allele of fmr1 are widely used to model FXS, although FXS-relevant phenotypes seen from morpholino antisense oligonucleotide (morpholino) suppression of fmr1 transcript translation were not observed when hu2787 was first described. The subsequent discovery of transcriptional adaptation (a form of genetic compensation), whereby mutations causing non-sense-mediated decay of transcripts can drive compensatory upregulation of homologous transcripts independent of protein feedback loops, suggested an explanation for the differences reported. We examined the whole-embryo transcriptome effects of homozygosity for fmr1 h u2787 at 2 days post fertilisation. We observed statistically significant changes in expression of a number of gene transcripts, but none from genes showing sequence homology to fmr1. Enrichment testing of differentially expressed genes implied effects on lysosome function and glycosphingolipid biosynthesis. The majority of the differentially expressed genes are located, like fmr1, on Chromosome 14. Quantitative PCR tests did not support that this was artefactual due to changes in relative chromosome abundance. Enrichment testing of the "leading edge" differentially expressed genes from Chromosome 14 revealed that their co-location on this chromosome may be associated with roles in brain development and function. The differential expression of functionally related genes due to mutation of fmr1, and located on the same chromosome as fmr1, is consistent with R.A. Fisher's assertion that the selective advantage of co-segregation of particular combinations of alleles of genes will favour, during evolution, chromosomal rearrangements that place them in linkage disequilibrium on the same chromosome. However, we cannot exclude that the apparent differential expression of genes on Chromosome 14 genes was, (if only in part), caused by differences between the expression of alleles of genes unrelated to the effects of the fmr1 h u2787 mutation and made manifest due to the limited, but non-zero, allelic diversity between the genotypes compared.

Selection signature analysis reveals genes underlying sheep milking performance

Sheep milk is the most important feed resource for newborn lambs and an important food resource for humans. Sheep milk production and ingredients are influenced by genetic and environmental factors. In this study, we implemented selection signature analysis using Illumina Ovine SNP50 BeadChip data of 78 meat Lacaune and 103 milk Lacaune sheep, which have similar genetic backgrounds, from the Sheep HapMap project to identify candidate genes related to ovine milk traits. Since different methods can detect different variation types and complement each other, we used a haplotype-based method (hapFLK) to implement selection signature analysis. The results revealed six selection signature regions showing signs of being selected ( P < 0.001 ): chromosomes 1, 2, 3, 6, 13 and 18. In addition, 38 quantitative trait loci (QTLs) related to sheep milk performance were identified in selection signature regions, which contain 334 candidate genes. Of those, SUCNR1 (succinate receptor 1) and PPARGC1A (PPARG coactivator 1 alpha) may be the most significant genes that affect sheep milking performance, which supply a significant indication for future studies to investigate candidate genes that play an important role in milk production and quality.

Lectin Sequence Distribution in QTLs from Rice (Oryza sativa) Suggest A Role in Morphological Traits and Stress Responses

Rice (Oryza sativa) is one of the main staple crops worldwide but suffers from important yield losses due to different abiotic and biotic stresses. Analysis of quantitative trait loci (QTL) is a classical genetic method which enables the creation of more resistant cultivars but does not yield information on the genes directly involved or responsible for the desired traits. Lectins are known as proteins with diverse functions in plants. Some of them are abundant proteins in seeds and are considered as storage/defense proteins while other lectins are known as stress-inducible proteins, implicated in stress perception and signal transduction as part of plant innate immunity. We investigated the distribution of lectin sequences in different QTL related to stress tolerance/resistance, morphology, and physiology through mapping of the lectin sequences and QTL regions on the chromosomes and subsequent statistical analysis. Furthermore, the domain structure and evolutionary relationships of the lectins in O. sativa spp. indica and japonica were investigated. Our results revealed that lectin sequences are statistically overrepresented in QTLs for (a)biotic resistance/tolerance as well as in QTLs related to economically important traits such as eating quality and sterility. These findings contribute to the characterization of the QTL sequences and can provide valuable information to the breeders.

Genomic regions associated with host response to porcine reproductive and respiratory syndrome vaccination and co-infection in nursery pigs

BACKGROUND

The WUR1000125 (WUR) single nucleotide polymorphism (SNP) can be used as a genetic marker for host response to porcine reproductive and respiratory syndrome (PRRS), PRRS vaccination, and co-infection with porcine circovirus type 2b (PCV2b). Objectives of this study were to identify genomic regions other than WUR associated with host response to PRRS vaccination and PRRSV/PCV2b co-infection and regions with a different effect on host response to co-infection, depending on previous vaccination for PRRS.

METHODS

Commercial crossbred nursery pigs were pre-selected for WUR genotype (n = 171 AA and 198 AB pigs) where B is the dominant and favorable allele. Half of the pigs were vaccinated for PRRS and 4 weeks later, all pigs were co-infected with PRRS virus and PCV2b. Average daily gain (ADG) and viral load (VL) were quantified post vaccination (Post Vx) and post co-infection (Post Co-X). Single-SNP genome-wide association analyses were then conducted to identify genomic regions associated with response to vaccination and co-infection.

RESULTS

Multiple SNPs near the major histocompatibility complex were significantly associated with PCV2b VL (-log ₁₀ P ≥ 5.5), regardless of prior vaccination for PRRS. Several SNPs were also significantly associated with ADG Post Vx and Post Co-X. SNPs with a different effect on ADG, depending on prior vaccination for PRRS, were identified Post Vx (-log ₁₀ P = 5.6) and Post Co-X (-log ₁₀ P = 5.5). No SNPs were significantly associated with vaccination VL (-log₁₀ P ≤ 4.7) or PRRS VL (-log₁₀ P ≤ 4.3). Genes near SNPs associated with vaccination VL, PRRS VL, and PCV2b VL were enriched (P ≤ 0.01) for immune-related pathways and genes near SNPs associated with ADG were enriched for metabolism pathways (P ≤ 0.04). SNPs associated with vaccination VL, PRRS VL, and PCV2b VL showed overrepresentation of health QTL identified in previous studies and SNPs associated with ADG Post Vx of Non-Vx pigs showed overrepresentation of growth QTL.

CONCLUSIONS

Multiple genomic regions were associated with PCV2b VL and ADG Post Vx and Post Co-X. Different SNPs were associated with ADG, depending on previous vaccination for PRRS. Results of functional annotation analyses and novel approaches of using previously-reported QTL support the identified regions.

Action of multiple intra-QTL genes concerted around a co-localized transcription factor underpins a large effect QTL

Sub-QTLs and multiple intra-QTL genes are hypothesized to underpin large-effect QTLs. Known QTLs over gene families, biosynthetic pathways or certain traits represent functional gene-clusters of genes of the same gene ontology (GO). Gene-clusters containing genes of different GO have not been elaborated, except in silico as coexpressed genes within QTLs. Here we demonstrate the requirement of multiple intra-QTL genes for the full impact of QTL qDTY12.1 on rice yield under drought. Multiple evidences are presented for the need of the transcription factor 'no apical meristem' (OsNAM12.1) and its co-localized target genes of separate GO categories for qDTY12.1 function, raising a regulon-like model of genetic architecture. The molecular underpinnings of qDTY12.1 support its effectiveness in further improving a drought tolerant genotype and for its validity in multiple genotypes/ecosystems/environments. Resolving the combinatorial value of OsNAM12.1 with individual intra-QTL genes notwithstanding, identification and analyses of qDTY12.1has fast-tracked rice improvement towards food security.

Genome-wide association analyses for boar taint components and testicular traits revealed regions having pleiotropic effects

Background

The aim of this study was to perform a genome-wide association analyses (GWAS) for androstenone, skatole and indole in different Pietrain sire lines and compare the results with previous findings in purebred populations. Furthermore, the genetic relationship of androstenone and skatole were investigated with respect to pleiotropy. In order to characterize the performance of intact boars, crossbred progenies of 136 Pietrain boars mated to crossbred sows from three different breeding companies were tested on four test stations. A total of 598 boars were performance tested according to the rules of stationary performance testing in Germany. Beside common fattening and carcass composition traits, the concentrations of the boar taint components and testicular size parameters were recorded. All boars were genotyped with the PorcineSNP60 Illumina BeadChip. The GWAS were performed using the whole data set as well as in sub groups according to the line of origin. Besides an univariate GWAS approach, principal component (PC) techniques were applied to identify common expression pattern affecting the biosynthesis and the metabolism of androstenone.

Results

In total, 33 SNPs were significantly associated with at least one of the boar taint components. Only one SNP was identified being significant in both subgroups. The analyses of the testes size parameters revealed 31 significant associations. The numbers of significant SNPs within the genetic groups evidenced the strong population specific effects. A multivariate approach using PC revealed 33 significant associations for five different PC.

Conclusions

Based on Pietrain sired cross bred boars, the mayor objective of our study was to identify QTL for boar taint components and to detect pleiotropy among boar taint and testes traits. The high number of identified QTL revealed that boar taint traits are influenced by a large number of loci. Analyzing pleiotropy allowed identifying a QTL affecting androstenone and the gonasomatic index. In this region, QTL for ovulation rate and age at puberty of sows have been described in literature. This supports the physiological findings that the androstenone level of boars and reproduction performance of sows might be linked by an antagonistic relationship.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0194-z) contains supplementary material, which is available to authorized users.

Genome-wide detection of selective signature in Chinese Holstein

Selective signatures in whole genome can help us understand the mechanisms of selection and target causal variants for breeding program. In present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions harboring such signals in Chinese Holstein, and then verified the biological significance of these identified regions based on commonly-used bioinformatics analyses. Results showed a total of 125 significant regions in entire genome containing some of important functional genes such as LEP, ABCG2, CSN1S1, CSN3 and TNF based on the Gene Ontology database. Some of these annotated genes involved in the core regions overlapped with those identified in our previous GWAS as well as those involved in a recently constructed candidate gene database for cattle, further indicating these genes under positive selection maybe underlie milk production traits and other important traits in Chinese Holstein. Furthermore, in the enrichment analyses for the second level GO terms and pathways, we observed some significant terms over represented in these identified regions as compared to the entire bovine genome. This indicates that some functional genes associated with milk production traits, as reflected by GO terms, could be clustered in core regions, which provided promising evidence for the exploitability of the core regions identified by EHH tests. Findings in our study could help detect functional candidate genes under positive selection for further genetic and breeding research in Chinese Holstein.

Phenotypic overlap in the contribution of individual genes to CNV pathogenicity revealed by cross-species computational analysis of single-gene mutations in humans, mice and zebrafish

Numerous disease syndromes are associated with regions of copy number variation (CNV) in the human genome and, in most cases, the pathogenicity of the CNV is thought to be related to altered dosage of the genes contained within the affected segment. However, establishing the contribution of individual genes to the overall pathogenicity of CNV syndromes is difficult and often relies on the identification of potential candidates through manual searches of the literature and online resources. We describe here the development of a computational framework to comprehensively search phenotypic information from model organisms and single-gene human hereditary disorders, and thus speed the interpretation of the complex phenotypes of CNV disorders. There are currently more than 5000 human genes about which nothing is known phenotypically but for which detailed phenotypic information for the mouse and/or zebrafish orthologs is available. Here, we present an ontology-based approach to identify similarities between human disease manifestations and the mutational phenotypes in characterized model organism genes; this approach can therefore be used even in cases where there is little or no information about the function of the human genes. We applied this algorithm to detect candidate genes for 27 recurrent CNV disorders and identified 802 gene-phenotype associations, approximately half of which involved genes that were previously reported to be associated with individual phenotypic features and half of which were novel candidates. A total of 431 associations were made solely on the basis of model organism phenotype data. Additionally, we observed a striking, statistically significant tendency for individual disease phenotypes to be associated with multiple genes located within a single CNV region, a phenomenon that we denote as pheno-clustering. Many of the clusters also display statistically significant similarities in protein function or vicinity within the protein-protein interaction network. Our results provide a basis for understanding previously un-interpretable genotype-phenotype correlations in pathogenic CNVs and for mobilizing the large amount of model organism phenotype data to provide insights into human genetic disorders.

Integrating genome annotation and QTL position to identify candidate genes for productivity, architecture and water-use efficiency in Populus spp

BACKGROUND

Hybrid poplars species are candidates for biomass production but breeding efforts are needed to combine productivity and water use efficiency in improved cultivars. The understanding of the genetic architecture of growth in poplar by a Quantitative Trait Loci (QTL) approach can help us to elucidate the molecular basis of such integrative traits but identifying candidate genes underlying these QTLs remains difficult. Nevertheless, the increase of genomic information together with the accessibility to a reference genome sequence (Populus trichocarpa Nisqually-1) allow to bridge QTL information on genetic maps and physical location of candidate genes on the genome. The objective of the study is to identify QTLs controlling productivity, architecture and leaf traits in a P. deltoides x P. trichocarpa F1 progeny and to identify candidate genes underlying QTLs based on the anchoring of genetic maps on the genome and the gene ontology information linked to genome annotation. The strategy to explore genome annotation was to use Gene Ontology enrichment tools to test if some functional categories are statistically over-represented in QTL regions.

RESULTS

Four leaf traits and 7 growth traits were measured on 330 F1 P. deltoides x P. trichocarpa progeny. A total of 77 QTLs controlling 11 traits were identified explaining from 1.8 to 17.2% of the variation of traits. For 58 QTLs, confidence intervals could be projected on the genome. An extended functional annotation was built based on data retrieved from the plant genome database Phytozome and from an inference of function using homology between Populus and the model plant Arabidopsis. Genes located within QTL confidence intervals were retrieved and enrichments in gene ontology (GO) terms were determined using different methods. Significant enrichments were found for all traits. Particularly relevant biological processes GO terms were identified for QTLs controlling number of sylleptic branches: intervals were enriched in GO terms of biological process like 'ripening' and 'adventitious roots development'.

CONCLUSION

Beyond the simple identification of QTLs, this study is the first to use a global approach of GO terms enrichment analysis to fully explore gene function under QTLs confidence intervals in plants. This global approach may lead to identification of new candidate genes for traits of interest.

Integrating genome annotation and QTL position to identify candidate genes for productivity, architecture and water-use efficiency in Populus spp

BACKGROUND

Hybrid poplars species are candidates for biomass production but breeding efforts are needed to combine productivity and water use efficiency in improved cultivars. The understanding of the genetic architecture of growth in poplar by a Quantitative Trait Loci (QTL) approach can help us to elucidate the molecular basis of such integrative traits but identifying candidate genes underlying these QTLs remains difficult. Nevertheless, the increase of genomic information together with the accessibility to a reference genome sequence (Populus trichocarpa Nisqually-1) allow to bridge QTL information on genetic maps and physical location of candidate genes on the genome. The objective of the study is to identify QTLs controlling productivity, architecture and leaf traits in a P. deltoides x P. trichocarpa F1 progeny and to identify candidate genes underlying QTLs based on the anchoring of genetic maps on the genome and the gene ontology information linked to genome annotation. The strategy to explore genome annotation was to use Gene Ontology enrichment tools to test if some functional categories are statistically over-represented in QTL regions.

RESULTS

Four leaf traits and 7 growth traits were measured on 330 F1 P. deltoides x P. trichocarpa progeny. A total of 77 QTLs controlling 11 traits were identified explaining from 1.8 to 17.2% of the variation of traits. For 58 QTLs, confidence intervals could be projected on the genome. An extended functional annotation was built based on data retrieved from the plant genome database Phytozome and from an inference of function using homology between Populus and the model plant Arabidopsis. Genes located within QTL confidence intervals were retrieved and enrichments in gene ontology (GO) terms were determined using different methods. Significant enrichments were found for all traits. Particularly relevant biological processes GO terms were identified for QTLs controlling number of sylleptic branches: intervals were enriched in GO terms of biological process like 'ripening' and 'adventitious roots development'.

CONCLUSION

Beyond the simple identification of QTLs, this study is the first to use a global approach of GO terms enrichment analysis to fully explore gene function under QTLs confidence intervals in plants. This global approach may lead to identification of new candidate genes for traits of interest.

Meta-analysis of QTL mapping experiments

Meta-analysis is an important tool for integrating information from multiple quantitative trait loci (QTLs) studies. Pooling of results from several studies allows greater statistical power for QTL detection and more precise estimation of their genetic effects. Hence, a meta-analysis can yield conclusions that are stronger than those of individual studies and can give greater insight into the genetic architecture of complex traits. In this chapter, we present basic theories and methods for meta-analysis of QTL mapping experiments. The meta-analytic procedures are described in a general context. The statistical methods cover both parametric and nonparametric statistical models. Finally, we illustrate the features of these statistical methods using simulated and real datasets.

Invited review: quantitative trait nucleotide determination in the era of genomic selection

Genome-wide association studies based on tens of thousands of single nucleotide polymorphisms have been completed for several dairy cattle populations. Methods have been proposed to directly incorporate genome scan data into breeding programs, chiefly by selection of young sires based on their genotypes for the genetic markers and pedigree without progeny test. Thus, the rate of genetic gain is increased by reduction of the mean generation interval. The methods developed so far for application of genomic selection do not require identification of the actual quantitative trait nucleotides (QTN) responsible for the observed variation of quantitative trait loci (QTL). To date, 2 QTN affecting milk production traits have been detected in dairy cattle: DGAT1 and ABCG2. This review will attempt to address the following questions based on the current state of bovine genomics and statistics. What are the pros and cons for QTN determination? How can data obtained from high-density, genome-wide scans be used most efficiently for QTN determination? Can the genome scan results already available and next-generation sequencing data be used to determine QTN? Should QTN be treated differently than markers at linkage disequilibrium with QTL in genetic evaluation programs? Data obtained by genome-wide association studies can be used to deduce QTL genotypes of sires via application of the a posteriori granddaughter design for concordance testing of putative QTN. This, together with next-generation sequencing technology, will dramatically reduce costs for QTN determination. By complete genome sequencing of 21 sires with many artificial insemination sons, it should be possible to determine concordance for all potential QTN, thus establishing the field of QTNomics.

Bovine Genome Database: integrated tools for genome annotation and discovery

The Bovine Genome Database (BGD; http://BovineGenome.org) strives to improve annotation of the bovine genome and to integrate the genome sequence with other genomics data. BGD includes GBrowse genome browsers, the Apollo Annotation Editor, a quantitative trait loci (QTL) viewer, BLAST databases and gene pages. Genome browsers, available for both scaffold and chromosome coordinate systems, display the bovine Official Gene Set (OGS), RefSeq and Ensembl gene models, non-coding RNA, repeats, pseudogenes, single-nucleotide polymorphism, markers, QTL and alignments to complementary DNAs, ESTs and protein homologs. The Bovine QTL viewer is connected to the BGD Chromosome GBrowse, allowing for the identification of candidate genes underlying QTL. The Apollo Annotation Editor connects directly to the BGD Chado database to provide researchers with remote access to gene evidence in a graphical interface that allows editing and creating new gene models. Researchers may upload their annotations to the BGD server for review and integration into the subsequent release of the OGS. Gene pages display information for individual OGS gene models, including gene structure, transcript variants, functional descriptions, gene symbols, Gene Ontology terms, annotator comments and links to National Center for Biotechnology Information and Ensembl. Each gene page is linked to a wiki page to allow input from the research community.

Mapping quantitative trait loci (QTL) in sheep. II. Meta-assembly and identification of novel QTL for milk production traits in sheep

An (Awassi × Merino) × Merino backcross family of 172 ewes was used to map quantitative trait loci (QTL) for different milk production traits on a framework map of 200 loci across all autosomes. From five previously proposed mathematical models describing lactation curves, the Wood model was considered the most appropriate due to its simplicity and its ability to determine ovine lactation curve characteristics. Derived milk traits for milk, fat, protein and lactose yield, as well as percentage composition and somatic cell score were used for single and two-QTL approaches using maximum likelihood estimation and regression analysis. A total of 15 significant (P < 0.01) and additional 25 suggestive (P < 0.05) QTL were detected across both single QTL methods and all traits. In preparation of a meta-analysis, all QTL results were compared with a meta-assembly of QTL for milk production traits in dairy ewes from various public domain sources and can be found on the ReproGen ovine gbrowser http://crcidp.vetsci.usyd.edu.au/cgi-bin/gbrowse/oaries_genome/. Many of the QTL for milk production traits have been reported on chromosomes 1, 3, 6, 16 and 20. Those on chromosomes 3 and 20 are in strong agreement with the results reported here. In addition, novel QTL were found on chromosomes 7, 8, 9, 14, 22 and 24. In a cross-species comparison, we extended the meta-assembly by comparing QTL regions of sheep and cattle, which provided strong evidence for synteny conservation of QTL regions for milk, fat, protein and somatic cell score data between cattle and sheep.

Unlocking the bovine genome

The draft genome sequence of cattle (Bos taurus) has now been analyzed by the Bovine Genome Sequencing and Analysis Consortium and the Bovine HapMap Consortium, which together represent an extensive collaboration involving more than 300 scientists from 25 different countries.