ezQTL: A Web Platform for Interactive Visualization and Colocalization of QTLs and GWAS Loci

Genome-wide association studies (GWAS) have identified thousands of genomic loci associated with complex diseases and traits, including cancer. The vast majority of common trait-associated variants identified via GWAS fall in non-coding regions of the genome, posing a challenge in elucidating the causal variants, genes, and mechanisms involved. Expression quantitative trait locus (eQTL) and other molecular QTL studies have been valuable resources in identifying candidate causal genes from GWAS loci through statistical colocalization methods. While QTL colocalization is becoming a standard analysis in post-GWAS investigation, an easy web tool for users to perform formal colocalization analyses with either user-provided or public GWAS and eQTL datasets has been lacking. Here, we present ezQTL, a web-based bioinformatic application to interactively visualize and analyze genetic association data such as GWAS loci and molecular QTLs under different linkage disequilibrium (LD) patterns (1000 Genomes Project, UK Biobank, or user-provided data). This application allows users to perform data quality control for variants matched between different datasets, LD visualization, and two-trait colocalization analyses using two state-of-the-art methodologies (eCAVIAR and HyPrColoc), including batch processing. ezQTL is a free and publicly available cross-platform web tool, which can be accessed online at https://analysistools.cancer.gov/ezqtl.     

Comparative mapping between<Emphasis Type="Italic"> Medicago sativa</Emphasis> and<Emphasis Type="Italic"> Pisum sativum</Emphasis>

Comparative genome analysis has been performed between alfalfa ( Medicago sativa) and pea ( Pisum sativum), species which represent two closely related tribes of the subfamily Papilionoideae with different basic chromosome numbers. The positions of genes on the most recent linkage map of diploid alfalfa were compared to those of homologous loci on the combined genetic map of pea to analyze the degree of co-linearity between their linkage groups. In addition to using unique genes, analysis of the map positions of multicopy (homologous) genes identified syntenic homologs (characterized by similar positions on the maps) and pinpointed the positions of non-syntenic homologs. The comparison revealed extensive conservation of gene order between alfalfa and pea. However, genetic rearrangements (due to breakage and reunion) were localized which can account for the difference in chromosome number (8 for alfalfa and 7 for pea). Based on these genetic events and our increasing knowledge of the genomic structure of pea, it was concluded that the difference in genome size between the two species (the pea genome is 5- to 10-fold larger than that of alfalfa) is not a consequence of genome duplication in pea. The high degree of synteny observed between pea and Medicago loci makes further map-based cloning of pea genes based on the genome resources now available for M. truncatula a promising strategy.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by W. R. McCombie     

Haplotype block and superblock structures of the alpha1-adrenergic receptor genes reveal echoes from the chromosomal past

A significant proportion of the human genome is contained within haplotype blocks across which pairwise linkage disequilibrium (LD) is very high. However, LD is also often high between markers at more remote distances, and within different haplotype blocks. Here, we evaluate the origins of haplotype block structure in the three genes for alpha1 adrenergic receptors (alpha1-AR) in the human genome ( ADRA1A, ADRA1B and ADRA1D) by genotyping dense single-nucleotide polymorphism (SNP) marker maps, and show that LD signals between distant markers are due to the presence of extended haplotype superblocks in individuals with ancient chromosomes which have escaped historic recombination. ARs mediate the physiological effects of epinephrine and norepinephrine, and are targets of many therapeutic drugs. This work has identified haplotype backgrounds of alpha1-AR missense variants, haplotype block structures in US Caucasians and African Americans, and haplotype tag SNPs for each block, and we present strong evidence for ancient haplotype block superstructure at these genes which has been partially disrupted by recombination, and evidence for reinstatement of linkage disequilibrium by subsequent recombination events. ADRA1A is comprised of four haplotype blocks in US Caucasians, while in African Americans Block 1 is split. ADRA1B has four blocks in US Caucasians, but in African Americans only the first two blocks are present. ADRA1D has two blocks in US Caucasians, and the first block is replaced by two smaller blocks in African Americans. For both ADRA1A and ADRA1B, haplotype superstructures may represent a novel, higher-level hierarchy in the human genome, which may reduce redundancy of testing by further aggregation of genotype data.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by W. R. McCombie     

