Evaluation of microsatellite markers in association studies: a search for an immune-related susceptibility gene in sarcoidosis

Association studies using linkage disequilibrium (LD) between candidate loci and nearby markers have been proposed to identify susceptibility genes for complex diseases. We analyzed polymorphisms of microsatellites (MSs) and LD patterns of the regions in which candidate genes related to the Th1 immune response have been annotated and attempted to identify a susceptibility gene for sarcoidosis in a marker-based association study. Nineteen MSs were identified in six Th1-related genes (IFNGR1, IFNGR2, IL12RB1, IL12RB2, STAT1 and STAT4) and then eight were further characterized as useful polymorphic markers. Most of these MSs showed LD with single nucleotide polymorphisms (SNPs) on both 5 and 3 ends of these candidate genes, in which r² values between at least one of the MS marker alleles and the SNPs were higher than 0.1. A significant association with one MS allele near STAT4 was shown and a cluster of SNPs in LD with the MS marker was associated with sarcoidosis. These results suggest that association studies using not only SNPs but also multi-allelic MS within or near candidate loci would be useful markers to search for a disease susceptibility gene, especially in populations with unknown LD structure.     

Genetic markers that influence feed efficiency phenotypes also affect cattle temperament as measured by flight speed

Flight speed is a predictive indicator of cattle temperament and is associated with feed efficiency phenotypes. Genetic markers associated with both traits may assist with selection of calmer animals with improved economic value. A preliminary genome‐wide association study determined chromosomal regions on BTA9, and 17 were associated with flight speed. The genes quaking (QKI), glutamate receptor, ionotropic, AMPA 2 (GRIA2) and glycine receptor β (GLRB) were identified in these regions as potential functional candidates. Beef steers (n = 1057) were genotyped with SNPs located within and flanking these genes. One SNP located near QKI and one near GRIA2 were nominally associated with flight speed (P ≤ 0.05) although neither was significant after Bonferroni correction. Several studies have shown a correlation between flight speed and feed intake or gain; therefore, we also analyzed SNPs on BTA6:38–39 Mb known to be associated with average daily gain (ADG) and average daily feed intake (ADFI) for association with flight speed. Several SNPs on BTA6 were associated with flight speed (P ≤ 0.005), and three were significant after Bonferroni correction. These results suggest that the genes tested are unlikely to contribute to flight speed variation for our cattle population, but SNPs on BTA6 associated with ADG and ADFI may influence temperament. Use of these markers to select for economically important feed efficiency phenotypes may produce cattle with more desirable temperaments.     

An association analysis for 14 candidate genes mapping to meat quality quantitative trait loci in a Duroc pig population reveals that the ATP1A2 genotype is highly associated with muscle electric conductivity

In previous GWAS carried out in a Duroc commercial line (Lipgen population), we detected on pig chromosomes 3, 4 and 14 several QTL for gluteus medius muscle redness (GM a*), electric conductivity in the longissimus dorsi muscle (LD CE) and vaccenic acid content in the LD muscle (LD C18:1 n − 7), respectively. We have genotyped, in the Lipgen population, 19 SNPs mapping to 14 genes located within these QTL. Subsequently, association analyses have been performed. After correction for multiple testing, two SNPs in the TGFBRAP1 (rs321173745) and SELENOI (rs330820437) genes were associated with GM a*, whereas ACADSB (rs81449951) and GPR26 (rs343087568) genotypes displayed significant associations with LD vaccenic content. Moreover, the polymorphisms located at the ATP1A2 (rs344748241), ATP8B2 (rs81382410) and CREB3L4 (rs321278469 and rs330133789) genes showed significant associations with LD CE. We made a second round of association analyses including the SNPs mentioned above as well as other SNPs located in the chromosomes to which they map. After performing a correction for multiple testing, the only association that remained significant at the chromosome-wide level was that between the ATP1A2 genotype and LD CE. From a functional point of view, this association is meaningful because this locus encodes a subunit of the Na⁺/K⁺-ATPase responsible for maintaining an electrochemical gradient across the plasma membrane.     

Genome-wide association analysis with selective genotyping identifies candidate loci for adult height at 8q21.13 and 15q22.33-q23 in Mongolians

