Shifting the trajectory of therapeutic development for neurological and psychiatric disorders

Clinical trials for central nervous system disorders often enroll patients with unrecognized heterogeneous diseases, leading to costly trials that have high failure rates. Here, we discuss the potential of emerging technologies and datasets to elucidate disease mechanisms and identify biomarkers to improve patient stratification and monitoring of disease progression in clinical trials for neuropsychiatric disorders. Greater efforts must be centered on rigorously standardizing data collection and sharing of methods, datasets, and analytical tools across sectors. To address health care disparities in clinical trials, diversity of genetic ancestries and environmental exposures of research participants and associated biological samples must be prioritized.

Addressing the challenges of polygenic scores in human genetic research

The genotyping of millions of human samples has made it possible to evaluate variants across the human genome for their possible association with risks for numerous diseases and other traits by using genome-wide association studies (GWASs). The associations between phenotype and genotype found in GWASs make possible the construction of polygenic scores (PGSs), which aim to predict a trait or disease outcome in an individual on the basis of their genotype (in the disease case, the term polygenic risk score [PRS] is often used). PGSs have shown promise for studying the biology of complex traits and as a tool for evaluating individual disease risks in clinical settings. Although the quantity and quality of data to compute PGSs are increasing, challenges remain in the technical aspects of developing PGSs and in the ethical and social issues that might arise from their use. This ASHG Guidance emphasizes three major themes for researchers working with or interested in the application of PGSs in their own research: (1) developing diverse research cohorts; (2) fostering robustness in the development, application, and interpretation of PGSs; and (3) improving the communication of PGS results and their implications to broad audiences.

Impact of natural selection on global patterns of genetic variation, and association with clinical phenotypes, at genes involved in SARS-CoV-2 infection

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (ACE2, TMPRSS2, DPP4, and LY6E). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2, we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

NNT mediates redox-dependent pigmentation via a UVB- and MITF-independent mechanism

Ultraviolet (UV) light and incompletely understood genetic and epigenetic variations determine skin color. Here we describe an UV- and microphthalmia-associated transcription factor (MITF)-independent mechanism of skin pigmentation. Targeting the mitochondrial redox-regulating enzyme nicotinamide nucleotide transhydrogenase (NNT) resulted in cellular redox changes that affect tyrosinase degradation. These changes regulate melanosome maturation and, consequently, eumelanin levels and pigmentation. Topical application of small-molecule inhibitors yielded skin darkening in human skin, and mice with decreased NNT function displayed increased pigmentation. Additionally, genetic modification of NNT in zebrafish alters melanocytic pigmentation. Analysis of four diverse human cohorts revealed significant associations of skin color, tanning, and sun protection use with various single-nucleotide polymorphisms within NNT. NNT levels were independent of UVB irradiation and redox modulation. Individuals with postinflammatory hyperpigmentation or lentigines displayed decreased skin NNT levels, suggesting an NNT-driven, redox-dependent pigmentation mechanism that can be targeted with NNT-modifying topical drugs for medical and cosmetic purposes.

Impact of natural selection on global patterns of genetic variation, and association with clinical phenotypes, at genes involved in SARS-CoV-2 infection

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection ( ACE2, TMPRSS2, DPP4 , and LY6E ). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2 , we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2 , we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

Meta-analysis of GWA studies provides new insights on the genetic architecture of skin pigmentation in recently admixed populations

BACKGROUND

Association studies in recently admixed populations are extremely useful to identify the genetic architecture of pigmentation, due to their high genotypic and phenotypic variation. However, to date only four Genome-Wide Association Studies (GWAS) have been carried out in these populations.

RESULTS

