Population genetic structure and demography analysis in palinurus elephas populations from atlantic and mediterranean waters

Report for the scientific sojourn at the University of Reading, United Kingdom, from January until May 2008. The main objectives have been firstly to infer population structure and parameters in demographic models using a total of 13 microsatellite loci for genotyping approximately 30 individuals per population in 10 Palinurus elephas populations both from Mediterranean and Atlantic waters. Secondly, developing statistical methods to identify discrepant loci, possibly under selection and implement those methods using the R software environment. It is important to consider that the calculation of the probability distribution of the demographic and mutational parameters for a full genetic data set is numerically difficult for complex demographic history (Stephens 2003). The Approximate Bayesian Computation (ABC), based on summary statistics to infer posterior distributions of variable parameters without explicit likelihood calculations, can surmount this difficulty. This would allow to gather information on different demographic prior values (i.e. effective population sizes, migration rate, microsatellite mutation rate, mutational processes) and assay the sensitivity of inferences to demographic priors by assuming different priors.

Quantification, self-renewal, and genetic tracing of FL1+ tumor-initiating cells in a large cohort of human gliomas.

Evidence has emerged that the initiation and growth of gliomas is sustained by a subpopulation of cancer-initiating cells (CICs). Because of the difficulty of using markers to tag CICs in gliomas, we have previously exploited more robust phenotypic characteristics, including a specific morphology and intrincic autofluorescence, to identify and isolate a subpopulation of glioma CICs, called FL1(+). The objective of this study was to further validate our method in a large cohort of human glioma and a mouse model of glioma. Seventy-four human gliomas of all grades and the GFAP-V(12)HA-ras B8 mouse model were analyzed for in vitro self-renewal capacity and their content of FL1(+). Nonneoplastic brain tissue and embryonic mouse brain were used as control. Genetic traceability along passages was assessed with microsatellite analysis. We found that FL1(+) cells from low-grade gliomas and from control nonneoplasic brain tissue show a lower level of autofluorescence and undergo a restricted number of cell divisions before dying in culture. In contrast, we found that FL1(+) cells derived from many but not all high-grade gliomas acquire high levels of autofluorescence and can be propagated in long-term cultures. Moreover, FL1(+) cells show a remarkable traceability over time in vitro and in vivo. Our results show that FL1(+) cells can be found in all specimens of a large cohort of human gliomas of different grades and in a model of genetically induced mouse glioma as well as nonneoplastic brain. However, their self-renewal capacity is variable and seems to be dependent on the tumor grade.

Sequencing of simple sequence repeat anchored polymerase chain reaction amplification products of Biomphalaria glabrata

Simple sequence repeat anchored polymerase chain reaction amplification (SSR-PCR) is a genetic typing technique based on primers anchored at the 5' or 3' ends of microsatellites, at high primer annealing temperatures. This technique has already been used in studies of genetic variability of several organisms, using different primer designs. In order to conduct a detailed study of the SSR-PCR genomic targets, we cloned and sequenced 20 unique amplification products of two commonly used primers, CAA(CT)6 and (CA)8RY, using Biomphalaria glabrata genomic DNA as template. The sequences obtained were novel B. glabrata genomic sequences. It was observed that 15 clones contained microsatellites between priming sites. Out of 40 clones, seven contained complex sequence repetitions. One of the repeats that appeared in six of the amplified fragments generated a single band in Southern analysis, indicating that the sequence was not widespread in the genome. Most of the annealing sites for the CAA(CT)6 primer contained only the six repeats found within the primer sequence. In conclusion, SSR-PCR is a useful genotyping technique. However, the premise of the SSR-PCR technique, verified with the CAA(CT)6 primer, could not be supported since the amplification products did not result necessarily from microsatellite loci amplification.

Extreme genetic differentiation among the remnant populations of marble trout (Salmo marmoratus) in Slovenia.

Populations of the marble trout (Salmo marmoratus) have declined critically due to introgression by brown trout (Salmo trutta) strains. In order to define strategies for long-term conservation, we examined the genetic structure of the 8 known pure populations using 15 microsatellite loci. The analyses reveal extraordinarily strong genetic differentiation among populations separated by < 15 km, and extremely low levels of intrapopulation genetic variability. As natural recolonization seems highly unlikely, appropriate management and conservation strategies should comprise the reintroduction of pure populations from mixed stocks (translocation) to avoid further loss of genetic diversity.

