Natural selection at genomic regions associated with obesity and type-2 diabetes: East Asians and sub-Saharan Africans exhibit high levels of differentiation at type-2 diabetes regions

Different populations suffer from different rates of obesity and type-2 diabetes (T2D). Little is known about the genetic or adaptive component, if any, that underlies these differences. Given the cultural, geographic, and dietary variation that accumulated among humans over the last 60,000 years, we examined whether loci identified by genome-wide association studies for these traits have been subject to recent selection pressures. Using genome-wide SNP data on 938 individuals in 53 populations from the Human Genome Diversity Panel, we compare population differentiation and haplotype patterns at these loci to the rest of the genome. Using an “expanding window” approach (100–1,600 kb) for the individual loci as well as the loci as ensembles, we find a high degree of differentiation for the ensemble of T2D loci. This differentiation is most pronounced for East Asians and sub-Saharan Africans, suggesting that these groups experienced natural selection at loci associated with T2D. Haplotype analysis suggests an excess of obesity loci with evidence of recent positive selection among South Asians and Europeans, compared to sub-Saharan Africans and Native Americans. We also identify individual loci that may have been subjected to natural selection, such as the T2D locus, HHEX, which displays both elevated differentiation and extended haplotype homozygosity in comparisons of East Asians with other groups. Our findings suggest that there is an evolutionary genetic basis for population differences in these traits, and we have identified potential group-specific genetic risk factors.     

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders

1. Download : Download high-res image (272KB)
2. Download : Download full-size image

     

A thesaurus of genetic variation for interrogation of repetitive genomic regions

Detecting genetic variation is one of the main applications of high-throughput sequencing, but is still challenging wherever aligning short reads poses ambiguities. Current state-of-the-art variant calling approaches avoid such regions, arguing that it is necessary to sacrifice detection sensitivity to limit false discovery. We developed a method that links candidate variant positions within repetitive genomic regions into clusters. The technique relies on a resource, a thesaurus of genetic variation, that enumerates genomic regions with similar sequence. The resource is computationally intensive to generate, but once compiled can be applied efficiently to annotate and prioritize variants in repetitive regions. We show that thesaurus annotation can reduce the rate of false variant calls due to mappability by up to three orders of magnitude. We apply the technique to whole genome datasets and establish that called variants in low mappability regions annotated using the thesaurus can be experimentally validated. We then extend the analysis to a large panel of exomes to show that the annotation technique opens possibilities to study variation in hereto hidden and under-studied parts of the genome.     

The subgingival microbiome of clinically healthy current and never smokers

Dysbiotic oral bacterial communities have a critical role in the etiology and progression of periodontal diseases. The goal of this study was to investigate the extent to which smoking increases risk for disease by influencing the composition of the subgingival microbiome in states of clinical health. Subgingival plaque samples were collected from 200 systemically and periodontally healthy smokers and nonsmokers. 16S pyrotag sequencing was preformed generating 1 623 713 classifiable sequences, which were compared with a curated version of the Greengenes database using the quantitative insights into microbial ecology pipeline. The subgingival microbial profiles of smokers and never-smokers were different at all taxonomic levels, and principal coordinate analysis revealed distinct clustering of the microbial communities based on smoking status. Smokers demonstrated a highly diverse, pathogen-rich, commensal-poor, anaerobic microbiome that is more closely aligned with a disease-associated community in clinically healthy individuals, suggesting that it creates an at-risk-for-harm environment that is primed for a future ecological catastrophe.     

Gastric cancers of Western European and African patients show different patterns of genomic instability

Background

Infection with H. pylori is important in the etiology of gastric cancer. Gastric cancer is infrequent in Africa, despite high frequencies of H. pylori infection, referred to as the African enigma. Variation in environmental and host factors influencing gastric cancer risk between different populations have been reported but little is known about the biological differences between gastric cancers from different geographic locations. We aim to study genomic instability patterns of gastric cancers obtained from patients from United Kingdom (UK) and South Africa (SA), in an attempt to support the African enigma hypothesis at the biological level.

Methods

DNA was isolated from 67 gastric adenocarcinomas, 33 UK patients, 9 Caucasian SA patients and 25 native SA patients. Microsatellite instability and chromosomal instability were analyzed by PCR and microarray comparative genomic hybridization, respectively. Data was analyzed by supervised univariate and multivariate analyses as well as unsupervised hierarchical cluster analysis.

Results