We present a GWAS of skin pigmentation in an admixed sample from Cuba (N = 762). Additionally, we conducted a meta-analysis including the Cuban sample, and admixed samples from Cape Verde, Puerto Rico and African-Americans from San Francisco. This meta-analysis is one of the largest efforts so far to characterize the genetic basis of skin pigmentation in admixed populations (N = 2,104). We identified five genome-wide significant regions in the meta-analysis, and explored if the markers observed in these regions are associated with the expression of relevant pigmentary genes in human melanocyte cultures. In three of the regions identified in the meta-analysis (SLC24A5, SLC45A2, and GRM5/TYR), the association seems to be driven by non-synonymous variants (rs1426654, rs16891982, and rs1042602, respectively). The rs16891982 polymorphism is strongly associated with the expression of the SLC45A2 gene. In the GRM5/TYR region, in addition to the rs1042602 non-synonymous SNP located on the TYR gene, variants located in the nearby GRM5 gene have an independent effect on pigmentation, possibly through regulation of gene expression of the TYR gene. We also replicated an association recently described near the MFSD12 gene on chromosome 19 (lead variant rs112332856). Additionally, our analyses support the presence of multiple signals in the OCA2/HERC2/APBA2 region on chromosome 15. A clear causal candidate is the HERC2 intronic variant rs12913832, which has a profound influence on OCA2 expression. This variant has pleiotropic effects on eye, hair, and skin pigmentation. However, conditional and haplotype-based analyses indicate the presence of other variants with independent effects on melanin levels in OCA2 and APBA2. Finally, a follow-up of genome-wide signals identified in a recent GWAS for tanning response indicates that there is a substantial overlap in the genetic factors influencing skin pigmentation and tanning response.

CONCLUSIONS

Our meta-analysis of skin pigmentation GWAS in recently admixed populations provides new insights about the genetic architecture of this complex trait.

Phylogeny Estimation by Integration over Isolation with Migration Models

Phylogeny estimation is difficult for closely related populations and species, especially if they have been exchanging genes. We present a hierarchical Bayesian, Markov-chain Monte Carlo method with a state space that includes all possible phylogenies in a full Isolation-with-Migration model framework. The method is based on a new type of genealogy augmentation called a "hidden genealogy" that enables efficient updating of the phylogeny. This is the first likelihood-based method to fully incorporate directional gene flow and genetic drift for estimation of a species or population phylogeny. Application to human hunter-gatherer populations from Africa revealed a clear phylogenetic history, with strong support for gene exchange with an unsampled ghost population, and relatively ancient divergence between a ghost population and modern human populations, consistent with human/archaic divergence. In contrast, a study of five chimpanzee populations reveals a clear phylogeny with several pairs of populations having exchanged DNA, but does not support a history with an unsampled ghost population.

The molecular epidemiology of Huntington disease is related to intermediate allele frequency and haplotype in the general population

Huntington disease (HD) is the most common monogenic neurodegenerative disorder in populations of European ancestry, but occurs at lower prevalence in populations of East Asian or black African descent. New mutations for HD result from CAG repeat expansions of intermediate alleles (IAs), usually of paternal origin. The differing prevalence of HD may be related to the rate of new mutations in a population, but no comparative estimates of IA frequency or the HD new mutation rate are available. In this study, we characterize IA frequency and the CAG repeat distribution in fifteen populations of diverse ethnic origin. We estimate the HD new mutation rate in a series of populations using molecular IA expansion rates. The frequency of IAs was highest in Hispanic Americans and Northern Europeans, and lowest in black Africans and East Asians. The prevalence of HD correlated with the frequency of IAs by population and with the proportion of IAs found on the HD-associated A1 haplotype. The HD new mutation rate was estimated to be highest in populations with the highest frequency of IAs. In European ancestry populations, one in 5,372 individuals from the general population and 7.1% of individuals with an expanded CAG repeat in the HD range are estimated to have a molecular new mutation. Our data suggest that the new mutation rate for HD varies substantially between populations, and that IA frequency and haplotype are closely linked to observed epidemiological differences in the prevalence of HD across major ancestry groups in different countries.

Abstract A33: Population and evolutionary genomics of prostate cancer-associated variants: implications for health disparities in men of African descent