Linkage disequilibrium and founder effect analysis of the NF1 gene in French Canadians from the Quebec population

We genotyped 19 neurofibromatosis type 1 (NF1) families from French Canadians of the Quebec population with four intragenic microsatellites (IVS26-2.3, IVS27AC28.4, IVS27AC33.1, and IVS38GT53.0). Linkage analysis of the four microsatellite markers among the 19 NF1 families indicates that the four microsatellites are strongly linked with NF1 disease (LOD = 2.76-3.64). The four markers are associated (P = 0-0.077) except marker pair IVS26-2.3/IVS27AC33.1 (P = 0.18 or 0.17). However, perhaps due to the high mutation rate of the NF1 gene, no founder effect for NF1 was detected in the Quebec French Canadians.     

Fine linkage mapping enables dissection of closely linked quantitative trait loci for seed dormancy and heading in rice

Two quantitative trait loci (QTLs) for seed dormancy (tentatively designated Sdr1) and heading date (Hd8) have been mapped to approximately the same region on chromosome 3 by interval mapping of backcross inbred lines derived from crosses between the rice cultivars Nipponbare (japonica) and Kasalath (indica). To clarify whether Sdr1 and Hd8 could be dissected genetically, we carried out fine-scale mapping with an advanced backcross progeny. We selected a BC₄F₁ plant, in which a small chromosomal region including Sdr1 and Hd8, on the short arm of chromosome 3, remained heterozygous, whereas all the other chromosomal regions were homozygous for Nipponbare. Days-to-heading and seed germination rate in the BC₄F₂ plants showed continuous variation. Ten BC₄F₂ plants with recombination in the vicinity of Sdr1 and Hd8 were selected on the basis of the genotypes of the restriction fragment length polymorphism (RFLP) markers flanking both QTLs. Genotypes of those plants for Sdr1 and Hd8 were determined by advanced progeny testing of BC₄F₄ families. Sdr1 was mapped between the RFLP markers R10942 and C2045, and co-segregated with C1488. Hd8 was also mapped between C12534S and R10942. Six recombination events were detected between Sdr1 and Hd8. These results clearly demonstrate that Sdr1 and Hd8 were tightly linked. Nearly isogenic lines for Sdr1 and Hd8 were selected by marker-assisted selection.Communicated by D. Mackill     

Patterns of polymorphism detected in the chloroplast and nuclear genomes of barley landraces sampled from Syria and Jordan

In order to examine how molecular polymorphism in barley landraces, sampled from five different ecogeographical regions of Syria and Jordan, is organised and partitioned, genetic variability at 21 nuclear and 10 chloroplast microsatellite loci were examined. Chloroplast polymorphism was detected, with most variation being ascribed to differences between the five regions (F_st 0.45) and to within sites within each region (F_st 0.44). Moreover, the distribution of chloroplast polymorphism is structured and not distributed randomly across the barley landraces sampled. From a total of 125 landrace accessions (five lines from each of five sites from each of five regions) genotyped with 21 SSRs a total of 244 alleles were detected, of which 38 were common to the five regions sampled. Most nuclear variation was detected within sites. Significant differentiation between sites (F_st 0.29) was detected with nuclear SSRs and this partially mirrored polymorphism in the chloroplast genome. Strong statistical associations/interaction was also detected between the chloroplast and nuclear SSRs, together with non-random association (linkage disequilibrium) of alleles at both linked and unlinked SSR loci. These results are discussed in the context of adaptation of landraces to the extreme environment, the concept of 'adapted gene complexes' and the exploitation of landraces in breeding programmes.Communicated by P. Langridge     

Identification and characterization of coding single-nucleotide polymorphisms within a human olfactory receptor gene cluster

