Population-specific common SNPs reflect demographic histories and highlight regions of genomic plasticity with functional relevance

A novel locus in CSMD1 gene is associated with increased susceptibility to severe malaria in Malian children

Background

Plasmodium falciparum malaria is still a leading cause of child mortality in sub-Saharan Africa. The clinical manifestations of malaria range from asymptomatic infection to severe disease. The variation in clinical presentation is partly attributed to host genetic factors with estimated narrow-sense heritability of 23%. Here, we investigate the associations between candidate gene polymorphisms and the likelihood of severe malaria (SM) in a cohort of Malian children.

Methods

Based on our previous genome-wide association studies (GWAS) analysis, candidate genes were selected for in-depth analysis using several criteria including gene-level GWAS scores, functional overlap with malaria pathogenesis, and evidence of association with protection or susceptibility to other infectious or inflammatory diseases. Single Nucleotide Polymorphisms (SNPs) residing within these genes were selected mainly based on p-values from previous severe malaria susceptibility GWAS studies and minor allele frequency (MAF) in West African populations.

Results

Of 182 candidate genes reported in our previous study, 11 genes and 22 SNPs residing in these genes were selected. The selected SNPs were genotyped using KASP technology in 477 DNA samples (87 SM and 390 controls). Logistic regression analysis revealed that a common intron variant, rs13340578 in CUB and Sushi Multi Domain (CSMD1) gene, is associated with increased odds of SM in recessive mode of inheritance (MAF = 0.42, OR = 1.8, 95% CI = [1.78, 1.84], p = 0.029). The SNP is in linkage disequilibrium (LD) with multiple variants with regulatory features.

Conclusion

Taken together, the current study showed that an intron variant rs13340578, residing in CSMD1 gene, is associated with increased susceptibility to malaria. This finding suggests that modified regulation of complement may contribute to malaria disease severity. Further studies are needed to identify the causal variants and the underlying molecular mechanisms.

PAPipe: A Pipeline for Comprehensive Population Genetic Analysis

Advancements in next-generation sequencing (NGS) technologies have led to a substantial increase in the availability of population genetic variant data, thus prompting the development of various population analysis tools to enhance our understanding of population structure and evolution. The tools that are currently used to analyze population genetic variant data generally require different environments, parameters, and formats of the input data, which can act as a barrier preventing the wide-spread usage of such tools by general researchers who may not be familiar with bioinformatics. To address this problem, we have developed an automated and comprehensive pipeline called PAPipe to perform nine widely used population genetic analyses using population NGS data. PAPipe seamlessly interconnects and serializes multiple steps, such as read trimming and mapping, genetic variant calling, data filtering, and format converting, along with nine population genetic analyses such as principal component analysis, phylogenetic analysis, population tree analysis, population structure analysis, linkage disequilibrium decay analysis, selective sweep analysis, population admixture analysis, sequentially Markovian coalescent analysis, and fixation index analysis. PAPipe also provides an easy-to-use web interface that allows for the parameters to be set and the analysis results to be browsed in intuitive manner. PAPipe can be used to generate extensive results that provide insights that can help enhance user convenience and data usability. PAPipe is freely available at https://github.com/jkimlab/PAPipe.

A combination of conserved and diverged responses underlies Theobroma cacao's defense response to Phytophthora palmivora

BACKGROUND

Plants have complex and dynamic immune systems that have evolved to resist pathogens. Humans have worked to enhance these defenses in crops through breeding. However, many crops harbor only a fraction of the genetic diversity present in wild relatives. Increased utilization of diverse germplasm to search for desirable traits, such as disease resistance, is therefore a valuable step towards breeding crops that are adapted to both current and emerging threats. Here, we examine diversity of defense responses across four populations of the long-generation tree crop Theobroma cacao L., as well as four non-cacao Theobroma species, with the goal of identifying genetic elements essential for protection against the oomycete pathogen Phytophthora palmivora.

RESULTS

We began by creating a new, highly contiguous genome assembly for the P. palmivora-resistant genotype SCA 6 (Additional file 1: Tables S1-S5), deposited in GenBank under accessions CP139290-CP139299. We then used this high-quality assembly to combine RNA and whole-genome sequencing data to discover several genes and pathways associated with resistance. Many of these are unique, i.e., differentially regulated in only one of the four populations (diverged 40 k-900 k generations). Among the pathways shared across all populations is phenylpropanoid biosynthesis, a metabolic pathway with well-documented roles in plant defense. One gene in this pathway, caffeoyl shikimate esterase (CSE), was upregulated across all four populations following pathogen treatment, indicating its broad importance for cacao's defense response. Further experimental evidence suggests this gene hydrolyzes caffeoyl shikimate to create caffeic acid, an antimicrobial compound and known inhibitor of Phytophthora spp.

CONCLUSIONS

Our results indicate most expression variation associated with resistance is unique to populations. Moreover, our findings demonstrate the value of using a broad sample of evolutionarily diverged populations for revealing the genetic bases of cacao resistance to P. palmivora. This approach has promise for further revealing and harnessing valuable genetic resources in this and other long-generation plants.

In-silico mining to glean SNPs of pharmaco-clinical importance: an investigation with reference to the Indian populated SNPs

