Genetic Differentiation in Populations with Different Classes of Individuals

I formulate and analyse a model of population structure with different classes of individuals. These different classes may be age classes, other demographic classes, or different types of habitats homogeneously distributed over a geographical area. The value of population differentiation under an island model of dispersal and the increase of differentiation with geographical distance in one- and two-dimensional "isolation by distance" models are then obtained for a generalization of the FST measure of population structure, as a function of "effective" mutation, migration, and population size parameters. The relevant effective subpopulation size is related to the "mutation effective population size" of a single isolated subpopulation and, in models of age-structured populations, to the inbreeding effective population size.     

Genetic Analysis of Hispanic Individuals with Cystic Fibrosis

We have performed molecular genetic analyses of Hispanic individuals with cystic fibrosis (CF) in the southwestern United States. Of 129 CF chromosomes analyzed, only 46% (59/129) carry ΔF508. The G542X mutation was found on 5% (7/129) of CF chromosomes. The 3849+10kbC→T mutation, detected primarily in Ashkenazi Jews, was present on 2% (3/129). R1162X and R334W, mutations identified in Spain and Italy, each occurred on 1.6% (2/129) of CF chromosomes. W1282X and R553X were each detected once. G551D and N1303K were not found. Overall, screening for 22 or more mutations resulted in detection of only 58% of CF transmembrane conductance regulator gene mutations among Hispanic individuals. Analysis of KM19/XV2c haplotypes revealed an unusual distribution. Although the majority of ΔF508 mutations are on chromosomes of B haplotypes, the other CF mutations are on A and C haplotypes at higher-than-expected frequencies. These genetic analyses demonstrate significant differences between Hispanic individuals with CF and those of the general North American population. Assessment of carrier/affected risk in Hispanic CF individuals cannot, therefore, be based on the mutation frequencies found through studies of the general population but must be adjusted to better reflect the genetic makeup of this ethnic group. Further studies are necessary to identify the causative mutation(s) in this population and to better delineate genotype/phenotype correlations. These will enable counselors to provide more accurate genetic counseling.     

Population-Based Resequencing of APOA1 in 10,330 Individuals: Spectrum of Genetic Variation,Phenotype, and Comparison with Extreme Phenotype Approach

Rare genetic variants, identified by in-detail resequencing of loci, may contribute to complex traits. We used the apolipoprotein A-I gene (APOA1), a major high-density lipoprotein (HDL) gene, and population-based resequencing to determine the spectrum of genetic variants, the phenotypic characteristics of these variants, and how these results compared with results based on resequencing only the extremes of the apolipoprotein A-I (apoA-I) distribution. First, we resequenced APOA1 in 10,330 population-based participants in the Copenhagen City Heart Study. The spectrum and distribution of genetic variants was determined as a function of the number of individuals resequenced. Second, apoA-I and HDL cholesterol phenotypes were determined for nonsynonymous (NS) and synonymous (S) variants and were validated in the Copenhagen General Population Study (n = 45,239). Third, observed phenotypes were compared with those predicted using an extreme phenotype approach based on the apoA-I distribution. Our results are as follows: First, population-based resequencing of APOA1 identified 40 variants of which only 7 (18%) had minor allele frequencies >1%, and most were exceedingly rare. Second, 0.27% of individuals in the general population were heterozygous for NS variants which were associated with substantial reductions in apoA-I (up to 39 mg/dL) and/or HDL cholesterol (up to 0.9 mmol/L) and, surprisingly, 0.41% were heterozygous for variants predisposing to amyloidosis. NS variants associated with a hazard ratio of 1.72 (1.09–2.70) for myocardial infarction (MI), largely driven by A164S, a variant not associated with apoA-I or HDL cholesterol levels. Third, using the extreme apoA-I phenotype approach, NS variants correctly predicted the apoA-I phenotype observed in the population-based resequencing. However, using the extreme approach, between 79% (screening 0–1^st percentile) and 21% (screening 0–20^th percentile) of all variants were not identified; among these were variants previously associated with amyloidosis. Population-based resequencing of APOA1 identified a majority of rare NS variants associated with reduced apoA-1 and HDL cholesterol levels and/or predisposing to amyloidosis. In addition, NS variants associated with increased risk of MI.     

