Genome-wide association studies: utility and limitations for research in physiology.

Physiological systems are subject to interindividual variation encoded by genetics. Genome-wide association studies (GWAS) operate by surveying thousands of genetic variants from a substantial number of individuals and assessing their association to a trait of interest, be it a physiological variable, a molecular phenotype (e.g. gene expression), or even a disease or condition. Through a myriad of methods, GWAS downstream analyses then explore the functional consequences of each variant and attempt to ascertain a causal relationship to the phenotype of interest, as well as to delve into its links to other traits. This type of investigation allows mechanistic insights into physiological functions, pathological disturbances and shared biological processes between traits (i.e. pleiotropy). An exciting example is the discovery of a new thyroid hormone transporter (SLC17A4) and hormone metabolising enzyme (AADAT) from a GWAS on free thyroxine levels. Therefore, GWAS have substantially contributed with insights into physiology and have been shown to be useful in unveiling the genetic control underlying complex traits and pathological conditions; they will continue to do so with global collaborations and advances in genotyping technology. Finally, the increasing number of trans-ancestry GWAS and initiatives to include ancestry diversity in genomics will boost the power for discoveries, making them also applicable to non-European populations.

Trans-ancestral fine-mapping of MHC reveals key amino acids associated with spontaneous clearance of hepatitis C in HLA-DQß1.

Spontaneous clearance of acute hepatitis C virus (HCV) infection is associated with single nucleotide polymorphisms (SNPs) on the MHC class II. We fine-mapped the MHC region in European (n = 1,600; 594 HCV clearance/1,006 HCV persistence) and African (n = 1,869; 340 HCV clearance/1,529 HCV persistence) ancestry individuals and evaluated HCV peptide binding affinity of classical alleles. In both populations, HLA-DQß1Leu26 (p valueMeta = 1.24 × 10-14) located in pocket 4 was negatively associated with HCV spontaneous clearance and HLA-DQß1Pro55 (p valueMeta = 8.23 × 10-11) located in the peptide binding region was positively associated, independently of HLA-DQß1Leu26. These two amino acids are not in linkage disequilibrium (r2 < 0.1) and explain the SNPs and classical allele associations represented by rs2647011, rs9274711, HLA-DQB1∗03:01, and HLA-DRB1∗01:01. Additionally, HCV persistence classical alleles tagged by HLA-DQß1Leu26 had fewer HCV binding epitopes and lower predicted binding affinities compared to clearance alleles (geometric mean of combined IC50 nM of persistence versus clearance; 2,321 nM versus 761.7 nM, p value = 1.35 × 10-38). In summary, MHC class II fine-mapping revealed key amino acids in HLA-DQß1 explaining allelic and SNP associations with HCV outcomes. This mechanistic advance in understanding of natural recovery and immunogenetics of HCV might set the stage for much needed enhancement and design of vaccine to promote spontaneous clearance of HCV infection.

Characterization and Use in Wheat Breeding of Leaf Rust Resistance Genes from Durable Varieties.

Leaf rust is one of the most significant diseases of wheat worldwide. In Argentina, it is one of the main reasons for variety replacement that becomes susceptible after large-scale use. Some varieties showed durable resistance to this disease, including Buck Manantial and Sinvalocho MA. RILs (Recombinant Inbred Lines) were developed for each of these varieties and used in genetics studies to identify components of resistance, both in greenhouse inoculations using leaf rust races, and in field evaluations under natural population infections. In Buck Manantial, the APR gene LrBMP1 was associated with resistance in field tests. In crosses involving Sinvalocho MA, four genes were previously identified and associated with resistance in field testing: APR (Adult Plant Resistance) gene LrSV1, the APR genetic system LrSV2 + LrcSV2 and the ASR (All Stage Resistance) gene LrG6. Using backcrosses, LrBMP1 was introgressed in four commercial susceptible varieties and LrSV1, LrSV2 + LrcSV2 and LrG6 were simultaneously introgressed in three susceptible commercial varieties. The use of molecular markers for recurrent parent background selection allowed us to select resistant lines with more than 80% similarity to commercial varieties. Additionally, progress towards positional cloning of the genetic system LrSV2 + LrcSV2 for leaf rust APR is reported.

A non-coding RNASEH1 gene variant associates with type 1 diabetes and interacts with HLA tagSNPs in families from Colombia.

