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.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Copy number variation and chromatin structure

The functional consequences of structural variation in the human genome range from adaptation, to phenotypic variation, to predisposition to diseases. Copy number variation (CNV) was shown to influence the phenotype by modifying, in a somewhat dose-dependent manner, the expression of genes that map within them, as well as that of genes located on their flanks. To assess the possible mechanism(s) behind this neighboring effect, we compared histone modification status of cell lines from patients affected by Williams-Beuren, Williams-Beuren region duplication, Smith-Magenis or DiGeorge Syndrome and control individuals using a high-throughput version of chromatin immuno-precipitation method (ChIP), called ChlP-seq. We monitored monomethylation of lysine K20 on histone H4 and trimethylation of lysine K27 on histone H3, as proxies for open and condensed chromatin, respectively. Consistent with the changes in expression levels observed for multiple genes mapping on the entire length of chromosomes affected by structural variants, we also detected regions with modified histone status between samples, up- and downstream from the critical regions, up to the end of the rearranged chromosome. We also gauged the intrachromosomal interactions of these cell lines utilizing chromosome conformation capture (4C-seq) technique. We observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together, possibly forming an interacting cluster with each other and the WBSCR. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We conclude, that large genomic rearrangements can lead to changes in the state of the chromatin spreading far away from the critical region, thus possibly affecting expression globally and as a result modifying the phenotype of the patients. - Les conséquences fonctionnelles des variations structurelles dans le génome humain sont vastes, allant de l'adaptation, en passant par les variations phénotypiques, aux prédispositions à certaines maladies. Il a été démontré que les variations du nombre de copies (CNV) influencent le phénotype en modifiant, d'une manière plus ou moins dose-dépendante, l'expression des gènes se situant à l'intérieur de ces régions, mais également celle des gènes se trouvant dans les régions flanquantes. Afin d'étudier les mécanismes possibles sous-jacents à cet effet de voisinage, nous avons comparé les états de modification des histones dans des lignées cellulaires dérivées de patients atteints du syndrome de Williams-Beuren, de la duplication de la région Williams-Beuren, du syndrome de Smith-Magenis ou du syndrome de Di- George et d'individus contrôles en utilisant une version haut-débit de la méthode d'immunoprécipitation de la chromatine (ChIP), appelée ChIP-seq. Nous avons suivi la mono-méthylation de la lysine K20 sur l'histone H4 et la tri-méthylation de la lysine K27 sur l'histone H3, marqueurs respectifs de la chromatine ouverte et fermée. En accord avec les changements de niveaux d'expression observés pour de multiples gènes tout le long des chromosomes affectés par les CNVs, nous avons aussi détecté des régions présentant des modifications d'histones entre les échantillons, situées de part et d'autre des régions critiques, jusqu'aux extrémités du chromosome réarrangé. Nous avons aussi évalué les interactions intra-chromosomiques ayant lieu dans ces cellules par l'utilisation de la technique de capture de conformation des chromosomes (4C-seq). Nous avons observé qu'un groupe de gènes flanquants la région critique du syndrome de Williams-Beuren (WBSCR) forment souvent une boucle, constituant un groupe d'interactions privilégiées entre ces gènes et la WBSCR. La délétion de la WBSCR perturbe l'expression de ce groupe de gènes flanquants, mais également les interactions à grande échelle entre eux et la région réarrangée. Nous en concluons que les larges réarrangements génomiques peuvent aboutir à des changements de l'état de la chromatine pouvant s'étendre bien plus loin que la région critique, affectant donc potentiellement l'expression de manière globale et ainsi modifiant le phénotype des patients.

Standardisation of 18F by a coincidence method using full solid angle detectors.