Low relatedness among cooperatively breeding workers of the greenhead ant Rhytidoponera metallica

The greenhead ant Rhytidoponera metallica has long been recognized as posing a potential challenge to kin selection theory because it has large queenless colonies where apparently many of the morphological workers are mated and reproducing. However this species has never been studied genetically and important elements of its breeding system and kin structure remain uncertain. We used microsatellite markers to measure the relatedness among nestmates unravel the fine-scale population genetic structure and infer the breeding system of R. metallica. The genetic relatedness among worker nestmates is very low but significantly greater than zero (r = 0.082 +/- 0.015) which demonstrates that nests contain many distantly related breeders. The inbreeding coefficient is very close to and not significantly different from zero indicating random mating and lack of microgeographic genetic differentiation. On average. closely located nests are not more similar genetically than distant nests which is surprising as new colonies form by budding and female dispersal is restricted. Lack of inbreeding and absence of population viscosity indicates high gene flow mediated by males. Overall the genetic pattern detected in R. metallica suggests that a high number of moderately related workers mate with unrelated males from distant nests. This breeding system results in the lowest relatedness among nestmates reported for social insect species where breeders and helpers are not morphologically differentiated. [References: 69]

Mutation screening of patients with Alzheimer disease identifies APP locus duplication in a Swedish patient.

BACKGROUND: Missense mutations in three different genes encoding amyloid-β precursor protein, presenilin 1 and presenilin 2 are recognized to cause familial early-onset Alzheimer disease. Also duplications of the amyloid precursor protein gene have been shown to cause the disease. At the Dept. of Geriatric Medicine, Karolinska University Hospital, Sweden, patients are referred for mutation screening for the identification of nucleotide variations and for determining copy-number of the APP locus. METHODS: We combined the method of microsatellite marker genotyping with a quantitative real-time PCR analysis to detect duplications in patients with Alzheimer disease. RESULTS: In 22 DNA samples from individuals diagnosed with clinical Alzheimer disease, we identified one patient carrying a duplication on chromosome 21 which included the APP locus. Further mapping of the chromosomal region by array-comparative genome hybridization showed that the duplication spanned a maximal region of 1.09 Mb. CONCLUSIONS: This is the first report of an APP duplication in a Swedish Alzheimer patient and describes the use of quantitative real-time PCR as a tool for determining copy-number of the APP locus.

Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations.

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.

Communally breeding Bechstein's bats have a stable social system that is independent from the postglacial history and location of the populations.

Investigating macro-geographical genetic structures of animal populations is crucial to reconstruct population histories and to identify significant units for conservation. This approach may also provide information about the intraspecific flexibility of social systems. We investigated the history and current structure of a large number of populations in the communally breeding Bechstein's bat (Myotis bechsteinii). Our aim was to understand which factors shape the species' social system over a large ecological and geographical range. Using sequence data from one coding and one noncoding mitochondrial DNA region, we identified the Balkan Peninsula as the main and probably only glacial refugium of the species in Europe. Sequence data also suggest the presence of a cryptic taxon in the Caucasus and Anatolia. In a second step, we used seven autosomal and two mitochondrial microsatellite loci to compare population structures inside and outside of the Balkan glacial refugium. Central European and Balkan populations both were more strongly differentiated for mitochondrial DNA than for nuclear DNA, had higher genetic diversities and lower levels of relatedness at swarming (mating) sites than in maternity (breeding) colonies, and showed more differentiation between colonies than between swarming sites. All these suggest that populations are shaped by strong female philopatry, male dispersal, and outbreeding throughout their European range. We conclude that Bechstein's bats have a stable social system that is independent from the postglacial history and location of the populations. Our findings have implications for the understanding of the benefits of sociality in female Bechstein's bats and for the conservation of this endangered species.

Sex-specific selective pressures on body mass in the greater white-toothed shrew, Crocidura russula.

