Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity.

Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations

Population History and Gene Divergence in Native Mexicans Inferred from 76 Human Exomes.

Native American genetic variation remains underrepresented in most catalogs of human genome sequencing data. Previous genotyping efforts have revealed that Mexico's Indigenous population is highly differentiated and substructured, thus potentially harboring higher proportions of private genetic variants of functional and biomedical relevance. Here we have targeted the coding fraction of the genome and characterized its full site frequency spectrum by sequencing 76 exomes from five Indigenous populations across Mexico. Using diffusion approximations, we modeled the demographic history of Indigenous populations from Mexico with northern and southern ethnic groups splitting 7.2 KYA and subsequently diverging locally 6.5 and 5.7 KYA, respectively. Selection scans for positive selection revealed BCL2L13 and KBTBD8 genes as potential candidates for adaptive evolution in Rarámuris and Triquis, respectively. BCL2L13 is highly expressed in skeletal muscle and could be related to physical endurance, a well-known phenotype of the northern Mexico Rarámuri. The KBTBD8 gene has been associated with idiopathic short stature and we found it to be highly differentiated in Triqui, a southern Indigenous group from Oaxaca whose height is extremely low compared to other Native populations

Towards Equity in Health: Researchers Take Stock.

For the 2016 end-of-the-year editorial, the PLOS Medicine editors asked 7 global health leaders to discuss developments relevant to the equitable provision of medical care to all populations. The result is a collection of expert views on ethical trial design, research during outbreaks, high-burden infectious diseases, diversity in research and protection of migrants

Population sequencing data reveal a compendium of mutational processes in human germline

Mechanistic processes underlying human germline mutations remain largely unknown.Variation in mutation rate and spectra along the genome is informative about the biological mechanisms. We statistically decompose this variation into separate processes using a blind source separation technique. The analysis of a large-scale whole genome sequencing dataset (TOPMed) reveals nine processes that explain the variation in mutation properties between loci. Seven of these processes lend themselves to a biological interpretation. One process is driven by bulky DNA lesions that resolve asymmetrically with respect to transcription and replication. Two processes independently track direction of replication fork and replication timing. We identify a mutagenic effect of active demethylation primarily acting in regulatory regions. We also demonstrate that a recently discovered mutagenic process specific to oocytes can be localized solely from population sequencing data. This process is spread across all chromosomes and is highly asymmetric with respect to the direction of transcription, suggesting a major role of DNA damage

PCA-Correlated SNPs for Structure Identification in Worldwide Human Populations

Existing methods to ascertain small sets of markers for the identification of human population structure require prior knowledge of individual ancestry. Based on Principal Components Analysis (PCA), and recent results in theoretical computer science, we present a novel algorithm that, applied on genomewide data, selects small subsets of SNPs (PCA-correlated SNPs) to reproduce the structure found by PCA on the complete dataset, without use of ancestry information. Evaluating our method on a previously described dataset (10,805 SNPs, 11 populations), we demonstrate that a very small set of PCA-correlated SNPs can be effectively employed to assign individuals to particular continents or populations, using a simple clustering algorithm. We validate our methods on the HapMap populations and achieve perfect intercontinental differentiation with 14 PCA-correlated SNPs. The Chinese and Japanese populations can be easily differentiated using less than 100 PCA-correlated SNPs ascertained after evaluating 1.7 million SNPs from HapMap. We show that, in general, structure informative SNPs are not portable across geographic regions. However, we manage to identify a general set of 50 PCA-correlated SNPs that effectively assigns individuals to one of nine different populations. Compared to analysis with the measure of informativeness, our methods, although unsupervised, achieved similar results. We proceed to demonstrate that our algorithm can be effectively used for the analysis of admixed populations without having to trace the origin of individuals. Analyzing a Puerto Rican dataset (192 individuals, 7,257 SNPs), we show that PCA-correlated SNPs can be used to successfully predict structure and ancestry proportions. We subsequently validate these SNPs for structure identification in an independent Puerto Rican dataset. The algorithm that we introduce runs in seconds and can be easily applied on large genome-wide datasets, facilitating the identification of population substructure, stratification assessment in multi-stage whole-genome association studies, and the study of demographic history in human populations

Expression of SMARCD1 interacts with age in association with asthma control on inhaled corticosteroid therapy.

BackgroundGlobal gene expression levels are known to be highly dependent upon gross demographic features including age, yet identification of age-related genomic indicators has yet to be comprehensively undertaken in a disease and treatment-specific context.MethodsWe used gene expression data from CD4+ lymphocytes in the Asthma BioRepository for Integrative Genomic Exploration (Asthma BRIDGE), an open-access collection of subjects participating in genetic studies of asthma with available gene expression data. Replication population participants were Puerto Rico islanders recruited as part of the ongoing Genes environments & Admixture in Latino Americans (GALA II), who provided nasal brushings for transcript sequencing. The main outcome measure was chronic asthma control as derived by questionnaires. Genomic associations were performed using regression of chronic asthma control score on gene expression with age in years as a covariate, including a multiplicative interaction term for gene expression times age.ResultsThe SMARCD1 gene (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily D member 1) interacted with age to influence chronic asthma control on inhaled corticosteroids, with a doubling of expression leading to an increase of 1.3 units of chronic asthma control per year (95% CI [0.86, 1.74], p = 6 × 10- 9), suggesting worsening asthma control with increasing age. This result replicated in GALA II (p = 3.8 × 10- 8). Cellular assays confirmed the role of SMARCD1 in glucocorticoid response in airway epithelial cells.ConclusionFocusing on age-dependent factors may help identify novel indicators of asthma medication response. Age appears to modulate the effect of SMARCD1 on asthma control with inhaled corticosteroids

