Neutrophil extracellular trap formation and gene programs distinguish TST/IGRA sensitization outcomes among Mycobacterium tuberculosis exposed persons living with HIV.

Persons living with HIV (PLWH) have an increased risk for tuberculosis (TB). After prolonged and repeated exposure, some PLWH never develop TB and show no evidence of immune sensitization to Mycobacterium tuberculosis (Mtb) as defined by persistently negative tuberculin skin tests (TST) and interferon gamma release assays (IGRA). This group has been identified and defined as HIV+ persistently TB, tuberculin and IGRA negative (HITTIN). To investigate potential innate mechanisms unique to individuals with the HITTIN phenotype we compared their neutrophil Mtb infection response to that of PLWH, with no TB history, but who test persistently IGRA positive, and tuberculin positive (HIT). Neutrophil samples from 17 HITTIN (PMNHITTIN) and 11 HIT (PMNHIT) were isolated and infected with Mtb H37Rv for 1h and 6h. RNA was extracted and used for RNAseq analysis. Since there was no significant differential transcriptional response at 1h between infected PMNHITTIN and PMNHIT, we focused on the 6h timepoint. When compared to uninfected PMN, PMNHITTIN displayed 3106 significantly upregulated and 3548 significantly downregulated differentially expressed genes (DEGs) (absolute cutoff of a log2FC of 0.2, FDR < 0.05) whereas PMNHIT demonstrated 3816 significantly upregulated and 3794 significantly downregulated DEGs following 6h Mtb infection. Contrasting the log2FC 6h infection response to Mtb from PMNHITTIN against PMNHIT, 2285 genes showed significant differential response between the two groups. Overall PMNHITTIN had a lower fold change response to Mtb infection compared to PMNHIT. According to pathway enrichment, Apoptosis and NETosis were differentially regulated between HITTIN and HIT PMN responses after 6h Mtb infection. To corroborate the blunted NETosis transcriptional response measured among HITTIN, fluorescence microscopy revealed relatively lower neutrophil extracellular trap formation and cell loss in PMNHITTIN compared to PMNHIT, showing that PMNHITTIN have a distinct response to Mtb.

Blood transcriptomic signatures for symptomatic tuberculosis in an African multicohort study.

BACKGROUND: Multiple host blood transcriptional signatures have been developed as non-sputum triage tests for tuberculosis (TB). We aimed to compare the diagnostic performance of 20 blood transcriptomic TB signatures for differentiating between symptomatic patients who have TB versus other respiratory diseases (ORD). METHODS: As part of a nested case-control study, individuals presenting with respiratory symptoms at primary healthcare clinics in Ethiopia, Malawi, Namibia, Uganda, South Africa and The Gambia were enrolled. TB was diagnosed based on clinical, microbiological and radiological findings. Transcriptomic signatures were measured in whole blood using microfluidic real-time quantitative PCR. Diagnostic performance was benchmarked against the World Health Organization Target Product Profile (TPP) for a non-sputum TB triage test. RESULTS: Among 579 participants, 158 had definite, microbiologically confirmed TB, 32 had probable TB, while 389 participants had ORD. Nine signatures differentiated between ORD and TB with equivalent performance (Satproedprai7: area under the curve 0.83 (95% CI 0.79-0.87); Jacobsen3: 0.83 (95% CI 0.79-0.86); Suliman2: 0.82 (95% CI 0.78-0.86); Roe1: 0.82 (95% CI 0.78-0.86); Kaforou22: 0.82 (95% CI 0.78-0.86); Sambarey10: 0.81 (95% CI 0.77-0.85); Duffy9: 0.81 (95% CI 0.76-0.86); Gliddon3: 0.8 (95% CI 0.75-0.85); Suliman4 0.79 (95% CI 0.75-0.84)). Benchmarked against a 90% sensitivity, these signatures achieved specificities between 44% (95% CI 38-49%) and 54% (95% CI 49-59%), not meeting the TPP criteria. Signature scores significantly varied by HIV status and country. In country-specific analyses, several signatures, such as Satproedprai7 and Penn-Nicholson6, met the minimal TPP criteria for a triage test in Ethiopia, Malawi and South Africa. CONCLUSION: No signatures met the TPP criteria in a pooled analysis of all countries, but several signatures met the minimum criteria for a non-sputum TB triage test in some countries.

