Bioinformatics-Based Identification of Candidate Genes from QTLs Associated with Cell Wall Traits in Populus

Quantitative trait locus (QTL) studies are an integral part of plant research and are used to characterize the genetic basis of phenotypic variation observed in structured populations and inform marker-assisted breeding efforts. These QTL intervals can span large physical regions on a chromosome comprising hundreds of genes, thereby hampering candidate gene identification. Genome history, evolution, and expression evidence can be used to narrow the genes in the interval to a smaller list that is manageable for detailed downstream functional genomics characterization. Our primary motivation for the present study was to address the need for a research methodology that identifies candidate genes within a broad QTL interval. Here we present a bioinformatics-based approach for subdividing candidate genes within QTL intervals into alternate groups of high probability candidates. Application of this approach in the context of studying cell wall traits, specifically lignin content and S/G ratios of stem and root in Populus plants, resulted in manageable sets of genes of both known and putative cell wall biosynthetic function. These results provide a roadmap for future experimental work leading to identification of new genes controlling cell wall recalcitrance and, ultimately, in the utility of plant biomass as an energy feedstock.     

Do large effect QTL fractionate? A case study at the maize domestication QTL teosinte branched1

Quantitative trait loci (QTL) mapping is a valuable tool for studying the genetic architecture of trait variation. Despite the large number of QTL studies reported in the literature, the identified QTL are rarely mapped to the underlying genes and it is usually unclear whether a QTL corresponds to one or multiple linked genes. Similarly, when QTL for several traits colocalize, it is usually unclear whether this is due to the pleiotropic action of a single gene or multiple linked genes, each affecting one trait. The domestication gene teosinte branched1 (tb1) was previously identified as a major domestication QTL with large effects on the differences in plant and ear architecture between maize and teosinte. Here we present the results of two experiments that were performed to determine whether the single gene tb1 explains all trait variation for its genomic region or whether the domestication QTL at tb1 fractionates into multiple linked QTL. For traits measuring plant architecture, we detected only one QTL per trait and these QTL all mapped to tb1. These results indicate that tb1 is the sole gene for plant architecture traits that segregates in our QTL mapping populations. For most traits related to ear morphology, we detected multiple QTL per trait in the tb1 genomic region, including a large effect QTL at tb1 itself plus one or two additional linked QTL. tb1 is epistatic to two of these additional QTL for ear traits. Overall, these results provide examples for both a major QTL that maps to a single gene, as well as a case in which a QTL fractionates into multiple linked QTL.     

基于元分析的抗玉米丝黑穗病QTL比较定位

以玉米遗传连锁图谱IBM2 2005 Neighbors为参考图谱,通过映射整合不同试验中的抗玉米丝黑穗病QTL,构建QTL综合图谱。在国内外种质中,共发现22个抗病QTL,分布在除第7染色体外的9条玉米染色体上。采用元分析技术,获得2个“一致性”抗病QTL,图距分别为8.79 cM和18.92cM。从MaizeGDB网站下载“一致性”QTL区间内基因和标记的原始序列;采用NCBI网站在线软件BLASTx通过同源比对在2个“一致性”QTL区间内初步获得4个抗病位置候选基因。借助比较基因电子定位策略,将69个水稻和玉米抗性基因定位于玉米IBM2图谱上,在2个“一致性”QTL区间内分别发现1个水稻抗性基因,初步推断为抗病位置候选基因。本文结果为抗玉米丝黑穗病QTL精细定位和分子育种提供了基础。     

Poor performance of bootstrap confidence intervals for the location of a quantitative trait locus