OBJECTIVES: RNASEH1 gene has recently been associated with type 1 diabetes (T1D) in Colombia. The purpose of this study was to fine mapping the putative functional variant in RNASEH1 and testing its interaction with HLA tagSNPs. METHODS: Two-hundred nuclear families with T1D were included in this study. Probands were tested for GAD65 and IA-2 autoantibodies. Genotyping was performed using 20 coding tagSNPs uncovered through Sanger sequencing (N = 96), in addition to 23 tagSNPs chosen from 1000genomes to cover the extent of the gene region. Also, 45 tagSNPs for classic HLA alleles associated with T1D were also genotyped. The transmission disequilibrium test (TDT) was used to test for association and a multiple testing correction was made using permutation. Interaction between RNASEH1 variants and HLA was evaluated by means of the M-TDT test. RESULTS: We identified 20 variants (15 were novel) in the 96 patients sequenced. None of these variants were in linkage disequilibrium. In total, 43 RNASEH1 variants were genotyped in the 200 families. Association between T1D and rs7607888 was identified (P = .002). Haplotype analysis involving rs7607888 variant revealed even stronger association with T1D (most significative P = .0003). HLA tagSNPs displayed stronger associations (OR = 6.39, 95% CI = 4.33-9.44, P-value = 9.74E-28). Finally, we found several statistically significant interactions of HLA variants with rs7607888 (P-value ranged from 8.77E-04 to 5.33E-12). CONCLUSION: Our results verify the association of rs7607888 in RNASEH1 gene with T1D. It is also shown in the interaction between RNASEH1 and HLA for conveying risk to T1D in Northwest Colombia. Work is underway aiming to identify the actual classic HLA alleles associated with the tagSNPs tested here.

Fine mapping under linkage peaks for symptomatic or asymptomatic outcomes of Leishmania infantum infection in Brazil.

Infection with the protozoan Leishmania infantum can lead to asymptomatic infection and protective immunity, or to the progressive and potentially fatal disease visceral leishmaniasis (VL). Published studies show host genetic background determines in part whether infected individuals will develop a symptomatic or asymptomatic outcome. The purpose of the current study was to fine map chromosome regions previously linked with risk for symptomatic (chromosome 9) or asymptomatic (chromosomes 15 and 19) manifestations of L. infantum infection. We conducted a family-based genetic study of VL and asymptomatic infection (detected by a DTH skin test) with a final post quality control sample of 961 individuals with full genotype and phenotype information from highly endemic neighborhoods of northeast Brazil. A total of 5485 SNPs under the linkage peaks on chromosomes 9, 15 and 19 were genotyped. No strong SNP associations were observed for the DTH phenotype. The most significant associations with the VL phenotype were with SNP rs1470217 (p=5.9e-05; pcorrected=0.057) on chromosome 9, and with SNP rs8107014 (p=1.4e-05; pcorrected=0.013) on chromosome 19. SNP rs1470217 is situated in a 180kb intergenic region between TMEM215 (Transmembrane protein 215) and APTX (Aprataxin). SNP rs8107014 lies in the intron between exons 26 and 27 of a 34 exon transcript (ENST00000204005) of LTBP4, (Latent transforming growth factor-beta-binding protein 4a). The latter supports growing evidence that the transforming growth factor-beta pathway is important in the immunopathogenesis of VL.

Genetic analysis and gene fine mapping of aroma in rice (Oryza sativa L. Cyperales, Poaceae)

We investigated inheritance and carried out gene fine mapping of aroma in crosses between the aromatic elite hybrid rice Oryza sativa indica variety Chuanxiang-29B (Ch-29B) and the non-aromatic rice O. sativa indica variety R2 and O. sativa japonica Lemont (Le). The F1 grains and leaves were non-aromatic while the F2 non-aroma to aroma segregation pattern was 3:1. The F3 segregation ratio was consistent with the expected 1:2:1 for a single recessive aroma gene in Ch-29B. Linkage analysis between simple sequence repeat (SSR) markers and the aroma locus for the aromatic F2 plants mapped the Ch-29B aroma gene to a chromosome 8 region flanked by SSR markers RM23120 at 0.52 cM and RM3459 at 1.23 cM, a replicate F2 population confirming these results. Three bacterial artificial chromosome (BAC) clones cover chromosome 8 markers RM23120 and RM3459. Our molecular mapping data from the two populations indicated that the aroma locus occurs in a 142.85 kb interval on BAC clones AP005301 or AP005537, implying that it might be the same gene reported by Bradbury et al (2005a; Plant Biotec J. 3:363-370). The flanking markers Aro7, RM23120 and RM3459 identified by us could greatly accelerate the efficiency and precision of aromatic rice breeding programs.