Obesity Has an Interactive Effect with Genetic Variation in the Activating Transcription Factor 6 Gene on the Risk of Pre-Diabetes in Individuals of Chinese Han Descent

Endoplasmic reticulum (ER) stress is one of the contributing factors to the development of β-cell failure in type 2 diabetes. ER stress response through ATF6 has been shown to play an important role in insulin resistance and pancreatic β-cell function. We investigated whether genetic polymorphisms in ATF6 were associated with the risk of pre-diabetes in a Chinese Han population, and whether they had a synergistic effect with obesity. Our samples included 828 individuals who were diagnosed as pre-diabetic, and 620 controls. The minor allele A at rs2340721 was associated with increased risk for pre-diabetes(p = 0.013), and this association was still significant after adjusting for gender, age, body mass index (BMI), and waist-hip ratio(p′ = 0.011). BMI, treated as a continuous variable, and rs2340721 had an interactive effect on pre-diabetic risk(p for interaction = 0.003, β = 0.106). Carriers of GG at rs7522210 were also at a higher risk compared to non-carriers (OR = 1.390, 95%CI:1.206–1.818, p = 0.013, adjusted OR′ = 1.516, 95%CI:1.101–2.006, p′ = 0.006). GG homozygotes had increased fasting blood glucose (FBG) levels(GG vs CX: 5.6±0.52 vs 5.5±0.57 mmol/L, p = 0.016), lower insulin levels (0,30,120 minutes after glucose load) (p<0.05), and reduced areas under the insulin curve than non-carriers(GG vs CX:67.3(44.2–102.3) vs 73.1(49.4–111.4), p = 0.014). rs10918270 was associated with FBG, and rs4657103 with 2 hour glucose levels after a 75 g glucose load. We also identified a haplotype of TTAG composed of rs4657103, rs2134697, rs2340721, and rs12079579, which was associated with pre-diabetes. The genetic variation in ATF6 is associated with pre-diabetes and has interactive effects with BMI on pre-diabetes in the Chinese Han population.     

Environmental Dependence of Genetic Constraint

The epistatic interactions that underlie evolutionary constraint have mainly been studied for constant external conditions. However, environmental changes may modulate epistasis and hence affect genetic constraints. Here we investigate genetic constraints in the adaptive evolution of a novel regulatory function in variable environments, using the lac repressor, LacI, as a model system. We have systematically reconstructed mutational trajectories from wild type LacI to three different variants that each exhibit an inverse response to the inducing ligand IPTG, and analyzed the higher-order interactions between genetic and environmental changes. We find epistasis to depend strongly on the environment. As a result, mutational steps essential to inversion but inaccessible by positive selection in one environment, become accessible in another. We present a graphical method to analyze the observed complex higher-order interactions between multiple mutations and environmental change, and show how the interactions can be explained by a combination of mutational effects on allostery and thermodynamic stability. This dependency of genetic constraint on the environment should fundamentally affect evolutionary dynamics and affects the interpretation of phylogenetic data.     

Genetic Variation in Healthy Oldest-Old