The aim of many genetic studies is to locate the genomic regions (called quantitative trait loci, QTL) that contribute to variation in a quantitative trait (such as body weight). Confidence intervals for the locations of QTL are particularly important for the design of further experiments to identify the gene or genes responsible for the effect. Likelihood support intervals are the most widely used method to obtain confidence intervals for QTL location, but the nonparametric bootstrap has also been recommended. Through extensive computer simulation, we show that bootstrap confidence intervals behave poorly and so should not be used in this context. The profile likelihood (or LOD curve) for QTL location has a tendency to peak at genetic markers, and so the distribution of the maximum-likelihood estimate (MLE) of QTL location has the unusual feature of point masses at genetic markers; this contributes to the poor behavior of the bootstrap. Likelihood support intervals and approximate Bayes credible intervals, on the other hand, are shown to behave appropriately.     

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.     

A genetical genomics approach to genome scans increases power for QTL mapping

We describe a method for integrating gene expression information into genome scans and show that this can substantially increase the statistical power of QTL mapping. The method has three stages. First, standard clustering methods identify small (size 5-20) groups of genes with similar expression patterns. Second, each gene group is tested for a causative genetic locus shared with the clinical trait of interest. This is done using an EM algorithm approach that treats genotype at the putative causative locus as an unobserved variable and combines expression information from all of the genes in the group to infer genotype information at the locus. Finally, expression QTL (eQTL) are mapped for each gene group that shares a causative locus with the clinical trait. Such eQTL are candidates for the causative locus. Simulation results show that this method has far superior power to standard QTL mapping techniques in many circumstances. We applied this method to existing data on mouse obesity. Our method identified 27 putative body weight QTL, whereas standard QTL mapping produced only one. Furthermore, most gene groups with body weight QTL included cis genes, so candidate genes could be immediately identified. Eleven body weight QTL produced 16 candidate genes that have been previously associated with body weight or body weight-related traits, thus validating our method. In addition, 15 of the 16 other loci produced 32 candidate genes that have not been associated with body weight. Thus, this method shows great promise for finding new causative loci for complex traits.     

Systems genetics approach reveals candidate genes for parasite resistance from quantitative trait loci studies in agricultural species

A systems genetics approach combining pathway analysis of quantitative trait loci (QTL) and gene expression information has provided strong evidence for common pathways associated with genetic resistance to internal parasites. Gene data, collected from published QTL regions in sheep, cattle, mice, rats and humans, and microarray data from sheep, were converted to human Entrez Gene IDs and compared to the KEGG pathway database. Selection of pathways from QTL data was based on a selection index that ensured that the selected pathways were in all species and the majority of the projects overall and within species. Pathways with either up- and down-regulated genes, primarily up-regulated genes or primarily down-regulated genes, were selected from gene expression data. After comparing the data sets independently, the pathways from each data set were compared and the common set of pathways and genes was identified. Comparisons within data sets identified 21 pathways from QTL data and 66 pathways from gene expression data. Both selected sets were enriched with pathways involved in immune functions, disease and cell responses to signals. The analysis identified 14 pathways that were common between QTL and gene expression data, and four directly associated with IFNγ or MHCII, with 31 common genes, including three MHCII genes. In conclusion, a systems genetics approach combining data from multiple QTL and gene expression projects led to the discovery of common pathways associated with genetic resistance to internal parasites. This systems genetics approach may prove significant for the discovery of candidate genes for many other multifactorial, economically important traits.     

Bioinformatics toolbox for narrowing rodent quantitative trait loci

Quantitative trait locus (QTL) analysis is a powerful method for localizing disease genes, but identifying the causal gene remains difficult. Rodent models of disease facilitate QTL gene identification, and causal genes underlying rodent QTL are often associated with the corresponding human diseases. Recently developed bioinformatics methods, including comparative genomics, combined cross analysis, interval-specific and genome-wide haplotype analysis, followed by sequence and expression analysis, each facilitated by public databases, provide new tools for narrowing rodent QTLs. Here we discuss each tool, illustrate its application and generate a bioinformatics strategy for narrowing QTLs. Combining these bioinformatics tools with classical experimental methods should accelerate QTL gene identification.     

Characterization of the porcine FABP5 gene and its association with the FAT1 QTL in an Iberian by Landrace cross

