High-fat-diet induced development of increased fasting glucose levels and impaired response to intraperitoneal glucose challenge in the collaborative cross mouse genetic reference population.

BACKGROUND: The prevalence of Type 2 Diabetes (T2D) mellitus in the past decades, has reached epidemic proportions. Several lines of evidence support the role of genetic variation in the pathogenesis of T2D and insulin resistance. Elucidating these factors could contribute to developing new medical treatments and tools to identify those most at risk. The aim of this study was to characterize the phenotypic response of the Collaborative Cross (CC) mouse genetic resource population to high-fat diet (HFD) induced T2D-like disease to evaluate its suitability for this purpose. RESULTS: We studied 683 mice of 21 different lines of the CC population. Of these, 265 mice (149 males and 116 females) were challenged by HFD (42% fat); and 384 mice (239 males and 145 females) of 17 of the 21 lines were reared as control group on standard Chow diet (18% fat). Briefly, 8 week old mice were maintained on HFD until 20 weeks of age, and subsequently assessed by intraperitoneal glucose tolerance test (IPGTT). Biweekly body weight (BW), body length (BL), waist circumstance (WC), and body mass index (BMI) were measured. On statistical analysis, trait measurements taken at 20 weeks of age showed significant sex by diet interaction across the different lines and traits. Consequently, males and females were analyzed, separately. Differences among lines were analyzed by ANOVA and shown to be significant (P <0.05), for BW, WC, BMI, fasting blood glucose, and IPGTT-AUC. We use these data to infer broad sense heritability adjusted for number of mice tested in each line; coefficient of genetic variation; genetic correlations between the same trait in the two sexes, and phenotypic correlations between different traits in the same sex. CONCLUSIONS: These results are consistent with the hypothesis that host susceptibility to HFD-induced T2D is a complex trait and controlled by multiple genetic factors and sex, and that the CC population can be a powerful tool for genetic dissection of this trait.

Targeted Recombinant Progeny: a design for ultra-high resolution mapping of Quantitative Trait Loci in crosses between inbred or pure lines.

BACKGROUND: High-resolution mapping of the loci (QTN) responsible for genetic variation in quantitative traits is essential for positional cloning of candidate genes, and for effective marker assisted selection. The confidence interval (QTL) flanking the point estimate of QTN-location is proportional to the number of individuals in the mapping population carrying chromosomes recombinant in the given interval. Consequently, many designs for high resolution QTN mapping are based on increasing the proportion of recombinants in the mapping population. The "Targeted Recombinant Progeny" (TRP) design is a new design for high resolution mapping of a target QTN in crosses between pure, or inbred lines. It is a three-generation procedure generating a large number of recombinant individuals within a QTL previously shown to contain a QTN. This is achieved by having individuals that carry chromosomes recombinant across the target QTL interval as parents of a large mapping population; most of whom will therefore carry recombinant chromosomes targeted to the given QTL. The TRP design is particularly useful for high resolution mapping of QTN that differentiate inbred or pure lines, and hence are not amenable to high resolution mapping by genome-wide association tests. RESULTS: In the absence of residual polygenic variation, population sizes required for achieving given mapping resolution by the TRP-F2 design relative to a standard F2 design ranged from 0.289 for a QTN with standardized allele substitution effect = 0.2, mapped to an initial QTL of 0.2 Morgan to 0.041 for equivalent QTN mapped to an initial QTL of 0.02 M. In the presence of residual polygenic variation, the relative effectiveness of the TRP design ranges from 1.068 to 0.151 for the same initial QTL intervals and QTN effect. Thus even in the presence of polygenic variation, the TRP can still provide major savings. Simulation showed that mapping by TRP should be based on 30-50 markers spanning the initial interval; and on at least 50 or more G2 families representing this number of recombination points,. CONCLUSIONS: The TRP design can be an effective procedure for achieving high and ultra-high mapping resolution of a target QTN previously mapped to a known confidence interval (QTL).

Heritability and coefficient of genetic variation analyses of phenotypic traits provide strong basis for high-resolution QTL mapping in the Collaborative Cross mouse genetic reference population.