We performed a genome-wide association study with 23,465 microsatellite markers to identify genes related to adult height. Selective genotyping was applied to extremely tall and extremely short individuals from the Khalkh-Mongolian population. Two loci, 8q21.13 and 15q22.33, which showed the strongest association with microsatellites were subjected to further analyses of SNPs in 782 tall and 773 short individuals. The most significant association was observed with SNP rs2220456 at 8q21.13 (P = 0.000016). In the LD block at 15q22.32, SNP rs8038652 located in intron 1 of IQCH was strongly associated (P = 0.0003), especially the AA genotype of the SNP under a recessive model was strongly associated with adult height (P = 0.000046).     

New plant binary vectors with selectable markers located proximal to the left T-DNA border

Five new binary vectors have been constructed which have the following features: (1) different plant selectable markers including neomycin phosphotransferase (nptII), hygromycin phosphotransferase (hpt), dihydrofolate reductase (dhfr), phosphinothricin acetyl transferase (bar), and bleomycin resistance (ble); (2) selectable markers are located near the T-DNA left border and; (3) selectable marker and -glucuronidase (uidA) reporter genes are divergently organized for efficient expression, and can easily be removed or replaced as needed.     

Molecular mapping of novel resistance genes against Barley Mild Mosaic Virus (BaMMV)

 In the present study three novel genes from barley accessions 10247 (ym8), Bulgarian 347 (ym9), and Russia 57 (ym11), which confer resistance to Barley Mild Mosaic Virus (BaMMV), were mapped using molecular markers. Bulked segregant analysis of four progenies segregating for resistance to BaMMV was followed by fine-scale mapping of the resistance genes using individual F₂ or BC₁F₂ plants. The resistance genes are inherited recessively and are located on the long arm of barley chromosome 4HL. A series of closely linked molecular markers are available for marker-assisted breeding programs. A marker (MWG2134) linked with resistance gene ym11 from Russia 57 was identified, which is diagnostic for the resistance gene. Received: 25 July 1997 / Accepted: 22 August 1997     

Genome-wide association analysis of meat quality traits in a porcine Large White × Minzhu intercross population

Pork quality is an economically important trait and one of the main selection criteria for breeding in the swine industry. In this genome-wide association study (GWAS), 455 pigs from a porcine Large White × Minzhu intercross population were genotyped using the Illumina PorcineSNP60K Beadchip, and phenotyped for intramuscular fat content (IMF), marbling, moisture, color L*, color a*, color b* and color score in the longissimus muscle (LM). Association tests between each trait and the SNPs were performed via the Genome Wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) approach. From the Ensembl porcine database, SNP annotation was implemented using Sus scrofa Build 9. A total of 45 SNPs showed significant association with one or multiple meat quality traits. Of the 45 SNPs, 36 were located on SSC12. These significantly associated SNPs aligned to or were in close approximation to previously reported quantitative trait loci (QTL) and some were located within introns of previously reported candidate genes. Two haplotype blocks ASGA0100525-ASGA0055225-ALGA0067099-MARC0004712-DIAS0000861, and ASGA0085522-H3GA0056170 were detected in the significant region. The first block contained the genes MYH1, MYH2 and MYH4. A SNP (ASGA0094812) within an intron of the USP43 gene was significantly associated with five meat quality traits. The present results effectively narrowed down the associated regions compared to previous QTL studies and revealed haplotypes and candidate genes on SSC12 for meat quality traits in pigs.     

Analysis of dyslexia candidate genes in the Raine cohort representing the general Australian population

Several genes have been suggested as dyslexia candidates. Some of these candidate genes have been recently shown to be associated with literacy measures in sample cohorts derived from the general population. Here, we have conducted an association study in a novel sample derived from the Australian population (the Raine cohort) to further investigate the role of dyslexia candidate genes. We analysed markers, previously reported to be associated with dyslexia, located within the MRPL19/C2ORF3, KIAA0319, DCDC2 and DYX1C1 genes in a sample of 520 individuals and tested them for association with reading and spelling measures. Association signals were detected for several single nucleotide polymorphisms (SNPs) within DYX1C1 with both the reading and spelling tests. The high linkage disequilibrium (LD) we observed across the DYX1C1 gene suggests that the association signal might not be refined by further genetic mapping.     

Loci of the immune system are implicated in diabetes frequency and age at onset of diabetes in BB rats