We have characterized and mapped the porcine fatty acid binding protein 5, epidermal (FABP5) gene. According to linkage and RH mapping, this gene is located close to the FABP4 (fatty acid binding protein 4, adipocyte) gene on swine chromosome 4. We resequenced 4.7 kb of the FABP5 gene in the parental population of an Iberian x Landrace cross (IBMAP), identifying seven SNPs arranged in two distinct FABP5 haplotypes. QTL and association analyses in the IBMAP population showed that this gene is strongly associated with fat deposition. QTL and haplotype analysis revealed that both FABP4 and FABP5 (clustered in mammals) are major candidate genes for the FAT1 QTL; the most likely position for the FAT1 QTL is between these two genes. Finally, our results suggest the presence of more than one QTL affecting fatness traits on porcine chromosome 4.     

Genetic and Molecular Basis of QTL of Diabetes in Mouse: Genes and Polymorphisms

A systematic study has been conducted of all available reports in PubMed and OMIM (Online Mendelian Inheritance in Man) to examine the genetic and molecular basis of quantitative genetic loci (QTL) of diabetes with the main focus on genes and polymorphisms. The major question is, What can the QTL tell us? Specifically, we want to know whether those genome regions differ from other regions in terms of genes relevant to diabetes. Which genes are within those QTL regions, and, among them, which genes have already been linked to diabetes? whether more polymorphisms have been associated with diabetes in the QTL regions than in the non-QTL regions.Our search revealed a total of 9038 genes from 26 type 1 diabetes QTL, which cover 667,096,006 bp of the mouse genomic sequence. On one hand, a large number of candidate genes are in each of these QTL; on the other hand, we found that some obvious candidate genes of QTL have not yet been investigated. Thus, the comprehensive search of candidate genes for known QTL may provide unexpected benefit for identifying QTL genes for diabetes. Key Words: Quantitative trait loci, type 1 diabetes, insulin-dependent diabetes mellitus (IDDM), candidate gene, polymorphism, mouse.     

Genome-Wide Association Studies Identifies Seven Major Regions Responsible for Iron Deficiency Chlorosis in Soybean (Glycine max)

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.     

作物QTL定位常用作图群体

作物大多重要农艺性状是数量性状,受多基因控制,基因之间及基因与环境之间都会发生互作,这为研究带来了很大的不便。因此,好的QTL作图群体,是研究QTL间的互作、QTL与环境的互作、QTL定位以及基因克隆的最根本保障。随着分子标记技术的发展,QTL定位的作图群体也在不断的发展并逐渐满足研究者对于QTL的精细定位及基因克隆等研究的进一步要求。文章主要综述了作物QTL定位常用作图群体的构建及优缺点。     

Evidence of Multiple Disease Resistance (MDR) and Implication of Meta-Analysis in Marker Assisted Selection

Meta-analysis was performed for three major foliar diseases with the aim to find out the total number of QTL responsible for these diseases and depict some real QTL for molecular breeding and marker assisted selection (MAS) in maize. Furthermore, we confirmed our results with some major known disease resistance genes and most well-known gene family of nucleotide binding site (NBS) encoding genes. Our analysis revealed that disease resistance QTL were randomly distributed in maize genome, but were clustered at different regions of the chromosomes. Totally 389 QTL were observed for these three major diseases in diverse maize germplasm, out of which 63 QTL were controlling more than one disease revealing the presence of multiple disease resistance (MDR). 44 real-QTLs were observed based on 4 QTL as standard in a specific region of genome. We also confirmed the Ht1 and Ht2 genes within the region of real QTL and 14 NBS-encoding genes. On chromosome 8 two NBS genes in one QTL were observed and on chromosome 3, several cluster and maximum MDR QTL were observed indicating that the apparent clustering could be due to genes exhibiting pleiotropic effect. Significant relationship was observed between the number of disease QTL and total genes per chromosome based on the reference genome B73. Therefore, we concluded that disease resistance genes are abundant in maize genome and these results can unleash the phenomenon of MDR. Furthermore, these results could be very handy to focus on hot spot on different chromosome for fine mapping of disease resistance genes and MAS.     