Most biological traits of human importance are complex in nature; their manifestation controlled by the cumulative effect of many genetic factors interacting with one another and with the individual's life history. Because of this, mouse genetic reference populations (GRPs) consisting of collections of inbred lines or recombinant inbred lines (RIL) derived from crosses between inbred lines are of particular value in analysis of complex traits, since massive amounts of data can be accumulated on the individual lines. However, existing mouse GRPs are derived from inbred lines that share a common history, resulting in limited genetic diversity, and reduced mapping precision due to long-range gametic disequilibrium. To overcome these limitations, the Collaborative Cross (CC) a genetically highly diverse collection of mouse RIL was established. The CC, now in advanced stages of development, will eventually consist of about 500 RIL derived from reciprocal crosses of eight divergent founder strains of mice, including three wild subspecies. Previous studies have shown that the CC indeed contains enormous diversity at the DNA level, that founder haplotypes are inherited in expected frequency, and that long-range gametic disequilibrium is not present. We here present data, primarily from our own laboratory, documenting extensive genetic variation among CC lines as expressed in broad-sense heritability (H(2)) and by the well-known "coefficient of genetic variation," demonstrating the ability of the CC resource to provide unprecedented mapping precision leading to identification of strong candidate genes.

Genome-wide association study for somatic cell score in Valdostana Red Pied cattle breed using pooled DNA.

BACKGROUND: Mastitis is a major disease of dairy cattle occurring in response to environmental exposure to infective agents with a great economic impact on dairy industry. Somatic cell count (SCC) and its log transformation in somatic cell score (SCS) are traits that have been used as indirect measures of resistance to mastitis for decades in selective breeding. A selective DNA pooling (SDP) approach was applied to identify Quantitative Trait Loci (QTL) for SCS in Valdostana Red Pied cattle using the Illumina Bovine HD BeadChip. RESULTS: A total of 171 SNPs reached the genome-wide significance for association with SCS. Fifty-two SNPs were annotated within genes, some of those involved in the immune response to mastitis. On BTAs 1, 2, 3, 4, 9, 13, 15, 17, 21 and 22 the largest number of markers in association to the trait was found. These regions identified novel genomic regions related to mastitis (1-Mb SNP windows) and confirmed those already mapped. The largest number of significant SNPs exceeding the threshold for genome-wide significant signal was found on BTA 15, located at 50.43-51.63 Mb. CONCLUSIONS: The genomic regions identified in this study contribute to a better understanding of the genetic control of the mastitis immune response in cattle and may allow the inclusion of more detailed QTL information in selection programs.

Genotype is an important determinant factor of host susceptibility to periodontitis in the Collaborative Cross and inbred mouse populations.

BACKGROUND: Periodontal infection (Periodontitis) is a chronic inflammatory disease, which results in the breakdown of the supporting tissues of the teeth. Previous epidemiological studies have suggested that resistance to chronic periodontitis is controlled to some extent by genetic factors of the host. The aim of this study was to determine the phenotypic response of inbred and Collaborative Cross (CC) mouse populations to periodontal bacterial challenge, using an experimental periodontitis model. In this model, mice are co-infected with Porphyromonas gingivalis and Fusobacterium nucleatum, bacterial strains associated with human periodontal disease. Six weeks following the infection, the maxillary jaws were harvested and analyzed for alveolar bone loss relative to uninfected controls, using computerized microtomography (microCT). Initially, four commercial inbred mouse strains were examined to calibrate the procedure and test for gender effects. Subsequently, we applied the same protocol to 23 lines (at inbreeding generations 10-18) from the newly developed mouse genetic reference population, the Collaborative Cross (CC) to determine heritability and genetic variation of control bone volume prior to infection (CBV, naïve bone volume around the teeth of uninfected mice), and residual bone volume (RBV, bone volume after infection) and loss of bone volume (LBV, the difference between CBV and RBV) following infection. RESULTS: BALB/CJ mice were highly susceptible (P<0.05) whereas DBA/2J, C57BL/6J and A/J mice were resistant. Six lines of the tested CC population were susceptible, whereas the remaining lines were resistant to alveolar bone loss. Gender effects on bone volume were tested across the four inbred and 23 CC lines, and found not to be significant. Based on ANOVA analyses, broad-sense heritabilities were statistically significant and equal to 0.4 for CBV and 0.2 for LBV. CONCLUSIONS: The moderate heritability values indicate that the variation in host susceptibility to the disease is controlled to an appreciable extent by genetic factors. These results strongly support the possibility of using the Collaborative Cross, as well as developing dedicated F2 (resistant x susceptible inbred strains) resource populations, for future dissection of genetic factors in periodontitis.

