Case studies and mathematical models of ecological speciation. 1. Cichlids in a crater lake.

A recent study of a pair of sympatric species of cichlids in Lake Apoyo in Nicaragua is viewed as providing probably one of the most convincing examples of sympatric speciation to date. Here, we describe and study a stochastic, individual-based, explicit genetic model tailored for this cichlid system. Our results show that relatively rapid (<20,000 generations) colonization of a new ecological niche and (sympatric or parapatric) speciation via local adaptation and divergence in habitat and mating preferences are theoretically plausible if: (i) the number of loci underlying the traits controlling local adaptation, and habitat and mating preferences is small; (ii) the strength of selection for local adaptation is intermediate; (iii) the carrying capacity of the population is intermediate; and (iv) the effects of the loci influencing nonrandom mating are strong. We discuss patterns and timescales of ecological speciation identified by our model, and we highlight important parameters and features that need to be studied empirically to provide information that can be used to improve the biological realism and power of mathematical models of ecological speciation.

A dual functional origin of transfer in the ICEclc genomic island of Pseudomonas knackmussii B13.

Genomic islands (GEIs) are large DNA segments, present in most bacterial genomes, that are most likely acquired via horizontal gene transfer. Here, we study the self-transfer system of the integrative and conjugative element ICEclc of Pseudomonas knackmussii B13, which stands model for a larger group of ICE/GEI with syntenic core gene organization. Functional screening revealed that unlike conjugative plasmids and other ICEs ICEclc carries two separate origins of transfer, with different sequence context but containing a similar repeat motif. Conjugation experiments with GFP-labelled ICEclc variants showed that both oriTs are used for transfer and with indistinguishable efficiencies, but that having two oriTs results in an estimated fourfold increase of ICEclc transfer rates in a population compared with having a single oriT. A gene for a relaxase essential for ICEclc transfer was also identified, but in vivo strand exchange assays suggested that the relaxase processes both oriTs in a different manner. This unique dual origin of transfer system might have provided an evolutionary advantage for distribution of ICE, a hypothesis that is supported by the fact that both oriT regions are conserved in several GEIs related to ICEclc.

Genome-wide association identifies nine common variants associated with fasting proinsulin levels and provides new insights into the pathophysiology of type 2 diabetes.