Tumors from Caucasian and native SA patients showed significantly more microsatellite instable tumors (p < 0.05). For the microsatellite stable tumors, geographical origin of the patients correlated with cluster membership, derived from unsupervised hierarchical cluster analysis (p = 0.001). Several chromosomal alterations showed significantly different frequencies in tumors from UK patients and native SA patients, but not between UK and Caucasian SA patients and between native and Caucasian SA patients.

Conclusions

Gastric cancers from SA and UK patients show differences in genetic instability patterns, indicating possible different biological mechanisms in patients from different geographical origin. This is of future clinical relevance for stratification of gastric cancer therapy.

     

Heterogeneity in genetic admixture across different regions of Argentina

The population of Argentina is the result of the intermixing between several groups, including Indigenous American, European and African populations. Despite the commonly held idea that the population of Argentina is of mostly European origin, multiple studies have shown that this process of admixture had an impact in the entire Argentine population. In the present study we characterized the distribution of Indigenous American, European and African ancestry among individuals from different regions of Argentina and evaluated the level of discrepancy between self-reported grandparental origin and genetic ancestry estimates. A set of 99 autosomal ancestry informative markers (AIMs) was genotyped in a sample of 441 Argentine individuals to estimate genetic ancestry. We used non-parametric tests to evaluate statistical significance. The average ancestry for the Argentine sample overall was 65% European (95%CI: 63–68%), 31% Indigenous American (28–33%) and 4% African (3–4%). We observed statistically significant differences in European ancestry across Argentine regions [Buenos Aires province (BA) 76%, 95%CI: 73–79%; Northeast (NEA) 54%, 95%CI: 49–58%; Northwest (NWA) 33%, 95%CI: 21–41%; South 54%, 95%CI: 49–59%; p<0.0001] as well as between the capital and immediate suburbs of Buenos Aires city compared to more distant suburbs [80% (95%CI: 75–86%) versus 68% (95%CI: 58–77%), p = 0.01]. European ancestry among individuals that declared all grandparents born in Europe was 91% (95%CI: 88–94%) compared to 54% (95%CI: 51–57%) among those with no European grandparents (p<0.001). Our results demonstrate the range of variation in genetic ancestry among Argentine individuals from different regions in the country, highlighting the importance of taking this variation into account in genetic association and admixture mapping studies in this population.     

Iodine levels in different regions of the human brain

BackgroundIodine is a key component of the thyroid hormones thyroxine (T4) and triiodothyronine (T3), which are crucial for proper growth and development of the human body. In particular, a great body of literature has been published on the link between thyroid hormones and brain development and functioning. However, there is a lack of knowledge on the iodine levels in the human brain. The aim of this work was to determine the brain iodine levels and to contribute to the establishment of “reference” levels for iodine in the different anatomical and functional regions of normal (i.e., subjects without neurological or psychiatric diseases) human brain.MethodsThe iodine levels were determined in 14 brain regions of 52 dead subjects without evidence of neurological or psychiatric disease (n = 728 samples). Iodine was extracted from brain samples using a standard procedure and determined by inductively coupled plasma – mass spectrometry (ICP-MS).ResultsFour subjects presented abnormally high brain iodine levels (26.0 ± 14.2 μg/g) and were excluded from the overall data analysis. The average brain iodine levels for the remaining 48 subjects was 0.14 ± 0.13 μg/g dry weight. Iodine showed very heterogeneous distribution across the different brain regions, with the frontal cortex, caudate nucleus and putamen showing the highest levels. Interestingly, these brain regions are closely related to cognitive function. Iodine levels also showed a tendency to increase with age. The high levels observed in four subjects seemed to be related to previous exposure to iodine-based contrast agents widely used in radiology and computed tomography exams.ConclusionsThis paper provides important data on iodine levels at different brain regions in “normal” people, which can be used to interpret eventual imbalances in subjects with mental disorders and neurodegenerative diseases.     

GIGANTEA and EARLY FLOWERING 4 in Arabidopsis exhibit differential phase-specific genetic influences over a diurnal cycle

The endogenous circadian clock regulates many physiological processes related to plant survival and adaptability. GIGANTEA (GI), a clock-associated protein, contributes to the maintenance of circadian period length and amplitude, and also regulates flowering time and hypocotyl growth in response to day length. Similarly, EARLY FLOWERING 4 (ELF4), another clock regulator, also contributes to these processes. However, little is known about either the genetic or molecular interactions between GI and ELF4 in Arabidopsis. In this study, we investigated the genetic interactions between GI and ELF4 in the regulation of circadian clock-controlled outputs. Our mutant analysis shows that GI is epistatic to ELF4 in flowering time determination, while ELF4 is epistatic to GI in hypocotyl growth regulation. Moreover, GI and ELF4 have a synergistic or additive effect on endogenous clock regulation. Gene expression profiling of gi, elf4, and gi elf4 mutants further established that GI and ELF4 have differentially dominant influences on circadian physiological outputs at dusk and dawn, respectively. This phasing of GI and ELF4 influences provides a potential means to achieve diversity in the regulation of circadian physiological outputs, including flowering time and hypocotyl growth.     