Identification of pleiotropic genes and gene sets underlying growth and immunity traits: a case study on Meishan pigs

Both growth and immune capacity are important traits in animal breeding. The animal quantitative trait loci (QTL) database is a valuable resource and can be used for interpreting the genetic mechanisms that underlie growth and immune traits. However, QTL intervals often involve too many candidate genes to find the true causal genes. Therefore, the aim of this study was to provide an effective annotation pipeline that can make full use of the information of Gene Ontology terms annotation, linkage gene blocks and pathways to further identify pleiotropic genes and gene sets in the overlapping intervals of growth-related and immunity-related QTLs. In total, 55 non-redundant QTL overlapping intervals were identified, 1893 growth-related genes and 713 immunity-related genes were further classified into overlapping intervals and 405 pleiotropic genes shared by the two gene sets were determined. In addition, 19 pleiotropic gene linkage blocks and 67 pathways related to immunity and growth traits were discovered. A total of 343 growth-related genes and 144 immunity-related genes involved in pleiotropic pathways were also identified, respectively. We also sequenced and genotyped 284 individuals from Chinese Meishan pigs and European pigs and mapped the single nucleotide polymorphisms (SNPs) to the pleiotropic genes and gene sets that we identified. A total of 971 high-confidence SNPs were mapped to the pleiotropic genes and gene sets that we identified, and among them 743 SNPs were statistically significant in allele frequency between Meishan and European pigs. This study explores the relationship between growth and immunity traits from the view of QTL overlapping intervals and can be generalized to explore the relationships between other traits.     

Investigation of Obesity Candidate Genes On Porcine Fat Deposition Quantitative Trait Loci Regions

Objectives: To investigate possible obesity candidate genes in regions of porcine quantitative trait loci (QTL) for fat deposition and obesity‐related phenotypes. Research Methods and Procedures: Chromosome mapping and QTL analyses of obesity candidate genes were performed using DNA panels from a reference pig family. Statistical association analyses of these genes were performed for fat deposition phenotypes in several other commercial pig populations. Results: Eight candidate genes were mapped to QTL regions of pig chromosomes in this study. These candidate genes also served as anchor loci to determine homologous human chromosomal locations of pig fat deposition QTL. Preliminary analyses of relationships among polymorphisms of individual candidate genes and a variety of phenotypic measurements in a large number of pigs were performed. On the basis of available data, gene‐gene interactions were also studied. Discussion: Comparative analysis of obesity‐related genes in the pig is not only important for development of marker‐assisted selection on growth and fat deposition traits in the pig but also provides for an understanding of their genetic roles in the development of human obesity.     

A potato molecular-function map for carbohydrate metabolism and transport

Molecular-linkage maps based on functional gene markers (molecular-function maps) are the prerequisite for a candidate-gene approach to identify genes responsible for quantitative traits at the molecular level. Genetic linkage between a quantitative trait locus (QTL) and a candidate-gene locus is observed when there is a causal relationship between alleles of the candidate gene and the QTL effect. Functional gene markers can also be used for marker-assisted selection and as anchors for structural and functional comparisons between distantly related plant species sharing the same metabolic pathways. A first molecular-function map with 85 loci was constructed in potato based on 69 genes. Priority was given to genes operating in carbohydrate metabolism and transport. Public databases were searched for genes of interest from potato, tomato, or other plant species. DNA sequence information was used to develop PCR-based marker assays that allowed the localization of corresponding potato genes on existing RFLP linkage maps. Comparing the molecular-function map for genes operating in carbohydrate metabolism and transport with a QTL map for tuber starch content indicates a number of putative candidate genes for this important agronomic trait. Received: 19 March 2000 / Accepted: 16 May 2000