Identification of CFTR variants in Latino patients with cystic fibrosis from the Dominican Republic and Puerto Rico

BackgroundIn cystic fibrosis (CF), the spectrum and frequency of CFTR variants differ by geography and race/ethnicity. CFTR variants in White patients are wellâ described compared with Latino patients. No studies of CFTR variants have been done in patients with CF in the Dominican Republic or Puerto Rico.MethodsCFTR was sequenced in 61 Dominican Republican patients and 21 Puerto Rican patients with CF and greater than â â â â 60â mmol/L sweat chloride. The spectrum of CFTR variants was identified and the proportion of patients with 0, 1, or 2 CFTR variants identified was determined. The functional effects of identified CFTR variants were investigated using clinical annotation databases and computational prediction tools.ResultsOur study found 10% of Dominican patients had two CFTR variants identified compared with 81% of Puerto Rican patients. No CFTR variants were identified in 69% of Dominican patients and 10% of Puerto Rican patients. In Dominican patients, there were 19 identified CFTR variants, accounting for 25 out of 122 disease alleles (20%). In Puerto Rican patients, there were 16 identified CFTR variants, accounting for 36 out of 42 disease alleles (86%) in Puerto Rican patients. Thirty CFTR variants were identified overall. The most frequent variants for Dominican patients were p.Phe508del and p.Ala559Thr and for Puerto Rican patients were p.Phe508del, p.Arg1066Cys, p.Arg334Trp, and p.I507del.ConclusionsIn this first description of the CFTR variants in patients with CF from the Dominican Republic and Puerto Rico, there was a low detection rate of two CFTR variants after full sequencing with the majority of patients from the Dominican Republic without identified variants.Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/153634/1/ppul24549.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/153634/2/ppul24549_am.pd

COMT Val 158 Met polymorphism is associated with post-traumatic stress disorder and functional outcome following mild traumatic brain injury

Mild traumatic brain injury (mTBI) results in variable clinical trajectories and outcomes. The source of variability remains unclear, but may involve genetic variations, such as single nucleotide polymorphisms (SNPs). A SNP in catechol-o-methyltransferase (COMT) is suggested to influence development of post-traumatic stress disorder (PTSD), but its role in TBI remains unclear. Here, we utilize the Transforming Research and Clinical Knowledge in Traumatic Brain Injury Pilot (TRACK-TBI Pilot) study to investigate whether the COMT Val158Met polymorphism is associated with PTSD and global functional outcome as measured by the PTSD Checklist - Civilian Version and Glasgow Outcome Scale Extended (GOSE), respectively. Results in 93 predominately Caucasian subjects with mTBI show that the COMT Met158 allele is associated with lower incidence of PTSD (univariate odds ratio (OR) of 0.25, 95% CI [0.09-0.69]) and higher GOSE scores (univariate OR 2.87, 95% CI [1.20-6.86]) 6-months following injury. The COMT Val158Met genotype and PTSD association persists after controlling for race (multivariable OR of 0.29, 95% CI [0.10-0.83]) and pre-existing psychiatric disorders/substance abuse (multivariable OR of 0.32, 95% CI [0.11-0.97]). PTSD emerged as a strong predictor of poorer outcome on GOSE (multivariable OR 0.09, 95% CI [0.03-0.26]), which persists after controlling for age, GCS, and race. When accounting for PTSD in multivariable analysis, the association of COMT genotype and GOSE did not remain significant (multivariable OR 1.73, 95% CI [0.69-4.35]). Whether COMT genotype indirectly influences global functional outcome through PTSD remains to be determined and larger studies in more diverse populations are needed to confirm these findings

Reconstructing native American migrations from whole-genome and whole-exome data.

There is great scientific and popular interest in understanding the genetic history of populations in the Americas. We wish to understand when different regions of the continent were inhabited, where settlers came from, and how current inhabitants relate genetically to earlier populations. Recent studies unraveled parts of the genetic history of the continent using genotyping arrays and uniparental markers. The 1000 Genomes Project provides a unique opportunity for improving our understanding of population genetic history by providing over a hundred sequenced low coverage genomes and exomes from Colombian (CLM), Mexican-American (MXL), and Puerto Rican (PUR) populations. Here, we explore the genomic contributions of African, European, and especially Native American ancestry to these populations. Estimated Native American ancestry is 48% in MXL, 25% in CLM, and 13% in PUR. Native American ancestry in PUR is most closely related to populations surrounding the Orinoco River basin, confirming the Southern American ancestry of the Taíno people of the Caribbean. We present new methods to estimate the allele frequencies in the Native American fraction of the populations, and model their distribution using a demographic model for three ancestral Native American populations. These ancestral populations likely split in close succession: the most likely scenario, based on a peopling of the Americas 16 thousand years ago (kya), supports that the MXL Ancestors split 12.2kya, with a subsequent split of the ancestors to CLM and PUR 11.7kya. The model also features effective populations of 62,000 in Mexico, 8,700 in Colombia, and 1,900 in Puerto Rico. Modeling Identity-by-descent (IBD) and ancestry tract length, we show that post-contact populations also differ markedly in their effective sizes and migration patterns, with Puerto Rico showing the smallest effective size and the earlier migration from Europe. Finally, we compare IBD and ancestry assignments to find evidence for relatedness among European founders to the three populations