Local ancestry inference in heterogeneous populations-Are recent recombination events more relevant?

To date, numerous software tools have been developed to infer recombination maps. Many of these software tools infer the recombination rate from linkage disequilibrium, and therefore they infer recombination many generations into the past. Other recently developed methods rely on the inference of recent recombination events to determine the recombination rate, such as identity by descent- and local ancestry inference (LAI)-based tools. Methods that mainly use recent recombination events to infer the recombination rate might be more relevant for certain analyses like LAI. We therefore describe a protocol for creating high-resolution, population-specific recombination maps using methods that mainly use recent recombination events and a method that uses recent and distant recombination events for recombination rate inference. Subsequently, we compared the effect of using maps inferred by these two paradigms on LAI accuracy.

Environmental and life-style risk factors for esophageal squamous cell carcinoma in Africa: a systematic review and meta-analysis.

BACKGROUND: The African Esophageal Squamous Cell Carcinoma (ESCC) corridor, which spans from Ethiopia down to South Africa, is an esophageal cancer hotspot. Disproportionately high incidence and mortality rates of esophageal cancer have been reported from this region. The aim of this study was to systematically assess the evidence on environmental and life-style risk factors associated with ESCC in African populations. METHODS: We followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines and carried out a comprehensive search of all African published studies up to March 2023 using PubMed, Embase, Scopus, and African Index Medicus databases. RESULTS: We identified 45 studies with measures of association [odds ratio (OR), relative risk (RR), and 95% confidence intervals (95%CI)], which reported on several environmental and lifestyle risk factors for ESCC in Africa. We performed a meta-analysis on 38 studies investigating tobacco, alcohol use, combined tobacco and alcohol use, polycyclic aromatic hydrocarbon exposure, hot food and beverages consumption (which served as a proxy for esophageal injury through exposure to high temperature), and poor oral health. We found significant associations between all the risk factors and ESCC development. Analysis of fruit and vegetable consumption showed a protective effect. Using population attributable fraction (PAF) analysis, we calculated the proportion of ESCC attributable to tobacco (18%), alcohol use (12%), combined tobacco and alcohol use (18%), polycyclic aromatic hydrocarbon exposure (12%), hot food and beverages intake (16%), poor oral health (37%), and fruit and vegetable consumption (-12%). CONCLUSIONS: Tobacco smoking and alcohol consumption were the most studied risk factors overall. Areas where there is an emerging body of evidence include hot food and beverages and oral health. Concurrently, new avenues of research are also emerging in PAH exposure, and diet as risk factors. Our results point to a multifactorial etiology of ESCC in African populations with further evidence on prevention potential.

Cascade Immune Mechanisms of Protection against Mycobacterium tuberculosis (IMPAc-TB): study protocol for the Household Contact Study in the Western Cape, South Africa.

BACKGROUND: Natural immunity against Mycobacterium tuberculosis exists, and > 90% of those infected remain disease-free. Innate and adaptive immune responses required to mediate such protection against tuberculosis (TB) are, however, poorly understood. METHODS: This is an analytical study exploring protective and non-protective pathways of immunity against Mycobacterium tuberculosis. Adults without HIV infection are recruited at community healthcare clinics in high TB incidence areas of the Western Cape Province, South Africa. Data regarding participants' medical, social and medication usage will be collected, and clinical examinations and point-of-care tests documented. Reference tests for TB (chest radiographs and sputum tests for GeneXpert MTB/RIF Ultra®, Auramine smear and liquid cultures) and investigations to classify infection states [interferon-gamma release assay (IGRA) and SARS-CoV-2 polymerase chain reaction (PCR) nasopharyngeal swab and IgG], are done on all participants who meet the inclusion criteria. 18F-Fluorodeoxyglucose positron emission tomography combined with computerized tomography will be done on all close contacts (contacts) and healthy control (controls) participants. Participants are divided into 12 study groups representing a spectrum of TB clinical phenotypes and prior SARS-CoV-2 infection based on their TB status, exposure history, results of IGRA test at baseline and 3 months, SARS-CoV-2 serology, and PCR results, and for contacts and controls, PET-CT imaging findings indicative of sub-clinical TB lesions. Samples for experimental assays include whole blood for isolation of peripheral blood mononuclear cells and blood in PAXgene® tubes for RNA isolation. All SARS-CoV-2 PCR negative study participants undergo bronchoscopy for collecting bronchoalveolar lavage samples. DISCUSSION: The paired blood and BAL samples will be used for comprehensive analyses of the tissue-specific and systemic immunity that will include e.g., cytometry by time-of-flight analyses, RNA-sequencing, multiplex immunoassays, epigenetic analysis, and mechanistic studies of control of infection by Mycobacterium tuberculosis. Results will be integrated with those from mice and non-human primate studies to provide a comprehensive analysis of protective pathways in natural and vaccine-induced immunity against Mycobacterium tuberculosis.