Trypanotolerance in N'Dama x Boran crosses under natural trypanosome challenge: effect of test-year environment, gender, and breed composition.

BACKGROUND: Trypanosomosis, a protozoal disease affecting livestock, transmitted by Glossina (tsetse) flies is a major constraint to agricultural production in Sub-Saharan Africa. It is accepted that utilization of the native trypanotolerance exhibited in some of the African cattle breeds to improve trypanotolerance of more productive but susceptible breeds, will offer a cost effective and sustainable solution to the problem. The success of this approach is based on the premise that quantitative trait loci previously identified under relatively controlled situations confer useful trypanotolerance under natural field situations. As part of a study to authenticate this hypothesis, a population of 192 cattle, consisting of six batches of N'Dama and Kenya-Boran backcross animals [(N'Dama x Kenya-Boran) x Kenya-Boran] born over the period 2002 to 2006 was constructed. Some of the batches also included pure Kenya-Boran cattle, or N'Dama x Kenya- Boran F1 animals. Each batch was exposed as yearlings to natural field trypanosomosis challenge over a period of about one year; the entire challenge period extending from December 2003 to June 2007. Performance of the animals was evaluated by weekly or biweekly measurements of body weight, packed blood cell volume (PCV), parasitemia score, and number of trypanocide treatments. From these basic data, 49 phenotypes were constructed reflecting dynamics of body weight, packed cell volume (PCV) and parasitemia under challenge. RESULTS: Females were distinctly more trypanotolerant than males. F1, backcross and pure Kenya- Boran animals ranked in that order with respect to trypanotolerance. Overall batch effects were highly significant (p<0.001) for most traits, and were generally more significant than the gender or genetic type effects. The superior trypanotolerance of the F1 animals was expressed in all three components of animal defense strategies against pathogens: Avoidance resistance, and tolerance. CONCLUSIONS: The results show that trypanotolerance derived from the N'Dama is expressed under field conditions; and that the trait is primarily additive in nature, being expressed in heterozygous condition and in a three-quarters Boran genetic background. The results further, underscore the complexity of the trait in the field manifesting all three host disease-control strategies, and show the importance of gender and local environmental conditions in determining response to challenge.

Sex Differences in Response to Marek's Disease: Mapping Quantitative Trait Loci Regions (QTLRs) to the Z Chromosome.

Marek's Disease (MD) has a significant impact on both the global poultry economy and animal welfare. The disease pathology can include neurological damage and tumour formation. Sexual dimorphism in immunity and known higher susceptibility of females to MD makes the chicken Z chromosome (GGZ) a particularly attractive target to study the chicken MD response. Previously, we used a Hy-Line F6 population from a full-sib advanced intercross line to map MD QTL regions (QTLRs) on all chicken autosomes. Here, we mapped MD QTLRs on GGZ in the previously utilized F6 population with individual genotypes and phenotypes, and in eight elite commercial egg production lines with daughter-tested sires and selective DNA pooling (SDP). Four MD QTLRs were found from each analysis. Some of these QTLRs overlap regions from previous reports. All QTLRs were tested by individuals from the same eight lines used in the SDP and genotyped with markers located within and around the QTLRs. All QTLRs were confirmed. The results exemplify the complexity of MD resistance in chickens and the complex distribution of p-values and Linkage Disequilibrium (LD) pattern and their effect on localization of the causative elements. Considering the fragments and interdigitated LD blocks while using LD to aid localization of causative elements, one must look beyond the non-significant markers, for possible distant markers and blocks in high LD with the significant block. The QTLRs found here may explain at least part of the gender differences in MD tolerance, and provide targets for mitigating the effects of MD.

High resolution mapping of trypanosomosis resistance loci Tir2 and Tir3 using F12 advanced intercross lines with major locus Tir1 fixed for the susceptible allele.

