Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

Recent studies identified MYH9 as a major susceptibility gene for common forms of non-diabetic end-stage kidney disease (ESKD). A set of African ancestry DNA sequence variants comprising the E-1 haplotype, was significantly associated with ESKD. In order to determine whether African ancestry variants are also associated with disease susceptibility in admixed populations with differing genomic backgrounds, we genotyped a total of 1425 African and Hispanic American subjects comprising dialysis patients with diabetic and non-diabetic ESKD and controls, using 42 single nucleotide polymorphisms (SNPs) within the MYH9 gene and 40 genome-wide and 38 chromosome 22 ancestry informative markers. Following ancestry correction, logistic regression demonstrated that three of the E-1 SNPs are also associated with non-diabetic ESKD in the new sample sets of both African and Hispanic Americans, with a stronger association in Hispanic Americans. We also identified MYH9 SNPs that are even more powerfully associated with the disease phenotype than the E-1 SNPs. These newly associated SNPs, could be divided into those comprising a haplotype termed S-1 whose association was significant under a recessive or additive inheritance mode (rs5750248, OR 4.21, P < 0.01, Hispanic Americans, recessive), and those comprising a haplotype termed F-1 whose association was significant under a dominant or additive inheritance mode (rs11912763, OR 4.59, P < 0.01, Hispanic Americans, dominant). These findings strengthen the contention that a sequence variant of MYH9, common in populations with varying degrees of African ancestry admixture, and in strong linkage disequilibrium with the associated SNPs and haplotypes reported herein, strongly predisposes to non-diabetic ESKD.

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aimed to identify the loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC) and hip circumference (HIP) adjusted for body mass index (adjBMI) by using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); determine if differences in associations differ by background group within HCHS/SOL and determine whether previously reported associations generalize to HCHS/SOL. Our analyses included 7472 women and 5200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban and Dominican) background residing in the USA. We performed genome-wide association analyses stratified and combined across sexes using linear mixed-model regression. We identified 16 variants for waist-to-hip ratio adjusted for body mass index (WHRadjBMI), 22 for waist circumference adjusted for body mass index (WCadjBMI) and 28 for hip circumference adjusted for body mass index (HIPadjBMI), which reached suggestive significance (P < 1 × 10-6). Many loci exhibited differences in strength of associations by ethnic background and sex. We brought a total of 66 variants forward for validation in cohorts (N = 34 161) with participants of Hispanic/Latino, African and European descent. We confirmed four novel loci (P < 0.05 and consistent direction of effect, and P < 5 × 10-8 after meta-analysis), including two for WHRadjBMI (rs13301996, rs79478137); one for WCadjBMI (rs3168072) and one for HIPadjBMI (rs28692724). Also, we generalized previously reported associations to HCHS/SOL, (8 for WHRadjBMI, 10 for WCadjBMI and 12 for HIPadjBMI). Our study highlights the importance of large-scale genomic studies in ancestrally diverse Hispanic/Latino populations for identifying and characterizing central obesity susceptibility that may be ancestry-specific.

Leukocyte telomere length (LTL) might be causal in cardiovascular disease and major cancers. To elucidate the roles of genetics and geography in LTL variability across humans, we compared LTL measured in 1295 sub-Saharan Africans (SSAs) with 559 African-Americans (AAms) and 2464 European-Americans (EAms). LTL differed significantly across SSAs (P = 0.003), with the San from Botswana (with the oldest genomic ancestry) having the longest LTL and populations from Ethiopia having the shortest LTL. SSAs had significantly longer LTL than AAms [P = 6.5(e-16)] whose LTL was significantly longer than EAms [P = 2.5(e-7)]. Genetic variation in SSAs explained 52% of LTL variance versus 27% in AAms and 34% in EAms. Adjustment for genetic variation removed the LTL differences among SSAs. LTL genetic variation among SSAs, with the longest LTL in the San, supports the hypothesis that longer LTL was ancestral in humans. Identifying factors driving LTL variation in Africa may have important ramifications for LTL-associated diseases.

Women of African ancestry have the highest mortality from triple-negative breast cancer (TNBC) of all racial groups. To understand the genomic basis of breast cancer in the populations, we previously conducted genome-wide association studies and identified single nucleotide polymorphisms (SNPs) associated with breast cancer in Black women. In this study, we investigated the functional significance of the top associated SNP rs13074711. We found the SNP served as an enhancer variant and regulated TNFSF10 (TRAIL) expression in TNBC cells, with a significant association between the SNP genotype and TNFSF10 expression in breast tumors. Mechanistically, rs13074711 modulated the binding activity of c-MYB at the motif and thereby controlled TNFSF10 expression. Interestingly, TNFSF10 expression in many cancers was consistently lower in African Americans compared with European Americans. Furthermore, TNFSF10 expression in TNBC was significantly correlated with the expression of antiviral immune genes and was regulated by type I interferons (IFNs). Accordingly, loss of TNFSF10 resulted in a profound decrease in apoptosis of TNBC cells in response to type I IFNs and poly(I:C), a synthetic analogue of double stranded virus. Lastly, in a syngeneic mouse model of breast cancer, TNFSF10-deficiency in breast tumors decreased tumor-infiltrated CD4+ and CD8+ T cell quantities. Collectively, our results suggested that TNFSF10 plays an important role in the regulation of antiviral immune responses in TNBC, and the expression is in part regulated by a genetic variant associated with breast cancer in Black women. Our results underscore the important contributions of genetic variants to immune defense mechanisms.

