Mapping Of Markers Research Articles

Targeted mapping of Cdu1, a major locus regulating grain cadmium concentration in durum wheat (Triticum turgidum L. var durum)

Some durum wheat (Triticum turgidum L. var durum) cultivars have the genetic propensity to accumulate cadmium (Cd) in the grain. A major gene controlling grain Cd concentration designated as Cdu1 has been reported on 5B, but the genetic factor(s) conferring the low Cd phenotype are currently unknown. The objectives of this study were to saturate the chromosomal region harboring Cdu1 with newly developed PCR-based markers and to investigate the colinearity of this wheat chromosomal region with rice (Oryza sativa L.) and Brachypodium distachyon genomes. Genetic mapping of markers linked to Cdu1 in a population of recombinant inbred substitution lines revealed that the gene(s) associated with variation in Cd concentration resides in wheat bin 5BL9 between fraction breakpoints 0.76 and 0.79. Genetic mapping and quantitative trait locus (QTL) analysis of grain Cd concentration was performed in 155 doubled haploid lines from the cross W9262-260D3 (low Cd) by Kofa (high Cd) revealed two expressed sequence tag markers (ESMs) and one sequence tagged site (STS) marker that co-segregated with Cdu1 and explained >80% of the phenotypic variation in grain Cd concentration. A second, minor QTL for grain Cd concentration was also identified on 5B, 67 cM proximal to Cdu1. The Cdu1 interval spans 286 kbp of rice chromosome 3 and 282 kbp of Brachypodium chromosome 1. The markers and rice and Brachypodium colinearity described here represent tools that will assist in the positional cloning of Cdu1 and can be used to select for low Cd accumulation in durum wheat breeding programs targeting this trait. The isolation of Cdu1 will further our knowledge of Cd accumulation in cereals as well as metal accumulation in general.

Significance Test and Genome Selection in Bayesian Shrinkage Analysis

Bayesian shrinkage analysis is the state-of-the-art method for whole genome analysis of quantitative traits. It can estimate the genetic effects for the entire genome using a dense marker map. The technique is now called genome selection. A nice property of the shrinkage analysis is that it can estimate effects of QTL as small as explaining 2% of the phenotypic variance in a typical sample size of 300–500 individuals. In most cases, QTL can be detected with simple visual inspection of the entire genome for the effect because the false positive rate is low. As a Bayesian method, no significance test is needed. However, it is still desirable to put some confidences on the estimated QTL effects. We proposed to use the permutation test to draw empirical thresholds to declare significance of QTL under a predetermined genome wide type I error. With the permutation test, Bayesian shrinkage analysis can be routinely used for QTL detection.

A resistance gene to Ustilago nuda in barley is located on chromosome 3H

Loose smut of barley is a common disease which can be controlled using resistant varieties. Information on the chromosome location of loci controlling loose smut resistance and the development of molecular markers to aid in selection for these genes can be beneficial in the resistant variety development process. The objectives of this work were to determine the resistance or susceptibility of doubled haploid barley lines arising from a cross of the varieties ‘Steptoe’ and ‘Morex’ to Ustilago nuda, the causal agent of loose smut of barley, and map the chromosome location of the loose smut resistance locus in ‘Morex’. The reaction to Ustilago nuda of the doubled-haploid barley plants was determined by inoculating spikelets of each line at anthesis by injection of a teliospore suspension using a needle inoculation method. Mature seeds from the inoculated spikelets were grown to determine the percentage of plants that developed with smutted heads. The lines were classified as susceptible if greater than 10% of the plants were smutted. The loose smut resistance locus from the resistant source ‘Morex’ was mapped using an existing DNA marker map of the ‘Steptoe’/‘Morex’ population. The distribution of the resistant and susceptible progeny from the loose smut testing fit a single gene model. The resistance gene was mapped to chromosome 3 (3H).

Bayesian Inference of Genetic Parameters Based on Conditional Decompositions of Multivariate Normal Distributions

