High resolution mapping and identification of new quantitative trait loci (QTL) affecting susceptibility to Marek's disease

Marek's disease (MD) is a lymphoproliferative disease of chickens that costs the poultry industry approximately $1 billion annually. Genetic resistance to MD is gaining increased attention to augment vaccinal control as disease outbreaks occur more frequently. Previously, analysis of a 272 F2 White Leghorn resource population measured for many MD traits and genotyped for 78 microsatellite markers revealed two and four quantitative trait loci (QTL) with significant and suggestive association, respectively, to one or more MD associated traits. Additional genetic markers have since been scored on the MD resource population to increase QTL resolution and genome coverage. Saturation of four of the QTL regions with 17 markers revealed five new QTL while 32 markers extended the genome coverage by 400 + CM and uncovered three more QTL. QTL analysis by single-point and interval mapping algorithms agreed well when marker saturation was approximately 20 CM or less. Currently 127 genetic markers cover approximately 68% of the genome that contain up to 14 MD QTL associated to one or more MD trait; seven at the significant level and seven at the suggestive level. Individually each QTL accounts for 2-10% of the variation and, in general, resistance was dominant although the resistant allele may come from either parental line. This study suggests that a limited number of genomic regions play a major role in the genetic control of MD resistance. Markers linked to these loci may be useful for selection of MD resistant stock by the poultry industry following verification of the association within their breeding populations.     

Mapping quantitative trait loci affecting susceptibility to Marek's disease virus in a backcross population of layer chickens

Marek's disease (MD), caused by the oncogenic MD avian herpes virus (MDV), is a major source of economic losses to the poultry industry. A reciprocal backcross (BC) population (total 2052 individuals) was generated by crossing two partially inbred commercial Leghorn layer lines known to differ in MDV resistance, measured as survival time after challenge with a (vv+) MDV. QTL affecting resistance were identified by selective DNA pooling using a panel of 198 microsatellite markers covering two-thirds of the chicken genome. Data for each BC were analyzed separately, and as a combined data set. Markers showing significant association with resistance generally appeared in blocks of two or three, separated by blocks of nonsignificant markers. Defined this way, 15 chromosomal regions (QTLR) affecting MDV resistance, distributed among 10 chromosomes (GGA 1, 2, 3, 4, 5, 7, 8, 9, 15, and Z), were identified. The identified QTLR include one gene and three QTL associated with resistance in previous studies of other lines, and three additional QTL associated with resistance in previous studies of the present lines. These QTL could be used in marker-assisted selection (MAS) programs for MDV resistance and as a platform for high-resolution mapping and positional cloning of the resistance genes.     

Identification of chicken lymphocyte antigen 6 complex, locus E (LY6E, alias SCA2) as a putative Marek's disease resistance gene via a virus-host protein interaction screen

Marek's disease virus (MDV) is a naturally occurring oncogenic avian herpesvirus that causes neurological disorders and T cell lymphoma disease in domestic chickens. Identification and functional characterization of the individual factors involved in Marek's disease (MD) resistance or pathogenesis will enhance our understanding of MDV pathogenesis and further genetic improvement of chickens. To study the genetic basis for resistance to MD, a strategy that combined protein-protein interaction screens followed by linkage analysis was performed. The MDV protein US10 was used as the bait in an E. COLI two-hybrid screening of a cDNA library derived from activated splenic T cells. The chicken LY6E, also known as SCA2 and TSA1, was found to specifically interact with US10. This interaction was confirmed by an in vitro protein-binding assay. Furthermore, LY6E was found to be significantly associated with MD traits in an MD resource population comprised of commercial chickens. Previously, LY6E was implicated in two independent DNA microarray experiments evaluating differential gene expression following MDV infection. Given that LY6E is involved in T cell differentiation and activation, we suggest that LY6E is a candidate gene for MD resistance and deserves further investigation on its role in MDV pathogenesis, especially with respect to the binding of US10.     