A multi-phenotype genome-wide association study of clades causing tuberculosis in a Ghanaian- and South African cohort.

Despite decades of research and advancements in diagnostics and treatment, tuberculosis remains a major public health concern. New computational methods are needed to interrogate the intersection of host- and bacterial genomes. Paired host genotype datum and infecting bacterial isolate information were analysed for associations using a multinomial logistic regression framework implemented in SNPTest. A cohort of 853 admixed South African participants and a Ghanaian cohort of 1359 participants were included. Two directly genotyped variants, namely rs529920 and rs41472447, were identified in the Ghanaian cohort as being statistically significantly associated with risk for infection with strains of different members of the MTBC. Thus, a multinomial logistic regression using paired host-pathogen data may prove valuable for investigating the complex relationships driving infectious disease.

Longitudinal telomere length profile does not reflect HIV and childhood trauma impacts on cognitive function in South African women.

HIV-associated neurocognitive disorders (HAND) present a challenge in South Africa where the burden of HIV infection is the highest. Identification of biological correlates of HAND is required to improve diagnosis and inform interventions. Telomeres maintain genomic integrity and their shortening is a marker of biological aging sensitive to environmental influences. This study examined relative telomere length (rTL) as a predictor of cognitive function in the context of HIV and childhood trauma (CT), a risk factor for HAND. Two hundred and eighty-six women completed a neurocognitive assessment battery and the Childhood Trauma Questionnaire-Short Form (CTQ). Quantitative polymerase chain reaction for amplification of telomeric repeats and the reference gene human beta-globin was used to calculate rTL. Neurocognitive and rTL assessments were repeated at 1 year in 110 participants. Cross-sectional and longitudinal data were assessed using linear and mixed models, respectively. Participants with HIV (n = 135 in cross-sectional and n = 62 in longitudinal study groups) reported more severe CT and had shorter baseline rTL compared to seronegative controls. Participants without HIV had a greater 1-year decline in rTL. Global cognitive and attention/working memory scores declined in participants with HIV. Our data indicate that baseline rTL in the context of CT and HIV did not predict decline in cognitive scores. HIV-associated pathophysiological processes driving cognitive decline may also engage mechanisms that protect against telomere shortening. The results highlight the importance of examining biological correlates in longitudinal studies.

Targeted next-generation sequencing identifies novel variants in candidate genes for Parkinson's disease in Black South African and Nigerian patients.

BACKGROUND: The prevalence of Parkinson's disease (PD) is increasing in sub-Saharan Africa, but little is known about the genetics of PD in these populations. Due to their unique ancestry and diversity, sub-Saharan African populations have the potential to reveal novel insights into the pathobiology of PD. In this study, we aimed to characterise the genetic variation in known and novel PD genes in a group of Black South African and Nigerian patients. METHODS: We recruited 33 Black South African and 14 Nigerian PD patients, and screened them for sequence variants in 751 genes using an Ion AmpliSeq™ Neurological Research panel. We used bcftools to filter variants and annovar software for the annotation. Rare variants were prioritised using MetaLR and MetaSVM prediction scores. The effect of a variant on ATP13A2's protein structure was investigated by molecular modelling. RESULTS: We identified 14,655 rare variants with a minor allele frequency ≤ 0.01, which included 2448 missense variants. Notably, no common pathogenic mutations were identified in these patients. Also, none of the known PD-associated mutations were found highlighting the need for more studies in African populations. Altogether, 54 rare variants in 42 genes were considered deleterious and were prioritized, based on MetaLR and MetaSVM scores, for follow-up studies. Protein modelling showed that the S1004R variant in ATP13A2 possibly alters the conformation of the protein. CONCLUSIONS: We identified several rare variants predicted to be deleterious in sub-Saharan Africa PD patients; however, further studies are required to determine the biological effects of these variants and their possible role in PD. Studies such as these are important to elucidate the genetic aetiology of this disorder in patients of African ancestry.