Extremely large chromosomal deletions are intimately involved in genetic instability and genomic rearrangements in Streptomyces glaucescens

Summary Genetic instability inStreptomyces glaucescens characteristically involves the occurrence of gross genomic rearrangements including high-level sequence amplification and extensive deletion. We investigated the relationship of the unstablemelC andstrS loci and a 100 kb region of the chromosome which frequently gives rise to intense heterogeneous DNA amplification. Standard chromosome walking using a cosmid bank in conjunction with a “reverse-blot” procedure enabled us to construct a contiguous genomicBamHI map of the unstable region exceeding 900 kb. The unstable genes and the amplifiable region (AUD locus) are physically linked within a 600 kb segment of the chromosome. The previously characterized deletions which affect these loci are merely components of much larger deletions ranging from 270 to over 800 kb which are polar in nature, effecting the sequential loss of thestrS andmelC loci. The more extensive deletions terminate either adjacent to, or in the vicinity of DNA reiterations at the AUD locus. Additionally, a deletion junction fragment and the corresponding deletion ends were cloned and analysed at the sequence level.     

Microsatellite instability in pediatric high grade glioma is associated with genomic profile and differential target gene inactivation

High grade gliomas (HGG) are one of the leading causes of cancer-related deaths in children, and there is increasing evidence that pediatric HGG may harbor distinct molecular characteristics compared to adult tumors. We have sought to clarify the role of microsatellite instability (MSI) in pediatric versus adult HGG. MSI status was determined in 144 patients (71 pediatric and 73 adults) using a well established panel of five quasimonomorphic mononucleotide repeat markers. Expression of MLH1, MSH2, MSH6 and PMS2 was determined by immunohistochemistry, MLH1 was assessed for mutations by direct sequencing and promoter methylation using MS-PCR. DNA copy number profiles were derived using array CGH, and mutations in eighteen MSI target genes studied by multiplex PCR and genotyping. MSI was found in 14/71 (19.7%) pediatric cases, significantly more than observed in adults (5/73, 6.8%; p = 0.02, Chi-square test). MLH1 expression was downregulated in 10/13 cases, however no mutations or promoter methylation were found. MSH6 was absent in one pediatric MSI-High tumor, consistent with an inherited mismatch repair deficiency associated with germline MSH6 mutation. MSI was classed as Type A, and associated with a remarkably stable genomic profile. Of the eighteen classic MSI target genes, we identified mutations only in MSH6 and DNAPKcs and described a polymorphism in MRE11 without apparent functional consequences in DNA double strand break detection and repair. This study thus provides evidence for a potential novel molecular pathway in a proportion of gliomas associated with the presence of MSI.     

The influence of landscape configuration and environment on population genetic structure in a sedentary passerine: insights from loci located in different genomic regions

The study of the factors structuring genetic variation can help to infer the neutral and adaptive processes shaping the demographic and evolutionary trajectories of natural populations. Here, we analyse the role of isolation by distance (IBD), isolation by resistance (IBR, defined by landscape composition) and isolation by environment (IBE, estimated as habitat and elevation dissimilarity) in structuring genetic variation in 25 blue tit (Cyanistes caeruleus) populations. We typed 1385 individuals at 26 microsatellite loci classified into two groups by considering whether they are located into genomic regions that are actively (TL; 12 loci) or not (NTL; 14 loci) transcribed to RNA. Population genetic differentiation was mostly detected using the panel of NTL. Landscape genetic analyses showed a pattern of IBD for all loci and the panel of NTL, but genetic differentiation estimated at TL was only explained by IBR models considering high resistance for natural vegetation and low resistance for agricultural lands. Finally, the absence for IBE suggests a lack of divergent selection pressures associated with differences in habitat and elevation. Overall, our study shows that markers located in different genomic regions can yield contrasting inferences on landscape‐level patterns of realized gene flow in natural populations.     

Phylogenetically closely related pseudomonads isolated from arthropods exhibit differential insect-killing abilities and genetic variations in insecticidal factors