BACKGROUND: Trypanosomosis is the most economically important disease constraint to livestock productivity in Africa. A number of trypanotolerant cattle breeds are found in West Africa, and identification of the genes conferring trypanotolerance could lead to effective means of genetic selection for trypanotolerance. In this context, high resolution mapping in mouse models are a promising approach to identifying the genes associated with trypanotolerance. In previous studies, using F2 C57BL/6J x A/J and C57BL/6J x BALB/cJ mouse resource populations, trypanotolerance QTL were mapped within a large genomic intervals of 20-40 cM to chromosomes MMU17, 5 and 1, and denoted Tir1, Tir2 and Tir3 respectively. Subsequently, using F6 C57BL/6J x A/J and C57BL/6J x BALB/cJ F6 advanced intercross lines (AIL), Tir1 was fine mapped to a confidence interval (CI) of less than 1 cM, while Tir2 and Tir3, were mapped within 5-12 cM. Tir1 represents the major trypanotolerance QTL. RESULTS: In order to improve map resolutions of Tir2 and Tir3, an F12 C57BL/6J x A/J AIL population fixed for the susceptible alleles at Tir1 QTL was generated. An F12 C57BL/6J x A/J AIL population, fixed for the resistant alleles at Tir1 QTL was also generated to provide an additional estimate of the gene effect of Tir1. The AIL populations homozygous for the resistant and susceptible Tir1 alleles and the parental controls were challenged with T. congolense and followed for survival times over 180 days. Mice from the two survival extremes of the F12 AIL population fixed for the susceptible alleles at Tir1 were genotyped with a dense panel of microsatellite markers spanning the Tir2 and Tir3 genomic regions and QTL mapping was performed. Tir2 was fine mapped to less than 1 cM CI while Tir3 was mapped to three intervals named Tir3a, Tir3b and Tir3c with 95% confidence intervals (CI) of 6, 7.2 and 2.2 cM, respectively. CONCLUSIONS: The mapped QTL regions encompass genes that are vital to innate immune response and can be potential candidate genes for the underlying QTL.

Mapping QTL Associated with Resistance to Avian Oncogenic Marek's Disease Virus (MDV) Reveals Major Candidate Genes and Variants.

Marek's disease (MD) represents a significant global economic and animal welfare issue. Marek's disease virus (MDV) is a highly contagious oncogenic and highly immune-suppressive α-herpes virus, which infects chickens, causing neurological effects and tumour formation. Though partially controlled by vaccination, MD continues to have a profound impact on animal health and on the poultry industry. Genetic selection provides an alternative and complementary method to vaccination. However, even after years of study, the genetic mechanisms underlying resistance to MDV remain poorly understood. The Major Histocompatability Complex (MHC) is known to play a role in disease resistance, along with a handful of other non-MHC genes. In this study, one of the largest to date, we used a multi-facetted approach to identify QTL regions (QTLR) influencing resistance to MDV, including an F6 population from a full-sib advanced intercross line (FSIL) between two elite commercial layer lines differing in resistance to MDV, RNA-seq information from virus challenged chicks, and genome wide association study (GWAS) from multiple commercial lines. Candidate genomic elements residing in the QTLR were further tested for association with offspring mortality in the face of MDV challenge in eight pure lines of elite egg-layer birds. Thirty-eight QTLR were found on 19 chicken chromosomes. Candidate genes, miRNAs, lncRNAs and potentially functional mutations were identified in these regions. Association tests were carried out in 26 of the QTLR, using eight pure lines of elite egg-layer birds. Numerous candidate genomic elements were strongly associated with MD resistance. Genomic regions significantly associated with resistance to MDV were mapped and candidate genes identified. Various QTLR elements were shown to have a strong genetic association with resistance. These results provide a large number of significant targets for mitigating the effects of MDV infection on both poultry health and the economy, whether by means of selective breeding, improved vaccine design, or gene-editing technologies.

Mapping QTL affecting resistance to Marek's disease in an F6 advanced intercross population of commercial layer chickens.

BACKGROUND: Marek's disease (MD) is a T-cell lymphoma of chickens caused by the Marek's disease virus (MDV), an oncogenic avian herpesvirus. MD is a major cause of economic loss to the poultry industry and the most serious and persistent infectious disease concern. A full-sib intercross population, consisting of five independent families was generated by crossing and repeated intercrossing of two partially inbred commercial White Leghorn layer lines known to differ in genetic resistance to MD. At the F6 generation, a total of 1615 chicks were produced (98 to 248 per family) and phenotyped for MD resistance measured as survival time in days after challenge with a very virulent plus (vv+) strain of MDV. RESULTS: QTL affecting MD resistance were identified by selective DNA pooling using a panel of 15 SNPs and 217 microsatellite markers. Since MHC blood type (BT) is known to affect MD resistance, a total of 18 independent pool pairs were constructed according to family x BT combination, with some combinations represented twice for technical reasons. Twenty-one QTL regions (QTLR) affecting post-challenge survival time were identified, distributed among 11 chromosomes (GGA1, 2, 3, 4, 5, 8, 9, 15, 18, 26 and Z), with about two-thirds of the MD resistance alleles derived from the more MD resistant parental line. Eight of the QTLR associated with MD resistance, were previously identified in a backcross (BC) mapping study with the same parental lines. Of these, 7 originated from the more resistant line, and one from the less resistant line. CONCLUSION: There was considerable evidence suggesting that MD resistance alleles tend to be recessive. The width of the QTLR for these QTL appeared to be reduced about two-fold in the F6 as compared to that found in the previous BC study. These results provide a firm basis for high-resolution linkage disequilibrium mapping and positional cloning of the resistance genes.