To determine why African men are more likely to suffer from prostate cancer (CaP), we integrated GWAS results and scans of selection with allele frequency data from 64 global populations. Despite substantial overlap in genetic risk scores across populations, we find that predicted CaP risk is highest in West Africans and that a small number of loci drive these differences in risk. There is a strong concordance between genetic risk scores and clinical estimates of CaP mortality. Although most CaP-associated loci are evolving neutrally, we find multiple instances where alleles have hitchhiked to high frequencies with linked locally adaptive alleles.

Citation Format: Joseph Lachance, Ali Berens, Matthew Hansen, Andrew Teng, Sarah Tishkoff, Timothy Rebbeck. Population and evolutionary genomics of prostate cancer-associated variants: implications for health disparities in men of African descent [abstract]. In: Proceedings of the AACR International Conference: New Frontiers in Cancer Research; 2017 Jan 18-22; Cape Town, South Africa. Philadelphia (PA): AACR; Cancer Res 2017;77(22 Suppl):Abstract nr A33.

Tracing the peopling of the world through genomics

Advances in the sequencing and the analysis of the genomes of both modern and ancient peoples have facilitated a number of breakthroughs in our understanding of human evolutionary history. These include the discovery of interbreeding between anatomically modern humans and extinct hominins; the development of an increasingly detailed description of the complex dispersal of modern humans out of Africa and their population expansion worldwide; and the characterization of many of the genetic adaptions of humans to local environmental conditions. Our interpretation of the evolutionary history and adaptation of humans is being transformed by analyses of these new genomic data.

A recent bottleneck of Y chromosome diversity coincides with a global change in culture

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50–100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192–307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47–52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

DCAF4, a novel gene associated with leucocyte telomere length

BACKGROUND

Leucocyte telomere length (LTL), which is fashioned by multiple genes, has been linked to a host of human diseases, including sporadic melanoma. A number of genes associated with LTL have already been identified through genome-wide association studies. The main aim of this study was to establish whether DCAF4 (DDB1 and CUL4-associated factor 4) is associated with LTL. In addition, using ingenuity pathway analysis (IPA), we examined whether LTL-associated genes in the general population might partially explain the inherently longer LTL in patients with sporadic melanoma, the risk for which is increased with ultraviolet radiation (UVR).

RESULTS

Genome-wide association (GWA) meta-analysis and de novo genotyping of 20 022 individuals revealed a novel association (p=6.4×10(-10)) between LTL and rs2535913, which lies within DCAF4. Notably, eQTL analysis showed that rs2535913 is associated with decline in DCAF4 expressions in both lymphoblastoid cells and sun-exposed skin (p=4.1×10(-3) and 2×10(-3), respectively). Moreover, IPA revealed that LTL-associated genes, derived from GWA meta-analysis (N=9190), are over-represented among genes engaged in melanoma pathways. Meeting increasingly stringent p value thresholds (p<0.05, <0.01, <0.005, <0.001) in the LTL-GWA meta-analysis, these genes were jointly over-represented for melanoma at p values ranging from 1.97×10(-169) to 3.42×10(-24).

CONCLUSIONS

We uncovered a new locus associated with LTL in the general population. We also provided preliminary findings that suggest a link of LTL through genetic mechanisms with UVR and melanoma in the general population.

Sequence diversity of Pan troglodytes subspecies and the impact of WFDC6 selective constraints in reproductive immunity