Individuals who live to 85 and beyond without developing major age-related diseases may achieve this, in part, by lacking disease susceptibility factors, or by possessing resistance factors that enhance their ability to avoid disease and prolong lifespan. Healthy aging is a complex phenotype likely to be affected by both genetic and environmental factors. We sequenced 24 candidate healthy aging genes in DNA samples from 47 healthy individuals aged eighty-five years or older (the ‘oldest-old’), to characterize genetic variation that is present in this exceptional group. These healthy seniors were never diagnosed with cancer, cardiovascular disease, pulmonary disease, diabetes, or Alzheimer disease. We re-sequenced all exons, intron-exon boundaries and selected conserved non-coding sequences of candidate genes involved in aging-related processes, including dietary restriction (PPARG, PPARGC1A, SIRT1, SIRT3, UCP2, UCP3), metabolism (IGF1R, APOB, SCD), autophagy (BECN1, FRAP1), stem cell activation (NOTCH1, DLL1), tumor suppression (TP53, CDKN2A, ING1), DNA methylation (TRDMT1, DNMT3A, DNMT3B) Progeria syndromes (LMNA, ZMPSTE24, KL) and stress response (CRYAB, HSPB2). We detected 935 variants, including 848 single nucleotide polymorphisms (SNPs) and 87 insertion or deletions; 41% (385) were not recorded in dbSNP. This study is the first to present a comprehensive analysis of genetic variation in aging-related candidate genes in healthy oldest-old. These variants and especially our novel polymorphisms are valuable resources to test for genetic association in models of disease susceptibility or resistance. In addition, we propose an innovative tagSNP selection strategy that combines variants identified through gene re-sequencing- and HapMap-derived SNPs.     

Genetic Variation in Heterogeneous Environments

A model of population structure in heterogeneous environments is described and conditions sufficient for maintaining a polymorphism are derived.

The absolute fitness of any genotype is regarded as a function of location in the niche space and the population density at that location. Two modes of habitat selection are examined: (1) organisms are distributed uniformly over the environment; and (2) each organism selects to occupy that habitat in which it is most fit ("optimal habitant selection").—Sufficient conditions for maintenance of genetic polymorphisms are derived for both models. In populations which do not practice habitat selection heterozygote superiority averaged over the environment is sufficient to guarantee the existence of polymorphisms. Comparable conditions for populations which practice optimal habitat selection are much less restrictive. If the heterozygotes are superior to one homozygote in any one part of the niche and to the other homozygote in any other part of the niche then a polymorphism will be defined.—A positive correlation between genetic and environmental variation follows from the model with habitat selection, but not from the other. The adaptive significance of polymorphisms thus depends on how animals behave.

     

Genetic Variation in Prosimian Species

Scheffrahn et al. (1998a,b) aimed to outline the possible contributions of genetics to conservation biology: the determination of genetic variation of wild groups and its use as a guideline in breeding programs. Unfortunately, the amount of genetic variation within wild groups of nonhuman primates is rarely known. In particular, this observation holds for prosimians. Average heterozygosity is widely considered to be a good indicator of the magnitude of intraspecific variation. We provide some values of average heterozygosity of wild populations of Otolemur garnettii, Eulemur macaco macaco and E. m. flavifrons.     

Association of Mitochondrial Genetic Variation with Carotid Atherosclerosis

In human pathology, several diseases are associated with somatic mutations in the mitochondrial genome (mtDNA). Even though mitochondrial dysfunction leads to increased oxidative stress, the role of mitochondrial mutations in atherosclerosis has not received much attention so far. In this study we analyzed the association of mitochondrial genetic variation with the severity of carotid atherosclerosis, as assessed by carotid intima-media thickness (cIMT) and the presence of coronary heart disease (CHD) in 190 subjects from Moscow, Russia, a population with high CHD occurrence. cIMT was measured by high-resolution B-mode ultrasonography and mtDNA heteroplasmies by a pyrosequencing-based method. We found that heteroplasmies for several mutations in the mtDNA in leukocytes, including C3256T, T3336C, G12315A, G13513A, G14459A, G14846A, and G15059A mutations, were significantly (p<0.001) associated with both the severity of carotid atherosclerosis and the presence of CHD. These findings indicate that somatic mitochondrial mutations have a role in the development of atherosclerosis.     

Variation in Heterozygosity Predicts Variation in Human Substitution Rates between Populations,Individuals and Genomic Regions

