Genetics and epidemiology of gallbladder disease in New World native peoples.

Native peoples of the New World, including Amerindians and admixed Latin Americans such as Mexican-Americans, are highly susceptible to diseases of the gallbladder. These include cholesterol cholelithiasis (gallstones) and its complications, as well as cancer of the gallbladder. Although there is clearly some necessary dietary or other environmental risk factor involved, the pattern of disease prevalence is geographically associated with the distribution of genes of aboriginal Amerindian origin, and levels of risk generally correspond to the degree of Amerindian admixture. This pattern differs from that generally associated with Westernization, which suggests a gene-environment interaction, and that within an admixed population there is a subset whose risk is underestimated when admixture is ignored. The risk that an individual of a susceptible New World genotype will undergo a cholecystectomy by age 85 can approach 40% in Mexican-American females, and their risk of gallbladder cancer can reach several percent. These are heretofore unrecognized levels of risk, especially of the latter, because previous studies have not accounted for admixture or for the loss of at-risk individuals due to cholecystectomy. A genetic susceptibility may, thus, be as "carcinogenic" in New World peoples as any known major environmental exposure; yet, while the risk has a genetic basis, its expression as gallbladder cancer is so delayed as to lead only very rarely to multiply-affected families. Estimates in this paper are derived in part from two studies of Mexican-Americans in Starr County and Laredo, Texas.

Recent speciation between sympatric Tanganyikan cichlid colour morphs

Lake Tanganyika, Africa’s oldest lake, harbours an impressive diversity of cichlid fishes. Although diversification in its radiating groups is thought to have been initially rapid, cichlids from Lake Tanganyika show little evidence for ongoing speciation. In contrast, examples of recent divergence among sympatric colour morphs are well known in haplochromine cichlids from Lakes Malawi and Victoria. Here, we report genetic evidence for recent divergence between two sympatric Tanganyikan cichlid colour morphs. These Petrochromis morphs share mitochondrial haplotypes, yet microsatellite loci reveal that their sympatric populations form distinct genetic groups. Nuclear divergence between the two morphs is equivalent to that which arises geographically within one of the morphs over short distances and is substantially smaller than that among other sympatric species in this genus. These patterns suggest that these morphs diverged only recently, yet that barriers to gene flow exist which prevent extensive admixture despite their sympatric distribution. The morphs studied here provide an unusual example of active diversification in Lake Tanganyika’s generally ancient cichlid fauna and enable comparisons of speciation processes between Lake Tanganyika and other African lakes.

Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

Imputation of sequence level genotypes in the Franches-Montagnes horse breed

BACKGROUND A cost-effective strategy to increase the density of available markers within a population is to sequence a small proportion of the population and impute whole-genome sequence data for the remaining population. Increased densities of typed markers are advantageous for genome-wide association studies (GWAS) and genomic predictions. METHODS We obtained genotypes for 54 602 SNPs (single nucleotide polymorphisms) in 1077 Franches-Montagnes (FM) horses and Illumina paired-end whole-genome sequencing data for 30 FM horses and 14 Warmblood horses. After variant calling, the sequence-derived SNP genotypes (~13 million SNPs) were used for genotype imputation with the software programs Beagle, Impute2 and FImpute. RESULTS The mean imputation accuracy of FM horses using Impute2 was 92.0%. Imputation accuracy using Beagle and FImpute was 74.3% and 77.2%, respectively. In addition, for Impute2 we determined the imputation accuracy of all individual horses in the validation population, which ranged from 85.7% to 99.8%. The subsequent inclusion of Warmblood sequence data further increased the correlation between true and imputed genotypes for most horses, especially for horses with a high level of admixture. The final imputation accuracy of the horses ranged from 91.2% to 99.5%. CONCLUSIONS Using Impute2, the imputation accuracy was higher than 91% for all horses in the validation population, which indicates that direct imputation of 50k SNP-chip data to sequence level genotypes is feasible in the FM population. The individual imputation accuracy depended mainly on the applied software and the level of admixture.

Late Quaternary history of North Eurasian Norway spruce ( Picea abies ) and Siberian spruce ( Picea obovata ) inferred from macrofossils, pollen and cytoplasmic DNA variation