It is widely recognized that the mixed linear model is an important tool for parameter estimation in the analysis of complex pedigrees, which includes both pedigree and genomic information, and where mutually dependent genetic factors are often assumed to follow multivariate normal distributions of high dimension. We have developed a Bayesian statistical method based on the decomposition of the multivariate normal prior distribution into products of conditional univariate distributions. This procedure permits computationally demanding genetic evaluations of complex pedigrees, within the user-friendly computer package WinBUGS. To demonstrate and evaluate the flexibility of the method, we analyzed two example pedigrees: a large noninbred pedigree of Scots pine (Pinus sylvestris L.) that includes additive and dominance polygenic relationships and a simulated pedigree where genomic relationships have been calculated on the basis of a dense marker map. The analysis showed that our method was fast and provided accurate estimates and that it should therefore be a helpful tool for estimating genetic parameters of complex pedigrees quickly and reliably.

Using eigenvalues as variance priors in the prediction of genomic breeding values by principal component analysis

Genome-wide selection aims to predict genetic merit of individuals by estimating the effect of chromosome segments on phenotypes using dense single nucleotide polymorphism (SNP) marker maps. In the present paper, principal component analysis was used to reduce the number of predictors in the estimation of genomic breeding values for a simulated population. Principal component extraction was carried out either using all markers available or separately for each chromosome. Priors of predictor variance were based on their contribution to the total SNP correlation structure. The principal component approach yielded the same accuracy of predicted genomic breeding values obtained with the regression using SNP genotypes directly, with a reduction in the number of predictors of about 96% and computation time of 99%. Although these accuracies are lower than those currently achieved with Bayesian methods, at least for simulated data, the improved calculation speed together with the possibility of extracting principal components directly on individual chromosomes may represent an interesting option for predicting genomic breeding values in real data with a large number of SNP. The use of phenotypes as dependent variable instead of conventional breeding values resulted in more reliable estimates, thus supporting the current strategies adopted in research programs of genomic selection in livestock.

Variations in C-reactive protein during a single haemodialysis session do not associate with mortality

An increase in C-reactive protein (CRP) levels during a single haemodialysis (HD) session has been associated with mortality. These associations, however, are difficult to understand from the current understanding of CRP metabolism. In 190 Swedish haemodialysis (HD) patients from the Mapping of Inflammatory Markers in Chronic Kidney Disease (MIMICK) cohort, CRP was measured before and after a HD session. During follow-up, events of death and censoring were recorded, and hazard ratios were calculated and analysed as a function of CRP variation. Results were replicated in 94 Dutch HD patients from the Netherlands Cooperative Study on the Adequacy of Dialysis (NECOSAD). In this cohort, also correlation and kappa statistics were calculated to assess concordance in CRP changes amid multiple dialysis sessions from the same individuals. In both cohorts, mean CRP values did not increase during a single HD session. In the MIMICK, median (interquartile range) dialysis vintage was 29.0 (14.8-57.0) months. In both crude [hazard ratio (95% confidence interval): 1.008 (0.971-1.047)] and multivariate Cox models [0.996 (0.949-1.046)], no association was observed with mortality. In the NECOSAD, individuals endured 6.0 (6.0-12.0) months on dialysis. No association was found with mortality neither in a crude [0.961 (0.908-1.018)] nor in an adjusted analysis [0.978 (0.923-1.037)]. Finally, the concordance between changes in different sessions was poor. CRP changes during a single HD session do not associate with mortality, thereby adding to the biological uncertainty concerning the ability of CRP to rise in such a short period.

Estimating genetic diversity across the neutral genome with the use of dense marker maps

BackgroundWith the advent of high throughput DNA typing, dense marker maps have become available to investigate genetic diversity on specific regions of the genome. The aim of this paper was to compare two marker based estimates of the genetic diversity in specific genomic regions lying in between markers: IBD-based genetic diversity and heterozygosity.MethodsA computer simulated population was set up with individuals containing a single 1-Morgan chromosome and 1665 SNP markers and from this one, an additional population was produced with a lower marker density i.e. 166 SNP markers. For each marker interval based on adjacent markers, the genetic diversity was estimated either by IBD probabilities or heterozygosity. Estimates were compared to each other and to the true genetic diversity. The latter was calculated for a marker in the middle of each marker interval that was not used to estimate genetic diversity.ResultsThe simulated population had an average minor allele frequency of 0.28 and an LD (r2) of 0.26, comparable to those of real livestock populations. Genetic diversities estimated by IBD probabilities and by heterozygosity were positively correlated, and correlations with the true genetic diversity were quite similar for the simulated population with a high marker density, both for specific regions (r = 0.19-0.20) and large regions (r = 0.61-0.64) over the genome. For the population with a lower marker density, the correlation with the true genetic diversity turned out to be higher for the IBD-based genetic diversity.ConclusionsGenetic diversities of ungenotyped regions of the genome (i.e. between markers) estimated by IBD-based methods and heterozygosity give similar results for the simulated population with a high marker density. However, for a population with a lower marker density, the IBD-based method gives a better prediction, since variation and recombination between markers are missed with heterozygosity.