A strategy to identify positional candidate genes conferring Marek''s disease resistance by integrating DNA microarrays and genetic mapping

Marker-assisted selection (MAS) to enhance genetic resistance to Marek's disease (MD), a herpesvirus-induced T cell cancer in chicken, is an attractive alternative to augment control with vaccines. Our earlier studies indicate that there are many quantitative trait loci (QTL) containing one or more genes that confer genetic resistance to MD. Unfortunately, it is difficult to sufficiently resolve these QTL to identify the causative gene and generate tightly linked markers. One possible solution is to identify positional candidate genes by virtue of gene expression differences between MD resistant and susceptible chicken using deoxyribonucleic acid (DNA) microarrays followed by genetic mapping of the differentially-expressed genes. In this preliminary study, we show that DNA microarrays containing approximately 1200 genes or expressed sequence tags (ESTs) are able to reproducibly detect differences in gene expression between the inbred ADOL lines 63 (MD resistant) and 72 (MD susceptible) of uninfected and Marek's disease virus (MDV)-infected peripheral blood lymphocytes. Microarray data were validated by quantitative polymerase chain reaction (PCR) and found to be consistent with previous literature on gene induction or immune response. Integration of the microarrays with genetic mapping data was achieved with a sample of 15 genes. Twelve of these genes had mapped human orthologues. Seven genes were located on the chicken linkage map as predicted by the human-chicken comparative map, while two other genes defined a new conserved syntenic group. More importantly, one of the genes with differential expression is known to confer genetic resistance to MD while another gene is a prime positional candidate for a QTL.     

Genetic mapping of quantitative trait loci affecting susceptibility to Marek's disease virus induced tumors in F2 intercross chickens.

Marek''s disease (MD) is a lymphoproliferative disease caused by the MD virus (MDV), which costs the poultry industry nearly $1 billion annually. To identify quantitative trait loci (QTL) affecting MD susceptibility, the inbred lines 6(3) (MD resistant) and 7(2) (MD susceptible) were mated to create more than 300 F2 chickens. The F2 chickens were challenged with MDV JM strain, moderately virulent) at 1 wk of age and assessed for MD susceptibility. The QTL analysis was divided into three stages. In stage 1, 65 DNA markers selected from the chicken genetic maps were typed on the 40 most MD-susceptible and the 40 most MD-resistant F2 chickens, and 21 markers residing near suggestive QTL were revealed by analysis of variance (ANOVA). In stage 2, the suggestive markers plus available flanking markers were typed on 272 F2 chickens, and three suggestive QTL were identified by ANOVA. In stage 3, using the interval mapping program Map Manager and permutation tests, two significant and two suggestive MD QTL were identified on four chromosomal subregions. Three to five loci collected explained between 11 and 23% of the phenotypic MD variation, or 32-68% of the genetic variance. This study constitutes the first report in the domestic chicken on the mapping of non-major histocompatibility complex QTL affecting MD susceptibility.     

玉米开花期相关性状的QTL分析

利用玉米强优势组合(Mo17×黄早四)自交衍生的191个F_2单株构建了由SSR和AFLP标记组成的分子连锁图谱,用F2进一步自交产生的184个F_(2:3)家系调查散粉期、吐丝期和开花-吐丝间隔期(ASI)的表型值,采用基于混合线性模型的复合区间作图法和相应的作图软件QTLmapper/V2.0,在两个生长环境下定位了与散粉期、吐丝期和ASI相关的QTL数目分别为13、7和5个,检测到3对控制散粉期、17对控制吐丝期和5对控制ASI的上位性效应位点;同时发现了与环境存在显著互作的3个散粉期、3个吐丝期和2个ASI单位点标记区域以及1对散粉期、3对吐丝期和2对ASI上位性效应区域.对玉米散粉期、吐丝期和ASI遗传基础中遗传因素相对作用大小分析表明,加性效应、部分显性效应和上位性效应是玉米开花期相关性状的重要遗传基础.     