The spontaneously diabetic BB rat is a well-established animal model in diabetes research developing an insulin-dependent type-1 diabetes mellitus closely resembling human diabetes. By several crossing studies using BB/OK rats it has been demonstrated that beside the MHC class-II genes of the RT1^u haplotype, Iddm1, and the lymphopenia, Iddm2, at least two additional non-MHC genes located on chromosomes 6 (Iddm4) and 18 (Iddm3) are involved in diabetes development. In addition, there are at least three genes located on chromosomes 6 (Dm1), 8 (Dm2) and 10 (Dm3) influencing the age at onset of diabetes. Comparing the homologous regions between rat and human, it is shown that most diabetogenic genes lie on human chromosomes near genes involved in immune processes providing human geneticist with new candidate regions for the analysis of diabetogenic non-MHC genes in human type-1 diabetes.     

A genome‐wide association study of heat stress‐associated SNPs in catfish

Heat tolerance is a complex and economically important trait for catfish genetic breeding programs. With global climate change, it is becoming an increasingly important trait. To better understand the molecular basis of heat stress, a genome‐wide association study (GWAS) was carried out using the 250 K catfish SNP array with interspecific backcross progenies, which derived from crossing female channel catfish with male F1 hybrid catfish (female channel catfish × male blue catfish). Three significant associated SNPs were detected by performing an EMMAX approach for GWAS. The SNP located on linkage group 14 explained 12.1% of phenotypical variation. The other two SNPs, located on linkage group 16, explained 11.3 and 11.5% of phenotypical variation respectively. A total of 14 genes with heat stress related functions were detected within the significant associated regions. Among them, five genes—TRAF2, FBXW5, ANAPC2, UBR1 and KLHL29— have known functions in the protein degradation process through the ubiquitination pathway. Other genes related to heat stress include genes involved in protein biosynthesis (PRPF4 and SYNCRIP), protein folding (DNAJC25), molecule and iron transport (SLC25A46 and CLIC5), cytoskeletal reorganization (COL12A1) and energy metabolism (COX7A2, PLCB1 and PLCB4) processes. The results provide fundamental information about genes and pathways that is useful for further investigation into the molecular mechanisms of heat stress. The associated SNPs could be promising candidates for selecting heat‐tolerant catfish lines after validating their effects on larger and various catfish populations.     

A genome-wide association study of seed protein and oil content in soybean

Background

Association analysis is an alternative to conventional family-based methods to detect the location of gene(s) or quantitative trait loci (QTL) and provides relatively high resolution in terms of defining the genome position of a gene or QTL. Seed protein and oil concentration are quantitative traits which are determined by the interaction among many genes with small to moderate genetic effects and their interaction with the environment. In this study, a genome-wide association study (GWAS) was performed to identify quantitative trait loci (QTL) controlling seed protein and oil concentration in 298 soybean germplasm accessions exhibiting a wide range of seed protein and oil content.

Results

A total of 55,159 single nucleotide polymorphisms (SNPs) were genotyped using various methods including Illumina Infinium and GoldenGate assays and 31,954 markers with minor allele frequency >0.10 were used to estimate linkage disequilibrium (LD) in heterochromatic and euchromatic regions. In euchromatic regions, the mean LD (r ² ) rapidly declined to 0.2 within 360 Kbp, whereas the mean LD declined to 0.2 at 9,600 Kbp in heterochromatic regions. The GWAS results identified 40 SNPs in 17 different genomic regions significantly associated with seed protein. Of these, the five SNPs with the highest associations and seven adjacent SNPs were located in the 27.6-30.0 Mbp region of Gm20. A major seed protein QTL has been previously mapped to the same location and potential candidate genes have recently been identified in this region. The GWAS results also detected 25 SNPs in 13 different genomic regions associated with seed oil. Of these markers, seven SNPs had a significant association with both protein and oil.

Conclusions

This research indicated that GWAS not only identified most of the previously reported QTL controlling seed protein and oil, but also resulted in narrower genomic regions than the regions reported as containing these QTL. The narrower GWAS-defined genome regions will allow more precise marker-assisted allele selection and will expedite positional cloning of the causal gene(s).     

A Scale-Corrected Comparison of Linkage Disequilibrium Levels between Genic and Non-Genic Regions