The X chromosome and sex-specific effects in infectious disease susceptibility.

The X chromosome and X-linked variants have largely been ignored in genome-wide and candidate association studies of infectious diseases due to the complexity of statistical analysis of the X chromosome. This exclusion is significant, since the X chromosome contains a high density of immune-related genes and regulatory elements that are extensively involved in both the innate and adaptive immune responses. Many diseases present with a clear sex bias, and apart from the influence of sex hormones and socioeconomic and behavioural factors, the X chromosome, X-linked genes and X chromosome inactivation mechanisms contribute to this difference. Females are functional mosaics for X-linked genes due to X chromosome inactivation and this, combined with other X chromosome inactivation mechanisms such as genes that escape silencing and skewed inactivation, could contribute to an immunological advantage for females in many infections. In this review, we discuss the involvement of the X chromosome and X inactivation in immunity and address its role in sexual dimorphism of infectious diseases using tuberculosis susceptibility as an example, in which male sex bias is clear, yet not fully explored.

Meta-Analysis of Genome-Wide Association Studies for Abdominal Aortic Aneurysm Identifies Four New Disease-Specific Risk Loci.

RATIONALE: Abdominal aortic aneurysm (AAA) is a complex disease with both genetic and environmental risk factors. Together, 6 previously identified risk loci only explain a small proportion of the heritability of AAA. OBJECTIVE: To identify additional AAA risk loci using data from all available genome-wide association studies. METHODS AND RESULTS: Through a meta-analysis of 6 genome-wide association study data sets and a validation study totaling 10 204 cases and 107 766 controls, we identified 4 new AAA risk loci: 1q32.3 (SMYD2), 13q12.11 (LINC00540), 20q13.12 (near PCIF1/MMP9/ZNF335), and 21q22.2 (ERG). In various database searches, we observed no new associations between the lead AAA single nucleotide polymorphisms and coronary artery disease, blood pressure, lipids, or diabetes mellitus. Network analyses identified ERG, IL6R, and LDLR as modifiers of MMP9, with a direct interaction between ERG and MMP9. CONCLUSIONS: The 4 new risk loci for AAA seem to be specific for AAA compared with other cardiovascular diseases and related traits suggesting that traditional cardiovascular risk factor management may only have limited value in preventing the progression of aneurysmal disease.

Four-Gene Pan-African Blood Signature Predicts Progression to Tuberculosis.

Rationale: Contacts of patients with tuberculosis (TB) constitute an important target population for preventive measures because they are at high risk of infection with Mycobacterium tuberculosis and progression to disease.Objectives: We investigated biosignatures with predictive ability for incident TB.Methods: In a case-control study nested within the Grand Challenges 6-74 longitudinal HIV-negative African cohort of exposed household contacts, we employed RNA sequencing, PCR, and the pair ratio algorithm in a training/test set approach. Overall, 79 progressors who developed TB between 3 and 24 months after diagnosis of index case and 328 matched nonprogressors who remained healthy during 24 months of follow-up were investigated.Measurements and Main Results: A four-transcript signature derived from samples in a South African and Gambian training set predicted progression up to two years before onset of disease in blinded test set samples from South Africa, the Gambia, and Ethiopia with little population-associated variability, and it was also validated in an external cohort of South African adolescents with latent M. tuberculosis infection. By contrast, published diagnostic or prognostic TB signatures were predicted in samples from some but not all three countries, indicating site-specific variability. Post hoc meta-analysis identified a single gene pair, C1QC/TRAV27 (complement C1q C-chain / T-cell receptor-α variable gene 27) that would consistently predict TB progression in household contacts from multiple African sites but not in infected adolescents without known recent exposure events.Conclusions: Collectively, we developed a simple whole blood-based PCR test to predict TB in recently exposed household contacts from diverse African populations. This test has potential for implementation in national TB contact investigation programs.

Epistatic Gene-Based Interaction Analyses for Glaucoma in eMERGE and NEIGHBOR Consortium.