OBJECTIVE: Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. RESEARCH DESIGN AND METHODS: We have conducted a meta-analysis of genome-wide association tests of ∼2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. RESULTS: Nine SNPs at eight loci were associated with proinsulin levels (P < 5 × 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC30A8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 × 10(-4)), improved β-cell function (P = 1.1 × 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 × 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. CONCLUSIONS: We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis.

Fast and simple epidemiological typing of Pseudomonas aeruginosa using the double-locus sequence typing (DLST) method.

Although the molecular typing of Pseudomonas aeruginosa is important to understand the local epidemiology of this opportunistic pathogen, it remains challenging. Our aim was to develop a simple typing method based on the sequencing of two highly variable loci. Single-strand sequencing of three highly variable loci (ms172, ms217, and oprD) was performed on a collection of 282 isolates recovered between 1994 and 2007 (from patients and the environment). As expected, the resolution of each locus alone [number of types (NT) = 35-64; index of discrimination (ID) = 0.816-0.964] was lower than the combination of two loci (NT = 78-97; ID = 0.966-0.971). As each pairwise combination of loci gave similar results, we selected the most robust combination with ms172 [reverse; R] and ms217 [R] to constitute the double-locus sequence typing (DLST) scheme for P. aeruginosa. This combination gave: (i) a complete genotype for 276/282 isolates (typability of 98%), (ii) 86 different types, and (iii) an ID of 0.968. Analysis of multiple isolates from the same patients or taps showed that DLST genotypes are generally stable over a period of several months. The high typability, discriminatory power, and ease of use of the proposed DLST scheme makes it a method of choice for local epidemiological analyses of P. aeruginosa. Moreover, the possibility to give unambiguous definition of types allowed to develop an Internet database ( http://www.dlst.org ) accessible by all.

A Novel Homozygous R764H Mutation in CRB1 Causes Autosomal Recessive Retinitis Pigmentosa in a Consanguineous Family

Purpose: Retinitis pigmentosa (RP; MIM 268000) is a hereditary disease characterized by poor night vision and progressive loss of photoreceptors, eventually leading to blindness. This degenerative process primarily affects peripheral vision due to the loss of rods. Autosomal recessive RP (arRP) is clinically and genetically heterogeneous. It has been associated with mutations in different genes, including CRB1 (Crumbs homolog 1). The aim of this study was to determine the causative gene in a Tunisian patient with arRP born to non consanguineous parents.Methods: Four accessible family members were included. They underwent full ophthalmic examination with best corrected Snellen visual acuity, fundus photography and fluoroangiography. Haplotype analyses were used to test linkage in the family to 20 arRP loci, including ABCA4, LRAT, USH2A, RP29, CERKL, CNGA1, CNGB1, CRB1, EYS, RP28, MERTK, NR2E3, PDE6A, PDE6B, RGR, RHO, RLBP1, TULP1. All exons and intron-exon junctions of candidate genes not excluded by haplotype analysis were PCR amplified and directly sequenced.Results: A 39 aged affected member was individualized. Best corrected visual acuity was OR: 20/63, OS: 20/80. Visual loss began at the third decade. Funduscopic examination and FA revealed typical advanced RP changes with bone spicule-shaped pigment deposits in the posterior pole and the mild periphery along with retinal atrophy, narrowing of the vessels and waxy optic discs. Haplotypes analysis revealed homozygosity with microsatellites markers D1S412 and D1S413 on chromosome 1q31.3. These markers flanked the CRB1 gene. Our results excluded linkage of all the other arRP loci/ genes tested. Sequencing of the 12 coding exons and splice sites of CRB1 gene disclosed a homozygous missense mutation in exon 7 at nucleotide c.(2291 G>A), resulting in an Arg to Hist substitution (p.R764H).Conclusions: R764H is a novel mutation associated with CRB1-related arRP. Previously, an R764C mutation was observed. Extending the mutation spectrum of CRB1 with additional families is important for genotype-phenotype correlations.

Clear detection of ADIPOQ locus as the major gene for plasma adiponectin: results of genome-wide association analyses including 4659 European individuals.

OBJECTIVE: Plasma adiponectin is strongly associated with various components of metabolic syndrome, type 2 diabetes and cardiovascular outcomes. Concentrations are highly heritable and differ between men and women. We therefore aimed to investigate the genetics of plasma adiponectin in men and women. METHODS: We combined genome-wide association scans of three population-based studies including 4659 persons. For the replication stage in 13795 subjects, we selected the 20 top signals of the combined analysis, as well as the 10 top signals with p-values less than 1.0 x 10(-4) for each the men- and the women-specific analyses. We further selected 73 SNPs that were consistently associated with metabolic syndrome parameters in previous genome-wide association studies to check for their association with plasma adiponectin. RESULTS: The ADIPOQ locus showed genome-wide significant p-values in the combined (p=4.3 x 10(-24)) as well as in both women- and men-specific analyses (p=8.7 x 10(-17) and p=2.5 x 10(-11), respectively). None of the other 39 top signal SNPs showed evidence for association in the replication analysis. None of 73 SNPs from metabolic syndrome loci exhibited association with plasma adiponectin (p>0.01). CONCLUSIONS: We demonstrated the ADIPOQ gene as the only major gene for plasma adiponectin, which explains 6.7% of the phenotypic variance. We further found that neither this gene nor any of the metabolic syndrome loci explained the sex differences observed for plasma adiponectin. Larger studies are needed to identify more moderate genetic determinants of plasma adiponectin.

Evolution of uni- and bifactorial sexual compatibility systems in fungi.

Mating systems, that is, whether organisms give rise to progeny by selfing, inbreeding or outcrossing, strongly affect important ecological and evolutionary processes. Large variations in mating systems exist in fungi, allowing the study of their origin and consequences. In fungi, sexual incompatibility is determined by molecular recognition mechanisms, controlled by a single mating-type locus in most unifactorial fungi. In Basidiomycete fungi, however, which include rusts, smuts and mushrooms, a system has evolved in which incompatibility is controlled by two unlinked loci. This bifactorial system probably evolved from a unifactorial system. Multiple independent transitions back to a unifactorial system occurred. It is still unclear what force drove evolution and maintenance of these contrasting inheritance patterns that determine mating compatibility. Here, we give an overview of the evolutionary factors that might have driven the evolution of bifactoriality from a unifactorial system and the transitions back to unifactoriality. Bifactoriality most likely evolved for selfing avoidance. Subsequently, multiallelism at mating-type loci evolved through negative frequency-dependent selection by increasing the chance to find a compatible mate. Unifactoriality then evolved back in some species, possibly because either selfing was favoured or for increasing the chance to find a compatible mate in species with few alleles. Owing to the existence of closely related unifactorial and bifactorial species and the increasing knowledge of the genetic systems of the different mechanisms, Basidiomycetes provide an excellent model for studying the different forces that shape breeding systems.

Evolutionary history of the greater white-toothed shrew (Crocidura russula) inferred from analysis of mtDNA, Y, and X chromosome markers.

We investigate the evolutionary history of the greater white-toothed shrew across its distribution in northern Africa and mainland Europe using sex-specific (mtDNA and Y chromosome) and biparental (X chromosome) markers. All three loci confirm a large divergence between eastern (Tunisia and Sardinia) and western (Morocco and mainland Europe) lineages, and application of a molecular clock to mtDNA divergence estimates indicates a more ancient separation (2.25 M yr ago) than described by some previous studies, supporting claims for taxonomic revision. Moroccan ancestry for the mainland European population is inconclusive from phylogenetic trees, but is supported by greater nucleotide diversity and a more ancient population expansion in Morocco than in Europe. Signatures of rapid population expansion in mtDNA, combined with low X and Y chromosome diversity, suggest a single colonization of mainland Europe by a small number of Moroccan shrews >38 K yr ago. This study illustrates that multilocus genetic analyses can facilitate the interpretation of species' evolutionary history but that phylogeographic inference using X and Y chromosomes is restricted by low levels of observed polymorphism.

Genetic polymorphisms of the main transcription factors for adiponectin gene promoter in regulation of adiponectin levels: association analysis in three European cohorts.

Adiponectin serum concentrations are an important biomarker in cardiovascular epidemiology with heritability etimates of 30-70%. However, known genetic variants in the adiponectin gene locus (ADIPOQ) account for only 2%-8% of its variance. As transcription factors are thought to play an under-acknowledged role in carrying functional variants, we hypothesized that genetic polymorphisms in genes coding for the main transcription factors for the ADIPOQ promoter influence adiponectin levels. Single nucleotide polymorphisms (SNPs) at these genes were selected based on the haplotype block structure and previously published evidence to be associated with adiponectin levels. We performed association analyses of the 24 selected SNPs at forkhead box O1 (FOXO1), sterol-regulatory-element-binding transcription factor 1 (SREBF1), sirtuin 1 (SIRT1), peroxisome-proliferator-activated receptor gamma (PPARG) and transcription factor activating enhancer binding protein 2 beta (TFAP2B) gene loci with adiponectin levels in three different European cohorts: SAPHIR (n = 1742), KORA F3 (n = 1636) and CoLaus (n = 5355). In each study population, the association of SNPs with adiponectin levels on log-scale was tested using linear regression adjusted for age, sex and body mass index, applying both an additive and a recessive genetic model. A pooled effect size was obtained by meta-analysis assuming a fixed effects model. We applied a significance threshold of 0.0033 accounting for the multiple testing situation. A significant association was only found for variants within SREBF1 applying an additive genetic model (smallest p-value for rs1889018 on log(adiponectin) = 0.002, β on original scale = -0.217 µg/ml), explaining ∼0.4% of variation of adiponectin levels. Recessive genetic models or haplotype analyses of the FOXO1, SREBF1, SIRT1, TFAPB2B genes or sex-stratified analyses did not reveal additional information on the regulation of adiponectin levels. The role of genetic variations at the SREBF1 gene in regulating adiponectin needs further investigation by functional studies.

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Eficiência de índices fenotípicos de comprimento de raiz seminal na avaliação de plantas individuais de milho quanto à tolerância ao alumínio

O objetivo deste trabalho foi testar a eficiência do comprimento relativo de raiz seminal (CRRS) e do comprimento líquido de raiz seminal (CLRS) como indicadores fenotípicos quanto à tolerância ao alumínio na avaliação de plantas individuais de milho. Plântulas de genótipos tolerantes e suscetíveis ao Al foram submetidas a solução nutritiva contendo nível tóxico deste elemento, por um período de sete dias, após o qual, determinaram-se os valores de CRRS e CLRS. Os resultados obtidos quando se utilizaram valores médios para CRRS e CLRS mostraram que ambos os índices foram capazes de discriminar com eficiência os materiais tolerantes dos suscetíveis. Entretanto, quando foram utilizados os valores de CRRS e CLRS obtidos a partir de plantas individuais, observou-se a existência de plantas tolerantes com valores típicos de plantas suscetíveis, o que indica que a avaliação fenotípica de plantas individuais pelos dois índices está sujeita a erros significativos, principalmente na caracterização de plantas suscetíveis. Portanto, em estudos para mapeamento de "quantitative trait loci" (QTLs) ligados à tolerância ao Al, nos quais utilizam-se estes índices fenotípicos, o mais apropriado é avaliar famílias F3, onde é possível obter valores médios para CRRS e CLRS e utilizar estas médias para representar os valores fenotípicos das respectivas plantas F2.

Inferring recent migration rates from individual genotypes.

We present a novel and straightforward method for estimating recent migration rates between discrete populations using multilocus genotype data. The approach builds upon a two-step sampling design, where individual genotypes are sampled before and after dispersal. We develop a model that estimates all pairwise backwards migration rates (m(ij), the probability that an individual sampled in population i is a migrant from population j) between a set of populations. The method is validated with simulated data and compared with the methods of BayesAss and Structure. First, we use data for an island model and then we consider more realistic data simulations for a metapopulation of the greater white-toothed shrew (Crocidura russula). We show that the precision and bias of estimates primarily depend upon the proportion of individuals sampled in each population. Weak sampling designs may particularly affect the quality of the coverage provided by 95% highest posterior density intervals. We further show that it is relatively insensitive to the number of loci sampled and the overall strength of genetic structure. The method can easily be extended and makes fewer assumptions about the underlying demographic and genetic processes than currently available methods. It allows backwards migration rates to be estimated across a wide range of realistic conditions.

Opportunity for sexual selection and effective population size in the lek-breeding European treefrog (Hyla arborea).

Sexual selection in lek-breeding species might drastically lower male effective population size, with potentially important consequences for evolutionary and conservation biology. Using field-monitoring and parental-assignment methods, we analyzed sex-specific variances in breeding success in a population of European treefrogs, to (1) help understanding the dynamics of genetic variance at sex-specific loci, and (2) better quantify the risk posed by genetic drift in this species locally endangered by habitat fragmentation. The variance in male mating success turned out to be markedly lower than values obtained from other amphibian species with polygamous mating systems. The ratio of effective breeding size to census breeding size was only slightly lower in males (0.44) than in females (0.57), in line with the patterns of genetic diversity previously reported from H. arborea sex chromosomes. Combining our results with data on age at maturity and adult survival, we show that the negative effect of the mating system is furthermore compensated by the effect of delayed maturity, so that the estimated instantaneous effective size broadly corresponded to census breeding size. We conclude that the lek-breeding system of treefrogs impacts only weakly the patterns of genetic diversity on sex-linked genes and the ability of natural populations to resist genetic drift.

Improved multiple-locus variable-number tandem-repeat assay for Staphylococcus aureus genotyping, providing a highly informative technique together with strong phylogenetic value.

We describe an improved multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) scheme for genotyping Staphylococcus aureus. We compare its performance to those of multilocus sequence typing (MLST) and spa typing in a survey of 309 strains. This collection includes 87 epidemic methicillin-resistant S. aureus (MRSA) strains of the Harmony collection, 75 clinical strains representing the major MLST clonal complexes (CCs) (50 methicillin-sensitive S. aureus [MSSA] and 25 MRSA), 135 nasal carriage strains (133 MSSA and 2 MRSA), and 13 published S. aureus genome sequences. The results show excellent concordance between the techniques' results and demonstrate that the discriminatory power of MLVA is higher than those of both MLST and spa typing. Two hundred forty-two genotypes are discriminated with 14 VNTR loci (diversity index, 0.9965; 95% confidence interval, 0.9947 to 0.9984). Using a cutoff value of 45%, 21 clusters are observed, corresponding to the CCs previously defined by MLST. The variability of the different tandem repeats allows epidemiological studies, as well as follow-up of the evolution of CCs and the identification of potential ancestors. The 14 loci can conveniently be analyzed in two steps, based upon a first-line simplified assay comprising a subset of 10 loci (panel 1) and a second subset of 4 loci (panel 2) that provides higher resolution when needed. In conclusion, the MLVA scheme proposed here, in combination with available on-line genotyping databases (including http://mlva.u-psud.fr/), multiplexing, and automatic sizing, can provide a basis for almost-real-time large-scale population monitoring of S. aureus.