Designing a QTL Mapping Study for Implementation in the Realized Collaborative Cross Genetic Reference Population.

The Collaborative Cross (CC) mouse resource is a next-generation mouse genetic reference population (GRP) designed for high-resolution mapping of quantitative trait loci (QTL) of large effect affecting complex traits during health and disease. The CC resource consists of a set of 72 recombinant inbred lines (RILs) generated by reciprocal crossing of five classical and three wild-derived mouse founder strains. Complex traits are controlled by variations within multiple genes and environmental factors, and their mutual interactions. These traits are observed at multiple levels of the animals' systems, including metabolism, body weight, immune profile, and susceptibility or resistance to the development and progress of infectious or chronic diseases. Herein, we present general guidelines for design of QTL mapping experiments using the CC resource-along with full step-by-step protocols and methods that were implemented in our lab for the phenotypic and genotypic characterization of the different CC lines-for mapping the genes underlying host response to infectious and chronic diseases. © 2019 by John Wiley & Sons, Inc. Basic Protocol 1: CC lines for whole body mass index (BMI) Basic Protocol 2: A detailed assessment of the power to detect effect sizes based on the number of lines used, and the number of replicates per line Basic Protocol 3: Obtaining power for QTL with given target effect by interpolating in Table 1 of Keele et al. (2019).

Mapping Ethanol Tolerance in Budding Yeast Reveals High Genetic Variation in a Wild Isolate.

Ethanol tolerance, a polygenic trait of the yeast Saccharomyces cerevisiae, is the primary factor determining industrial bioethanol productivity. Until now, genomic elements affecting ethanol tolerance have been mapped only at low resolution, hindering their identification. Here, we explore the genetic architecture of ethanol tolerance, in the F6 generation of an Advanced Intercrossed Line (AIL) mapping population between two phylogenetically distinct, but phenotypically similar, S. cerevisiae strains (a common laboratory strain and a wild strain isolated from nature). Under ethanol stress, 51 quantitative trait loci (QTLs) affecting growth and 96 QTLs affecting survival, most of them novel, were identified, with high resolution, in some cases to single genes, using a High-Resolution Mapping Package of methodologies that provided high power and high resolution. We confirmed our results experimentally by showing the effects of the novel mapped genes: MOG1, MGS1, and YJR154W. The mapped QTLs explained 34% of phenotypic variation for growth and 72% for survival. High statistical power provided by our analysis allowed detection of many loci with small, but mappable effects, uncovering a novel "quasi-infinitesimal" genetic architecture. These results are striking demonstration of tremendous amounts of hidden genetic variation exposed in crosses between phylogenetically separated strains with similar phenotypes; as opposed to the more common design where strains with distinct phenotypes are crossed. Our findings suggest that ethanol tolerance is under natural evolutionary fitness-selection for an optimum phenotype that would tend to eliminate alleles of large effect. The study provides a platform for development of superior ethanol-tolerant strains using genome editing or selection.

Hormonal levels of estradiol, testosterone, and progesterone at entry into lay of year 1980 vs. 2000 broiler breeder females under fast and slow release from feed restriction.