The “heterozygote instability” (HI) hypothesis suggests that gene conversion events focused on heterozygous sites during meiosis locally increase the mutation rate, but this hypothesis remains largely untested. As humans left Africa they lost variability, which, if HI operates, should have reduced the mutation rate in non-Africans. Relative substitution rates were quantified in diverse humans using aligned whole genome sequences from the 1,000 genomes project. Substitution rate is consistently greater in Africans than in non-Africans, but only in diploid regions of the genome, consistent with a role for heterozygosity. Analysing the same data partitioned into a series of non-overlapping 2 Mb windows reveals a strong, non-linear correlation between the amount of heterozygosity lost “out of Africa” and the difference in substitution rate between Africans and non-Africans. Putative recent mutations, derived variants that occur only once among the 80 human chromosomes sampled, occur preferentially at the centre of 2 Kb windows that have elevated heterozygosity compared both with the same region in a closely related population and with an immediately adjacent region in the same population. More than half of all substitutions appear attributable to variation in heterozygosity. This observation provides strong support for HI with implications for many branches of evolutionary biology.     

Genetic Variation on Chromosome 6 Influences F Cell Levels in Healthy Individuals of African Descent and HbF Levels in Sickle Cell Patients

Fetal haemoglobin (HbF) is a major ameliorating factor in sickle cell disease. We investigated if a quantitative trait locus on chromosome 6q23 was significantly associated with HbF and F cell levels in individuals of African descent. Single nucleotide polymorphisms (SNPs) in a 24-kb intergenic region, 33-kb upstream of the HBS1L gene and 80-kb upstream of the MYB gene, were typed in 177 healthy Afro-Caribbean subjects (AC) of approximately 7% European admixture, 631 healthy Afro-Germans (AG, a group of African and German descendents located in rural Jamaica with about 20% European admixture), 87 West African and Afro-Caribbean individuals with sickle cell anaemia (HbSS), as well as 75 Northern Europeans, which served as a contrasting population. Association with a tag SNP for the locus was detected in all four groups (AC, P = 0.005, AG, P = 0.002, HbSS patients, P = 0.019, Europeans, P = 1.5×10⁻⁷). The association signal varied across the interval in the African-descended groups, while it is more uniform in Europeans. The 6q QTL for HbF traits is present in populations of African origin and is also acting in sickle cell anaemia patients. We have started to distinguish effects originating from European and African ancestral populations in our admixed study populations.     

Selection Transforms the Landscape of Genetic Variation Interacting with Hsp90

The protein-folding chaperone Hsp90 has been proposed to buffer the phenotypic effects of mutations. The potential for Hsp90 and other putative buffers to increase robustness to mutation has had major impact on disease models, quantitative genetics, and evolutionary theory. But Hsp90 sometimes contradicts expectations for a buffer by potentiating rapid phenotypic changes that would otherwise not occur. Here, we quantify Hsp90’s ability to buffer or potentiate (i.e., diminish or enhance) the effects of genetic variation on single-cell morphological features in budding yeast. We corroborate reports that Hsp90 tends to buffer the effects of standing genetic variation in natural populations. However, we demonstrate that Hsp90 tends to have the opposite effect on genetic variation that has experienced reduced selection pressure. Specifically, Hsp90 tends to enhance, rather than diminish, the effects of spontaneous mutations and recombinations. This result implies that Hsp90 does not make phenotypes more robust to the effects of genetic perturbation. Instead, natural selection preferentially allows buffered alleles to persist and thereby creates the false impression that Hsp90 confers greater robustness.     

Association of Genetic Variants Influencing Lipid Levels with Coronary Artery Disease in Japanese Individuals

Background/Objective

In Japanese populations, we performed a replication study of genetic loci previously identified in European-descent populations as being associated with lipid levels and risk of coronary artery disease (CAD).

Methods

We genotyped 48 single nucleotide polymorphisms (SNPs) from 22 candidate loci that had previously been identified by genome-wide association (GWA) meta-analyses for low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and/or triglycerides in Europeans. We selected 22 loci with 2 parallel tracks from 95 reported loci: 16 significant loci (p<1×10⁻³⁰ in Europeans) and 6 other loci including those with suggestive evidence of lipid associations in 1292 GWA-scanned Japanese samples. Genotyping was done in 4990 general population samples, and 1347 CAD cases and 1337 controls. For 9 SNPs, we further examined CAD associations in an additional panel of 3052 CAD cases and 6335 controls.