Aim We used combined palaeobotanical and genetic data to assess whether Norway spruce (Picea abies) and Siberian spruce (Picea obovata), two major components of the Eurasian boreal forests, occupied separate glacial refugia, and to test previous hypotheses on their distinction, geographical delimitation and introgression. Location The range of Norway spruce in northern Europe and Siberian spruce in northern Asia. Methods Pollen data and recently compiled macrofossil records were summarized for the Last Glacial Maximum (LGM), late glacial and Holocene. Genetic variation was assessed in 50 populations using one maternally (mitochondrial nad1) and one paternally (chloroplast trnT–trnL) inherited marker and analysed using spatial analyses of molecular variance (SAMOVA). Results Macrofossils showed that spruce was present in both northern Europe and Siberia at the LGM. Congruent macrofossil and pollen data from the late glacial suggested widespread expansions of spruce in the East European Plain, West Siberian Plain, southern Siberian mountains and the Baikal region. Colonization was largely completed during the early Holocene, except in the formerly glaciated area of northern Europe. Both DNA markers distinguished two highly differentiated groups that correspond to Norway spruce and Siberian spruce and coincide spatially with separate LGM spruce occurrences. The division of the mtDNA variation was geographically well defined and occurred to the east of the Ural Mountains along the Ob River, whereas the cpDNA variation showed widespread admixture. Genetic diversity of both DNA markers was higher in western than in eastern populations. Main conclusions North Eurasian Norway spruce and Siberian spruce are genetically distinct and occupied separate LGM refugia, Norway spruce on the East European Plain and Siberian spruce in southern Siberia, where they were already widespread during the late glacial. They came into contact in the basin of the Ob River and probably hybridized. The lower genetic diversity in the eastern populations may indicate that Siberian spruce suffered more from past climatic fluctuations than Norway spruce.

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y-chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least-divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi-absence of Y-chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non-neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.

Models of hybridization during range expansions and their application to recent human evolution

Several lines of genetic, archeological and paleontological evidence suggest that anatomically modern humans (Homo sapiens) colonized the world in the last 60,000 years by a series of migrations originating from Africa (e.g. Liu et al., 2006; Handley et al., 2007; Prugnolle, Manica, and Balloux, 2005; Ramachandran et al. 2005; Li et al. 2008; Deshpande et al. 2009; Mellars, 2006a, b; Lahr and Foley, 1998; Gravel et al., 2011; Rasmussen et al., 2011). With the progress of ancient DNA analysis, it has been shown that archaic humans hybridized with modern humans outside Africa. Recent direct analyses of fossil nuclear DNA have revealed that 1–4 percent of the genome of Eurasian has been likely introgressed by Neanderthal genes (Green et al., 2010; Reich et al., 2010; Vernot and Akey, 2014; Sankararaman et al., 2014; Prufer et al., 2014; Wall et al., 2013), with Papua New Guineans and Australians showing even larger levels of admixture with Denisovans (Reich et al., 2010; Skoglund and Jakobsson, 2011; Reich et al., 2011; Rasmussen et al., 2011). It thus appears that the past history of our species has been more complex than previously anticipated (Alves et al., 2012), and that modern humans hybridized several times with local hominins during their expansion out of Africa, but the exact mode, time and location of these hybridizations remain to be clarifi ed (Ibid.; Wall et al., 2013). In this context, we review here a general model of admixture during range expansion, which lead to some predictions about expected patterns of introgression that are relevant to modern human evolution.

Hybrid ‘superswarm’ leads to rapid divergence and establishment of populations during a biological invasion

Understanding the genetic background of invading species can be crucial information clarifying why they become invasive. Intraspecific genetic admixture among lineages separated in the native ranges may promote the rate and extent of an invasion by substantially increasing standing genetic variation. Here we examine the genetic relationships among threespine stickleback that recently colonized Switzerland. This invasion results from several distinct genetic lineages that colonized multiple locations and have since undergone range expansions, where they coexist and admix in parts of their range. Using 17 microsatellites genotyped for 634 individuals collected from 17 Swiss and two non-Swiss European sites, we reconstruct the invasion of stickleback and investigate the potential and extent of admixture and hybridization among the colonizing lineages from a population genetic perspective. Specifically we test for an increase in standing genetic variation in populations where multiple lineages coexist. We find strong evidence of massive hybridization early on, followed by what appears to be recent increased genetic isolation and the formation of several new genetically distinguishable populations, consistent with a hybrid ‘superswarm’. This massive hybridization and population formation event(s) occurred over approximately 140 years and likely fuelled the successful invasion of a diverse range of habitats. The implications are that multiple colonizations coupled with hybridization can lead to the formation of new stable genetic populations potentially kick-starting speciation and adaptive radiation over a very short time.