Drugs pharmacology is defined by pharmacokinetics and pharmacodynamics and both of them are affected by genetic variability. Genetic variability varies from population to population, and sometimes even within the population, it exists. Single nucleotide polymorphisms (SNPs) are one of the major genetic variability factors which are found to be associated with the pharmacokinetics and pharmacodynamics process of a drug and are responsible for variable drug response and clinical phenotypes. Studies of SNPs can help to perform genome-wide association studies for their association with pharmacological and clinical events, at the same time; their information can direct genome-wide association studies for their use as biomarkers. With the aim to mine and characterize Indian populated SNPs of pharmacological and clinical importance. Two hundred six candidate SNPs belonging to 43 genes were retrieved from Indian Genome Variation Database. The distribution pattern of considered SNPs was observed against all five world super-populations (AFR, AMR, EAS, EUR, and SAS). Further, their annotation was done through SNP-nexus by considering Human genome reference builds - hg38, pharmacological and clinical information was supplemented by PharmGKB and ClinVar database. At last, to find out the association between SNPs linkage disequilibrium was observed in terms of r2. Overall, the study reported 53 pharmaco-clinical active SNPs and found 24 SNP-pairs as potential markers, and recommended their clinical and experimental validation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-023-00154-4.

A stimulus-contingent positive feedback loop enables IFN-β dose-dependent activation of pro-inflammatory genes

Type I interferons (IFN) induce powerful antiviral and innate immune responses via the transcription factor, IFN-stimulated gene factor (ISGF3). However, in some pathological contexts, type I IFNs are responsible for exacerbating inflammation. Here, we show that a high dose of IFN-β also activates an inflammatory gene expression program in contrast to IFN-λ3, a type III IFN, which elicits only the common antiviral gene program. We show that the inflammatory gene program depends on a second, potentiated phase in ISGF3 activation. Iterating between mathematical modeling and experimental analysis, we show that the ISGF3 activation network may engage a positive feedback loop with its subunits IRF9 and STAT2. This network motif mediates stimulus-specific ISGF3 dynamics that are dependent on ligand, dose, and duration of exposure, and when engaged activates the inflammatory gene expression program. Our results reveal a previously underappreciated dynamical control of the JAK-STAT/IRF signaling network that may produce distinct biological responses and suggest that studies of type I IFN dysregulation, and in turn therapeutic remedies, may focus on feedback regulators within it.

Quantitative trait locus (xQTL) approaches identify risk genes and drug targets from human non-coding genomes

Advances and reduction of costs in various sequencing technologies allow for a closer look at variations present in the non-coding regions of the human genome. Correlating non-coding variants with large-scale multi-omic data holds the promise not only of a better understanding of likely causal connections between non-coding DNA and expression of traits but also identifying potential disease-modifying medicines. Genome-phenome association studies have created large datasets of DNA variants that are associated with multiple traits or diseases, such as Alzheimer's disease; yet, the functional consequences of variants, in particular of non-coding variants, remain largely unknown. Recent advances in functional genomics and computational approaches have led to the identification of potential roles of DNA variants, such as various quantitative trait locus (xQTL) techniques. Multi-omic assays and analytic approaches toward xQTL have identified links between genetic loci and human transcriptomic, epigenomic, proteomic and metabolomic data. In this review, we first discuss the recent development of xQTL from multi-omic findings. We then highlight multimodal analysis of xQTL and genetic data for identification of risk genes and drug targets using Alzheimer's disease as an example. We finally discuss challenges and future research directions (e.g. artificial intelligence) for annotation of non-coding variants in complex diseases.

Analysis of Common SNPs across Continents Reveals Major Genomic Differences between Human Populations

Common alleles tend to be more ancient than rare alleles. These common SNPs appeared thousands of years ago and reflect intricate human evolution including various adaptations, admixtures, and migration events. Eighty-four thousand abundant region-specific alleles (ARSAs) that are common in one continent but absent in the rest of the world have been characterized by processing 3100 genomes from 230 populations. Also computed were 17,446 polymorphic sites with regional absence of common alleles (RACAs), which are widespread globally but absent in one region. A majority of these region-specific SNPs were found in Africa. America has the second greatest number of ARSAs (3348) and is even ahead of Europe (1911). Surprisingly, East Asia has the highest number of RACAs (10,524) and the lowest number of ARSAs (362). ARSAs and RACAs have distinct compositions of ancestral versus derived alleles in different geographical regions, reflecting their unique evolution. Genes associated with ARSA and RACA SNPs were identified and their functions were analyzed. The core 100 genes shared by multiple populations and associated with region-specific natural selection were examined. The largest part of them (42%) are related to the nervous system. ARSA and RACA SNPs are important for both association and human evolution studies.

Association of the 5HTR2C gene Ser23 variation with childhood allergic asthma

Objective: Allergic asthma is the most frequently observed subtype of individuals with asthma. The effects of serotonin plays in the pathophysiology of asthma has not been clearly determined. Thus, this study aimed to investigate the association between the 5HTR2C gene rs 6318 G>C polymorphism and allergic asthma in pediatric patients in Cyprus.

Methods: This study included total number of 177 individuals with 118 control and 59 pediatric patients (43 atopic and 16 non-atopic asthma patient). A skin prick test was performed for each patient to confirm asthma diagnosis and to evaluate atopic status. Genotyping for the 5HTR2C was completed by Real Time-PCR analysis.