A genome wide association study for QTL affecting direct and maternal effects of stillbirth and dystocia in cattle

Dystocia and stillbirth are significant causes of female and neonatal death in many species and there is evidence for a genetic component to both traits. Identifying causal mutations affecting these traits through genome wide association studies could reveal the genetic pathways involved and will be a step towards targeted interventions. Norwegian Red cattle are an ideal model breed for such studies as very large numbers of records are available. We conducted a genome wide association study for direct and maternal effects of dystocia and stillbirth using almost 1 million records of these traits. Genotyping costs were minimized by genotyping the sires of the recorded cows, and using daughter averages as phenotypes. A dense marker map containing 17,343 single nucleotide polymorphisms covering all autosomal chromosomes was utilized. The genotyped sires were assigned to one of two groups in an attempt to ensure independence between the groups. Associations were only considered validated if they occurred in both groups. Strong associations were found and validated on chromosomes 4, 5, 6, 9, 12, 20, 22 and 28. The QTL region on chromosome 6 was refined using LDLA analysis. The results showed that this chromosome most probably contains two QTL for direct effect on dystocia and one for direct effect on stillbirth. Several candidate genes may be identified close to these QTL. Of these, a cluster of genes expected to affect bone and cartilage formation (i.e. SPP1, IBSP and MEPE) are of particular interest and we suggest that these genes are screened in candidate gene studies for dystocia and stillbirth in cattle as well as other species.

Abstract 4256: Neuroendocrine, stem, and stromal cell interactions involved in the organization of the normal colonic stem cell niche which become disorganized during colon tumorigenesis

Abstract In development of colorectal cancer (CRC), stem cell (SC) overpopulation underlies tumor initiation and progression, but it is incompletely understood which dysregulated mechanisms cause SC overpopulation. Because neuroendocrine cell (NEC) and SC populations both reside in the SC niche at the normal crypt bottom, we hypothesized that i) NECs play a role in maintaining the SC population, and ii) changes in the NEC population contributes to SC overpopulation and altered epithelial: stromal interactions during tumorigenesis. Immunofluorescence mapping and flow cytometry were used to analyze SC niche structure and cell-cell interactions in normal and neoplastic colonic tissues from Caucasian and African-American patients. In normal colonic crypts, most cells (77%) staining positively for the colonic SC marker ALDH1, also co-stained for the NE cell marker chromogranin-A (CGA). In the progression from normal to normal-appearing FAP to adenomatous to malignant colonic tissue, the proportion of ALDH1+/CGA+ cells progressively diminished while the proportion of ALDH1+/CGA- cells increased. The epithelial phenotype of tumor cells was confirmed by co-staining for pancytokeratin. These data suggest that colonic SC express a NE phenotype, and loss of this phenotype is associated with an increased number of SC during colon tumorigenesis. Our studies using flow cytometry and CRC cell lines (COLO320, DiFi, LIM1863, LoVo, SW480, HCT116 & HT29) revealed that each line is characterized by the proportion of cells having a SC phenotype (ALDEFLUOR+) and NEC phenotype (SSTR+ or GLPR2+). have a low proportion. These findings suggest that these proportions are maintained constant through feedback mechanisms involving SC and NEC subpopulations. Our studies of epithelial-stromal interactions used immunofluorescence mapping of basement membrane proteins (Perlecan; collagen IV), inflammation markers (Cox2; Caveolin 1), pericryptal fibroblasts (smooth muscle actin), endothelial cells (CD34) and macrophages (CD68). Our results show that the stroma undergoes progressive disorganization (particularly pericryptal fibroblasts & basement membrane) in parallel with increasing size of the SC population during colon tumorigenesis. These results suggest that genetic changes occurring during colon tumorigenesis lead to breakdown of stromal-epithelial interactions that are crucial to the maintenance of the stem cell niche, which contributes to colonic tissue disorganization that is a cardinal feature of cancer. Taken together, our investigations yield insights into: (i) the control of normal SC populations (including how the SC niche is maintained), (ii) molecular mechanisms underlying SC overpopulation in tumorigenesis, and (iii) potential targets for novel anti-cancer therapeutic strategies. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 4256.