The direction, intensity and shape of viability-, sexual- and fecundity selection on body mass were investigated in a natural population of the greater white-toothed shrew (Crocidura russula), combining parentage assignment through molecular techniques and mark-recapture data over several generations. A highly significant stabilizing viability selection was found in both sexes, presumably stemming from the constraints imposed by their insectivorous habits and high metabolic costs. Sexual selection, directional in both sexes, was twice as large in males than in females. Our results suggest that body mass matters in this context by facilitating the acquisition and defense of a breeding territory. No fecundity selection could be detected. The direction of sexual size dimorphism (SSD) was in agreement with the observed pattern of selective pressures: males were heavier than females, because of stronger sexual selection. SSD intensity, however, was low compared with other mammals, because of the low level of polygyny, the active role of females in territory defense and the intensity of stabilizing viability selection.

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia.

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F(ST)=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5-30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.

Testing advances in molecular discrimination among Chinook salmon life histories: evidence from a blind test.

The application of DNA-based markers toward the task of discriminating among alternate salmon runs has evolved in accordance with ongoing genomic developments and increasingly has enabled resolution of which genetic markers associate with important life-history differences. Accurate and efficient identification of the most likely origin for salmon encountered during ocean fisheries, or at salvage from fresh water diversion and monitoring facilities, has far-reaching consequences for improving measures for management, restoration and conservation. Near-real-time provision of high-resolution identity information enables prompt response to changes in encounter rates. We thus continue to develop new tools to provide the greatest statistical power for run identification. As a proof of concept for genetic identification improvements, we conducted simulation and blind tests for 623 known-origin Chinook salmon (Oncorhynchus tshawytscha) to compare and contrast the accuracy of different population sampling baselines and microsatellite loci panels. This test included 35 microsatellite loci (1266 alleles), some known to be associated with specific coding regions of functional significance, such as the circadian rhythm cryptochrome genes, and others not known to be associated with any functional importance. The identification of fall run with unprecedented accuracy was demonstrated. Overall, the top performing panel and baseline (HMSC21) were predicted to have a success rate of 98%, but the blind-test success rate was 84%. Findings for bias or non-bias are discussed to target primary areas for further research and resolution.

Higher genetic variation estimated by microsatellites compared to isoenzyme markers in Aedes aegypti from Rio de Janeiro

Aedes aegypti populations from five districts in Rio de Janeiro were analyzed using five microsatellites and six isoenzyme markers, to assess the amount of variation and patterns of gene flow at local levels. Microsatellite loci were polymorphic enough to detect genetic differentiation of populations collected at small geographic scales (e.g. within a city). Ae. aegypti populations were highly differentiated as well in the city center as in the outskirt. Thus, dengue virus propagation by mosquitoes could be as efficient in the urban area as in the outskirt of Rio de Janeiro, the main entry point of dengue in Brazil.

A large study of androgen receptor germline variants and their relation to sex hormone levels and prostate cancer risk. Results from the National Cancer Institute Breast and Prostate Cancer Cohort Consortium

Background: Androgens are key regulators of prostate gland maintenance and prostate cancer growth, and androgen deprivation therapy has been the mainstay of treatment for advanced prostate cancer for many years. A long-standing hypothesis has been that inherited variation in the androgen receptor (AR) gene plays a role in prostate cancer initiation. However, studies to date have been inconclusive and often suffered from small sample sizes. Objective and Methods: We investigated the association of AR sequence variants with circulating sex hormone levels and prostate cancer risk in 6058 prostate cancer cases and 6725 controls of Caucasian origin within the Breast and Prostate Cancer Cohort Consortium. We genotyped a highly polymorphic CAG microsatellite in exon 1 and six haplotype tagging single nucleotide polymorphisms and tested each genetic variant for association with prostate cancer risk and with sex steroid levels. Results: We observed no association between AR genetic variants and prostate cancer risk. However, there was a strong association between longer CAG repeats and higher levels of testosterone (P = 4.73 × 10−5) and estradiol (P = 0.0002), although the amount of variance explained was small (0.4 and 0.7%, respectively). Conclusions: This study is the largest to date investigating AR sequence variants, sex steroid levels, and prostate cancer risk. Although we observed no association between AR sequence variants and prostate cancer risk, our results support earlier findings of a relation between the number of CAG repeats and circulating levels of testosterone and estradiol.