A solution of (18)F was standardised with a 4pibeta-4pigamma coincidence counting system in which the beta detector is a one-inch diameter cylindrical UPS89 plastic scintillator, positioned at the bottom of a well-type 5''x5'' NaI(Tl) gamma-ray detector. Almost full detection efficiency-which was varied downwards electronically-was achieved in the beta-channel. Aliquots of this (18)F solution were also measured using 4pigamma NaI(Tl) integral counting and Monte Carlo calculated efficiencies as well as the CIEMAT-NIST method. Secondary measurements of the same solution were also performed with an IG11 ionisation chamber whose equivalent activity is traceable to the Système International de Référence through the contribution IRA-METAS made to it in 2001; IRA's degree of equivalence was found to be close to the key comparison reference value (KCRV). The (18)F activity predicted by this coincidence system agrees closely with the ionisation chamber measurement and is compatible within one standard deviation of the other primary measurements. This work demonstrates that our new coincidence system can standardise short-lived radionuclides used in nuclear medicine.

Number of transferred embryos: how to reduce multiple pregnancies.

Because the diagnostic tools for predicting whether an early cleavage stage embryo can lead to a viable pregnancy are still elusive, transfer of more than one embryo remains quite common. However, the only way to reduce multiple pregnancies, considered as the main adverse effect of assisted reproductive technology, is to transfer a single embryo. In countries such as Switzerland and Germany, the law allows cryopreservation only at the 2-pronuclear stage. This restricts considerably the possibility of selecting the embryos to be transferred. Therefore, a good cryopreservation program at the 2-pronuclear stage is an essential tool to optimize the efficiency of in vitro fertilization (IVF). We therefore recommend the Cumulated Singleton Delivery Rate (CUSIDERA) as a measure of standard IVF efficiency. This rate averages approximately 23.5% when calculated over the last 10 years in our unit and reaches a value above 35% for patients with more than 10 zygotes. Elective single-embryo transfers and the decrease of iatrogenic multiple pregnancies in IVF remain dependent on better prognostic tools for the appropriate selection of patients, gametes, and zygotes.

Chromosomal number aberrations and transformation in adult mouse retinal stem cells in vitro.

PURPOSE: The potential of stem cells (SCs) as a source for cell-based therapy on a wide range of degenerative diseases and damaged tissues such as retinal degeneration has been recognized. Generation of a high number of retinal stem cells (RSCs) in vitro would thus be beneficial for transplantation in the retina. However, as cells in prolonged cultivation may be unstable and thus have a risk of transformation, it is important to assess the stability of these cells. METHODS: Chromosomal aberrations were analyzed in mouse RSC lines isolated from adult and from postnatal day (PN)1 mouse retinas. Moreover, selected cell lines were tested for anchorage-dependent proliferation, and SCs were transplanted into immunocompromised mice to assess the possibility of transformation. RESULTS: Marked aneuploidy occurred in all adult cell lines, albeit to different degrees, and neonatal RSCs were the most stable and displayed a normal karyotype until at least passage 9. Of interest, the level of aneuploidy of adult RSCs did not necessarily correlate with cell transformation. Only the adult RSC lines passaged for longer periods and with a higher dilution ratio underwent transformation. Furthermore, we identified several cell cycle proteins that might support the continuous proliferation and transformation of the cells. CONCLUSIONS: Adult RSCs rapidly accumulated severe chromosomal aberrations during cultivation, which led to cell transformation in some cell lines. The culture condition plays an important role in supporting the selection and growth of transformed cells.

Gene copy number polymorphisms in an arbuscular mycorrhizal fungal population.

Gene copy number polymorphism was studied in a population of the arbuscular mycorrhizal fungus Glomus intraradices by using a quantitative PCR approach on four different genomic regions. Variation in gene copy number was found for a pseudogene and for three ribosomal genes, providing conclusive evidence for a widespread occurrence of macromutational events in the population.

Plastic and evolutionary responses of cell size and number to larval malnutrition in Drosophila melanogaster.