Intersystem crossing rates of S1 state keto-amino cytosine at low excess energy

The amino-keto tautomer of supersonic jet-cooled cytosine undergoes intersystem crossing (ISC) from the v = 0 and low-lying vibronic levels of its S1(¹ππ*) state. We investigate these ISC rates experimentally and theoretically as a function of S1 state vibrational excess energy Eexc. The S1 vibronic levels are pumped with a ~5 ns UV laser, the S1 and triplet state ion signals are separated by prompt or delayed ionization with a second UV laser pulse. After correcting the raw ISC yields for the relative S1 and T1ionization cross sections, we obtain energy dependent ISC quantum yields Q corr ISC =1% –5%. These are combined with previously measured vibronic state-specific decay rates, giving ISC rates kISC = 0.4–1.5 ⋅ 10⁹ s⁻¹, the corresponding S1⇝S0internal conversion (IC) rates are 30–100 times larger. Theoretical ISC rates are computed using SCS-CC2 methods, which predict rapid ISC from the S1; v = 0 state with kISC = 3 ⋅ 10⁹ s⁻¹ to the T1(³ππ*) triplet state. The surprisingly high rate of this El Sayed-forbidden transition is caused by a substantial admixture of ¹nOπ* character into the S1(¹ππ*) wave function at its non-planar minimum geometry. The combination of experiment and theory implies that (1) below Eexc = 550 cm⁻¹ in the S1 state, S1⇝S0internal conversion dominates the nonradiative decay with kIC ≥ 2 ⋅ 10¹⁰ s⁻¹, (2) the calculated S1⇝T1 (¹ππ*⇝³ππ*) ISC rate is in good agreement with experiment, (3) being El-Sayed forbidden, the S1⇝T1 ISC is moderately fast (kISC = 3 ⋅ 10⁹ s⁻¹), and not ultrafast, as claimed by other calculations, and (4) at Eexc ~ 550 cm⁻¹ the IC rate increases by ~50 times, probably by accessing the lowest conical intersection (the C5-twist CI) and thereby effectively switching off the ISC decay channels.

Volcanic terrain and the possible periglacial formation of "excess ice" at the mid-latitudes of Utopia Planitia, Mars

At the mid-latitudes of Utopia Planitia (UP), Mars, a suite of spatially-associated landforms exhibit geomorphological traits that, on Earth, would be consistent with periglacial processes and the possible freeze-thaw cycling of water. The suite comprises small-sized polygonally-patterned ground, polygon-junction and -margin pits, and scalloped, rimless depressions. Typically, the landforms incise a dark-toned terrain that is thought to be ice-rich. Here, we investigate the dark-toned terrain by using high resolution images from the HiRISE as well as near-infrared spectral-data from the OMEGA and CRISM. The terrain displays erosional characteristics consistent with a sedimentary nature and near-infrared spectra characterised by a blue slope similar to that of weathered basaltic-tephra. We also describe volcanic terrain that is dark-toned and periglacially-modified in the Kamchatka mountain-range of eastern Russia. The terrain is characterised by weathered tephra inter-bedded with snow, ice-wedge polygons and near-surface excess ice. The excess ice forms in the pore space of the tephra as the result of snow-melt infiltration and, subsequently, in-situ freezing. Based on this possible analogue, we construct a three-stage mechanism that explains the possible ice-enrichment of a broad expanse of dark-toned terrain at the mid-latitudes of UP: (1) the dark-toned terrain accumulates and forms via the regional deposition of sediments sourced from explosive volcanism; (2) the volcanic sediments are blanketed by atmospherically-precipitated (H2O) snow, ice or an admixture of the two, either concurrent with the volcanic-events or between discrete events; and, (3) under the influence of high obliquity or explosive volcanism, boundary conditions tolerant of thaw evolve and this, in turn, permits the migration, cycling and eventual formation of excess ice in the volcanic sediments. Over time, and through episodic iterations of this scenario, excess ice forms to decametres of depth. (C) 2015 Elsevier B.V. All rights reserved.

Traditional linkage analysis in admixed families