Results: The genotype distribution frequencies were not in agreement with the Hardy-Weinberg Equilibrium in the patients’ group (p<0.00001). The frequency of the risk allele (allele C) was not significantly different between the patient and control groups (p=0.255). The genotypic distribution between atopic asthma and non-atopic asthma within the patientsts groups was not in agreement with the Hardy-Weinberg Equilibrium (p=0.006). However, risk allele presence showed a statistically significant association with atopy-related asthma (p=0.037).

Discussion & Conclusion: Overall, despite the finding of no association between the 5HTR2C rs6318 C allele and childhood asthma, the current results indicated that there is a strong association between the 5HTR2C rs6318 C variant and childhood atopic asthma.

High-Throughput Variant Detection Using a Color-Mixing Strategy

Many diseases are related to multiple genetic alterations along a single gene. Probing for highly multiple (>10) variants in a single quantitative PCR tube is not possible because of a limited number of fluorescence channels and one variant per channel, so many more tubes are needed. Herein, a novel color-mixing strategy that uses fluorescence combinations as digital color codes to probe multiple variants simultaneously was experimentally validated. The color-mixing strategy relies on a simple intratube assay that can probe for 15 variants as part of an intertube assay that can probe for an exponentially increased number of variants. This strategy is achieved by using multiplex double-stranded toehold probes modified with fluorophores and quenchers; the probes are designed to be quenched or luminous after binding to wild-type or variant templates. The color-mixing strategy was used to probe for 21 pathogenic variants in thalassemia and to distinguish between heterozygous and homozygous variants in six tubes, with a specificity of 99% and a sensitivity of 94%. To support tuberculosis diagnosis, the same strategy was applied to simultaneously probe in Mycobacterium tuberculosis for rifampicin-resistance mutations occurring within one 81-bp region and one 48-bp region in the rpoB gene, plus five isoniazid-resistance mutations in the inhA and katG genes.

Patagonian sheepdog: Genomic analyses trace the footprints of extinct UK herding dogs to South America

Most modern dog breeds were developed within the last two hundred years, following strong and recent human selection based predominantly on aesthetics, with few modern breeds constructed solely to maximize their work potential. In many cases, these working breeds represent the last remnants of now lost populations. The Patagonian sheepdog (PGOD), a rare herding breed, is a remarkable example of such a population. Maintained as an isolated population for over 130 years, the PGOD offers a unique opportunity to understand the genetic relationship amongst modern herding breeds, determine key genomic structure of the founder PGOD populations, and investigate how canine genomic data can mirror human migration patterns. We thus analyzed the population structure of 159 PGOD, comparing them with 1514 dogs representing 175 established breeds. Using 150,069 SNPs from a high-density SNP genotyping array, we establish the genomic composition, ancestry, and genetic diversity of the population, complementing genomic data with the PGOD's migratory history to South America. Our phylogenetic analysis reveals that PGODs are most closely related to modern herding breeds hailing from the United Kingdom. Admixture models illustrate a greater degree of diversity and genetic heterogeneity within the very small PGOD population than in Western European herding breeds, suggesting the PGOD predates the 200-year-old construction of most pure breeds known today. We thus propose that PGODs originated from the foundational herding dogs of the UK, prior to the Victorian explosion of breeds, and that they are the closest link to a now-extinct population of herding dogs from which modern herding breeds descended.

Whole-genome resequencing reveals genomic footprints of Italian sweet and hot pepper heirlooms giving insight into genes underlying key agronomic and qualitative traits

Background

Pepper is a major crop species of the Solanaceae family, largely appreciated for its high nutritional and healthy contribution to human diets. In the Mediterranean basin, the favorable pedoclimatic conditions enhanced the selection of several diversified landraces cultivated pepper (Capsicum annuum), for whom Italy can be considered a main pole of diversification. Hence, a survey of traditional C. annuum genetic resources is essential for deep understanding of such diversity and for applications in genomics assisted breeding. Here, we report whole-genome resequencing analyses of two sweet and two pungent genotypes highly diffused in South Italy and representative of the variability for shape, colour and nutritional properties.

Results

The four genomes were reconstructed at a chromosomal scale using a reference-guided approach, based on a dataset of 2.6 billion paired-end reads, corresponding to 20× genome coverage and a mapping rate above 99% for a final genomes size of approximately 3 Gb. After five iterations of variant calling, a total of 29,258,818 single nucleotide polymorphisms (SNPs) and 1,879,112 InDels, were identified. Substantial differences were observed among the four genomes based on geographical origin, with chromosomes 9 and 11 showing more polymorphisms in the accessions with higher fruit weight and absence of pungency. Among the identified variants, a small private indel (T - > TA) shared between sweet and big fruits accessions induces a frameshift with the generation of a new stop codon in a gene annotated as extensin, whereas two private SNPs within hot types were identified in 1-aminocyclopropane-1-carboxylate oxidase (ACO), a key gene involved in fruit ripening. The estimation of repetitive elements highlights a preponderant presence of Long Terminal Repeats (LTRs), the majority of which belonged to Gypsy superfamily. By comparing the four genomes with publicly available references including ‘CM334’ and Zunla-1 highlight the presence of 49,475 shared gene families.

Conclusions

The new genomic sequences aim to enrich the whole genome information of pepper local varieties, providing a valuable tool for precision gene mapping, marker discovery, comparative studies. Such knowledge widens the frontiers to understand the selection history of Italian pepper landraces toward the recognition of specificity local agri-food products marks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01039-9.