Recent efforts have attempted to describe the population structure of common chimpanzee, focusing on four subspecies: Pan troglodytes verus, P. t. ellioti, P. t. troglodytes, and P. t. schweinfurthii. However, few studies have pursued the effects of natural selection in shaping their response to pathogens and reproduction. Whey acidic protein (WAP) four-disulfide core domain (WFDC) genes and neighboring semenogelin (SEMG) genes encode proteins with combined roles in immunity and fertility. They display a strikingly high rate of amino acid replacement (d_N/d_S), indicative of adaptive pressures during primate evolution. In human populations, three signals of selection at the WFDC locus were described, possibly influencing the proteolytic profile and antimicrobial activities of the male reproductive tract. To evaluate the patterns of genomic variation and selection at the WFDC locus in chimpanzees, we sequenced 17 WFDC genes and 47 autosomal pseudogenes in 68 chimpanzees (15 P. t. troglodytes, 22 P. t. verus, and 31 P. t. ellioti). We found a clear differentiation of P. t. verus and estimated the divergence of P. t. troglodytes and P. t. ellioti subspecies in 0.173 Myr; further, at the WFDC locus we identified a signature of strong selective constraints common to the three subspecies in WFDC6—a recent paralog of the epididymal protease inhibitor EPPIN. Overall, chimpanzees and humans do not display similar footprints of selection across the WFDC locus, possibly due to different selective pressures between the two species related to immune response and reproductive biology.

Genome-wide meta-analysis points to CTC1 and ZNF676 as genes regulating telomere homeostasis in humans

Leukocyte telomere length (LTL) is associated with a number of common age-related diseases and is a heritable trait. Previous genome-wide association studies (GWASs) identified two loci on chromosomes 3q26.2 (TERC) and 10q24.33 (OBFC1) that are associated with the inter-individual LTL variation. We performed a meta-analysis of 9190 individuals from six independent GWAS and validated our findings in 2226 individuals from four additional studies. We confirmed previously reported associations with OBFC1 (rs9419958 P = 9.1 × 10⁻¹¹) and with the telomerase RNA component TERC (rs1317082, P = 1.1 × 10⁻⁸). We also identified two novel genomic regions associated with LTL variation that map near a conserved telomere maintenance complex component 1 (CTC1; rs3027234, P = 3.6 × 10⁻⁸) on chromosome17p13.1 and zinc finger protein 676 (ZNF676; rs412658, P = 3.3 × 10⁻⁸) on 19p12. The minor allele of rs3027234 was associated with both shorter LTL and lower expression of CTC1. Our findings are consistent with the recent observations that point mutations in CTC1 cause short telomeres in both Arabidopsis and humans affected by a rare Mendelian syndrome. Overall, our results provide novel insights into the genetic architecture of inter-individual LTL variation in the general population.

Global distribution of a novel trinucleotide microsatellite polymorphism (ATA)n in intron 8 of the SLC11A1 gene and susceptibility to pulmonary tuberculosis

We identified a novel trinucleotide (ATA)n repeat polymorphism in intron 8 of SLC11AI, a candidate gene for susceptibility to tuberculosis (TB) infection. We characterized the frequency of this polymorphism in 485 individuals originating from eight globally diverse human populations and compared the distribution of (ATA)n alleles in 146 adults and in 80 cord blood samples from newborns in the Gambian population. Lastly, we tested for association of this microsatellite with pulmonary TB in 318 TB cases and 146 controls in the Gambian population. We found no significant difference in frequency or distribution of alleles in adult and cord blood samples, and we found no significant association between this marker and pulmonary TB.

Estimating European admixture in African Americans by using microsatellites and a microsatellite haplotype (CD4/Alu)

We have analyzed 10 unlinked microsatellites and a linked Alu deletion polymorphism at the CD4 locus in an African American population sample from Chicago (USA). Heterozygosity estimates at the microsatellite loci range from 0.727+/-0.025 (D3S1358) to 0.873+/-0.017 (D18S51), with an average of 0.794+/-0.016. These values are comparable to or higher than those reported for Europeans, with only one exception (D3S1358). The CD4/Alu haplotypic diversity (0.887+/-0.012) is comparable to diversity levels observed in sub-Saharan African populations and is higher than the diversity levels reported in European populations. No consistent pattern of within, between, or multi-locus deviations from Hardy-Weinberg expectations is observed, suggesting a low sub-heterogeneity within the sampled population. We have applied a maximum likelihood method and estimated the proportion of European admixture to the African American gene pool to be 0.26+/-0.02. The narrow confidence interval indicates that allele frequency data from multiple microsatellite loci, whether analyzed independently or as haplotypes, are particularly useful for estimating genetic admixture.