The understanding of non-random association between loci, termed linkage disequilibrium (LD), plays a central role in genomic research. Since causal mutations are generally not included in genomic marker data, LD between those and available markers is essential for capturing the effects of causal loci on localizing genes responsible for traits. Thus, the interpretation of association studies requires a detailed knowledge of LD patterns. It is well known that most LD measures depend on minor allele frequencies (MAF) of the considered loci and the magnitude of LD is influenced by the physical distances between loci. In the present study, a procedure to compare the LD structure between genomic regions comprising several markers each is suggested. The approach accounts for different scaling factors, namely the distribution of MAF, the distribution of pair-wise differences in MAF, and the physical extent of compared regions, reflected by the distribution of pair-wise physical distances. In the first step, genomic regions are matched based on similarity in these scaling factors. In the second step, chromosome- and genome-wide significance tests for differences in medians of LD measures in each pair are performed. The proposed framework was applied to test the hypothesis that the average LD is different in genic and non-genic regions. This was tested with a genome-wide approach with data sets for humans (Homo sapiens), a highly selected chicken line (Gallus gallus domesticus) and the model plant Arabidopsis thaliana. In all three data sets we found a significantly higher level of LD in genic regions compared to non-genic regions. About 31% more LD was detected genome-wide in genic compared to non-genic regions in Arabidopsis thaliana, followed by 13.6% in human and 6% chicken. Chromosome-wide comparison discovered significant differences on all 5 chromosomes in Arabidopsis thaliana and on one third of the human and of the chicken chromosomes.     

A genome‐wide association study for the incidence of persistent bovine viral diarrhea virus infection in cattle

Bovine viral diarrhea viruses (BVDV) comprise a diverse group of viruses that cause disease in cattle. BVDV may establish both transient and persistent infections depending on the developmental stage of the animal at exposure. The objective was to determine whether genomic regions harboring single nucleotide polymorphisms (SNPs) could be associated with the presence or absence of persistent BVDV infection. A genome‐wide association approach based on 777 000 SNP markers was used. Samples of animals identified as positive (n = 1200) or negative (n = 1200) for the presence of BVDV in skin samples (n = 1200) were used. DNA samples were combined in 24 pools (100 animals per pool). One SNP, significant at the 5 percent genome‐wide level (P = 9.41 × 10^?8), was detected on chromosome 14, located at position 80 675 176 bp. Fifteen SNPs, residing on chromosomes 1, 2, 6, 8, 10, 15 and 18, were moderately associated (P < 1 × 10^?5) with persistent BVDV infection. Results show that genes harboring or neighboring significant SNPs are involved in leucopenia, signal transduction, RNA splicing and DNA methylation processes.     

Soybean aphid intrabiotype variability based on colonization of specific soybean genotypes

The soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae), is one of the most destructive insect pests on soybeans in the United States. One method for managing this pest is through host plant resistance. Since its arrival in 2000, 4 aphid biotypes have been identified that are able to overcome soybean aphid resistance (Rag) genes. A soybean aphid isolate collected from Moline, Illinois readily colonized soybean plants with the soybean aphid resistance gene Rag2, unlike biotypes 1 and 2, but similar to soybean aphid biotype 3. Two no‐choice experiments compared the virulence of the Moline isolate with biotype 3. In both experiments, differences in aphid population counts were not significant (P > 0.05) on soybean genotypes LD08–12957a (Rag2) and LD11–5413a (Rag2), but the aphid counts for the Moline isolate were significantly (P < 0.05) lower than the aphid counts for the biotype 3 isolate on the soybean genotypes Dowling (Rag1), LD05–16611 (Rag1), LD11–4576a (Rag1), and PI 567598B (rag1b and rag3). The Moline isolate was a variant of aphid biotype 3, which is the first report showing that soybean aphid isolates classified as the same biotype, based on virulence against specific Rag genes, can differ in aggressiveness or ability to colonize specific host genotypes.     