AFA: Ancestry-specific allele frequency estimation in admixed populations: The Hispanic Community Health Study/Study of Latinos

Allele frequency estimates in admixed populations, such as Hispanics and Latinos, rely on the sample's specific admixture composition and thus may differ between two seemingly similar populations. However, ancestry-specific allele frequencies, i.e., pertaining to the ancestral populations of an admixed group, may be particularly useful for prioritizing genetic variants for genetic discovery and personalized genomic health. We developed a method, ancestry-specific allele frequency estimation in admixed populations (AFA), to estimate the frequencies of biallelic variants in admixed populations with an unlimited number of ancestries. AFA uses maximum-likelihood estimation by modeling the conditional probability of having an allele given proportions of genetic ancestries. It can be applied using either local ancestry interval proportions encompassing the variant (local-ancestry-specific allele frequency estimations in admixed populations [LAFAs]) or global proportions of genetic ancestries (global-ancestry-specific allele frequency estimations in admixed populations [GAFAs]), which are easier to compute and are more widely available. Simulations and comparisons to existing software demonstrated the high accuracy of the method. We implemented AFA on high-quality imputed data of ∼9,000 Hispanics and Latinos from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL), an understudied, admixed population with three predominant continental ancestries: Amerindian, European, and African. Comparison of the European and African estimated frequencies to the respective gnomAD frequencies demonstrated high correlations (Pearson R² = 0.97-0.99). We provide a genome-wide dataset of the estimated ancestry-specific allele frequencies for available variants with allele frequency between 5% and 95% in at least one of the three ancestral populations. Association analysis of Amerindian-enriched variants with cardiometabolic traits identified five loci associated with lipid traits in Hispanics and Latinos, demonstrating the utility of ancestry-specific allele frequencies in admixed populations.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

An integrated personal and population-based Egyptian genome reference

A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

Spinal Muscular Atrophy in the Black South African Population: A Matter of Rearrangement?

Spinal muscular atrophy (SMA) is a neuromuscular disorder, characterized by muscle atrophy and impaired mobility. A homozygous deletion of survival motor neuron 1 (SMN1), exon 7 is the main cause of SMA in ~94% of patients worldwide, but only accounts for 51% of South African (SA) black patients. SMN1 and its highly homologous centromeric copy, survival motor neuron 2 (SMN2), are located in a complex duplicated region. Unusual copy number variations (CNVs) have been reported in black patients, suggesting the presence of complex pathogenic rearrangements. The aim of this study was to further investigate the genetic cause of SMA in the black SA population. Multiplex ligation-dependent probe amplification (MLPA) testing was performed on 197 unrelated black patients referred for SMA testing (75 with a homozygous deletion of SMN1, exon 7; 50 with a homozygous deletion of SMN2, exon 7; and 72 clinically suggestive patients with no homozygous deletions). Furthermore, 122 black negative controls were tested. For comparison, 68 white individuals (30 with a homozygous deletion of SMN1, exon 7; 8 with a homozygous deletion of SMN2, exon 7 and 30 negative controls) were tested. Multiple copies (>2) of SMN1, exon 7 were observed in 50.8% (62/122) of black negative controls which could mask heterozygous SMN1 deletions and potential pathogenic CNVs. MLPA is not a reliable technique for detecting carriers in the black SA population. Large deletions extending into the rest of SMN1 and neighboring genes were more frequently observed in black patients with homozygous SMN1, exon 7 deletions when compared to white patients. Homozygous SMN2, exon 7 deletions were commonly observed in black individuals. No clear pathogenic CNVs were identified in black patients but discordant copy numbers of exons suggest complex rearrangements, which may potentially interrupt the SMN1 gene. Only 8.3% (6/72) of clinically suggestive patients had heterozygous deletions of SMN1, exon 7 (1:0) which is lower than previous SA reports of 69.5%. This study emphasizes the lack of understanding of the architecture of the SMN region as well as the cause of SMA in the black SA population. These factors need to be taken into account when counseling and performing diagnostic testing in black populations.

Genome-Wide SNP Discovery in Indigenous Cattle Breeds of South Africa

Single nucleotide polymorphism arrays have created new possibilities for performing genome-wide studies to detect genomic regions harboring sequence variants that affect complex traits. However, the majority of validated SNPs for which allele frequencies have been estimated are limited primarily to European breeds. The objective of this study was to perform SNP discovery in three South African indigenous breeds (Afrikaner, Drakensberger, and Nguni) using whole genome sequencing. DNA was extracted from blood and hair samples, quantified and prepared at 50 ng/μl concentration for sequencing at the Agricultural Research Council Biotechnology Platform using an Illumina HiSeq 2500. The fastq files were used to call the variants using the Genome Analysis Tool Kit. A total of 1,678,360 were identified as novel using Run 6 of 1000 Bull Genomes Project. Annotation of the identified variants classified them into functional categories. Within the coding regions, about 30% of the SNPs were non-synonymous substitutions that encode for alternate amino acids. The study of distribution of SNP across the genome identified regions showing notable differences in the densities of SNPs among the breeds and highlighted many regions of functional significance. Gene ontology terms identified genes such as MLANA, SYT10, and CDC42EP5 that have been associated with coat color in mouse, and ADAMS3, DNAJC3, and PAG5 genes have been associated with fertility in cattle. Further analysis of the variants detected 688 candidate selective sweeps (ZHp Z-scores ≤ -4) across all three breeds, of which 223 regions were assigned as being putative selective sweeps (ZHp scores ≤-5). We also identified 96 regions with extremely low ZHp Z-scores (≤-6) in Afrikaner and Nguni. Genes such as KIT and MITF that have been associated with skin pigmentation in cattle and CACNA1C, which has been associated with biopolar disorder in human, were identified in these regions. This study provides the first analysis of sequence data to discover SNPs in indigenous South African cattle breeds. The information will play an important role in our efforts to understand the genetic history of our cattle and in designing appropriate breed improvement programmes.

