Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Recombination between heterologous human acrocentric chromosomes

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome-the Telomere-to-Telomere Consortium's CHM13 assembly (T2T-CHM13)-provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9.

Human-SARS-CoV-2 interactome and human genetic diversity: TMPRSS2-rs2070788, associated with severe influenza, and its population genetics caveats in Native Americans

For human/SARS-CoV-2 interactome genes ACE2, TMPRSS2 and BSG, there is a convincing evidence of association in Asians with influenza-induced SARS for TMPRSS2-rs2070788, tag-SNP of the eQTL rs383510. This case illustrates the importance of population genetics and of sequencing data in the design of genetic association studies in different human populations: the high linkage disequilibrium (LD) between rs2070788 and rs383510 is Asian-specific. Leveraging on a combination of genotyping and sequencing data for Native Americans (neglected in genetic studies), we show that while their frequencies of the Asian tag-SNP rs2070788 is, surprisingly, the highest worldwide, it is not in LD with the eQTL rs383510, that therefore, should be directly genotyped in genetic association studies of SARS in populations with Native American ancestry.

Adaptive selection drives TRPP3 loss-of-function in an Ethiopian population

TRPP3 (also called PKD2L1) is a nonselective, cation-permeable channel activated by multiple stimuli, including extracellular pH changes. TRPP3 had been considered a candidate for sour sensor in humans, due to its high expression in a subset of tongue receptor cells detecting sour, along with its membership to the TRP channel family known to function as sensory receptors. Here, we describe the functional consequences of two non-synonymous genetic variants (R278Q and R378W) found to be under strong positive selection in an Ethiopian population, the Gumuz. Electrophysiological studies and 3D modelling reveal TRPP3 loss-of-functions produced by both substitutions. R278Q impairs TRPP3 activation after alkalinisation by mislocation of H⁺ binding residues at the extracellular polycystin mucolipin domain. R378W dramatically reduces channel activity by altering conformation of the voltage sensor domain and hampering channel transition from closed to open state. Sour sensitivity tests in R278Q/R378W carriers argue against both any involvement of TRPP3 in sour detection and the role of such physiological process in the reported evolutionary positive selection past event.

Linkage-based genome assembly improvement of oil palm (Elaeis guineensis)

Meiotic crossovers in outbred species, such as oil palm (Elaeis guineensis Jacq., 2n = 32) contribute to allelic re-assortment in the genome. Such genetic variation is usually exploited in breeding to combine positive alleles for trait superiority. A good quality reference genome is essential for identifying the genetic factors underlying traits of interest through linkage or association studies. At the moment, an AVROS pisifera genome is publicly available for oil palm. Distribution and frequency of crossovers throughout chromosomes in different origins of oil palm are still unclear. Hence, an ultrahigh-density genomic linkage map of a commercial Deli dura x AVROS pisifera family was constructed using the OP200K SNP array, to evaluate the genetic alignment with the genome assembly. A total of 27,890 linked SNP markers generated a total map length of 1,151.7 cM and an average mapping interval of 0.04 cM. Nineteen linkage groups represented 16 pseudo-chromosomes of oil palm, with 61.7% of the mapped SNPs present in the published genome. Meanwhile, the physical map was also successfully extended from 658 Mb to 969 Mb by assigning unplaced scaffolds to the pseudo-chromosomes. A genic linkage map with major representation of sugar and lipid biosynthesis pathways was subsequently built for future studies on oil related quantitative trait loci (QTL). This study improves the current physical genome of the commercial oil palm, and provides important insights into its recombination landscape, eventually unlocking the full potential genome sequence-enabled biology for oil palm.

A neutrality test for detecting selection on DNA methylation using single methylation polymorphism frequency spectrum

Inheritable epigenetic mutations (epimutations) can contribute to transmittable phenotypic variation. Thus, epimutations can be subject to natural selection and impact the fitness and evolution of organisms. Based on the framework of the modified Tajima’s D test for DNA mutations, we developed a neutrality test with the statistic “D^m” to detect selection forces on DNA methylation mutations using single methylation polymorphisms. With computer simulation and empirical data analysis, we compared the D^m test with the original and modified Tajima’s D tests and demonstrated that the D^m test is suitable for detecting selection on epimutations and outperforms original/modified Tajima’s D tests. Due to the higher resetting rate of epimutations, the interpretation of D^m on epimutations and Tajima’s D test on DNA mutations could be different in inferring natural selection. Analyses using simulated and empirical genome-wide polymorphism data suggested that genes under genetic and epigenetic selections behaved differently. We applied the D^m test to recently originated Arabidopsis and human genes, and showed that newly evolved genes contain higher level of rare epialleles, suggesting that epimutation may play a role in origination and evolution of genes and genomes. Overall, we demonstrate the utility of the D^m test to detect whether the loci are under selection regarding DNA methylation. Our analytical metrics and methodology could contribute to our understanding of evolutionary processes of genes and genomes in the field of epigenetics. The Perl script for the “D^m” test is available at http://fanlab.wayne.edu/ (last accessed December 18, 2014).

Population and genomic lessons from genetic analysis of two Indian populations