Primary open angle glaucoma (POAG) is a complex disease and is one of the major leading causes of blindness worldwide. Genome-wide association studies have successfully identified several common variants associated with glaucoma; however, most of these variants only explain a small proportion of the genetic risk. Apart from the standard approach to identify main effects of variants across the genome, it is believed that gene-gene interactions can help elucidate part of the missing heritability by allowing for the test of interactions between genetic variants to mimic the complex nature of biology. To explain the etiology of glaucoma, we first performed a genome-wide association study (GWAS) on glaucoma case-control samples obtained from electronic medical records (EMR) to establish the utility of EMR data in detecting non-spurious and relevant associations; this analysis was aimed at confirming already known associations with glaucoma and validating the EMR derived glaucoma phenotype. Our findings from GWAS suggest consistent evidence of several known associations in POAG. We then performed an interaction analysis for variants found to be marginally associated with glaucoma (SNPs with main effect p-value <0.01) and observed interesting findings in the electronic MEdical Records and GEnomics Network (eMERGE) network dataset. Genes from the top epistatic interactions from eMERGE data (Likelihood Ratio Test i.e. LRT p-value <1e-05) were then tested for replication in the NEIGHBOR consortium dataset. To replicate our findings, we performed a gene-based SNP-SNP interaction analysis in NEIGHBOR and observed significant gene-gene interactions (p-value <0.001) among the top 17 gene-gene models identified in the discovery phase. Variants from gene-gene interaction analysis that we found to be associated with POAG explain 3.5% of additional genetic variance in eMERGE dataset above what is explained by the SNPs in genes that are replicated from previous GWAS studies (which was only 2.1% variance explained in eMERGE dataset); in the NEIGHBOR dataset, adding replicated SNPs from gene-gene interaction analysis explain 3.4% of total variance whereas GWAS SNPs alone explain only 2.8% of variance. Exploring gene-gene interactions may provide additional insights into many complex traits when explored in properly designed and powered association studies.

Genetic risk models: Influence of model size on risk estimates and precision.

Disease risk estimation plays an important role in disease prevention. Many studies have found that the ability to predict risk improves as the number of risk single-nucleotide polymorphisms (SNPs) in the risk model increases. However, the width of the confidence interval of the risk estimate is often not considered in the evaluation of the risk model. Here, we explore how the risk and the confidence interval width change as more SNPs are added to the model in the order of decreasing effect size, using both simulated data and real data from studies of abdominal aortic aneurysms and age-related macular degeneration. Our results show that confidence interval width is positively correlated with model size and the majority of the bigger models have wider confidence interval widths than smaller models. Once the model size is bigger than a certain level, the risk does not shift markedly, as 100% of the risk estimates of the one-SNP-bigger models lie inside the confidence interval of the one-SNP-smaller models. We also created a confidence interval-augmented reclassification table. It shows that both more effective SNPs with larger odds ratios and less effective SNPs with smaller odds ratios contribute to the correct decision of whom to screen. The best screening strategy is selected and evaluated by the net benefit quantity and the reclassification rate. We suggest that individuals whose upper bound of their risk confidence interval is above the screening threshold, which corresponds to the population prevalence of the disease, should be screened.

Penetrance of Hemochromatosis in HFE Genotypes Resulting in p.Cys282Tyr and p.[Cys282Tyr];[His63Asp] in the eMERGE Network.

Hereditary hemochromatosis (HH) is a common autosomal-recessive disorder associated with pathogenic HFE variants, most commonly those resulting in p.Cys282Tyr and p.His63Asp. Recommendations on returning incidental findings of HFE variants in individuals undergoing genome-scale sequencing should be informed by penetrance estimates of HH in unselected samples. We used the eMERGE Network, a multicenter cohort with genotype data linked to electronic medical records, to estimate the diagnostic rate and clinical penetrance of HH in 98 individuals homozygous for the variant coding for HFE p.Cys282Tyr and 397 compound heterozygotes with variants resulting in p.[His63Asp];[Cys282Tyr]. The diagnostic rate of HH in males was 24.4% for p.Cys282Tyr homozygotes and 3.5% for compound heterozygotes (p < 0.001); in females, it was 14.0% for p.Cys282Tyr homozygotes and 2.3% for compound heterozygotes (p < 0.001). Only males showed differences across genotypes in transferrin saturation levels (100% of homozygotes versus 37.5% of compound heterozygotes with transferrin saturation > 50%; p = 0.003), serum ferritin levels (77.8% versus 33.3% with serum ferritin > 300 ng/ml; p = 0.006), and diabetes (44.7% versus 28.0%; p = 0.03). No differences were found in the prevalence of heart disease, arthritis, or liver disease, except for the rate of liver biopsy (10.9% versus 1.8% [p = 0.013] in males; 9.1% versus 2% [p = 0.035] in females). Given the higher rate of HH diagnosis than in prior studies, the high penetrance of iron overload, and the frequency of at-risk genotypes, in addition to other suggested actionable adult-onset genetic conditions, opportunistic screening should be considered for p.[Cys282Tyr];[Cys282Tyr] individuals with existing genomic data.