Japaneseplex: A forensic SNP assay for identification of Japanese people using Japanese-specific alleles

It is sometimes necessary to determine whether a forensic biological sample came from a Japanese person. In this study, we developed a 60-locus SNP assay designed for the differentiation of Japanese people from other East Asians using entirely and nearly Japanese-specific alleles. This multiplex assay consisted of 6 independent PCR reactions followed by single nucleotide extension. The average number and standard deviation of Japanese-specific alleles possessed by an individual were 0.81 ± 0.93 in 108 Koreans from Seoul, 8.87 ± 2.89 in 103 Japanese from Tottori, 17.20 ± 3.80 in 88 Japanese from Okinawa, and 0 in 220 Han Chinese from Wuxi and Changsha. The Koreans had 0-4 Japanese-specific alleles per individual, whereas the Japanese had 4-26 Japanese-specific alleles. Almost all Japanese were distinguished from the Koreans and other people by the factorial correspondence and principal component analyses. The Snipper program was also useful to estimate the degree of Japaneseness. The method described here was successfully applied to the differentiation of Japanese from non-Japanese people in forensic cases. This Japanese-specific SNP assay was named Japaneseplex.

Population Genetics and Natural Selection in Rheumatic Disease

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Natural selection is an important influence on human genetic variation. Because immune and inflammatory function genes are enriched for signals of positive selection, the prevalence of rheumatic disease-risk alleles seen in different populations is partially the result of differing selective pressures (eg, due to pathogens). This review summarizes the genetic regions associated with susceptibility to different rheumatic diseases and concomitant evidence for natural selection, including known agents of selection exerting selective pressure in these regions.

Towards the application of precision medicine in Age-Related Macular Degeneration

The review essentially describes genetic and non-genetic variables contributing to the onset and progression of exudative Age-related Macular Degeneration (AMD) in Italian population. In particular, AMD susceptibility within Italian population is contributed to by genetic variants, accounting for 23% of disease and non-genetic variants, accounting for 10% of AMD. Our data highlighted prominent differences concerning genetic and non-genetic contributors to AMD in our cohort with respect to worldwide populations. Among genetic variables, SNPs of CFH, ARMS2, IL-8, TIMP3, SLC16A8, RAD51B, VEGFA and COL8A1 were significantly associated with the risk of AMD in the Italian cohort. Surprisingly, other susceptibility variants described in European, American and Asiatic populations, did not reach the significance threshold in our cohort. As expected, advanced age, smoking and dietary habits were associated with the disease. In addition, we also describe a number of gene-gene and gene-phenotype interactions. In fact, AMD-associated genes may be involved in the alteration of Bruch's membrane and induction of angiogenesis, contributing to exacerbate the damage caused by aging and environmental factors. Our review provides an overview of genetic and non-genetic factors characterizing AMD susceptibility in Italian population, outlining the differences with respect to the worldwide populations. Altogether, these data reflect historical, geographic, demographic and lifestyle peculiarities of Italian population. The role of epigenetics, pharmacogenetics, comorbities and genetic counseling in the management of AMD patients have been described, in the perspective of the application of a "population-specific precision medicine" approach addressed to prevent AMD onset and improve patients' quality of life.

SLC2A3 single-nucleotide polymorphism and duplication influence cognitive processing and population-specific risk for attention-deficit/hyperactivity disorder

BACKGROUND

Attention-deficit/hyperactivity disorder (ADHD) is a common, highly heritable neurodevelopmental disorder with profound cognitive, behavioral, and psychosocial impairments with persistence across the life cycle. Our initial genome-wide screening approach for copy number variants (CNVs) in ADHD implicated a duplication of SLC2A3, encoding glucose transporter-3 (GLUT3). GLUT3 plays a critical role in cerebral glucose metabolism, providing energy for the activity of neurons, which, in turn, moderates the excitatory-inhibitory balance impacting both brain development and activity-dependent neural plasticity. We therefore aimed to provide additional genetic and functional evidence for GLUT3 dysfunction in ADHD.

METHODS

Case-control association analyses of SLC2A3 single-nucleotide polymorphisms (SNPs) and CNVs were conducted in several European cohorts of patients with childhood and adult ADHD (SNP, n = 1,886 vs. 1,988; CNV, n = 1,692 vs. 1,721). These studies were complemented by SLC2A3 expression analyses in peripheral cells, functional EEG recordings during neurocognitive tasks, and ratings of food energy content.

RESULTS