Indian demographic history includes special features such as founder effects, interpopulation segregation, complex social structure with a caste system and elevated frequency of consanguineous marriages. It also presents a higher frequency for some rare mendelian disorders and in the last two decades increased prevalence of some complex disorders. Despite the fact that India represents about one-sixth of the human population, deep genetic studies from this terrain have been scarce. In this study, we analyzed high-density genotyping and whole-exome sequencing data of a North and a South Indian population. Indian populations show higher differentiation levels than those reported between populations of other continents. In this work, we have analyzed its consequences, by specifically assessing the transferability of genetic markers from or to Indian populations. We show that there is limited genetic marker portability from available genetic resources such as HapMap or the 1,000 Genomes Project to Indian populations, which also present an excess of private rare variants. Conversely, tagSNPs show a high level of portability between the two Indian populations, in contrast to the common belief that North and South Indian populations are genetically very different. By estimating kinship from mates and consanguinity in our data from trios, we also describe different patterns of assortative mating and inbreeding in the two populations, in agreement with distinct mating preferences and social structures. In addition, this analysis has allowed us to describe genomic regions under recent adaptive selection, indicating differential adaptive histories for North and South Indian populations. Our findings highlight the importance of considering demography for design and analysis of genetic studies, as well as the need for extending human genetic variation catalogs to new populations and particularly to those with particular demographic histories.

Worldwide population distribution of the common LCE3C-LCE3B deletion associated with psoriasis and other autoimmune disorders

Background

There is increasing evidence of the importance of copy number variants (CNV) in genetic diversity among individuals and populations, as well as in some common genetic diseases. We previously characterized a common 32-kb insertion/deletion variant of the PSORS4 locus at chromosome 1q21 that harbours the LCE3C and LCE3B genes. This variant allele (LCE3C_LCE3B-del) is common in patients with psoriasis and other autoimmune disorders from certain ethnic groups.

Results

Using array-CGH (Agilent 244 K) in samples from the HapMap and Human Genome Diversity Panel (HGDP) collections, we identified 54 regions showing population differences in comparison to Africans. We provided here a comprehensive population-genetic analysis of one of these regions, which involves the 32-kb deletion of the PSORS4 locus. By a PCR-based genotyping assay we characterised the profiles of the LCE3C_LCE3B-del and the linkage disequilibrium (LD) pattern between the variant allele and the tag SNP rs4112788. Our results show that most populations tend to have a higher frequency of the deleted allele than Sub-Saharan Africans. Furthermore, we found strong LD between rs4112788G and LCE3C_LCE3B-del in most non-African populations (r2 >0.8), in contrast to the low concordance between loci (r2 <0.3) in the African populations.

Conclusions

These results are another example of population variability in terms of biomedical interesting CNV. The frequency distribution of the LCE3C_LCE3B-del allele and the LD pattern across populations suggest that the differences between ethnic groups might not be due to natural selection, but the consequence of genetic drift caused by the strong bottleneck that occurred during “out of Africa” expansion.

Detecting genetic isolation in human populations: a study of European language minorities

The identification of isolation signatures is fundamental to better understand the genetic structure of human populations and to test the relations between cultural factors and genetic variation. However, with current approaches, it is not possible to distinguish between the consequences of long-term isolation and the effects of reduced sample size, selection and differential gene flow. To overcome these limitations, we have integrated the analysis of classical genetic diversity measures with a Bayesian method to estimate gene flow and have carried out simulations based on the coalescent. Combining these approaches, we first tested whether the relatively short history of cultural and geographical isolation of four “linguistic islands” of the Eastern Alps (Lessinia, Sauris, Sappada and Timau) had left detectable signatures in their genetic structure. We then compared our findings to previous studies of European population isolates. Finally, we explored the importance of demographic and cultural factors in shaping genetic diversity among the groups under study. A combination of small initial effective size and continued genetic isolation from surrounding populations seems to provide a coherent explanation for the diversity observed among Sauris, Sappada and Timau, which was found to be substantially greater than in other groups of European isolated populations. Simulations of micro-evolutionary scenarios indicate that ethnicity might have been important in increasing genetic diversity among these culturally related and spatially close populations.

Perspectives on human population structure at the cusp of the sequencing era

Human groups show structured levels of genetic similarity as a consequence of factors such as geographical subdivision and genetic drift. Surveying this structure gives us a scientific perspective on human origins, sheds light on evolutionary processes that shape both human adaptation and disease, and is integral to effectively carrying out the mission of global medical genetics and personalized medicine. Surveys of population structure have been ongoing for decades, but in the past three years, single-nucleotide-polymorphism (SNP) array technology has provided unprecedented detail on human population structure at global and regional scales. These studies have confirmed well-known relationships between distantly related populations and uncovered previously unresolvable relationships among closely related human groups. SNPs represent the first dense genome-wide markers, and as such, their analysis has raised many challenges and insights relevant to the study of population genetics with whole-genome sequences. Here we draw on the lessons from these studies to anticipate the directions that will be most fruitful to pursue during the emerging whole-genome sequencing era.

Population-genetic comparison of the Sorbian isolate population in Germany with the German KORA population using genome-wide SNP arrays