Both development and evolution under chronic malnutrition lead to reduced adult size in Drosophila. We studied the contribution of changes in size vs. number of epidermal cells to plastic and evolutionary reduction of wing size in response to poor larval food. We used flies from six populations selected for tolerance to larval malnutrition and from six unselected control populations, raised either under standard conditions or under larval malnutrition. In the control populations, phenotypic plasticity of wing size was mediated by both cell size and cell number. In contrast, evolutionary change in wing size, which was only observed as a correlated response expressed on standard food, was mediated entirely by reduction in cell number. Plasticity of cell number had been lost in the selected populations, and cell number did not differ between the sexes despite males having smaller wings. Results of this and other experimental evolution studies are consistent with the hypothesis that alleles which increase body size through prolonged growth affect wing size mostly via cell number, whereas alleles which increase size through higher growth rate do so via cell size.

Inference about the number of contributors to a DNA mixture: Comparative analyses of a Bayesian network approach and the maximum allele count method.

In the forensic examination of DNA mixtures, the question of how to set the total number of contributors (N) presents a topic of ongoing interest. Part of the discussion gravitates around issues of bias, in particular when assessments of the number of contributors are not made prior to considering the genotypic configuration of potential donors. Further complication may stem from the observation that, in some cases, there may be numbers of contributors that are incompatible with the set of alleles seen in the profile of a mixed crime stain, given the genotype of a potential contributor. In such situations, procedures that take a single and fixed number contributors as their output can lead to inferential impasses. Assessing the number of contributors within a probabilistic framework can help avoiding such complication. Using elements of decision theory, this paper analyses two strategies for inference on the number of contributors. One procedure is deterministic and focuses on the minimum number of contributors required to 'explain' an observed set of alleles. The other procedure is probabilistic using Bayes' theorem and provides a probability distribution for a set of numbers of contributors, based on the set of observed alleles as well as their respective rates of occurrence. The discussion concentrates on mixed stains of varying quality (i.e., different numbers of loci for which genotyping information is available). A so-called qualitative interpretation is pursued since quantitative information such as peak area and height data are not taken into account. The competing procedures are compared using a standard scoring rule that penalizes the degree of divergence between a given agreed value for N, that is the number of contributors, and the actual value taken by N. Using only modest assumptions and a discussion with reference to a casework example, this paper reports on analyses using simulation techniques and graphical models (i.e., Bayesian networks) to point out that setting the number of contributors to a mixed crime stain in probabilistic terms is, for the conditions assumed in this study, preferable to a decision policy that uses categoric assumptions about N.

Quantitative analysis of postnatal neurogenesis and neuron number in the macaque monkey dentate gyrus

The dentate gyrus is one of only two regions of the mammalian brain where substantial neurogenesis occurs postnatally. However, detailed quantitative information about the postnatal structural maturation of the primate dentate gyrus is meager. We performed design-based, stereological studies of neuron number and size, and volume of the dentate gyrus layers in rhesus macaque monkeys (Macaca mulatta) of different postnatal ages. We found that about 40% of the total number of granule cells observed in mature 5-10-year-old macaque monkeys are added to the granule cell layer postnatally; 25% of these neurons are added within the first three postnatal months. Accordingly, cell proliferation and neurogenesis within the dentate gyrus peak within the first 3 months after birth and remain at an intermediate level between 3 months and at least 1 year of age. Although granule cell bodies undergo their largest increase in size during the first year of life, cell size and the volume of the three layers of the dentate gyrus (i.e. the molecular, granule cell and polymorphic layers) continue to increase beyond 1 year of age. Moreover, the different layers of the dentate gyrus exhibit distinct volumetric changes during postnatal development. Finally, we observe significant levels of cell proliferation, neurogenesis and cell death in the context of an overall stable number of granule cells in mature 5-10-year-old monkeys. These data identify an extended developmental period during which neurogenesis might be modulated to significantly impact the structure and function of the dentate gyrus in adulthood.

Prevalence of methylphenidate prescription among school-aged children in a Swiss population: increase in the number of prescriptions in the Swiss Canton of Vaud, from 2002 to 2005, and changes in patient demographics.