Viruses as indicators of contemporary host dispersal and phylogeography: an example of feline immunodeficiency virus (FIVPle ) in free-ranging African lion (Panthera leo).

Measuring contemporary dispersal in highly mobile terrestrial species is challenging, especially when species are characterized by low levels of population differentiation. Directly transmitted viruses can be used as a surrogate for traditional methods of tracking host movement. Feline immunodeficiency virus (FIV) is a species-specific lentivirus, which has an exceptionally high mutation rate and circulates naturally in wild felids. Using samples derived from 35 lion (Panthera leo) prides, we tested the prediction that FIV in lions (FIVPle ) can be used to track the dispersal of individuals between prides. As FIVPle subtypes are geographically structured throughout Africa, we predicted that this marker could be used to detect phylogeographic structure of lions at smaller spatial scales. Phylogenetic analyses of FIVPle pol-RT sequences showed that core pride members (females and subadults) shared evolutionary close viral lineages which differed from neighbouring core prides, whereas sequences from sexually mature males associated with the same pride were always the most divergent. In six instances, natal pride associations of divergent male lions could be inferred, on the assumption that FIVPle infections are acquired during early life stages. Congruence between the genetic pattern of FIV and pride structure suggests that vertical transmission plays an important role in lion FIV dynamics. At a fine spatial scale, significant viral geographic structuring was also detected between lions occurring north of the Olifants River within the Kruger National Park (KNP) and those occupying the southern and central regions. This pattern was further supported by phylogenetic analyses and the confinement of FIVPle subtype E to the northern region of KNP. The study provides new insights into the use of retroviral sequences to predict host dispersal and fine-scale contemporary geographic structure in a social felid species.

