Genome-wide linkage analysis of 972 bipolar pedigrees using single-nucleotide polymorphisms.

Because of the high costs associated with ascertainment of families, most linkage studies of Bipolar I disorder (BPI) have used relatively small samples. Moreover, the genetic information content reported in most studies has been less than 0.6. Although microsatellite markers spaced every 10 cM typically extract most of the genetic information content for larger multiplex families, they can be less informative for smaller pedigrees especially for affected sib pair kindreds. For these reasons we collaborated to pool family resources and carried out higher density genotyping. Approximately 1100 pedigrees of European ancestry were initially selected for study and were genotyped by the Center for Inherited Disease Research using the Illumina Linkage Panel 12 set of 6090 single-nucleotide polymorphisms. Of the ~1100 families, 972 were informative for further analyses, and mean information content was 0.86 after pruning for linkage disequilibrium. The 972 kindreds include 2284 cases of BPI disorder, 498 individuals with bipolar II disorder (BPII) and 702 subjects with recurrent major depression. Three affection status models (ASMs) were considered: ASM1 (BPI and schizoaffective disorder, BP cases (SABP) only), ASM2 (ASM1 cases plus BPII) and ASM3 (ASM2 cases plus recurrent major depression). Both parametric and non-parametric linkage methods were carried out. The strongest findings occurred at 6q21 (non-parametric pairs LOD 3.4 for rs1046943 at 119 cM) and 9q21 (non-parametric pairs logarithm of odds (LOD) 3.4 for rs722642 at 78 cM) using only BPI and schizoaffective (SA), BP cases. Both results met genome-wide significant criteria, although neither was significant after correction for multiple analyses. We also inspected parametric scores for the larger multiplex families to identify possible rare susceptibility loci. In this analysis, we observed 59 parametric LODs of 2 or greater, many of which are likely to be close to maximum possible scores. Although some linkage findings may be false positives, the results could help prioritize the search for rare variants using whole exome or genome sequencing

Reconstructing Druze population history

The Druze are an aggregate of communities in the Levant and Near East living almost exclusively in the mountains of Syria, Lebanon and Israel whose ~1000 year old religion formally opposes mixed marriages and conversions. Despite increasing interest in genetics of the population structure of the Druze, their population history remains unknown. We investigated the genetic relationships between Israeli Druze and both modern and ancient populations. We evaluated our findings in light of three hypotheses purporting to explain Druze history that posit Arabian, Persian or mixed Near Eastern-Levantine roots. The biogeographical analysis localised proto-Druze to the mountainous regions of southeastern Turkey, northern Iraq and southeast Syria and their descendants clustered along a trajectory between these two regions. The mixed Near Eastern-Middle Eastern localisation of the Druze, shown using both modern and ancient DNA data, is distinct from that of neighbouring Syrians, Palestinians and most of the Lebanese, who exhibit a high affinity to the Levant. Druze biogeographic affinity, migration patterns, time of emergence and genetic similarity to Near Eastern populations are highly suggestive of Armenian-Turkish ancestries for the proto-Druze

Reconstructing Roma History from Genome-Wide Data