In the mid-1960s egg production, fertility, and hatchability of broiler breeder females dropped precipitously. Due to disrupted follicle hierarchies and development of the erratic oviposition and defective eggs (EODES) syndrome. EODES was controlled by restricting feed. In the 1990s, another set of problems arose at entry of broiler breeders into lay and characterized by high mortality followed by lower peak lay and reduction in egg and chick production. These problems are induced by even slight over-feeding, and hence we termed it the "Over Feeding Complex" (OFC). We have speculated that OFC is a quasi-EODES condition, induced by the intense selection for increased breast proportion. To test this, we compared, under fast (FF) and slow (SF) release from feed restriction, body composition and reproductive performance of a broiler breeder from year 1980 (B1980) and kept without selection for performance traits since then, to a line hatched in 2000 (B2000). During the first 16 d of lay, feeding treatment had little effect on egg mass or Laying % for the B1980 birds, while for the B2000 birds, SF treatment resulted in significantly greater egg mass and Laying % compared to FF, showing that the OFC indeed manifested in this experiment. However, contrary to hypothesis, follicle hierarchies were normal for both lines under both feeding treatments. To gain further insight into the OFC syndrome, we here report levels of estradiol, testosterone, and progesterone for these line and treatment groups in the time period leading up to and into lay. A significant line × feeding treatment interaction effect was found for estradiol and testosterone, to a lesser extent for progesterone. For all 3 hormones, for B1980 levels 2 to 3 wk post entry into lay were similar and intermediate under FF and SF, but differed significantly for B2000, being much greater under SF than under FF. Thus, the hormonal effects were parallel and may explain the egg mass and Laying % effects of FF and SF in the 2 genetic types.

Mapping a syntenic modifier on mouse chromosome 1 influencing the expressivity of the compact phenotype in the myostatin mutant (MstnCmpt-dl1Abc) compact mouse.

A novel method for mapping a modifier gene that is syntenic to its major gene was used to map a male-sex-limited modifier of the expressivity of the Compact phenotype in the myostatin mutant (Mstn(Cmpt-dl1Abc)) Compact mouse. The modifier was mapped to the general region of D1Mit262, 40 cM distal to Mstn on chromosome 1. Myogenin, a postulated downstream target of myostatin, maps to the same region.

The Use of Kosher Phenotyping for Mapping QTL Affecting Susceptibility to Bovine Respiratory Disease.

Bovine respiratory disease (BRD) is the leading cause of morbidity and mortality in feedlot cattle, caused by multiple pathogens that become more virulent in response to stress. As clinical signs often go undetected and various preventive strategies failed, identification of genes affecting BRD is essential for selection for resistance. Selective DNA pooling (SDP) was applied in a genome wide association study (GWAS) to map BRD QTLs in Israeli Holstein male calves. Kosher scoring of lung adhesions was used to allocate 122 and 62 animals to High (Glatt Kosher) and Low (Non-Kosher) resistant groups, respectively. Genotyping was performed using the Illumina BovineHD BeadChip according to the Infinium protocol. Moving average of -logP was used to map QTLs and Log drop was used to define their boundaries (QTLRs). The combined procedure was efficient for high resolution mapping. Nineteen QTLRs distributed over 13 autosomes were found, some overlapping previous studies. The QTLRs contain polymorphic functional and expression candidate genes to affect kosher status, with putative immunological and wound healing activities. Kosher phenotyping was shown to be a reliable means to map QTLs affecting BRD morbidity.

Mapping modifiers affecting muscularity of the myostatin mutant (Mstn(Cmpt-dl1Abc)) compact mouse.

The hypermuscular Compact phenotype was first noted in a line of mice selected for high body weight and protein content. A new line, based on mice showing the Compact phenotype, was formed and selected for maximum expression of the Compact phenotype. Previously we mapped and identified a 12-bp deletion in the myostatin gene, denoted Mstn(Cmpt-dl1Abc), which can be considered as a major gene responsible for the hypermuscular phenotype. Genetic analysis revealed that full expression of the hypermuscular phenotype requires the action of modifier loci in addition to Mstn(Cmpt-dl1Abc). To map these modifier loci, an interspecific F(2) population was generated between Comp9, an inbred line homozygous for Mstn(Cmpt-dl1Abc), and CAST/Ei, an inbred line generated from Mus musculus castaneus. Selective DNA pooling and genotyping, separately by gender, was carried out within a subpopulation of the F(2) consisting of individuals homozygous for Mstn(Cmpt-dl1Abc). Significant association with hypermuscularity at a false discovery rate (FDR) of 0.05 was found for markers on chromosomes 3, 5, 7, 11, 16, and X. In all cases, the marker allele derived from the Comp9 parent showed a higher frequency in the hypermuscular group and the CAST/Ei allele in the normal group. The modifier loci apparently exerted their effects on muscularity only in the presence of Mstn(Cmpt-dl1Abc).

If a bull were a cow, how much milk would he give?