Currently there is no general method to study the impact of population admixture within families on the assumptions of random mating and consequently, Hardy-Weinberg equilibrium (HWE) and linkage equilibrium (LE) and on the inference obtained from traditional linkage analysis. ^ First, through simulation, the effect of admixture of two populations on the log of the odds (LOD) score was assessed, using Prostate Cancer as the typical disease model. Comparisons between simulated mixed and homogeneous families were performed. LOD scores under both models of admixture (within families and within a data set of homogeneous families) were closest to the homogeneous family scores of the population having the highest mixing proportion. Random sampling of families or ascertainment of families with disease affection status did not affect this observation, nor did the mode of inheritance (dominant/recessive) or sample size. ^ Second, after establishing the effect of admixture on the LOD score and inference for linkage, the presence of induced disequilibria by population admixture within families was studied and an adjustment procedure was developed. The adjustment did not force all disequilibria to disappear but because the families were adjusted for the population admixture, those replicates where the disequilibria exist are no longer affected by the disequilibria in terms of maximization for linkage. Furthermore, the adjustment was able to exclude uninformative families or families that had such a high departure from HWE and/or LE that their LOD scores were not reliable. ^ Together these observations imply that the presence of families of mixed population ancestry impacts linkage analysis in terms of the LOD score and the estimate of the recombination fraction. ^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

Planktonic foraminifera of the Mediterranean Sea

A high-resolution biochronology is presented for the Late Quaternary of the central Mediterranean. In the Late Pleistocene-Holocene successions three assemblage zones are distinguished on the basis of frequency patterns of planktic foraminifera. The age of these zones is determined by Accelerator Mass Spectrometry (AMS)14C dating. The zonal boundaries are dated at 12,700 yr B.P. (the end of Termination Ia) and 9600 yr B.P. (the start of Termination Ib), respectively. The AMS dates show that major changes in the planktic and benthic realms occurred synchronously over wide areas, although records of individual species may show important regional differences. In the studied areas, resedimentation processes revealed by anomalous successions of14C dates, play a far more important role than indicated by the sedimentological and micropaleontological data. Possibly these processes contribute to the very high accumulation rates in the glacial Zone III. Although the AMS technique has increased the accuracy of14C-measurements, admixture of older carbonate may still lead to substantial age differences between areas with different sedimentary regimes.

(Table 3) Chemical composition of phosphatized rocks from the Kammu Seamount

Phosphatized biogenic limestones and phosphorites with initial Fe-Mn mineralization dredged from the summit surface of the Kammu Seamount (Milwaukee Seamounts, northwestern Pacific) are studied. The rocks are largely composed of nannofossils and planktonic foraminifers with an admixture of benthic foraminifers, bryozoans, and other organic remains, presumably including bacterial ones. The nannofosssil and foraminiferal assemblages indicate Quaternary age of sediments, and their phosphatization is consistent with the phosphatization age determined previously based on nonequilibrium uranium (within the limits of 1 My). The age of phosphatization and the Fe-Mn mineralization in the sediments from Pacific seamounts that young implies dependence of these ore-forming processes on oceanic environments favorable for ore accumulation rather than on their age.

Major and trace element analyses from DSDP Hole 26-253 at the Ninetyeast Ridge in the Indian Ocean

A basaltic sequence of Eocene submarine-erupted pyroclastic sediments totals at least 388 m at DSDP Site 253 on the Ninetyeast Ridge. These fossiliferous hyaloclastic sediments have been erupted and fragmented by explosive volcanism (hydroexplosions) in shallow water. The occurrence of interbedded basaltic ash-fall tuffs within the younger horizons of the hyaloclastic sequence marks the emergence of some Ninetyeast Ridge volcanic vents above sea level. Considerable textural variation allows subdivision of the sequence into six informal lithostratigraphic units. Hydrothermal and diagenetic alteration has caused the complete replacement of all original glass by smectites, and the introduction of abundant zeolite and calcite cements. The major and trace element contents of the hyaloclastites vary due to the alteration, and the admixture of biogenous calcite. On a calcium carbonate-free basis systematic variations are recognisable. Mg, Ni, Cr and Cu are enriched, and Li and Zn depleted in the three older units relative to the younger three. The chemical variability is reflected by the development of saponite in the older part of the sequence and montmorillonite in the younger; and by the presence of a quartz-normative basalt flow occurring in Unit II, in contrast to the Mg-rich highly olivine-normative basalt at the base of the sequence. The younger and older parts of the sequence therefore appear to have been derived from magmas of different chemistry. The sequence, like other basaltic rocks recovered from the Ninetyeast Ridge, is enriched in the light relative to the heavy rare earth elements (REE) although the REE contents vary unsystematically with depth, probably because of the high-temperature subaqueous alteration and the presence of biogenous calcite. This REE data indicates that the Ninetyeast Ridge volcanism was different from that which produces mid-ocean ridge basalts.