The grape microvine - a model system for rapid forward and reverse genetics of grapevines

A grapevine model system is described that is suitable for rapid forward and reverse genetic studies in small controlled environments. It is based on the Vvgai1 mutant allele that confers a dwarf stature, short generation cycles and continuous flowering ('microvine'). Black and white berry microvine genotypes were developed that can be transformed by Agrobacterium tumefaciens. Near-homozygous lines were created for efficient bi-allelic single nucleotide polymorphism (SNP) marker mapping and mutagenesis studies. A genetic mapping strategy based on picovine-derived microvine progeny populations was used to rapidly phenotype and map the flower sex and fleshless berry loci and identify a new lethal recessive locus, Vvlrl1. The microvine provides a unique model system for rapid genetic studies of grapevine by changing the perennial long life cycle of the plant to one with features and advantages similar to an annual plant.

Genetic mapping of microsatellite markers around the arcelin bruchid resistance locus in common bean

The deployment in common beans (Phaseolus vulgaris L.) of arcelin-based bruchid resistance could help reduce post-harvest storage losses to the Mexican bean weevil [(Zabrotes subfasciatus (Boheman)]. Arcelin is a member of the arcelin-phytohemagglutinin-α-amylase inhibitor (APA) family of seed proteins, which has been extensively studied but not widely used in bean breeding programs. The purpose of this study was to evaluate microsatellite markers for genetic analysis of arcelin-based bruchid resistance and to determine the orientation of markers and the rate of recombination around the APA locus. A total of 10 previously developed microsatellites and 22 newly developed markers based on a sequenced BAC from the APA locus were screened for polymorphism and of these 15 were mapped with an F2 population of 157 individuals resulting from a susceptible × resistant cross of SEQ1006 × RAZ106 that segregated for both the arcelin 1 allele and resistance to the bruchid, Z. subfasciatus. Microsatellites derived from APA gene sequences were linked within 0.8 cM of each other and were placed relative to the rest of the b04 linkage group. In a comparison of genetic to physical distance on the BAC sequence, recombination was found to be moderate with a ratio of 125 kb/cM, but repressed within the APA locus itself. Several markers were predicted to be very effective for genetic studies or marker-assisted selection, based on their significant associations with bruchid resistance and on low adult insect emergence and positions flanking the arcelin and phytohemagglutinin genes.

Fine-mapping quantitative trait loci for twinning rate on Bos taurus chromosome 14 in North American Holsteins1

A previous genome-wide search with a moderate-density 10,000-SNP set identified many marker associations with twinning rate on BTA14 through either single-marker analysis or combined linkage-linkage disequilibrium (LLD; haplotype) analysis. The objective of the current study was to fine-map putative QTL using a more densely populated marker map and both a larger and an independent set of phenotypic data. Holstein bulls (n = 921) from 100 paternal half-sib families were genotyped for 129 SNP markers that included both original and additionally selected markers for increasing marker density in the targeted 34 megabase region. Twinning rate predicted transmitting abilities were calculated using calving records from 1994 to 1998 (data I) and 1999 to 2006 (data II), and the underlying liability scores from threshold model analysis were used as the trait in marker association analyses. The previous analysis used 201 bulls with daughter records in data I. In the current analysis, this was increased to 434, providing a revised estimate of effect and significance. Bulls with daughter records in data II totaled 851, and analysis of these data provided an opportunity for an independent analysis separate from data I. Single-marker association and LLD analyses were performed. Fifteen significant single-marker associations were found (minimally exceeding P < 8.74 x 10(-3)) to concur between data sets. Three and 12 regions in data I and data II, respectively, showed positive results for the presence of QTL from LLD analysis (P < 0.001) within the respective data sets. After combining results from single-marker association, LLD analysis, and model-building strategies, 3 QTL were identified on BTA14. Based on single-marker results from data II, BTA14 harbors QTL responsible for approximately 24% of the variation in twinning rate predicted transmitting ability.

Mapping and validation of simple sequence repeat markers linked to a major gene controlling seed cadmium accumulation in soybean [Glycine max (L.) Merr

Daily consumption of cadmium (Cd) contaminated foods poses a risk to human health. Cultivar selection is an important method to limit Cd uptake and accumulation; however, analyzing grain Cd concentration is costly and time-consuming. Developing markers for low Cd accumulation will facilitate marker assisted selection (MAS). Inheritance studies using a threshold value of 0.2 mg kg(-1) for low and high and an F(2:3) population showed that low Cd accumulation in soybean seed is under the control of a major gene (Cda1, proposed name) with the allele for low accumulation being dominant. A recombinant inbred line (RIL) population (F(6:8)) derived from the cross AC Hime (high Cd accumulation) and Westag-97 (low Cd accumulation) was used to identify the DNA markers linked to Cda gene(s) or quantitative trait loci (QTLs) controlling low Cd accumulation. We screened 171 simple sequence repeat (SSR) primers that showed polymorphism between parents on the 166 RILs. Of these, 40 primers were newly developed from the soybean genomic DNA sequence. Seven SSR markers, SatK138, SatK139, SatK140 (0.5 cM), SatK147, SacK149, SaatK150 and SattK152 (0.3 cM), were linked to Cda1 in soybean seed. All the linked markers were mapped to the same linkage group (LG) K. The closest flanking SSR markers linked to Cda1 were validated using a parallel population (RILs) involving Leo x Westag-97. Linked markers were also validated with diverse soybean genotypes differing in their seed Cd concentration and showed that SSR markers SatK147, SacK149, and SattK152 clearly differentiated the high and low Cd accumulating genotypes tested. To treat Cd uptake as a quantitative trait, QTL analysis using a linkage map constructed with 161 markers identified a major QTL associated with low Cd concentration in the seeds. The QTL was also mapped to the same location as Cda1 on LG-K. This QTL accounted for 57.3% of the phenotypic variation. Potential candidate genes (genes with known or predicted function that could influence the seed Cd concentration) like protein kinase, putative Adagio-like protein, and plasma membrane H(+)-ATPase were found to be located in the locus of interest. Of the four SSR markers located in the region, SattK152 was localized in the plasma membrane H(+)-ATPase gene. SSR markers closely linked to Cda1 in seeds of soybean were identified and have potential to be used for MAS to develop low Cd accumulating cultivars in a breeding program.

Bayesian estimation of marker dosage in sugarcane and other autopolyploids

In sugarcane or other autopolyploids, after generating the data, the first step in constructing molecular marker maps is to determine marker dosage. Improved methods for correctly allocating marker dosage will result in more accurate maps and increased efficiency of QTL linkage detection. When employing dominant markers like AFLPs, single-dose markers represent alleles present as one copy in one parent and null in the other parent, double-dose markers are those present as two copies in one parent and null in the other parent and so on. Observed segregation ratios in the offspring are employed to infer marker dosage in the parent from which the marker was inherited. Commonly, for each marker, a chi (2) test is used to assign dosage. Such an approach does not address important practical considerations such as multiple testing and departures from theoretical assumptions. In particular, extra-binomial variation or overdispersion has been observed in sugarcane studies and standard methods may result in fewer correct dosage allocations than the data warrant. To address these shortcomings, a Bayesian mixture model is proposed where all markers are considered simultaneously. Since analytic solutions are not available, Markov chain Monte Carlo methods are employed. Marker dosage allocation for each individual marker employs the estimated posterior probability of each dosage. For a sugarcane study these methods resulted in more markers being allocated a dosage than by standard approaches. Simulation studies demonstrated that, in general, not only are more markers classified but that more markers are also correctly classified, particularly if overdispersion is present.

Development, characterization, segregation, and mapping of microsatellite markers for European hazelnut (Corylus avellana L.) from enriched genomic libraries and usefulness in genetic diversity studies

Eighty-six new microsatellite loci were developed for European hazelnut (Corylus avellana L.) by screening two genomic libraries enriched for dinucleotide repeats. The loci, 73 developed from a GA-enriched library and 13 from a CA-enriched library, showed a high level of polymorphism in 50 accessions. The number of alleles per locus ranged from five to 21, with a mean of 10.55. Mean values for expected heterozygosity, observed heterozygosity, and polymorphism information content were high, averaging 0.76, 0.69, and 0.73, respectively. In a mapping population, loci segregated 1:1, 1:2:1, and 1:1:1:1. Segregation distortion and null alleles were observed at some loci. Eighty-one of the 86 loci were assigned to linkage groups. A neighbor-joining dendrogram reflected great diversity among the 50 accessions and showed clustering by geographic origin.

Association between genes on chromosome 4p16 and non-syndromic oral clefts in four populations

Isolated cleft lip with or without cleft palate and cleft palate are among the most common human birth defects. Several candidate gene studies on MSX1 have shown significant association between markers in MSX1 and risk of oral clefts, and re-sequencing studies have identified multiple mutations in MSX1 in a small minority of cases, which may account for 1-2% of all isolated oral clefts cases. We explored the 2-Mb region around MSX1, using a marker map of 393 single nucleotide polymorphisms (SNPs) in 297 cleft lip, with or without cleft palate, case-parent trios and 84 cleft palate trios from Maryland, Taiwan, Singapore, and Korea. Both individual markers and haplotypes of two to five SNPs showed several regions yielding statistical evidence for linkage and disequilibrium. Two genes (STK32B and EVC) yielded consistent evidence from cleft lip, with or without cleft palate, trios in all four populations. These two genes plus EVC2 also yielded suggestive evidence for linkage and disequilibrium among cleft palate trios. This analysis suggests that several genes, not just MSX1, in this region may influence risk of oral clefts.

Genetic dissection of scent metabolic profiles in diploid rose populations

The scent of flowers is a very important trait in ornamental roses in terms of both quantity and quality. In cut roses, scented varieties are a rare exception. Although metabolic profiling has identified more than 500 scent volatiles from rose flowers so far, nothing is known about the inheritance of scent in roses. Therefore, we analysed scent volatiles and molecular markers in diploid segregating populations. We resolved the patterns of inheritance of three volatiles (nerol, neryl acetate and geranyl acetate) into single Mendelian traits, and we mapped these as single or oligogenic traits in the rose genome. Three other volatiles (geraniol, beta-citronellol and 2-phenylethanol) displayed quantitative variation in the progeny, and we mapped a total of six QTLs influencing the amounts of these volatiles onto the rose marker map. Because we included known scent related genes and newly generated ESTs for scent volatiles as markers, we were able to link scent related QTLs with putative candidate genes. Our results serve as a starting point for both more detailed analyses of complex scent biosynthetic pathways and the development of markers for marker-assisted breeding of scented rose varieties.

Molecular Mapping of Wheat Leaf Rust Resistance Gene Lr42

ABSTRACTLeaf rust, caused by Puccinia triticina Eriks., is an important foliar disease of wheat (Triticum aestivum L.) worldwide. Leaf rust resistance gene Lr42 from Aegilops tauschii Coss. has been used as a source of rust resistance in breeding programs. To identify molecular markers closely linked to Lr42, a segregating population of near‐isogenic lines contrasting for the presence of Lr42 was developed in the hard winter wheat cultivar Century background and evaluated for rust infection type at both seedling and adult‐plant stages. Simple sequence repeat (SSR) markers were screened using bulked‐segregant analysis. Two markers closely linked to Lr42 were identified on chromosome 1DS. The closest marker, Xwmc432, is about 0.8 cM from Lr42. Physical mapping of both SSR markers using Chinese Spring nullitetrasomic and ditelosomic genetic stocks confirmed that the markers linked to Lr42 were on 1DS. Markers for Lr42 were highly polymorphic between parents and among a diverse set of wheat germplasm collected from several countries, indicating that these markers are useful for marker‐assisted selection for Lr42.

Increase in accuracy using multi-trait genomic breeding value estimation

Introduction Genomic selection is becoming common practice in animal breeding. It uses genome-wide dense marker maps, to accurately predict the genetic ability of animals, without the need to record phenotypic performance from the animal itself or from close relatives. Presented applications of genomic selection have mainly been limited to implementations where genomic breeding values are estimated using single trait models. A major breakthrough in traditional breeding value estimation was the application of multi-trait breeding value estimation, for instance to combine mastitis and somatic cell count information. Therefore, our objective was to develop multi-trait genomic breeding value estimation methods.

Human and livestock genetics: parallel evolution and horizontal exchange

Human and livestock genetics: parallel evolution and horizontal exchange