with the goal of discovering the genes that contribute to the risk of common diseases , numerous susceptibility loci have been identified using linkage analysis . despite replication of many of these linkages in a second sample ,
some exceptions include tbc1d1 , identified as the gene responsible for the obesity linkage on chromosome 4p15 - 14 , and hoxb13 , identified as the gene responsible for the prostate cancer linkage on chromosome 17q21 - 22 .
however , for many other common disease linkage signals , the underlying causal genes and variants await discovery .
the expectation that a single gene accounts for a linkage peak may contribute to the difficulty in identifying causal genes . on the contrary
attributed a triglyceride linkage on chromosome 7q36 to variants in five genes and christians et al .
found that fine - mapping caused a single quantitative trait locus ( qtl ) for body size in mice to resolve first into three qtls , then one of those to split into two .
furthermore , the clustering of causal genes may have stymied the multi - group effort to identify type 2 diabetes ( t2d ) genes on chromosome1q ; two strong associations present in populations of european ancestry failed replication and confirmation in other ethnic groups .
these examples suggest that abandoned linkage findings might yet reveal susceptibility genes if reappraised while considering the possibility of multiple causal genes . as for other common diseases ,
the challenge is even greater in african american ( aa ) populations where both prevalence and genetic diversity are higher , and pathophysiology may differ .
genome wide case - control association studies ( gwass ) have identified single nucleotide polymorphisms ( snps ) with small effects on t2d risk , but few of the associated snps , identified primarily in european ancestry populations , are replicated in aas and few snps have been identified in aa populations . in general , the widespread rejection of family studies in favor of gwass has failed to produce the promised prognostic and diagnostic variants for t2d . as a resource for the discovery of genes related to t2d and its complications , the american diabetes association established the genetics of niddm ( gennid ) study . from 1993 to 2003
we have used this resource to increase understanding of the genetics of t2d by applying linkage and family - based association analysis to the african american ( aa ) subset of the gennid sample .
genome - wide linkage analysis using genotypes on 5,914 snps identified chromosomal regions that potentially harbor risk genes for t2d and age of t2d diagnosis ( aod ) .
the strongest signal for t2d occurred on chromosome 2 at 95121 megabases ( mb ) , with weaker support for aod in the same region and for t2d at 6895 mb .
both t2d and aod also showed linkage on chromosome 13 at 1930 mb , but linkage was limited to t2d on chromosome 7 at 5079 mb and to aod on chromosome 18 at 3165 mb .
other analyses inferred pleiotropy with triglyceride in the chromosome 2p region and pleiotropy with obesity in the chromosome 13 region . herein we present linkage and association analyses on the aa subset of the gennid sample using genotypes on 9,203 fine - mapping snps added to the genome scan snps used previously .
our first goal was to more precisely localize the t2d susceptibility genes in each of the five chromosomal regions identified in the genome scan ( chromosome 2 at 6895 and 95121 mb , chromosome 7 at 5079 mb , chromosome 13 at 1930 mb , and chromosome 18 at 3165 mb ) .
our second goal was to test for association of the gene - based fine - mapping snps with t2d and aod in order to identify genes related to t2d in aas .
the gennid study ascertained families through sibling pairs each with a t2d diagnosis . during phase 1 ,
extended family members were also studied ; one site ascertained aas . during phase 2 , data collection beyond the sibling pair was limited to parents , or , if parents were unavailable , unaffected siblings ; five sites ascertained aas . during phase 3 , only affected sibling pairs and trios were studied ; 5 additional sites ascertained aas . in total , 1,496 aa members of 580 pedigrees were studied at 10 sites .
data cleaning and re - evaluation of t2d diagnoses by current criteria reduced the analysis sample to 1,344 members of 524 pedigrees ; for this study , we selected an informative subset of 1,077 members of 415 pedigrees .
t2d diagnoses , originally following national diabetes data group criteria , were re - evaluated using current criteria prior to all analyses ; 84 affected and 13 unaffected sample members were re - assigned as unknown .
body mass index ( bmi ) was computed from height and weight obtained from physical examination .
we selected a subset of the sample for fine - mapping by excluding unaffected individuals below age 53 , the age by which onset had occurred in 75% of cases , and then any individual , either affected or unaffected , who consequently had no family members in the sample .
this sample subset maintained significant lod scores in each of the linkage regions and produced power over 80% to detect association .
power to produce a p - value of 0.00001 was estimated from simulation of 1000 replicates of a snp with minor allele frequency 0.4 and heterozygous effect a 10% increase in penetrance at age 50 for t2d or 3 year increase in aod . for fine - map genotyping , we used snagger ( http://snagger.sourceforge.net/ ) to select tagsnps within genes in the linkage regions on chromosomes 2 , 7 , 13 , and 18 . using hapmap phase 1 and 2 ( http://hapmap.ncbi.nlm.nih.gov ) ,
variants in the region were tagged at linkage disequilibrium ( ld ) r  0.7 .
each snp selected had minor allele frequency > 0.1 and illumina design score > 0.4 ; snp pairs had minimum spacing of 60 base pairs ( bp ) .
the genetic map locations of the snps in centimorgans ( cm ) were obtained from the rutgers combined linkage - physical map of the human genome ( http://compgen.rutgers.edu/maps ) ; locations of snps not mapped directly were estimated using physical positions from ncbi dbsnp build 123 in a smoothing calculation ( conway institute bioinformatics service ; http://integrin.ucd.ie ) .
genotyping errors were identified using pedcheck and merlin ; we zeroed 4978 genotypes of 710 snps , 74 identified by pedcheck and 4904 by merlin .
genotypes for 24 snps , available on less than 95% of the sample following data cleaning , were retained after confirmation of minimal effects on the results .
multi - point identity by descent ( ibd ) probabilities were computed at each cm using merlin , treating as haplotypes snp sets with pairwise ld r > 0.7 .
we identified 279 , 175 , 133 , and 259 snp sets on chromosomes 2 , 7 , 13 , and 18 , respectively .
likelihood analysis , as implemented in jpap , was used for univariate linkage analysis of unadjusted t2d ( ut2d ) , t2d adjusted for bmi ( bt2d ) , and aod and for bivariate association analysis of ut2d and aod .
aod was adjusted for gender and modeled as a normal density with mean , standard deviation , and gender effect as parameters .
t2d risk was modeled to account for aod in affected pedigree members , while allowing for censored observations with age , gender , and bmi ( for bt2d only ) effects and penetrance as parameters . for each trait ,
additional parameters included heritability , a quantitative trait locus ( qtl ) effect in linkage analysis , and an additive snp effect in association analysis .
we applied variance components linkage analysis using the univariate model for each trait in conjunction with the ibd probabilities . at each cm across the 5 regions , all parameters were estimated for aod while only heritability and qtl effect were estimated for ut2d and bt2d with all other parameters for those traits fixed at estimates obtained upon maximizing the likelihood while correcting the likelihood for the ascertainment of each pedigree through an affected sib pair .
the lod score at each cm was computed as the common logarithm of the ratio of the maximized likelihoods with the qtl effect estimated to the maximized likelihood with the qtl effect set to zero .
we assessed each fine - mapping snp for association with ut2d and aod by coding the snp genotype as an additive covariate and testing its effect on both traits in a bivariate model .
p values were determined using a 2 degree of freedom  statistic computed as twice the natural logarithm of the ratio of the maximized likelihood with covariate effects estimated for both t2d and aod to the maximized likelihood with both effects set equal to zero .
we controlled for multiple testing separately within each of the 5 linkage regions by specifying a false discovery rate of 0.05 accounting for 1291 , 1780 , 1819 , 1336 , and 2977 fine - mapping snps for chromosomes 2p , 2q , 7 , 13 , and 18 , respectively .
we tested for an independent effect of each secondary associated snp , conditional on the snp within the same gene that attained the highest significance , by comparing the maximized likelihood estimating t2d and aod covariate effects for both snps simultaneously to the maximized likelihood estimating t2d and aod covariate effects for the most significant snp alone .
p < 0.05 after bonferroni correction for the number of secondary snps in the gene supported the independence of the secondary snp .
we genotyped 9,203 snps within five regions encompassing a total of 133 cm and 128 mb ( table 1 ) and merged them with 244 ( 36 duplicates ) genome scan snps .
spacing between snps averaged 0.0145 cm ( 13,910 bp ) and ranged from 0.0 to 5.3 cm and from 26 bp to 6 mb . to comply with the sample size limitations of the cidr genotyping platform , we selected a subset of the original linkage sample to be informative for both linkage and association .
the genotyped subset included 1077 members of 415 pedigrees that ranged from 2 to 11 members ( 1 to 8 were genotyped ) . in this subset ,
aod and bmi differed little from the complete sample , but the unaffected sample members , selected for a minimum age of 53 years , were older ( table 2 ) .
we previously identified 5 linkage peaks within four chromosomal regions from autosomal scans of t2d and aod ; the chromosome 2 region contained 2 peaks separated by 30 mb and the centromere .
using updated genetic map positions for the genome scan snps , we repeated the linkage analyses on the smaller sample : the genome scan lod scores for all 5 peaks ( table 3 ) generally agreed with our published lod scores for bt2d and aod and for ut2d . upon adding the fine - mapping snps , linkage evidence remained in all 5 regions , but strengthened only on chromosome 2 ( table 3 ) . on chromosome 18 , only aod supported linkage , although ut2d provided weak support ( lod = 2.62 ) upon elimination of sample members with bmi > 45
. the most remarkable effect of fine - mapping was the splitting of the lod score curves into multiple peaks in all of the regions except on chromosome 18 .
additional evidence that multiple susceptibility genes contribute to the linkage signals derived from the identification of associated snps that reside in multiple genes within each of the 5 regions , including on chromosome 18 ( table 4 ) .
the number of associated snps ranged from 4 on chromosome 13 to 17 on chromosome 2q ; the number of associated genes ranged from 2 on chromosome 13 to 9 on chromosome 7 .
although 20 of the 27 associated genes were identified through a single snp , 5 genes were identified through 2 or 3 snps and arhgap25 and dpp10 were identified by 8 and 10 , respectively
. however , six associated snp pairs exceeded ld r  0.7 : arhgap25 ( rs6714065 and rs7605681 ) , ctnna2 ( rs968820 and rs1368915 ) , dpp10 ( rs843417 and rs1823267 , rs10204212 and rs13432035 , rs4848376 and rs11694256 ) , and hip1 ( rs1179625 and rs1179622 ) . nevertheless , conditional association analysis supported the independence of two or more snps in arhgap25 , polr1b , dpp10 , and mtus2 ( supplementary table s1 ) .
the fine - mapping of five chromosomal regions allowed us not only to confirm our genome scan linkages , but also to infer the presence of multiple t2d susceptibility genes in each of the regions .
two types of evidence supported the presence of multiple susceptibility genes underlying the linkage peaks .
first , except chromosome 18 , each lod score peak split into at least two peaks when fine - mapping snps were added to the linkage analysis .
second , in every region including chromosome 18 , two or more potential susceptibility genes were identified through the association with t2d and aod of fine - mapping snps residing in those genes .
in fact , martin et al . attributed many failures to identify causal variants to the incorrect assumption that a single or a limited number of variants are responsible for a linkage signal .
as evidence , they identified variants in five genes that fully accounted for a linkage peak for plasma triglyceride level
. if small effects are typical of common disease susceptibility genes , their detection through linkage analysis may be limited to locations where the genes cluster .
however , fine - mapping using next generation sequencing may provide sufficient power necessary to detect even isolated genes . although chromosome 18 failed to show t2d linkage and chromosome 7 failed to show aod linkage
, we nevertheless expect that the susceptibility genes residing in all regions both increase t2d risk and decrease aod .
variation between regions in the information available undoubtedly affects the relative strength of the corresponding lod scores for t2d or aod . as evidence , the elimination of t2d cases with bmi > 45 kg / m revealed a weak t2d linkage on chromosome 18 and fine - mapping revealed aod linkages on chromosomes 2p and 2q .
the lod scores are generally lower for aod than for t2d as expected since aod reflects onset age imprecisely and the familial correlation of aod may partially result from temporal clustering of t2d diagnoses among relatives . nevertheless , aod proved sufficiently accurate to produce linkage evidence in 4 of the chromosomal regions .
few of the associated genes have previously been reported to associate with t2d or related traits .
the exceptions include grb10 with t2d , 
 nedd4l with diabetic nephropathy , and lipg with lipid metabolism .
novel candidates reported herein include a cytokine ( il36b ) and genes involved in lipid metabolism ( acoxl ) and cell - cell and cell - matrix adhesions ( magi2 , cldn4 , ctnna2 ) .
interestingly , the candidates also included genes involved in williams - beuren syndrome ( wbscr28 , wbscr17 , cldn4 ) , whose sufferers have a high prevalence of diabetes and pre - diabetes .
dpp10 is related to dpp4 , whose inhibitors are oral anti - hyperglycemics used in t2d therapy , and dpp6 , which contributes to a triglyceride linkage , however , unlike dpp4 , ddp10 lacks serine protease / dipeptidyl peptidase activity .
support for the functionality of these genes derived from evidence that their expression levels are affected by our associated snps or proxies with high ld .
we tested 14 snps in the 5 genes for which expression levels from transformed lymphocytes were measured on a subset of 160 sample members ; nominal significance was obtained for one : p = 0.0148 for an effect of rs7579103 on arhgap2 expression levels .
in other populations and various tissues , expression levels of arhgap2 , aff3 , polr1b , hip1 , lipg , and nedd4l showed evidence of an effect of an associated snp or its proxy ( p<1105 ) . for further insight into the nature of the associations , we haplotyped sets of 1530 fine - mapping snps , each encompassing two or more of the associated snps reported in table 4 .
for arhgap25 , two haplotypes associated with t2d : the first extended across 4 of the associated snps and contained the risk alleles for all 4 ; the second extended across 6 associated snps and contained the risk alleles for all 6 , but differed from the first haplotype for alleles at other snps .
we propose that each haplotype harbors a distinct causal variant , despite sharing associations with the same risk alleles , demonstrating the complexity of disease associations .
the advantages include the opportunity to exploit information on local ancestry to localize disease susceptibility genes .
as with any such analysis , confirmation of these findings awaits replication in another sample . in summary , linkage and association analysis using genotypes on an additional 9,203 fine - mapping snps added to five chromosomal regions confirmed each linkage and identified potential susceptibility genes in each of the regions .
 we previously localized type 2 diabetes ( t2d ) susceptibility genes to five chromosomal regions through a genome - wide linkage scan of t2d and age of diagnosis ( aod ) in the african american subset of the gennid sample . 
 to follow up these findings , we repeated the linkage and association analysis using genotypes on an additional 9,203 fine - mapping snps selected to tag genes under the linkage peaks . in each of the five regions , we confirmed linkage and inferred the presence two or more susceptibility genes . the evidence of multiple susceptibility genes comprised : 1 ) multiple linkage peaks in four of the five regions ; 2 ) association of t2d and aod with snps within 2 or more genes in every region . 
 the associated genes included 3 previously reported to associate with t2d or related traits ( grb10 , nedd4l , lipg ) , and 24 novel candidate genes including genes in lipid metabolism ( acoxl ) and cell - cell and cell - matrix adhesion ( magi2 , cldn4 , ctnna2 ) . 
Introduction
Subjects and Methods
Results
Discussion
Supplementary Material