Simultaneous mapping of epistatic QTL in DU6i × DBA/2 mice

We have mapped epistatic quantitative trait loci (QTL) in an F₂ cross between DU6i × DBA/2 mice. By including these epistatic QTL and their interaction parameters in the genetic model, we were able to increase the genetic variance explained substantially (8.8%–128.3%) for several growth and body composition traits. We used an analysis method based on a simultaneous search for epistatic QTL pairs without assuming that the QTL had any effect individually. We were able to detect several QTL that could not be detected in a search for marginal QTL effects because the epistasis cancelled out the individual effects of the QTL. In total, 23 genomic regions were found to contain QTL affecting one or several of the traits and eight of these QTL did not have significant individual effects. We identified 44 QTL pairs with significant effects on the traits, and, for 28 of the pairs, an epistatic QTL model fit the data significantly better than a model without interactions. The epistatic pairs were classified by the significance of the epistatic parameters in the genetic model, and visual inspection of the two-locus genotype means identified six types of related genotype–phenotype patterns among the pairs. Five of these patterns resembled previously published patterns of QTL interactions.     

Epistatic interaction is an important genetic basis of grain yield and its components in maize

A population of 294 recombinant inbred lines (RIL) derived from Yuyu22, an elite maize hybrid extending broadly in China, has been constructed to investigate the genetic basis of grain yield, and associated yield components in maize. The main-effect quantitative trait loci (QTL), digenic epistatic interactions, and their interactions with the environment for grain yield and its three components were identified by using the mixed linear model approach. Thirty-two main-effect QTL and forty-four pairs of digenic epistatic interactions were detected for the four measured traits in four environments. Our results suggest that both additive effects and epistasis (additive × additive) effects are important genetic bases of grain yield and its components in the RIL population. Only 30.4% of main-effect QTL for ear length were involved in epistatic interactions. This implies that many loci in epistatic interactions may not have significant effects for traits alone but may affect trait expression by epistatic interaction with the other loci.     

Major and minor QTL and epistasis contribute to fatty acid compositions and oil concentration in high-oil maize

High-oil maize is a useful genetic resource for genomic investigation in plants. To determine the genetic basis of oil concentration and composition in maize grain, a recombinant inbred population derived from a cross between normal line B73 and high-oil line By804 was phenotyped using gas chromatography, and genotyped with 228 molecular markers. A total of 42 individual QTL, associated with fatty acid compositions and oil concentration, were detected in 21 genomic regions. Five major QTL were identified for measured traits, one each of which explained 42.0% of phenotypic variance for palmitic acid, 15.0% for stearic acid, 27.7% for oleic acid, 48.3% for linoleic acid, and 15.7% for oil concentration in the RIL population. Thirty-six loci were involved in 24 molecular marker pairs of epistatic interactions across all traits, which explained phenotypic variances ranging from 0.4 to 6.1%. Seven of 18 mapping candidate genes related to lipid metabolism were localized within or were close to identified individual QTL, explaining 0.7–13.2% of the population variance. These results demonstrated that a few major QTL with large additive effects could play an important role in attending fatty acid compositions and increasing oil concentration in used germplasm. A larger number of minor QTL and a certain number of epistatic QTL, both with additive effects, also contributed to fatty acid compositions and oil concentration.     

Epistatic QTL pairs associated with meat quality and carcass composition traits in a porcine Duroc × Pietrain population

Background

Quantitative trait loci (QTL) analyses in pig have revealed numerous individual QTL affecting growth, carcass composition, reproduction and meat quality, indicating a complex genetic architecture. In general, statistical QTL models consider only additive and dominance effects and identification of epistatic effects in livestock is not yet widespread. The aim of this study was to identify and characterize epistatic effects between common and novel QTL regions for carcass composition and meat quality traits in pig.

Methods