SOS2 and ACP1 Loci Identified through Large-Scale Exome Chip Analysis Regulate Kidney Development and Function.

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium (nStage1: 111,666; nStage2: 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea (PPM1J, EDEM3, ACP1, SPEG, EYA4, CYP1A1, and ATXN2L; PStage1<3.7×10-7), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 (P=5.4×10-8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2-knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

Identifying gene-gene interactions that are highly associated with four quantitative lipid traits across multiple cohorts.

Genetic loci explain only 25-30 % of the heritability observed in plasma lipid traits. Epistasis, or gene-gene interactions may contribute to a portion of this missing heritability. Using the genetic data from five NHLBI cohorts of 24,837 individuals, we combined the use of the quantitative multifactor dimensionality reduction (QMDR) algorithm with two SNP-filtering methods to exhaustively search for SNP-SNP interactions that are associated with HDL cholesterol (HDL-C), LDL cholesterol (LDL-C), total cholesterol (TC) and triglycerides (TG). SNPs were filtered either on the strength of their independent effects (main effect filter) or the prior knowledge supporting a given interaction (Biofilter). After the main effect filter, QMDR identified 20 SNP-SNP models associated with HDL-C, 6 associated with LDL-C, 3 associated with TC, and 10 associated with TG (permutation P value <0.05). With the use of Biofilter, we identified 2 SNP-SNP models associated with HDL-C, 3 associated with LDL-C, 1 associated with TC and 8 associated with TG (permutation P value <0.05). In an independent dataset of 7502 individuals from the eMERGE network, we replicated 14 of the interactions identified after main effect filtering: 11 for HDL-C, 1 for LDL-C and 2 for TG. We also replicated 23 of the interactions found to be associated with TG after applying Biofilter. Prior knowledge supports the possible role of these interactions in the genetic etiology of lipid traits. This study also presents a computationally efficient pipeline for analyzing data from large genotyping arrays and detecting SNP-SNP interactions that are not primarily driven by strong main effects.

Biology-Driven Gene-Gene Interaction Analysis of Age-Related Cataract in the eMERGE Network.

Bioinformatics approaches to examine gene-gene models provide a means to discover interactions between multiple genes that underlie complex disease. Extensive computational demands and adjusting for multiple testing make uncovering genetic interactions a challenge. Here, we address these issues using our knowledge-driven filtering method, Biofilter, to identify putative single nucleotide polymorphism (SNP) interaction models for cataract susceptibility, thereby reducing the number of models for analysis. Models were evaluated in 3,377 European Americans (1,185 controls, 2,192 cases) from the Marshfield Clinic, a study site of the Electronic Medical Records and Genomics (eMERGE) Network, using logistic regression. All statistically significant models from the Marshfield Clinic were then evaluated in an independent dataset of 4,311 individuals (742 controls, 3,569 cases), using independent samples from additional study sites in the eMERGE Network: Mayo Clinic, Group Health/University of Washington, Vanderbilt University Medical Center, and Geisinger Health System. Eighty-three SNP-SNP models replicated in the independent dataset at likelihood ratio test P < 0.05. Among the most significant replicating models was rs12597188 (intron of CDH1)-rs11564445 (intron of CTNNB1). These genes are known to be involved in processes that include: cell-to-cell adhesion signaling, cell-cell junction organization, and cell-cell communication. Further Biofilter analysis of all replicating models revealed a number of common functions among the genes harboring the 83 replicating SNP-SNP models, which included signal transduction and PI3K-Akt signaling pathway. These findings demonstrate the utility of Biofilter as a biology-driven method, applicable for any genome-wide association study dataset.

The complete genome sequence of the African buffalo (Syncerus caffer).

BACKGROUND: The African buffalo (Syncerus caffer) is an important role player in the savannah ecosystem. It has become a species of relevance because of its role as a wildlife maintenance host for an array of infectious and zoonotic diseases some of which include corridor disease, foot-and-mouth disease and bovine tuberculosis. To date, no complete genome sequence for S. caffer had been available for study and the genomes of other species such as the domestic cow (Bos taurus) had been used as a proxy for any genetics analysis conducted on this species. Here, the high coverage genome sequence of the African buffalo (S. caffer) is presented. RESULTS: A total of 19,765 genes were predicted and 19,296 genes could be successfully annotated to S. caffer while 469 genes remained unannotated. Moreover, in order to extend a detailed annotation of S. caffer, gene clusters were constructed using twelve additional mammalian genomes. The S. caffer genome contains 10,988 gene clusters, of which 62 are shared exclusively between B. taurus and S. caffer. CONCLUSIONS: This study provides a unique genomic perspective for the S. caffer, allowing for the identification of novel variants that may play a role in the natural history and physiological adaptations.

Regulatory polymorphisms in human DBH affect peripheral gene expression and sympathetic activity.

RATIONALE: Dopamine ß-hydroxylase (DBH) catalyzes the conversion of dopamine to norepinephrine in the central nervous system and peripherally. DBH variants are associated with large changes in circulating DBH and implicated in multiple disorders; yet causal relationships and tissue-specific effects remain unresolved. OBJECTIVE: To characterize regulatory variants in DBH, effect on mRNA expression, and role in modulating sympathetic tone and disease risk. METHODS AND RESULTS: Analysis of DBH mRNA in human tissues confirmed high expression in the locus coeruleus and adrenal gland, but also in sympathetically innervated organs (liver>lung>heart). Allele-specific mRNA assays revealed pronounced allelic expression differences in the liver (2- to 11-fold) attributable to promoter rs1611115 and exon 2 rs1108580, but only small differences in locus coeruleus and adrenals. These alleles were also associated with significantly reduced mRNA expression in liver and lung. Although DBH protein is expressed in other sympathetically innervated organs, mRNA levels were too low for analysis. In mice, hepatic Dbh mRNA levels correlated with cardiovascular risk phenotypes. The minor alleles of rs1611115 and rs1108580 were associated with sympathetic phenotypes, including angina pectoris. Testing combined effects of these variants suggested protection against myocardial infarction in 3 separate clinical cohorts. CONCLUSIONS: We demonstrate profound effects of DBH variants on expression in 2 sympathetically innervated organs, liver and lung, but not in adrenals and brain. Preliminary results demonstrate an association of these variants with clinical phenotypes responsive to peripheral sympathetic tone. We hypothesize that in addition to endocrine effects via circulating DBH and norepinephrine, the variants act in sympathetically innervated target organs.