Estrogen modifies human emotion and cognition and impacts symptoms of schizophrenia. We hypothesized that the variation in the estrogen receptor alpha (ESR1) gene and cortical ESR1 mRNA is associated with schizophrenia. In a small case-control genetic association analysis of postmortem brain tissue, genotype CC (rs2234693) and haplotypes containing the C allele of a single-nucleotide polymorphism (SNP) in intron1 (PvuII) were more frequent in African American schizophrenics (P = 0.01-0.001). In a follow-up family-based association analysis, we found overtransmission of PvuII allele C and a PvuII C-containing haplotype (P = 0.01-0.03) to African American and Caucasian patients with schizophrenia. Schizophrenics with the 'at risk' PvuII genotype had lower ESR1 mRNA levels in the frontal cortex. Eighteen ESR1 splice variants and decreased frequencies of the wild-type ESR1 mRNA were detected in schizophrenia. In one patient, a unique ESR1 transcript with a genomic insert encoding a premature stop codon and a truncated ESR1 protein lacking most of the estrogen binding domain was the only transcript detected. Using a luciferase assay, we found that mRNA encoding a truncated ESR1 significantly attenuates gene expression at estrogen-response elements demonstrating a dominant negative function. An intron 6 SNP [rs2273207(G)] was associated with an ESR1 splice variant missing exon seven. The T allele of another intron 6 SNP was part of a 3' haplotype less common in schizophrenia [rs2273206(T), rs2273207(G), rs2228480(G)]. Thus, the variation in the ESR1 gene is associated with schizophrenia and the mechanism of this association may involve alternative gene regulation and transcript processing.

The children of related parents show increased risk of early mortality. The Native American genome typically exhibits long stretches of homozygosity, and Latin Americans are highly heterogeneous regarding the individual burden of homozygosity, the proportion and the type of Native American ancestry. We analysed nationwide mortality and genome-wide genotype data from admixed Chileans to investigate the relationship between common causes of child mortality, homozygosity and Native American ancestry. Results from two-stage linear-Poisson regression revealed a strong association between the sum length of runs of homozygosity (SROH) above 1.5 Megabases (Mb) in each genome and mortality due to intracranial non-traumatic haemorrhage of foetus and newborn (5% increased risk of death per Mb in SROH, P = 1 × 10-3) and disorders related to short gestation and low birth weight (P = 3 × 10-4). The major indigenous populations in Chile are Aymara-Quechua in the north of the country and the Mapuche-Huilliche in the south. The individual proportion of Aymara-Quechua ancestry was associated with an increased risk of death due to anencephaly and similar malformations (P = 4 × 10-5), and the risk of death due to Edwards and Patau trisomy syndromes decreased 4% per 1% Aymara-Quechua ancestry proportion (P = 4 × 10-4) and 5% per 1% Mapuche-Huilliche ancestry proportion (P = 2 × 10-3). The present results suggest that short gestation, low birth weight and intracranial non-traumatic haemorrhage mediate the negative effect of inbreeding on human selection. Independent validation of the identified associations between common causes of child death, homozygosity and fine-scale ancestry proportions may inform paediatric medicine.

ITGAM was recently found to be associated with systemic lupus erythematosus (SLE) in populations of not only European ancestry, but also in Hispanic- and African-Americans, Mexicans and Colombians. The risk alleles in the gene, however, were found to be monomorphic in two Asian populations examined: Japanese and Korean. In this study, using a collection of 910 SLE patients and 2360 controls from Chinese living in Hong Kong, analyzed by both genome-wide association and direct sequencing, we confirmed the association of the same risk alleles in ITGAM with the disease. These findings were further replicated in the Thai population with 278 patients and 383 ethnicity- and geography-matched controls. Subphenotype stratification analyses showed significantly more involvement of the gene in patients with renal nephritis and neurological disorders. Although our results support a pivotal role by rs1143679 (R77H) in disease association, our data also suggests an additional contribution from rs1143683, another non-synonymous polymorphism in this gene (A858V). Therefore, despite the low-allele frequencies of the risk alleles of the gene in our two Asian populations, ITGAM was confirmed to be a risk factor related to disease susceptibility and probably severe manifestations of SLE.