Background

The Sorbs are an ethnic minority in Germany with putative genetic isolation, making the population interesting for disease mapping. A sample of N = 977 Sorbs is currently analysed in several genome-wide meta-analyses. Since genetic differences between populations are a major confounding factor in genetic meta-analyses, we compare the Sorbs with the German outbred population of the KORA F3 study (N = 1644) and other publically available European HapMap populations by population genetic means. We also aim to separate effects of over-sampling of families in the Sorbs sample from effects of genetic isolation and compare the power of genetic association studies between the samples.

Results

The degree of relatedness was significantly higher in the Sorbs. Principal components analysis revealed a west to east clustering of KORA individuals born in Germany, KORA individuals born in Poland or Czech Republic, Half-Sorbs (less than four Sorbian grandparents) and Full-Sorbs. The Sorbs cluster is nearest to the cluster of KORA individuals born in Poland. The number of rare SNPs is significantly higher in the Sorbs sample. FST between KORA and Sorbs is an order of magnitude higher than between different regions in Germany. Compared to the other populations, Sorbs show a higher proportion of individuals with runs of homozygosity between 2.5 Mb and 5 Mb. Linkage disequilibrium (LD) at longer range is also slightly increased but this has no effect on the power of association studies.

Oversampling of families in the Sorbs sample causes detectable bias regarding higher FST values and higher LD but the effect is an order of magnitude smaller than the observed differences between KORA and Sorbs. Relatedness in the Sorbs also influenced the power of uncorrected association analyses.

Conclusions

Sorbs show signs of genetic isolation which cannot be explained by over-sampling of relatives, but the effects are moderate in size. The Slavonic origin of the Sorbs is still genetically detectable.

Regarding LD structure, a clear advantage for genome-wide association studies cannot be deduced. The significant amount of cryptic relatedness in the Sorbs sample results in inflated variances of Beta-estimators which should be considered in genetic association analyses.

Genetic diversity analysis of abiotic stress response gene TaSnRK2.7-A in common wheat

Sucrose non-fermenting1-related protein kinase 2 (SnRK2) plays a key role in plant stress signaling transduction pathways. In this study, one copy of TaSnRK2.7, a SnRK2 member of common wheat, was isolated and characterized for nucleotide diversity among 45 wheat accessions with different stress-response features. Most of the accessions were elite wheat cultivars, which had been subject to population bottlenecks and intensive selection during breeding. Nucleotide and haplotype diversity across the entire TaSnRK2.7-A region was 0.00076 and 0.590, respectively, and diversity in non-coding regions was higher than that in coding regions. Sliding-window analysis showed variable levels of nucleotide variation along the entire TaSnRK2.7-A region; the sixth intron and ninth exon represented variation-enriched regions. As predicted, neutrality tests revealed that population bottlenecks or purifying selection had acted on the TaSnRK2.7-A gene, a relatively conserved gene. Furthermore, strong linkage disequilibrium between SNP loci extends across the entire TaSnRK2.7-A region. These findings demonstrate that the TaSnRK2.7-A genomic region has evolved under extensive selection pressure during crop breeding.

Genetic variation in the Sorbs of eastern Germany in the context of broader European genetic diversity

Population isolates have long been of interest to genetic epidemiologists because of their potential to increase power to detect disease-causing genetic variants. The Sorbs of Germany are considered as cultural and linguistic isolates and have recently been the focus of disease association mapping efforts. They are thought to have settled in their present location in eastern Germany after a westward migration from a largely Slavic-speaking territory during the Middle Ages. To examine Sorbian genetic diversity within the context of other European populations, we analyzed genotype data for over 30 000 autosomal single-nucleotide polymorphisms from over 200 Sorbs individuals. We compare the Sorbs with other European individuals, including samples from population isolates. Despite their geographical proximity to German speakers, the Sorbs showed greatest genetic similarity to Polish and Czech individuals, consistent with the linguistic proximity of Sorbian to other West Slavic languages. The Sorbs also showed evidence of subtle levels of genetic isolation in comparison with samples from non-isolated European populations. The level of genetic isolation was less than that observed for the Sardinians and French Basque, who were clear outliers on multiple measures of isolation. The finding of the Sorbs as only a minor genetic isolate demonstrates the need to genetically characterize putative population isolates, as they possess a wide range of levels of isolation because of their different demographic histories.

Brief communication: patterns of linkage disequilibrium and haplotype diversity at Xq13 in six Native American populations

Comparative studies of linkage disequilibrium (LD) can provide insights into human demographic history. Here, we characterize LD in six Native American populations using seven microsatellite markers in Xq13, a region of the genome extensively studied in populations around the world. Native Americans show relatively low diversity and high LD, in agreement with recent genome-wide survey and a scenario of sequential founder effects accompanying human population dispersal around the globe. LD in Native Americans is similar to that observed in some recently described small population isolates and higher than in large European isolates (e.g., Finns), which have been extensively analyzed in medical genetics studies. Haplotype analyses are consistent with a colonization of the New World by a differentiated East Asian population, followed by extensive genetic drift in the Americas.