The Roma people, living throughout Europe and West Asia, are a diverse population linked by the Romani language and culture. Previous linguistic and genetic studies have suggested that the Roma migrated into Europe from South Asia about 1,000–1,500 years ago. Genetic inferences about Roma history have mostly focused on the Y chromosome and mitochondrial DNA. To explore what additional information can be learned from genome-wide data, we analyzed data from six Roma groups that we genotyped at hundreds of thousands of single nucleotide polymorphisms (SNPs). We estimate that the Roma harbor about 80% West Eurasian ancestry–derived from a combination of European and South Asian sources–and that the date of admixture of South Asian and European ancestry was about 850 years before present. We provide evidence for Eastern Europe being a major source of European ancestry, and North-west India being a major source of the South Asian ancestry in the Roma. By computing allele sharing as a measure of linkage disequilibrium, we estimate that the migration of Roma out of the Indian subcontinent was accompanied by a severe founder event, which appears to have been followed by a major demographic expansion after the arrival in Europe.Országos Tudományos Kutatási Alapprogramok (OTKA K 103983)Országos Tudományos Kutatási Alapprogramok (OTKA 73430)National Science Foundation (U.S.) (HOMINID grant 1032255)National Institutes of Health (U.S.) (grant GM100233

Extensive Copy-Number Variation of Young Genes across Stickleback Populations

MM received funding from the Max Planck innovation funds for this project. PGDF was supported by a Marie Curie European Reintegration Grant (proposal nr 270891). CE was supported by German Science Foundation grants (DFG, EI 841/4-1 and EI 841/6-1). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

Reconstructing Native American Population History

The peopling of the Americas has been the subject of extensive genetic, archaeological and linguistic research; however, central questions remain unresolved1–5. One contentious issue is whether the settlement occurred via a single6–8 or multiple streams of migration from Siberia9–15. The pattern of dispersals within the Americas is also poorly understood. To address these questions at higher resolution than was previously possible, we assembled data from 52 Native American and 17 Siberian groups genotyped at 364,470 single nucleotide polymorphisms. We show that Native Americans descend from at least three streams of Asian gene flow. Most descend entirely from a single ancestral population that we call “First American”. However, speakers of Eskimo-Aleut languages from the Arctic inherit almost half their ancestry from a second stream of Asian gene flow, and the Na-Dene-speaking Chipewyan from Canada inherit roughly one-tenth of their ancestry from a third stream. We show that the initial peopling followed a southward expansion facilitated by the coast, with sequential population splits and little gene flow after divergence, especially in South America. A major exception is in Chibchan-speakers on both sides of the Panama Isthmus, who have ancestry from both North and South America

Whole-genome analysis of introgressive hybridization and characterization of the bovine legacy of Mongolian yaks

The yak is remarkable for its adaptation to high altitude and occupies a central place in the economies of the mountainous regions of Asia. At lower elevations, it is common to hybridize yaks with cattle to combine the yak’s hardiness with the productivity of cattle. Hybrid males are sterile, however, preventing the establishment of stable hybrid populations, but not a limited introgression after backcrossing several generations of female hybrids to male yaks. Here we inferred bovine haplotypes in the genomes of 76 Mongolian yaks using high-density SNP genotyping and whole-genome sequencing. These yaks inherited ~1.3% of their genome from bovine ancestors after nearly continuous admixture over at least the last 1,500 years. The introgressed regions are enriched in genes involved in nervous system development and function, and particularly in glutamate metabolism and neurotransmission. We also identified a novel mutation associated with a polled (hornless) phenotype originating from Mongolian Turano cattle. Our results suggest that introgressive hybridization contributed to the improvement of yak management and breeding

Lymph node hemophagocytosis in rickettsial diseases: a pathogenetic role for CD8 T lymphocytes in human monocytic ehrlichiosis (HME)?

BACKGROUND: Human monocytic ehrlichiosis (HME) and Rocky Mountain spotted fever (RMSF) are caused by Ehrlichia chaffeensis and Rickettsia rickettsii, respectively. The pathogenesis of RMSF relates to rickettsia-mediated vascular injury, but it is unclear in HME. METHODS: To study histopathologic responses in the lymphatic system for correlates of immune injury, lymph nodes from patients with HME (n = 6) and RMSF (n = 5) were examined. H&E-stained lymph node tissues were examined for five histopathologic features, including hemophagocytosis, cellularity, necrosis, and vascular congestion and edema. The relative proportions of CD68 macrophages, CD8 and CD4 T lymphocytes, and CD20 B lymphocytes were evaluated by immunohistochemical staining. RESULTS: Hemophagocytosis was similar in HME and RMSF, and was greater than in control cases (p = .015). Cellularity in HME was not different from controls, whereas RMSF lymph nodes were markedly less cellular (p < 0.002). E. chaffeensis-infected mononuclear phagocytes were infrequent compared to R. rickettsii-infected endothelial cells. More CD8 cells in lymph nodes were observed with HME (p < .001), but no quantitative differences in CD4 lymphocytes, macrophages, or B lymphocytes were identified. CONCLUSION: Hemophagocytosis, CD8 T cell expansion, and the paucity of infected cells in HME, suggest that E. chaffeensis infection leads to macrophage activation and immune-mediated injury

Analyses and Comparison of Imputation-Based Association Methods

Genotype imputation methods have become increasingly popular for recovering untyped genotype data. An important application with imputed genotypes is to test genetic association for diseases. Imputation-based association test can provide additional insight beyond what is provided by testing on typed tagging SNPs only. A variety of effective imputation-based association tests have been proposed. However, their performances are affected by a variety of genetic factors, which have not been well studied. In this study, using both simulated and real data sets, we investigated the effects of LD, MAF of untyped causal SNP and imputation accuracy rate on the performances of seven popular imputation-based association methods, including MACH2qtl/dat, SNPTEST, ProbABEL, Beagle, Plink, BIMBAM and SNPMStat. We also aimed to provide a comprehensive comparison among methods. Results show that: 1). imputation-based association tests can boost signals and improve power under medium and high LD levels, with the power improvement increasing with strengthening LD level; 2) the power increases with higher MAF of untyped causal SNPs under medium to high LD level; 3). under low LD level, a high imputation accuracy rate cannot guarantee an improvement of power; 4). among methods, MACH2qtl/dat, ProbABEL and SNPTEST perform similarly and they consistently outperform other methods. Our results are helpful in guiding the choice of imputation-based association test in practical application