Five hundred and eighty five F₂ pigs from a Duroc × Pietrain resource population were genotyped using 131 genetic markers (microsatellites and SNP) spread over the 18 pig autosomes. Phenotypic information for 26 carcass composition and meat quality traits was available for all F₂ animals. Linkage analysis was performed in a two-step procedure using a maximum likelihood approach implemented in the QxPak program.

Results

A number of interacting QTL was observed for different traits, leading to the identification of a variety of networks among chromosomal regions throughout the porcine genome. We distinguished 17 epistatic QTL pairs for carcass composition and 39 for meat quality traits. These interacting QTL pairs explained up to 8% of the phenotypic variance.

Conclusions

Our findings demonstrate the significance of epistasis in pigs. We have revealed evidence for epistatic relationships between different chromosomal regions, confirmed known QTL loci and connected regions reported in other studies. Considering interactions between loci allowed us to identify several novel QTL and trait-specific relationships of loci within and across chromosomes.     

Identification of main effect and epistatic quantitative trait loci for morphological and yield-related traits in peanut (<Emphasis Type="Italic">Arachis hypogaea</Emphasis> L.)

An effort was made in the present study to identify the main effect and epistatic quantitative trait locus (QTL) for the morphological and yield-related traits in peanut. A recombinant inbred line (RIL) population derived from TAG 24 × GPBD 4 was phenotyped in seven environments at two locations. QTL analysis with available genetic map identified 62 main-effect QTLs (M-QTLs) for ten morphological and yield-related traits with the phenotypic variance explained (PVE) of 3.84–15.06%. Six major QTLs (PVE >?10%) were detected for PLHT, PPP, YPP, and SLNG. Stable M-QTLs appearing in at least two environments were detected for PLHT, LLN, YPP, YKGH, and HSW. Five M-QTLs governed two traits each, and 16 genomic regions showed co-localization of two to four M-QTLs. Intriguingly, a major QTL reported to be linked to rust resistance showed pleiotropic effect for yield-attributing traits like YPP (15.06%, PVE) and SLNG (13.40%, PVE). Of the 24 epistatic interactions identified across the traits, five interactions involved six M-QTLs. Three interactions were additive × additive and remaining two involved QTL × environment (QE) interactions. Only one major M-QTL governing PLHT showed epistatic interaction. Overall, this study identified the major M-QTLs for the important productivity traits and also described the lack of epistatic interactions for majority of them so that they can be conveniently employed in peanut breeding.     

Exploring the quantitative resistance to <Emphasis Type="Italic">Pseudomonas syringae</Emphasis> pv. <Emphasis Type="Italic">phaseolicola</Emphasis> in common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.)

Pseudomonas syringae pv. phaseolicola is an important disease that causes halo blight in common bean. The genetic mechanisms underlying quantitative halo blight resistance are poorly understood in this species, as most disease studies have focused on qualitative resistance. The present work examines the genetic basis of quantitative resistance to the nine halo blight races in different organs (primary and trifoliate leaf, stem and pod) of an Andean recombinant inbred line (RIL) progeny. Using a multi-environment quantitative trait locus (QTL) mapping approach, 76 and 101 main-effect and epistatic QTLs were identified, respectively. Most of the epistatic interactions detected were due to loci without detectable QTL additive main effects. Main and epistatic QTLs detected were mainly consistent across the environment conditions. The homologous genomic regions corresponding to 26 of the 76 main-effect detected QTLs were positive for the presence of resistance-associated gene cluster encoding nucleotide-binding and leucine-rich repeat (NL) proteins and known defence genes. Main-effect QTLs for resistance to races 3, 4 and 5 in leaf, stem and pod were located on chromosome 2 within a 3.01-Mb region, where a cluster of nine NL genes was detected. The NL gene Phvul.002G323300 is located in this region, which can be considered an important putative candidate gene for the non-organ-specific QTL identified here. The present research provides essential information not only for the better understanding of the plant-pathogen interaction but also for the application of genomic assisted breeding for halo blight resistance in common bean.     