Principal Findings

Significant lipid associations (one-tailed p<0.05) were replicated for 18 of 22 loci in Japanese samples, with significant inter-ethnic heterogeneity at 4 loci–APOB, APOE-C1, CETP, and APOA5–and allelic heterogeneity. The strongest association was detected at APOE rs7412 for LDL-C (p = 1.3×10⁻⁴¹), CETP rs3764261 for HDL-C (p = 5.2×10⁻²⁴), and APOA5 rs662799 for triglycerides (p = 5.8×10⁻⁵⁴). CAD association was replicated and/or verified for 4 loci: SORT1 rs611917 (p = 1.7×10⁻⁸), APOA5 rs662799 (p = 0.0014), LDLR rs1433099 (p = 2.1×10⁻⁷), and APOE rs7412 (p = 6.1×10⁻¹³).

Conclusions

Our results confirm that most of the tested lipid loci are associated with lipid traits in the Japanese, further indicating that in genetic susceptibility to lipid levels and CAD, the related metabolic pathways are largely common across the populations, while causal variants at individual loci can be population-specific.     

Genetic Variation in the Free-Living Amoeba Naegleria fowleri

In this study, 30 strains of the pathogenic free-living amoeba Naegleria fowleri were investigated by using the randomly amplified polymorphic DNA (RAPD) method. The present study confirmed our previous finding that RAPD variation is not correlated with geographical origin. In particular, Mexican strains belong to the variant previously detected in Asia, Europe, and the United States. In France, surprisingly, strains from Cattenom gave RAPD patterns identical to those of the Japanese strains. In addition, all of these strains, together with an additional French strain from Chooz, exhibited similarities to South Pacific strains. The results also confirmed the presence of numerous variants in Europe, whereas only two variants were detected in the United States. The two variants found in the United States were different from the South Pacific variants. These findings do not support the previous hypothesis concerning the origin and modes of dispersal of N. fowleri.     

Multivariate Analysis of Genetic Variation

Multiple factor analysis was used to interpret geographical variation of gene frequencies. Allelic frequencies at three loci (two esterase loci, Esr and Esh, and a malic dehydrogenase locus) from ants (Pogonomyrmex barbatus) collected throughout Texas and reported by Johnsonet al. (1969) were reexamined for patterns of correlation with five environmental measurements: mean January temperature, mean July temperature, mean annual precipitation, elevation, and annual growing season. These measurements and the associated gene frequencies at each sampling location were subjected to factor analysis. Variables highly correlated with the same factor were hypothesized to be causally related.Both orthogonal and oblique rotations of the factor solution provided four factors with essentially the same interpretation. Variation at the Esh locus was associated with a continuum from cold-wet to hot-dry. Variation at the Mdh locus and some of that at the Esr locus was related to the severity of winter months. Other allelic frequencies at the Esr locus had high correlations with a third factor which seemed to be independent of the environmental measurements. None of the allelic frequencies had high correlations with the fourth factor which was interpretable as an altitudinal gradient.     

Genetic Variation of Hypocotyl Length in the Brussels Sprout

A study was made of variation in the hypocotyl length of twoinbred lines of Brassica oleracea gemmifera (Brussels sprout).The pattern of means and variances obtained from parental, F₁,F₂ and backcross generations suggested polygenic control ofhypocotyl length. Measurements of parental seed weights, embryocell numbers, and the growth of seedlings from seeds of knownweight support the view that variation between the lines inhypocotyl length, when the seedlings were exposed to light,was determined by differential extension rather than by initialdifferences in seed weight. The effect of reduction in lightintensity on hypocotyl extension was studied and the differencesbetween the lines was maintained under all light intensities.Differential reaction of the lines was only observed when theperiod allowed for growth was short. The inbred line with thelonger hypocotyl was shown to suppress the other under competitiveconditions.