Strains belonging to the Pseudomonas protegens and Pseudomonas chlororaphis species are able to control soilborne plant pathogens and to kill pest insects by producing virulence factors such as toxins, chitinases, antimicrobials or two-partner secretion systems. Most insecticidal Pseudomonas described so far were isolated from roots or soil. It is unknown whether these bacteria naturally occur in arthropods and how they interact with them. Therefore, we isolated P. protegens and P. chlororaphis from various healthy insects and myriapods, roots and soil collected in an agricultural field and a neighbouring grassland. The isolates were compared for insect killing, pathogen suppression and host colonization abilities. Our results indicate that neither the origin of isolation nor the phylogenetic position mirror the degree of insecticidal activity. Pseudomonas protegens strains appeared homogeneous regarding phylogeny, biocontrol and insecticidal capabilities, whereas P. chlororaphis strains were phylogenetically and phenotypically more heterogenous. A phenotypic and genomic analysis of five closely related P. chlororaphis isolates displaying varying levels of insecticidal activity revealed variations in genes encoding insecticidal factors that may account for the reduced insecticidal activity of certain isolates. Our findings point towards an adaption to insects within closely related pseudomonads and contribute to understand the ecology of insecticidal Pseudomonas.     

The consequences of structural genomic alterations in humans: genomic disorders, genomic instability and cancer

Over the last decade or so, sophisticated technological advances in array-based genomics have firmly established the contribution of structural alterations in the human genome to a variety of complex developmental disorders, and also to diseases such as cancer. In fact, multiple 'novel' disorders have been identified as a direct consequence of these advances. Our understanding of the molecular events leading to the generation of these structural alterations is also expanding. Many of the models proposed to explain these complex rearrangements involve DNA breakage and the coordinated action of DNA replication, repair and recombination machinery. Here, and within the context of Genomic Disorders, we will briefly overview the principal models currently invoked to explain these chromosomal rearrangements, including Non-Allelic Homologous Recombination (NAHR), Fork Stalling Template Switching (FoSTeS), Microhomology Mediated Break-Induced Repair (MMBIR) and Breakage-fusion-bridge cycle (BFB). We will also discuss an unanticipated consequence of certain copy number variations (CNVs) whereby the CNVs potentially compromise fundamental processes controlling genomic stability including DNA replication and the DNA damage response. We will illustrate these using specific examples including Genomic Disorders (DiGeorge/Veleocardiofacial syndrome, HSA21 segmental aneuploidy and rec (3) syndrome) and cell-based model systems. Finally, we will review some of the recent exciting developments surrounding specific CNVs and their contribution to cancer development as well as the latest model for cancer genome rearrangement; 'chromothripsis'.     

Human genomic regions with exceptionally high levels of population differentiation identified from 911 whole-genome sequences

Background

Population differentiation has proved to be effective for identifying loci under geographically localized positive selection, and has the potential to identify loci subject to balancing selection. We have previously investigated the pattern of genetic differentiation among human populations at 36.8 million genomic variants to identify sites in the genome showing high frequency differences. Here, we extend this dataset to include additional variants, survey sites with low levels of differentiation, and evaluate the extent to which highly differentiated sites are likely to result from selective or other processes.

Results

We demonstrate that while sites with low differentiation represent sampling effects rather than balancing selection, sites showing extremely high population differentiation are enriched for positive selection events and that one half may be the result of classic selective sweeps. Among these, we rediscover known examples, where we actually identify the established functional SNP, and discover novel examples including the genes ABCA12, CALD1 and ZNF804, which we speculate may be linked to adaptations in skin, calcium metabolism and defense, respectively.

Conclusions

We identify known and many novel candidate regions for geographically restricted positive selection, and suggest several directions for further research.     

Genomic instability and mobile genetic elements in regions surrounding two discoidin I genes of Dictyostelium discoideum.

We have found that the genomic regions surrounding the linked discoidin I genes of various Dictyostelium discoideum strains have undergone rapid changes. Wild-type strain NC-4 has three complete discoidin I genes; its axenic derivative strain Ax-3L has duplicated a region starting approximately 1 kilobase upstream from the two linked genes and extending for at least 8 kilobases past the genes. A separately maintained stock, strain Ax-3K, does not have this duplication but has undergone a different rearrangement approximately 3 kilobases farther upstream. We show that there are repeat elements in these rapidly changing regions. At least two of these elements, Tdd-2 and Tdd-3, have characteristics associated with mobile genetic elements. The Tdd-3 element is found in different locations in related strains and causes a 9- to 10-base-pair duplication of the target site DNA. The Tdd-2 and Tdd-3 elements do not cross-hybridize, but they share a 22-base-pair homology near one end. At two separate sites, the Tdd-3 element has transposed into the Tdd-2 element, directly adjacent to the 22-base-pair homology. The Tdd-3 element may use this 22-base-pair region as a preferential site of insertion.