A Novel Tandem Mass Spectrometry Method for Rapid Confirmation of Medium- and Very Long-Chain acyl-CoA Dehydrogenase Deficiency in Newborns

BACKGROUND:Newborn screening for medium- and very long-chain acyl-CoA dehydrogenase (MCAD and VLCAD, respectively) deficiency, using acylcarnitine profiling with tandem mass spectrometry, has increased the number of patients with fatty acid oxidation disorders due to the identification of additional milder, and so far silent, phenotypes. However, especially for VLCADD, the acylcarnitine profile can not constitute the sole parameter in order to reliably confirm disease. Therefore, we developed a new liquid chromatography tandem mass spectrometry (LC-MS/MS) method to rapidly determine both MCAD- and/or VLCAD-activity in human lymphocytes in order to confirm diagnosis. METHODOLOGY:LC-MS/MS was used to measure MCAD- or VLCAD-catalyzed production of enoyl-CoA and hydroxyacyl-CoA, in human lymphocytes. PRINCIPAL FINDINGS:VLCAD activity in controls was 6.95+/-0.42 mU/mg (range 1.95 to 11.91 mU/mg). Residual VLCAD activity of 4 patients with confirmed VLCAD-deficiency was between 0.3 and 1.1%. Heterozygous ACADVL mutation carriers showed residual VLCAD activities of 23.7 to 54.2%. MCAD activity in controls was 2.38+/-0.18 mU/mg. In total, 28 patients with suspected MCAD-deficiency were assayed. Nearly all patients with residual MCAD activities below 2.5% were homozygous 985A>G carriers. MCAD-deficient patients with one other than the 985A>G mutation had higher MCAD residual activities, ranging from 5.7 to 13.9%. All patients with the 199T>C mutation had residual activities above 10%. CONCLUSIONS:Our newly developed LC-MS/MS method is able to provide ample sensitivity to correctly and rapidly determine MCAD and VLCAD residual activity in human lymphocytes. Importantly, based on measured MCAD residual activities in correlation with genotype, new insights were obtained on the expected clinical phenotype