Horizontal transmission of Marek's disease virus requires US2, the UL13 protein kinase, and gC

Marek's disease virus (MDV) causes a general malaise in chickens that is mostly characterized by the development of lymphoblastoid tumors in multiple organs. The use of bacterial artificial chromosomes (BACs) for cloning and manipulation of the MDV genome has facilitated characterization of specific genes and genomic regions. The development of most MDV BACs, including pRB-1B-5, derived from a very virulent MDV strain, involved replacement of the US2 gene with mini-F vector sequences. However, when reconstituted viruses based on pRB-1B were used in pathogenicity studies, it was discovered that contact chickens housed together with experimentally infected chickens did not contract Marek's disease (MD), indicating a lack of horizontal transmission. Staining of feather follicle epithelial cells in the skins of infected chickens showed that virus was present but was unable to be released and/or infect susceptible chickens. Restoration of US2 and removal of mini-F sequences within viral RB-1B did not alter this characteristic, although in vivo viremia levels were increased significantly. Sequence analyses of pRB-1B revealed that the UL13, UL44, and US6 genes encoding the UL13 serine/threonine protein kinase, glycoprotein C (gC), and gD, respectively, harbored frameshift mutations. These mutations were repaired individually, or in combination, using two-step Red mutagenesis. Reconstituted viruses were tested for replication, MD incidence, and their abilities to horizontally spread to contact chickens. The experiments clearly showed that US2, UL13, and gC in combination are essential for horizontal transmission of MDV and that none of the genes alone is able to restore this phenotype.     

Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse

The extent to which epistasis contributes to adaptation, population differentiation, and speciation is a long-standing and important problem in evolutionary genetics. Using recombinant inbred (RI) lines of Arabidopsis thaliana grown under natural field conditions, we have examined the genetic architecture of fitness-correlated traits with respect to epistasis; we identified both single-locus additive and two-locus epistatic QTL for natural variation in fruit number, germination, and seed length and width. For fruit number, we found seven significant epistatic interactions, but only two additive QTL. For seed germination, length, and width, there were from two to four additive QTL and from five to eight epistatic interactions. The epistatic interactions were both positive and negative. In each case, the magnitude of the epistatic effects was roughly double that of the effects of the additive QTL, varying from -41% to +29% for fruit number and from -5% to +4% for seed germination, length, and width. A number of the QTL that we describe participate in more than one epistatic interaction, and some loci identified as additive also may participate in an epistatic interaction; the genetic architecture for fitness traits may be a network of additive and epistatic effects. We compared the map positions of the additive and epistatic QTL for germination, seed width, and seed length from plants grown in both the field and the greenhouse. While the total number of significant additive and epistatic QTL was similar under the two growth conditions, the map locations were largely different. We found a small number of significant epistatic QTL x environment effects when we tested directly for them. Our results support the idea that epistatic interactions are an important part of natural genetic variation and reinforce the need for caution in comparing results from greenhouse-grown and field-grown plants.     

The genetic basis of the interspecific differences in wing size in Nasonia (Hymenoptera; Pteromalidae): major quantitative trait loci and epistasis

There is a 2.5-fold difference in male wing size between two haplodiploid insect species, Nasonia vitripennis and N. giraulti. The haploidy of males facilitated a full genomic screen for quantitative trait loci (QTL) affecting wing size and the detection of epistatic interactions. A QTL analysis of the interspecific wing-size difference revealed QTL with major effects and epistatic interactions among loci affecting the trait. We analyzed 178 hybrid males and initially found two major QTL for wing length, one for wing width, three for a normalized wing-size variable, and five for wing seta density. One QTL for wing width explains 38.1% of the phenotypic variance, and the same QTL explains 22% of the phenotypic variance in normalized wing size. This corresponds to a region previously introgressed from N. giraulti into N. vitripennis that accounts for 44% of the normalized wing-size difference between the species. Significant epistatic interactions were also found that affect wing size and density of setae on the wing. Screening for pairwise epistatic interactions between loci on different linkage groups revealed four additional loci for wing length and four loci for normalized wing size that were not detected in the original QTL analysis. We propose that the evolution of smaller wings in N. vitripennis males is primarily the result of major mutations at few genomic regions and involves epistatic interactions among some loci.     