Single-nucleotide polymorphisms (SNPs) were studied in 15 olfactory receptor (OR) coding regions, one control region and two noncoding sequences all residing within a 412 kb OR gene cluster on human chromosome 17p13.3, as well as in other G-protein coupled receptors (GPCRs). A total of 26 SNPs were identified in ORs, 21 of which are coding SNPs (cSNPs). The mean nucleotide diversity of OR coding regions was 0.078% (ranging from 0 to 0.16%), which is about twice higher than that of other GPCRs, and similar to the nucleotide diversity levels of noncoding regions along the human genome. The high polymorphism level in the OR coding regions might be due to a weak positive selection pressure acting on the OR genes. In two cases, OR genes have been found to share the same cSNP. This could be explained by recent gene conversion events, which might be a part of a concerted evolution mechanism acting on the OR superfamily. Using the genotype data of 85 unrelated individuals in 15 SNPs, we found linkage disequilibrium (LD) between pairs of SNPs located on the centromeric part of the cluster. On the other hand, no LD was found between SNPs located on the telomeric part of the cluster, suggesting the presence of several hot-spots for recombination within this cluster. Thus, different regions of this gene cluster may have been subject to different recombination rates.     

Footprints of intragenic recombination at HLA loci

 To evaluate the effect of balancing selection and intragenic recombination (or gene conversion) at six individual HLA loci, synonymous nucleotide diversity in different exon groups is examined within (π_w) and between (π_b) allelic lineages that may be defined by either serological or DNA sequence differences. Both π values are high in exons which encode for the peptide binding region (PBR) and tend to decrease in other exons. The value of π_w is significantly smaller than that of π_b in any exon of any locus. However, even π_w is much greater than nucleotide diversity at non-HLA loci. These observations provide additional strong evidence for the operation of balancing selection in PBR-encoding exons and its indirect effects on polymorphism at linked neighboring regions. It appears that allelic lineages have generally evolved in isolation but the linkage relationships within and between exons are incomplete throughout the long evolutionary history. To quantify intragenic recombination and account for the large discrepancy between the HLA and non-HLA diversity, a population genetics model is analyzed with special reference to the evolution of modern humans. The analysis suggests that the recombination rate between two sites 1000 base pairs apart is about 10^–5 per generation and that the effective size of human populations (equivalent roughly to the number of breeding individuals in a randomly mating population) has dropped from 10⁵ to 10⁴ in most of the Quaternary. One possibility for this reduction is discussed. Received: 11 August 1997 / Revised: 8 October 1997     

AFLP genetic maps of Eucalyptus globulus and E. tereticornis

 Amplified fragment length polymorphism (AFLP) analysis is a rapid and efficient technique for detecting large numbers of DNA markers in eucalypts. We have used AFLP markers in a two-way pseudo-testcross strategy to generate genetic maps of two clones of different Eucalyptus species (E. tereticornis and E. globulus). Of 606 polymorphic fragments scored, 487 segregated in a 1 : 1 ratio, corresponding to DNA polymorphisms heterozygous in one parent and null in the other. In the maternal E. tereticornis map, 268 markers were ordered in 14 linkage groups (919 cM); the paternal E. globulus map had 200 markers in 16 linkage groups (967 cM). Results from PGRI software were compared with MAPMAKER. The average density of markers was approximately 1 per 3.9 cM. Framework markers were ordered with an average confidence level of 90%, covering 80–100% of the estimated Eucalyptus genome size. In order to investigate the homologies between the E. tereticornis and the E. globulus genetic linkage maps, we included 19 markers segregating 3 : 1 in the analysis. Some homeologous linkage groups were recognized. The linkage data developed in these maps will be used to detect loci controlling commercially important traits. Received: 17 July 1997 / Accepted: 13 October 1997     

陆地棉枯萎病抗性基因的等位性测定及连锁分析

1995-1996年,对我国育成的有代表性的5个抗病品种进行抗枯萎病基因的等位性测定。结果表明:在所选用的5个抗病品种中至少存在两个不同的抗病基因(暂定名为Fwl和Fw2)。连锁分析显示:Fwl与T586的8个标志性状间、Fw2与T582、T586的13个标志性状间无连锁关系。 Abstract:Allelism in vestigation of genes resistante to Fusarium wilt in cotton suggested that there were 2 genes(assigned symbols Fw1 and Fw2)in 5 cultivars used.No linkage was found between Fw1 and the marker genes in T586 and between Fw2 and those marker genes in T582 and T586.