Meta-analysis of all cohorts detected an association of SNP rs12842 with ADHD. While CNV analysis detected a population-specific enrichment of SLC2A3 duplications only in German ADHD patients, the CNV + rs12842 haplotype influenced ADHD risk in both the German and Spanish cohorts. Duplication carriers displayed elevated SLC2A3 mRNA expression in peripheral blood cells and altered event-related potentials reflecting deficits in working memory and cognitive response control, both endophenotypic traits of ADHD, and an underestimation of energy units of high-caloric food.

CONCLUSIONS

Taken together, our results indicate that both common and rare SLC2A3 variation impacting regulation of neuronal glucose utilization and energy homeostasis may result in neurocognitive deficits known to contribute to ADHD risk.

Intricacies in arrangement of SNP haplotypes suggest "Great Admixture" that created modern humans

Background

Inferring history from genomic sequences is challenging and problematic because chromosomes are mosaics of thousands of small Identicalby-descent (IBD) fragments, each of them having their own unique story. However, the main events in recent evolution might be deciphered from comparative analysis of numerous loci. A paradox of why humans, whose effective population size is only 10⁴, have nearly three million frequent SNPs is formulated and examined.

Results

We studied 5398 loci evenly covering all human autosomes. Common haplotypes built from frequent SNPs that are present in people from various populations have been examined. We demonstrated highly non-random arrangement of alleles in common haplotypes. Abundance of mutually exclusive pairs of common haplotypes that have different alleles at every polymorphic position (so-called Yin/Yang haplotypes) was found in 56% of loci. A novel widely spread category of common haplotypes named Mosaic has been described. Mosaic consists of numerous pieces of Yin/Yang haplotypes and represents an ancestral stage of one of them. Scenarios of possible appearance of large number of frequent human SNPs and their habitual arrangement in Yin/Yang common haplotypes have been evaluated with an advanced genomic simulation algorithm.

Conclusions

Computer modeling demonstrated that the observed arrangement of 2.9 million frequent SNPs could not originate from a sole stand-alone population. A “Great Admixture” event has been proposed that can explain peculiarities with frequent SNP distributions. This Great Admixture presumably occurred 100–300 thousand years ago between two ancestral populations that had been separated from each other about a million years ago. Our programs and algorithms can be applied to other species to perform evolutionary and comparative genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3776-5) contains supplementary material, which is available to authorized users.

Fast diffusion of domesticated maize to temperate zones

Adaptation to a temperate climate was a prerequisite for the spread of maize across a broad geographical range. To explicitly explore the demographic process underlying maize adaptation, we used a diffusion-based method to model the differentiation between temperate and tropical populations using the Non-Stiff Stalk group as a proxy for temperate maize. Based on multiple sequential Markovian coalescent approaches, we estimate that tropical and temperate maize diverged approximately 3'000 to 5'000 years ago and the population size shrank after the split. Using composite likelihood approaches, we identified a distinct tropical-temperate divergence event initiated 4'958 years ago (95% confidence interval (CI): 4'877-5'039) from an ancestral population whose effective size was 24,162 (95% CI: 23,914-24,409). We found that continuous gene flow between tropical and temperate maize accompanied the differentiation of temperate maize. Long identical-by-descent tracts shared by tropical and temperate inbred lines have been identified, which might be the result of gene flow between tropical and temperate maize or artificial selection during domestication and crop improvement. Understanding the demographic history of maize diffusion not only provides evidence for population dynamics of maize, but will also assist the identification of regions under selection and the genetic basis of complex traits of agronomic importance.

The screening of the functional microRNA binding site SNPs in sporadic colorectal cancer genes

Sporadic colorectal cancer (sCRC) is one of the most commonly diagnosed cancers worldwide, but few genetic markers have been identified and used for its early detection. MicroRNAs are diverse cellular regulators in cancer pathogenesis that bind to the 3'-untranslated region (3'-UTR) of their target mRNAs, and variants within the miRNA target sites on sCRC-related genes may influence its pathogenesis. To investigate this possibility, we used a bioinformatical method to screen SNPs for putative changes in miRNA recognition sites within the 3'-UTR of sCRC-related genes. The rs11466537 single nucleotide polymorphism was predicted to modify the regulation of hsa-miR-1193 on the Transforming Growth Factor β Receptor II (TGFBR2) gene. Additionally, luciferase reporter assays indicated that hsa-miR-1193 bound the T allele more strongly than the A allele of rs11466537 (with A being the less frequent variant), and real time-polymerase chain reaction and western blot analysis showed that TGFBR2 is significantly repressed by hsa-miR-1193. Furthermore, overexpression of hsa-miR-1193 promoted HT-29 cell proliferation, while the loss of hsa-miR-1193 inhibited the process. Finally, the rs11466537 genotyping result revealed that the frequency of A allele carriers was 1.5% in the control blood samples, but 0 in the sCRC patients' normal colon tissue samples. Our results demonstrated that hsa-miR-1193 may be involved in sCRC tumourigenesis at least in part by suppression of TGFBR2, and the A allele of rs11466537 disturbed the regulation of hsa-miR-1193 on TGFBR2.

PREMIX: PRivacy-preserving EstiMation of Individual admiXture