An epistatic genetic basis for physical activity traits in mice

We recently identified several (4-8) quantitative trait loci (QTL) for 3 physical activity traits (daily distance, duration, and speed voluntarily run) in an F(2) population of mice derived from an original intercross of 2 strains that exhibited large differences in activity. These QTL cumulatively explained from 11% to 34% of the variation in these traits, but this was considerably less than their total genetic variability estimated from differences among inbred strains. We therefore decided to test whether epistatic interactions might account for additional genetic variation in these traits in this same population of mice. We conducted a full genome epistasis scan for all possible interactions of QTL between each pair of 20 chromosomes. The results of this scan revealed an abundance of epistasis, with QTL throughout the genome being involved in significant interactions. Overall, epistatic effects contributed an average of 26% of the total variation among the 3 activity traits. These results suggest that epistatic interactions of genes may play as important a role in the genetic architecture of physical activity traits as single-locus effects and need to be considered in future candidate gene identification studies.     

The main effects, epistatic effects and environmental interactions of QTLs on the cooking and eating quality of rice in a doubled-haploid line population

Amylose content (AC), gel consistency (GC) and gelatinazation temperature (GT) are three important traits that influence the cooking and eating quality of rice. The objective of this study was to characterize the genetic components, including main-effect quantitative trait loci (QTLs), epistatic QTLs and QTL-by-environment interactions (QEs), that are involved in the control of these three traits. A population of doubled haploid (DH) lines derived from a cross between two indica varieties Zhenshan 97 and H94 was used, and data were collected from a field experiment conducted in two different environments. A genetic linkage map consisting of 218 simple sequence repeat (SSR) loci was constructed, and QTL analysis performed using qtlmapper 1.6 resolved the genetic components into main-effect QTLs, epistatic QTLs and QEs. The analysis detected a total of 12 main-effect QTLs for the three traits, with a QTL corresponding to the Wx locus showing a major effect on AC and GC, and a QTL corresponding to the Alk locus having a major effect on GT. Ten digenic interactions involving 19 loci were detected for the three traits, and six main-effect QTLs and two pairs of epistatic QTLs were involved in QEs. While the main-effect QTLs, especially the ones corresponding to known major loci, apparently played predominant roles in the genetic basis of the traits, under certain conditions epistatic effects and QEs also played important roles in controlling the traits. The implications of the findings for rice quality improvement are discussed.     

The role of epistasis in the manifestation of heterosis: a systems-oriented approach

Heterosis is widely used in breeding, but the genetic basis of this biological phenomenon has not been elucidated. We postulate that additive and dominance genetic effects as well as two-locus interactions estimated in classical QTL analyses are not sufficient for quantifying the contributions of QTL to heterosis. A general theoretical framework for determining the contributions of different types of genetic effects to heterosis was developed. Additive x additive epistatic interactions of individual loci with the entire genetic background were identified as a major component of midparent heterosis. On the basis of these findings we defined a new type of heterotic effect denoted as augmented dominance effect di* that comprises the dominance effect at each QTL minus half the sum of additive x additive interactions with all other QTL. We demonstrate that genotypic expectations of QTL effects obtained from analyses with the design III using testcrosses of recombinant inbred lines and composite-interval mapping precisely equal genotypic expectations of midparent heterosis, thus identifying genomic regions relevant for expression of heterosis. The theory for QTL mapping of multiple traits is extended to the simultaneous mapping of newly defined genetic effects to improve the power of QTL detection and distinguish between dominance and overdominance.