Genome-wide association studies (GWAS) of kidney function have uncovered hundreds of loci, primarily in populations of European ancestry. We have undertaken the first continental African GWAS of estimated glomerular filtration rate (eGFR), a measure of kidney function used to define chronic kidney disease (CKD). We conducted GWAS of eGFR in 3288 East Africans from the Uganda General Population Cohort (GPC) and replicated in 8224 African Americans from the Women's Health Initiative. Loci attaining genome-wide significant evidence for association (P < 5 × 10-8) were followed up with Bayesian fine-mapping to localize potential causal variants. The predictive power of a genetic risk score (GRS) constructed from previously reported trans-ancestry eGFR lead single nucleotide polymorphism (SNPs) was evaluated in the Uganda GPC. We identified and validated two eGFR loci. At the glycine amidinotransferase (GATM) locus, the association signal (lead SNP rs2433603, P = 1.0 × 10-8) in the Uganda GPC GWAS was distinct from previously reported signals at this locus. At the haemoglobin beta (HBB) locus, the association signal (lead SNP rs141845179, P = 3.0 × 10-8) has been previously reported. The lead SNP at the HBB locus accounted for 88% of the posterior probability of causality after fine-mapping, but did not colocalise with kidney expression quantitative trait loci. The trans-ancestry GRS of eGFR was not significantly predictive into the Ugandan population. In the first GWAS of eGFR in continental Africa, we validated two previously reported loci at GATM and HBB. At the GATM locus, the association signal was distinct from that previously reported. These results demonstrate the value of performing GWAS in continental Africans, providing a rich genomic resource to larger consortia for further discovery and fine-mapping. The study emphasizes that additional large-scale efforts in Africa are warranted to gain further insight into the genetic architecture of CKD.

Obscurin is a giant muscle protein that connects the sarcomere with the sarcoplasmic reticulum, and has poorly understood structural and signalling functions. Increasingly, obscurin variants are implicated in the pathophysiology of cardiovascular diseases. The Arg4344Gln variant (R4344Q) in obscurin domain Ig58, initially discovered in a patient with hypertrophic cardiomyopathy, has been reported to reduce binding to titin domains Z8-Z9, impairing obscurin's Z-disc localization. An R4344Q knock-in mouse developed a cardiomyopathy-like phenotype with abnormal Ca2+-handling and arrhythmias, which were attributed to an enhanced affinity of a putative interaction between obscurin Ig58 and phospholamban (PLN) due to the R4344Q variant. However, the R4344Q variant is found in 15% of African Americans, arguing against its pathogenicity. To resolve this apparent paradox, we quantified the influence of the R4344Q variant (alongside another potentially pathogenic variant: Arg4444Trp (R4444W)) on binding to titin Z8-Z9, novex-3 and PLN using pull-down assays and microscale thermophoresis and characterized the influence on domain stability using differential scanning fluorimetry. We found no changes in titin binding and thermostability for both variants and modestly increased affinities of PLN for R4344Q and R4444W. While we could not confirm the novex-3/obscurin interaction, the PLN/obscurin interaction relies on the transmembrane region of PLN and is not reproducible in mammalian cells, suggesting it is an in vitro artefact. Without clear clinical evidence for disease involvement, we advise against classifying these obscurin variants as pathogenic.

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP), and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the major histocompatibility complex region [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66-0.83, lowest P-value = 4.53 × 10-8]. The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95% CI = 0.90-0.97, P-value = 3.54 × 10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the human leukocyte antigen genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

Admixture mapping recently identified MYH9 as a susceptibility gene for idiopathic focal segmental glomerulosclerosis (FSGS), HIV-associated nephropathy (HIVAN) and end-stage kidney disease attributed to hypertension (H-ESKD) in African Americans (AA). MYH9 encodes the heavy chain of non-muscle myosin IIA, a cellular motor involved in motility. A haplotype and its tagging SNPs spanning introns 12-23 were most strongly associated with kidney disease (OR 2-7; P < 10(-8), recessive). To narrow the region of association and identify potential causal variation, we performed a dense-mapping study using 79 MYH9 SNPs in AA populations with FSGS, HIVAN and H-ESKD (typed for a subset of 46 SNPs), for a total of 2496 cases and controls. The strongest associations were for correlated SNPs rs5750250, rs2413396 and rs5750248 in introns 13, 14 and 15, a region of 5.6 kb. Rs5750250 showed OR 5.0, 8.0 and 2.8; P = 2 x 10(-17), 2 x 10(-10) and 3 x 10(-22), respectively, for FSGS, HIVAN and H-ESKD; OR 5.7; P = 9 x 10(-27) for combined FSGS and HIVAN, recessive. An independent association was observed for rs11912763 in intron 33. Neither the highly associated SNPs nor the results of resequencing MYH9 in 40 HIVAN or FSGS cases and controls revealed non-synonymous changes that could account for the disease associations. Rs2413396 and one of the highly associated SNPs in intron 23, rs4821480, are predicted splicing motif modifiers. Rs5750250 combined with rs11912763 had receiver operator characteristic (ROC) C statistics of 0.80, 0.73 and 0.65 for HIVAN, FSGS and H-ESKD, respectively, allowing prediction of genetic risk by typing two SNPs.