In this paper we proposed a framework: PRivacy-preserving EstiMation of Individual admiXture (PREMIX) using Intel software guard extensions (SGX). SGX is a suite of software and hardware architectures to enable efficient and secure computation over confidential data. PREMIX enables multiple sites to securely collaborate on estimating individual admixture within a secure enclave inside Intel SGX. We implemented a feature selection module to identify most discriminative Single Nucleotide Polymorphism (SNP) based on informativeness and an Expectation Maximization (EM)-based Maximum Likelihood estimator to identify the individual admixture. Experimental results based on both simulation and 1000 genome data demonstrated the efficiency and accuracy of the proposed framework. PREMIX ensures a high level of security as all operations on sensitive genomic data are conducted within a secure enclave using SGX.

Evidence of Polygenic Adaptation in the Systems Genetics of Anthropometric Traits

Many signals of natural selection have been identified in the human genome. However, except for some single-locus mechanisms, most molecular processes generating these adaptation signals are still unknown. We developed an approach that integrates datasets related to genome-wide association studies (GWAS) with information about systems biology and genetic signatures of natural selection to identify evidence of polygenic adaptation. Specifically, we focused on five anthropometric measurements: body mass index (BMI), height, waist-to-hip ratio adjusted for BMI (WHR), and waist circumference adjusted for BMI (WC), and sex differences for WHR and WC. We performed an enrichment analysis for signals of natural selection in protein interaction networks associated with anthropometric traits in European populations. The adaptation signals-enriched gene networks associated highlighted epistatic interactions in the context of polygenic selection for the investigated traits. These polygenic mechanisms indicated intriguing selective mechanisms related to the anthropometric traits: adult locomotory behavior for BMI, infection resistance for height, interplay between lipid transport and immune systems for WHR, and female-specific polygenic adaptation for WHR and WC. In conclusion, we observed evidence of polygenic adaptation in the context of systems genetics of anthropometric traits that indicates polygenic mechanisms related to the natural selection in European populations.

High light intensity plays a major role in emergence of population level variation in Arabidopsis thaliana along an altitudinal gradient

Environmental conditions play an important role in the emergence of genetic variations in natural populations. We identified genome-wide patterns of nucleotide variations in the coding regions of natural Arabidopsis thaliana populations. These populations originated from 700 m to 3400 m a.m.s.l. in the Western Himalaya. Using a pooled RNA-Seq approach, we identified the local and global level population-specific SNPs. The biological functions of the SNP-containing genes were primarily related to the high light intensity prevalent at high-altitude regions. The novel SNPs identified in these genes might have arisen de novo in these populations. In another approach, the F_STs of SNP-containing genes were correlated with the corresponding climatic factors. ‘Radiation in the growing season’ was the only environmental factor found to be strongly correlated with the gene-level F_STs. In both the approaches, the high light intensity was identified as the primary abiotic stress associated with the variations in these populations. The differential gene expression analysis between field and controlled condition grown plants also showed high light intensity as the primary abiotic stress, particularly for the high altitude populations. Our results provide a genome-wide perspective of nucleotide variations in populations along altitudinal gradient and their putative role in emergence of these variations.

Genetic variants involved in oxidative stress, base excision repair, DNA methylation, and folate metabolism pathways influence myeloid neoplasias susceptibility and prognosis

Myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) share common features: elevated oxidative stress, DNA repair deficiency, and aberrant DNA methylation. We performed a hospital-based case-control study to evaluate the association in variants of genes involved in oxidative stress, folate metabolism, DNA repair, and DNA methylation with susceptibility and prognosis of these malignancies. To that end, 16 SNPs (one per gene: CAT, CYBA, DNMT1, DNMT3A, DNMT3B, GPX1, KEAP1, MPO, MTRR, NEIL1, NFE2F2, OGG1, SLC19A1, SOD1, SOD2, and XRCC1) were genotyped in 191 patients (101 MDS and 90 AML) and 261 controls. We also measured oxidative stress (reactive oxygen species/total antioxidant status ratio), DNA damage (8-hydroxy-2'-deoxyguanosine), and DNA methylation (5-methylcytosine) in 50 subjects (40 MDS and 10 controls). Results showed that five genes (GPX1, NEIL1, NFE2L2, OGG1, and SOD2) were associated with MDS, two (DNMT3B and SLC19A1) with AML, and two (CYBA and DNMT1) with both diseases. We observed a correlation of CYBA TT, GPX1 TT, and SOD2 CC genotypes with increased oxidative stress levels, as well as NEIL1 TT and OGG1 GG genotypes with higher DNA damage. The 5-methylcytosine levels were negatively associated with DNMT1 CC, DNMT3A CC, and MTRR AA genotypes, and positively with DNMT3B CC genotype. Furthermore, DNMT3A, MTRR, NEIL1, and OGG1 variants modulated AML transformation in MDS patients. Additionally, DNMT3A, OGG1, GPX1, and KEAP1 variants influenced survival of MDS and AML patients. Altogether, data suggest that genetic variability influence predisposition and prognosis of MDS and AML patients, as well AML transformation rate in MDS patients. © 2016 Wiley Periodicals, Inc.

Genetics of autoimmune diseases: insights from population genetics

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Autoimmune diseases (ADs) are a family of complex heterogeneous disorders with similar underlying mechanisms characterized by immune responses against self. Collectively, ADs are common, exhibit gender and ethnic disparities, and increasing incidence. As natural selection is an important influence on human genetic variation, and immune function genes are enriched for signals of positive selection, it is thought that the prevalence of AD risk alleles seen in different population is partially the result of differing selective pressures (for example, due to pathogens). With the advent of high-throughput technologies, new analytical methodologies and large-scale projects, evidence for the role of natural selection in contributing to the heritable component of ADs keeps growing. This review summarizes the genetic regions associated with susceptibility to different ADs and concomitant evidence for selection, including known agents of selection exerting selective pressure in these regions. Examples of specific adaptive variants with phenotypic effects are included as an evidence of natural selection increasing AD susceptibility. Many of the complexities of gene effects in different ADs can be explained by population genetics phenomena. Integrating AD susceptibility studies with population genetics to investigate how natural selection has contributed to genetic variation that influences disease risk will help to identify functional variants and elucidate biological mechanisms. As such, the study of population genetics in human population holds untapped potential for elucidating the genetic causes of human disease and more rapidly focusing to personalized medicine.

Genotypic and allelic variability in CYP19A1 among populations of African and European ancestry

CYP19A1 facilitates the bioconversion of estrogens from androgens. CYP19A1 intron single nucleotide polymorphisms (SNPs) may alter mRNA splicing, resulting in altered CYP19A1 activity, and potentially influencing disease susceptibility. Genetic studies of CYP19A1 SNPs have been well documented in populations of European ancestry; however, studies in populations of African ancestry are limited. In the present study, ten 'candidate' intronic SNPs in CYP19A1 from 125 African Americans (AA) and 277 European Americans (EA) were genotyped and their frequencies compared. Allele frequencies were also compared with HapMap and ASW 1000 Genomes populations. We observed significant differences in the minor allele frequencies between AA and EA in six of the ten SNPs including rs10459592 (p<0.0001), rs12908960 (p<0.0001), rs1902584 (p = 0.016), rs2470144 (p<0.0001), rs1961177 (p<0.0001), and rs6493497 (p = 0.003). While there were no significant differences in allele frequencies between EA and CEU in the HapMap population, a 1.2- to 19-fold difference in allele frequency for rs10459592 (p = 0.004), rs12908960 (p = 0.0006), rs1902584 (p<0.0001), rs2470144 (p = 0.0006), rs1961177 (p<0.0001), and rs6493497 (p = 0.0092) was observed between AA and the Yoruba (YRI) population. Linkage disequilibrium (LD) blocks and haplotype clusters that is unique to the EA population but not AA was also observed. In summary, we demonstrate that differences in the allele frequencies of CYP19A1 intron SNPs are not consistent between populations of African and European ancestry. Thus, investigations into whether CYP19A1 intron SNPs contribute to variations in cancer incidence, outcomes and pharmacological response seen in populations of different ancestry may prove beneficial.

Genomic medicine and risk prediction across the disease spectrum

Genomic medicine is based on the knowledge that virtually every medical condition, disease susceptibility or response to treatment is caused, regulated or influenced by genes. Genetic testing may therefore add value across the disease spectrum, ranging from single-gene disorders with a Mendelian inheritance pattern to complex multi-factorial diseases. The critical factors for genomic risk prediction are to determine: (1) where the genomic footprint of a particular susceptibility or dysfunction resides within this continuum, and (2) to what extent the genetic determinants are modified by environmental exposures. Regarding the small subset of highly penetrant monogenic disorders, a positive family history and early disease onset are mostly sufficient to determine the appropriateness of genetic testing in the index case and to inform pre-symptomatic diagnosis in at-risk family members. In more prevalent polygenic non-communicable diseases (NCDs), the use of appropriate eligibility criteria is required to ensure a balance between benefit and risk. An additional screening step may therefore be necessary to identify individuals most likely to benefit from genetic testing. This need provided the stimulus for the development of a pathology-supported genetic testing (PSGT) service as a new model for the translational implementation of genomic medicine in clinical practice. PSGT is linked to the establishment of a research database proven to be an invaluable resource for the validation of novel and previously described gene-disease associations replicated in the South African population for a broad range of NCDs associated with increased cardio-metabolic risk. The clinical importance of inquiry concerning family history in determining eligibility for personalized genotyping was supported beyond its current limited role in diagnosing or screening for monogenic subtypes of NCDs. With the recent introduction of advanced microarray-based breast cancer subtyping, genetic testing has extended beyond the genome of the host to also include tumor gene expression profiling for chemotherapy selection. The decreasing cost of next generation sequencing over recent years, together with improvement of both laboratory and computational protocols, enables the mapping of rare genetic disorders and discovery of shared genetic risk factors as novel therapeutic targets across diagnostic boundaries. This article reviews the challenges, successes, increasing inter-disciplinary integration and evolving strategies for extending PSGT towards exome and whole genome sequencing (WGS) within a dynamic framework. Specific points of overlap are highlighted between the application of PSGT and exome or WGS, as the next logical step in genetically uncharacterized patients for whom a particular disease pattern and/or therapeutic failure are not adequately accounted for during the PSGT pre-screen. Discrepancies between different next generation sequencing platforms and low concordance among variant-calling pipelines caution against offering exome or WGS as a stand-alone diagnostic approach. The public reference human genome sequence (hg19) contains minor alleles at more than 1 million loci and variant calling using an advanced major allele reference genome sequence is crucial to ensure data integrity. Understanding that genomic risk prediction is not deterministic but rather probabilistic provides the opportunity for disease prevention and targeted treatment in a way that is unique to each individual patient.