OBJECTIVE: Methylphenidate is prescribed for children and adolescents to treat ADHD. As in many Western countries, the increase in methylphenidate consumption is a public concern in Switzerland. The article discusses the authors' assessment of prescription prevalence in 2002 and 2005 for school-aged children in the canton of Vaud. METHOD: Pharmacy prescription information is available from the regional public health authority. Descriptive analyses are conducted on an anonymized database of the years 2002 and 2005. Data for each year are compared to assess trends in methylphenidate prescription prevalence. RESULTS: The findings show an increase from 0.74% to 1.02% in the number of prescriptions for 5- to 14-year-old children, particularly in prescriptions for girls. Data also show important geographical differences in prescription. CONCLUSION: The prevalence of methylphenidate prescription is lower in Switzerland than other Western countries, particularly the United States. However, some aspects of prevalence are similar, including the increase per year, demographics, and geographic characteristics.

The average inter-crossing number of equilateral random walks and polygons

In this paper, we study the average inter-crossing number between two random walks and two random polygons in the three-dimensional space. The random walks and polygons in this paper are the so-called equilateral random walks and polygons in which each segment of the walk or polygon is of unit length. We show that the mean average inter-crossing number ICN between two equilateral random walks of the same length n is approximately linear in terms of n and we were able to determine the prefactor of the linear term, which is a = (3 In 2)/(8) approximate to 0.2599. In the case of two random polygons of length n, the mean average inter-crossing number ICN is also linear, but the prefactor of the linear term is different from that of the random walks. These approximations apply when the starting points of the random walks and polygons are of a distance p apart and p is small compared to n. We propose a fitting model that would capture the theoretical asymptotic behaviour of the mean average ICN for large values of p. Our simulation result shows that the model in fact works very well for the entire range of p. We also study the mean ICN between two equilateral random walks and polygons of different lengths. An interesting result is that even if one random walk (polygon) has a fixed length, the mean average ICN between the two random walks (polygons) would still approach infinity if the length of the other random walk (polygon) approached infinity. The data provided by our simulations match our theoretical predictions very well.

NADPH-diaphorase-positive ganglion cells of the rat adrenal gland: age- and sex-related changes in their number, size, and distribution.

The rat adrenal gland contains ganglion cells able to synthesize nitric oxide (NO). This messenger molecule controls and modulates adrenal secretory activity and blood flow. The present study analyzed the number, size, and distribution of NO-producing adrenal neurons in adulthood and during postnatal development by means of beta-nicotinamide adenine dinucleotide phosphate-diaphorase (NADPH-d) histochemistry. This method reliably visualizes the enzyme responsible for NO generation. The reactive neurons per adrenal gland were 350-400 in both male and female adult rats. The positive nerve cell bodies were mostly located in the medulla, few being detected within the cortex and the subcapsular region. Dual labeling with anti-microtubule-associated protein 2 antibody, specific for neuronal elements, confirmed this distribution. Anti-microtubule-associated protein 1b antibody identified a subset of NADPH-d-positive neurons, displaying different degrees of maturation according to their position within the adrenal gland. At birth, there were about 220 NADPH-d-labeled neurons per adrenal gland in both sexes. As confirmed by dual immunocytochemical labeling, their great majority was evenly distributed between the cortex and the subcapsular region, the medulla being practically devoid of stained neurons. After birth, the number of adrenal NADPH-d-positive ganglion cells displayed a strong postnatal increase and reached the adult-like distribution after 1-2 months. During the period of increase, there was a transient difference in the numbers of these cells in the two sexes. Thus we present here evidence of plasticity in the number, size, and distribution of NADPH-d-positive adrenal neurons between birth and adulthood; in addition, we describe transient sex-related differences in their number and distribution during the 2nd postnatal week, which are possibly related to the epigenetic action of gonadal hormones during this period.