I became enamored of genetics at an early age. The desire to participate in the rebuilding of my people in our ancient homeland led to dairy science at Rutgers University, to Animal Breeding Plans by J.L. Lush, and to the realization that I could combine genetics and dairy science in animal breeding. It is to my mother-in-law that I owe the felicitous phrasing of the titular scientific question that has occupied my professional life: If a bull were a cow, how much milk would he give? Following my PhD (in 1956), I joined the Volcani Institute in Israel and, later (in 1972), the Applied Genetics group at the Hebrew University. The Applied Genetics group had an active marker lab, and this and a paper by Spickett & Thoday led me to explore genetic markers for quantitative trait loci mapping and marker-assisted selection. A chance encounter with Jacques Beckmann in 1980 opened my eyes to the potential of DNA-level markers for these purposes, and the rest followed.

Characterization and comparison of the leukocyte transcriptomes of three cattle breeds.

In this study, mRNA-Seq was used to characterize and compare the leukocyte transcriptomes from two taurine breeds (Holstein and Jersey), and one indicine breed (Cholistani). At the genomic level, we identified breed-specific base changes in protein coding regions. Among 7,793,425 coding bases, only 165 differed between Holstein and Jersey, and 3,383 (0.04%) differed between Holstein and Cholistani, 817 (25%) of which resulted in amino acid changes in 627 genes. At the transcriptional level, we assembled transcripts and estimated their abundances including those from more than 3,000 unannotated intergeneic regions. Differential gene expression analysis showed a high similarity between Holstein and Jersey, and a much greater difference between the taurine breeds and the indicine breed. We identified gene ontology pathways that were systematically altered, including the electron transport chain and immune response pathways that may contribute to different levels of heat tolerance and disease resistance in taurine and indicine breeds. At the post-transcriptional level, sequencing mRNA allowed us to identify a number of genes undergoing differential alternative splicing among different breeds. This study provided a high-resolution survey of the variation between bovine transcriptomes at different levels and may provide important biological insights into the phenotypic differentiation among cattle breeds.

Expression of trypanotolerance in N'Dama x Boran crosses under field challenge in relation to N'Dama genome content.

BACKGROUND: Animal trypanosomosis in sub-Saharan Africa is a major obstacle to livestock based agriculture. Control relies on drugs with increasing incidence of multiple-drug resistance. A previous mapping experiment in an F2 population derived from the indigenous trypanotolerant N'Dama cattle crossed to susceptible (Kenya)-Boran cattle under controlled challenge, uncovered a number of trypanotolerance QTL (T-QTL). The present study was to determine expression of N'Dama trypanotolerance in a backcross to the Boran under conditions of field challenge, and whether chromosomal regions associated with trypanotolerance in the F2 experiment showed similar effects in the BC population. METHODS: 192 backcross animals to the Boran were produced in six batches from June 2001 to December 2006. At one year of age animals were moved to the field and exposed to natural challenge over about one year in Southwest Kenya (Narok). The animals were individually recorded weekly for body weight, packed cell volume, parasitaemia score, and drug treatments, and were genotyped using 35 microsatellite markers spanning 5 chromosomes found in the F2 study to harbour T-QTL. RESULTS: The F1 were most trypanotolerant, Boran least, and BC intermediate. Females showed distinctly higher trypanotolerance than males. There was a positive correlation in the BC population between trypanotolerance and number of N'Dama origin marker alleles. QTL mapping revealed T-QTL distributed among all five targeted chromosomes, corresponding in part to the results obtained in the F2 experiment. CONCLUSIONS: N'Dama origin trypanotolerance is expressed in a BC population under field conditions in proportion to N'Dama origin marker alleles. Consequently, marker assisted selection in such populations may be a means of increasing trypanotolerance, while retaining the desirable productive qualities of the recurrent parent.

An analytical formula to estimate confidence interval of QTL location with a saturated genetic map as a function of experimental design.

Analytical formulae are derived for the confidence interval for location of a quantitative trait locus (QTL) using a saturated genetic map, as a function of the experimental design, the QTL allele substitution effect, and the number of individuals genotyped and phenotyped. The formulae are derived assuming evenly spaced recombination events, rather than the actual unevenly spaced distribution. The formulae are useful for determining desired sample size when designing a wide variety of QTL mapping experiments, and for evaluating a priori the potential of a given mapping population for defining the location of a QTL. The formulae do not take into account the finite number of recombination events in a given sample.