Genetic structure of the ancestral population of modern humans

Neutral DNA polymorphisms from an 8-kb segment of the dystrophin gene, previously ascertained in a worldwide sample (n = 250 chromosomes), were used to characterize the population ancestral to the present-day human groups. The ancestral state of each polymorphic site was determined by comparing human variants with their orthologous sites in the great apes. The "age before fixation" of the underlying mutations was estimated from the frequencies of the new alleles and analyzed in the context of these polymorphisms' distribution among 13 populations from Africa, Europe, Asia, New Guinea, and the Americas (n = 860 chromosomes in total). Seventeen polymorphisms older tan 100,000-200,000 years, which contributed approximately 90% to the overall nucleotide diversity, were common to all human groups. Polymorphisms endemic to human groups or continentally restricted were younger than 100,000-200,000 years. Africans (six populations) with 13 such sites stood out from the rest of the world (seven populations), where only 2 population-specific variants were observed. The similarity of the frequencies of the old polymorphisms in Africans and non-Africans suggested a similar profile of genetic variability in the population before the modern human's divergence. This ancestral population was characterized by an effective size of about 10,000 as estimated from the nucleotide diversity; this size may describe the number of breeding individuals over a long time during the Middle Pleistocene or reflect a speciation bottleneck from an initially larger population at the end of this period.

Nuclear DNA diversity in worldwide distributed human populations

Nucleotide variation was examined in an 8 kb intronic DNA bordering exon 44 of the human dystrophin gene on Xp21. Thirty-six polymorphisms (substitutions, small insertions/deletions and one (T)n microsatellite) were found using SSCP/heteroduplex analysis of DNA samples from mixed Europeans, Papua New Guineans as well as from six African, three Asian and two Amerindian populations. In this way the European bias in the nuclear polymorphism ascertainment has been avoided. In a maximum likelihood tree constructed from the frequency data, Africans clustered separately from the non-African populations. Fifteen polymorphisms were shared among most of the populations compared, whereas 13 sites were found to be endemic to Africans and four to non-Africans. The common sites contributed most to the average heterozygosity (Hn=0.101%+/-0.023), whereas the endemic ones, being rare, had little effect on this estimate. The F(ST) values were lower for Africans (0.072) than for non-Africans (0.158), suggesting a higher level of gene exchange within Africa, corroborating the observation of a greater number of segregating sites on this continent than elsewhere. The data suggest a recent common origin of the African and non-African populations, where a greater geographical isolation of the latter resulted in a smaller number of newly acquired polymorphisms.

Detection of DNA sequence variation via deliberate heteroduplex formation from genomic DNAs amplified en masse in "population tubes"

We have developed the population tube (poptube) system for sensitive detection and large-scale sampling of DNA sequence variation in several human populations of wide geographic distribution. In this methodology, genomic DNAs from five individuals in a population are PCR amplified en masse to maximize deliberately the chances of forming heteroduplexes among allelic variants. Interpopulation mixing is performed in a separate set of tubes containing one individual from each of five populations as well as a reference chimpanzee sample deliberately chosen to be different from all humans. Mismatches at sites of allelic variation retard the electrophoresis and reduce the stability of heteroduplex molecules. The products are electrophoresed on denaturing gradient gels where detection of heteroduplexes is accomplished readily. Using poptubes, we have discovered a rare variant in an otherwise highly conserved 440-bp segment in the long intron of the glucose-6-phosphate dehydrogenase (G6PD) gene. The polymorphism at this X-chromosome locus could be only detected in males by mixing samples, as homoduplexes for both alleles co-focus on denaturing gradient electrophoresis.