Dissection of the genetic architecture of rice resistance to Xanthomonas oryzae pv. oryzae using a genomewide association study

Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), usually causes serious rice yield loss in many countries. Rice breeders have used resistance (R) genes to control the disease but many of the resistant cultivars become susceptible few years after releasing. Identification of new R genes to Xoo is one of the main objectives in rice breeding programs. In this study, we used a genomewide association study (GWAS) to analyse the resistance against the Xoo race C1 using the Rice Diversity Panel 1 (RDP1). Disease evaluation of the RDP1 population to C1 indicated that the AUS subgroup conferred a higher level of resistance to C1 than other subgroups. Genomewide association mapping identified 15 QTLs that are distributed on chromosomes 1, 2, 3, 4, 5, 6, 8, 9, 10 and 12. Some of them are located in the regions without known resistance loci or QTLs. This study demonstrated the effectiveness of GWAS on the genetic dissection of rice resistance to Xoo and provided many Xoo resistance‐associated SNP markers for rice breeding.     

Genome-wide identification and analysis of sulfatase and sulfatase modifying factor genes in Bemisia tabaci (Hemiptera: Aleyrodidae)

The invasive pest whitefly (Bemisia tabaci) is a complex species, of which Middle East-Minor Asia 1 (MEAM1) and Mediterranean (MED) are the two most damaging members. Previous research showed that cabbage is frequently infested with MEAM1 but seldomly with MED, and this difference in performance is associated with glucosinolate (GS) content. Some insects can modify GS using glucosinolate sulfatase (SULF), the activity of which is regulated by sulfatase modifying factor 1 (SUMF1); therefore, to increase our understanding of different performances of MEAM1 and MED on cabbage plants, we identified and compared nine putative SULFs and one SUMF in MEAM1 and MED. We found that the lengths of two genes, BtSulf2 and BtSulf4, differed between MEAM1 and MED. The messenger RNA levels of BtSulf4 increased more than 20-fold after MEAM1 and MED adults were exposed to GS, but BtSulf2 expression was only induced by GS in MEAM1. Knockdown of BtSulf2 and BtSulf4 in MEAM1 resulted in a substantial increase in the mortality of GS-treated adults but not in MED. These results indicate that differences in BtSulf2 and BtSulf4 sequences and/or expression may explain why MEAM1 performs better than MED on cabbage. Our results provide a basis for future functional research on SULF and SUMF in B. tabaci.     

Linkage Disequilibrium in Salt Tolerant Genotypes of Rice (Oryza sativa L)

Genome wide linkage disequilibrium (LD) was investigated in a set of 32 genotypes representing salt tolerant improved varieties and landraces and six salt sensitive genotypes of rice with 64 microsatellite markers to identify the genomic regions that are associated with salt tolerance in rice. Out of 64 markers analyzed, 36% SSR pairs exhibited significant LD at 0.05. A few regions were identified as targets of selection in 10 chromosomes with high r ² values. The model-based groups from Bayesian clustering analysis are largely consistent with known pedigrees of the lines. The increased percentage of association of SSR loci in the improved varieties indicated the role of selection in linkage disequilibrium especially for salt tolerance. LD was extended as far as 100 cM in the present study. Most of the markers (43.8%) with significant LD values were observed in the genomic regions of reported QTL for salt tolerance in rice.     

QTL mapping ofFusarium moniliforme ear rot resistance in maize. 1. Map construction with microsatellite and AFLP markers

To map the QTLsof Fusarium moniliforme ear rot resistance inZea mays L., a total of 230 F₂ individuals, derived from a single cross between inbred maize lines R15 (resistant) and Ye478 (susceptible), were genotyped for genetic map construction using simple sequence repeat (SSR) markers and amplified fragment length polymorphism (AFLP) markers. We used 778 pairs of SSR primers and 63 combinations of AFLP primers to detect the polymorphisms between parents, R15 and Ye478. From the polymorphic 30 AFLP primer combinations and 159 SSR primers, we scored 260 loci in the F₂ population, among which 8 SSR and 13 AFLP loci could not be assigned to any of the linkage groups. An integrated molecular genetic linkage map was constructed by the remaining 151 SSR and 88 AFLP markers, which distributed throughout the 10 linkage groups of maize and spanned the genome of about 3463.5 cM with an average of 14.5 cM between two markers. On 4 chromosomes, we detected 5 putative segregation distortion regions (SDR_s), including 2 new ones (SDR₂ and SDR₇). The other 3 SDR_s were located near the regions where gametophyte genes were mapped, indicating that segregation distortion could be partially caused by gametophytic factors.