The fixation of locally beneficial alleles in a metapopulation.

Extinction, recolonization, and local adaptation are common in natural spatially structured populations. Understanding their effect upon genetic variation is important for systems such as genetically modified organism management or avoidance of drug resistance. Theoretical studies on the effect of extinction and recolonization upon genetic variance started appearing in the 1970s, but the role of local adaptation still has no good theoretical basis. Here we develop a model of a haploid species in a metapopulation in which a locally adapted beneficial allele is introduced. We study the effect of different spatial patterns of local adaptation, and different metapopulation dynamics, upon the fixation probability of the beneficial allele. Controlling for the average selection pressure, we find that a small area of positive selection can significantly increase the global probability of fixation. However, local adaptation becomes less important as extinction rate increases. Deme extinction and recolonization have a spatial smoothing effect that effectively reduces spatial variation in fitness.

Importance of IL28B gene polymorphisms in hepatitis C virus genotype 2 and 3 infected patients.

BACKGROUND & AIMS: Genetic variation in the interleukin 28B (IL28B) gene has been associated with the response to interferon-alfa/ribavirin therapy in hepatitis C virus (HCV) genotype 1-infected patients. The importance of three IL28B single nucleotide polymorphisms (rs8099917, rs12980275 and rs12979860) for HCV genotype 2/3-infected patients is unknown. METHODS: In patients with chronic hepatitis C genotype 2/3 (n=267), IL28B host genotypes (rs8099917, rs12980275 and rs12979860) were analyzed for associations with sustained virologic response (SVR) to antiviral therapy with (pegylated) interferon-alfa and ribavirin and with respect to epidemiological, biochemical, and virological parameters. For comparison, hepatitis C genotype 1 patients (n=378) and healthy controls (n=200) were included. RESULTS: The rs12979860 CC genotype, lower age, and genotype 2 were significantly associated with SVR in HCV genotype 2/3-infected patients (p=0.01, p=0.03 and p=0.03, respectively). No association was observed for rs8099917 and rs12980275. In addition, an SVR in patients with rapid virologic response (RVR) was associated with the rs12979860 CC genotype (p=0.05), while for non-RVR no association was found. Furthermore, a significant association with a higher baseline viral load was observed for all three IL28B genotypes in genotype 1/2/3-infected patients. Finally, increasing frequencies of the rs12979860 CC genotypes were observed in genotype 1- (33.9%), genotype 3- (38.9%), and genotype 2-infected (51.9%) patients in comparison with healthy controls (49.0%) (p<0.01). CONCLUSIONS: In genotype 2/3-infected patients, rs12979860 was significantly associated with SVR. The frequency of the rs12979860 CC genotype is lower in HCV genotype 1 vs. genotype 2/3 patients. All major IL28B genotypes are associated with HCV-RNA concentration.

Licorice-induced hypertension and common variants of genes regulating renal sodium reabsorption.

AIM: To study if gene alterations affecting renal sodium reabsorption associate with susceptibility to licorice-induced hypertension.METHODS: Finnish subjects (n = 30) with a previously documented incident of licorice-induced hypertension were recruited for the study using a newspaper announcement. Their previous clinical and family histories as well as serum electrolyte levels were examined. DNA samples from all individuals were screened for variants of the genes encoding 11beta-hydroxysteroid dehydrogenase type 2 (11betaHSD2) and alpha-, beta-, and gamma-subunits of the epithelial sodium channel (ENaC).RESULTS: Upon licorice predisposition, the patients had a mean blood pressure of 201/118 mmHg. Circulating potassium, renin, and aldosterone levels were low. No significant DNA variations were identified in the 11betaHSD2 gene. Four subjects were heterozygous for beta- and gammaENaC variants previously shown to be associated with hypertension. Furthermore, a novel G insertion (2004-2005insG) in the SCNN1A gene encoding the alphaENaC was identified in two subjects. The frequency of these ENaC variants was significantly higher in subjects with licorice-induced hypertension (6/30 i.e. 20%) than in blood donors (11/301 i.e. 3.7%, P = 0.002).CONCLUSIONS: Defects of the 11betaHSD2 gene do not constitute a likely cause for licorice-induced hypertension. Variants of the ENaC subunits may render some individuals sensitive to licorice-induced metabolic alterations and hypertension.

The population genetics of clonal and partially clonal diploids.

The consequences of variable rates of clonal reproduction on the population genetics of neutral markers are explored in diploid organisms within a subdivided population (island model). We use both analytical and stochastic simulation approaches. High rates of clonal reproduction will positively affect heterozygosity. As a consequence, nearly twice as many alleles per locus can be maintained and population differentiation estimated as F(ST) value is strongly decreased in purely clonal populations as compared to purely sexual ones. With increasing clonal reproduction, effective population size first slowly increases and then points toward extreme values when the reproductive system tends toward strict clonality. This reflects the fact that polymorphism is protected within individuals due to fixed heterozygosity. Contrarily, genotypic diversity smoothly decreases with increasing rates of clonal reproduction. Asexual populations thus maintain higher genetic diversity at each single locus but a lower number of different genotypes. Mixed clonal/sexual reproduction is nearly indistinguishable from strict sexual reproduction as long as the proportion of clonal reproduction is not strongly predominant for all quantities investigated, except for genotypic diversities (both at individual loci and over multiple loci).

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.

'Good-genes' and 'compatible-genes' effects in an Alpine whitefish and the information content of breeding tubercles over the course of the spawning season.

Some models of sexual selection predict that individuals vary in their genetic quality and reveal some of this variation in their secondary sexual characteristics. Alpine whitefish (Coregonus sp.) develop breeding tubercles shortly before their spawning season. These tubercles are epidermal structures that are distributed regularly along the body sides of both males and females. There is still much unexplained variation in the size of breeding tubercles within both sexes and with much overlap between the sexes. It has been suggested that breeding tubercles function to maintain body contact between the mating partners during spawning, act as weapons for defence of spawning territories, or are sexual signals that reveal aspects of genetic quality. We took two samples of whitefish from their spawning place, one at the beginning and one around the peak of spawning season. We found that females have on average smaller breeding tubercles than males, and that tubercle size partly reveals the stage of gonad maturation. Two independent full-factorial breeding experiments revealed that embryo mortality was significantly influenced by male and female effects. This finding demonstrates that the males differed in their genetic quality (because offspring get nothing but genes from their fathers). Tubercle size was negatively linked to some aspects of embryo mortality in the first breeding experiment but not significantly so in the second. This lack of consistency adds to inconsistent results that were reported before and suggests that (i) some aspects of genetic quality are not revealed in breeding tubercles while others are, or (ii) individuals vary in their signaling strategies and the information content of breeding tubercles is not always reliable. Moreover, the fact that female whitefish have breeding tubercles of significant size while males seem to have few reasons to be choosy suggests that the tubercles might also serve some functions that are not linked to sexual signaling.

Reduced IFNλ4 activity is associated with improved HCV clearance and reduced expression of interferon-stimulated genes

Hepatitis C virus (HCV) infections are the major cause of chronic liver disease, cirrhosis and hepatocellular carcinoma worldwide. Both spontaneous and treatment-induced clearance of HCV depend on genetic variation within the interferon-lambda locus, but until now no clear causal relationship has been established. Here we demonstrate that an amino-acid substitution in the IFNλ4 protein changing a proline at position 70 to a serine (P70S) substantially alters its antiviral activity. Patients harbouring the impaired IFNλ4-S70 variant display lower interferon-stimulated gene (ISG) expression levels, better treatment response rates and better spontaneous clearance rates, compared with patients coding for the fully active IFNλ4-P70 variant. Altogether, these data provide evidence supporting a role for the active IFNλ4 protein as the driver of high hepatic ISG expression as well as the cause of poor HCV clearance.

Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. Results: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays.

A Role for cis Interaction between the Inhibitory Ly49A receptor and MHC class I for natural killer cell education.

Natural killer (NK) cells show enhanced functional competence when they express inhibitory receptors specific for inherited major histocompatibility complex class I (MHC-I) molecules. Current models imply that NK cell education requires an interaction of inhibitory receptors with MHC-I expressed on other cells. However, the inhibitory Ly49A receptor can also bind MHC-I ligand on the NK cell itself (in cis). Here we describe a Ly49A variant, which can engage MHC-I expressed on other cells but not in cis. Even though this variant inhibited NK cell effector function, it failed to educate NK cells. The association with MHC-I in cis sequestered wild-type Ly49A, and this was found to relieve NK cells from a suppressive effect of unengaged Ly49A. These data explain how inhibitory MHC-I receptors can facilitate NK cell activation. They dissociate classical inhibitory from educating functions of Ly49A and suggest that cis interaction of Ly49A is necessary for NK cell education.

Inbreeding and sex-biased gene flow in the ant Formica exsecta.

The objective of this study was to assess breeding and dispersal patterns of both males and females in a monogyne (a single queen per colony) population of ants. Monogyny is commonly associated with extensive nuptial flights, presumably leading to considerable gene flow over large areas. Opposite to these expectations we found evidence of both inbreeding and sex-biased gene flow in a monogyne population of Formica exsecta. We found a significant degree of population subdivision at a local scale (within islands) for queens (females heading established colonies) and workers, but not for colony fathers (the males mated to the colony queens). However, we found little evidence of population subdivision at a larger scale (among islands). More conclusive support for sex-biased gene flow comes from the analysis of isolation by distance on the largest island, and from assignment tests revealing differences in female and male philopatry. The genetic similarity between pairs of queens decreased significantly when geographical distance increased, demonstrating limited dispersal and isolation by distance in queens. By contrast, we found no such pattern for colony fathers. Furthermore, a significantly greater fraction of colony queens were assigned as having originated from the population of residence, as compared to colony fathers. Inbreeding coefficients were significantly positive for workers, but not for mother queens. The queen-male relatedness coefficient of 0.23 (regression relatedness) indicates that mating occurs between fairly close relatives. These results suggest that some monogyne species of ants have complex dispersal and mating systems that can result in genetic isolation by distance over small geographical scales. More generally, this study also highlights the importance of identifying the relevant scale in analyses of population structure and dispersal.

Mobile gibberellin directly stimulates Arabidopsis hypocotyl xylem expansion.

Secondary growth of the vasculature results in the thickening of plant structures and continuously produces xylem tissue, the major biological carbon sink. Little is known about the developmental control of this quantitative trait, which displays two distinct phases in Arabidopsis thaliana hypocotyls. The later phase of accelerated xylem expansion resembles the secondary growth of trees and is triggered upon flowering by an unknown, shoot-derived signal. We found that flowering-dependent hypocotyl xylem expansion is a general feature of herbaceous plants with a rosette growth habit. Flowering induction is sufficient to trigger xylem expansion in Arabidopsis. By contrast, neither flower formation nor elongation of the main inflorescence is required. Xylem expansion also does not depend on any particular flowering time pathway or absolute age. Through analyses of natural genetic variation, we found that ERECTA acts locally to restrict xylem expansion downstream of the gibberellin (GA) pathway. Investigations of mutant and transgenic plants indicate that GA and its signaling pathway are both necessary and sufficient to directly trigger enhanced xylogenesis. Impaired GA signaling did not affect xylem expansion systemically, suggesting that it acts downstream of the mobile cue. By contrast, the GA effect was graft transmissible, suggesting that GA itself is the mobile shoot-derived signal.

Genome-wide association study to identify common variants associated with brachial circumference: a meta-analysis of 14 cohorts.

Brachial circumference (BC), also known as upper arm or mid arm circumference, can be used as an indicator of muscle mass and fat tissue, which are distributed differently in men and women. Analysis of anthropometric measures of peripheral fat distribution such as BC could help in understanding the complex pathophysiology behind overweight and obesity. The purpose of this study is to identify genetic variants associated with BC through a large-scale genome-wide association scan (GWAS) meta-analysis. We used fixed-effects meta-analysis to synthesise summary results across 14 GWAS discovery and 4 replication cohorts comprising overall 22,376 individuals (12,031 women and 10,345 men) of European ancestry. Individual analyses were carried out for men, women, and combined across sexes using linear regression and an additive genetic model: adjusted for age and adjusted for age and BMI. We prioritised signals for follow-up in two-stages. We did not detect any signals reaching genome-wide significance. The FTO rs9939609 SNP showed nominal evidence for association (p<0.05) in the age-adjusted strata for men and across both sexes. In this first GWAS meta-analysis for BC to date, we have not identified any genome-wide significant signals and do not observe robust association of previously established obesity loci with BC. Large-scale collaborations will be necessary to achieve higher power to detect loci underlying BC.

Typing Candida albicans oral isolates from human immunodeficiency virus-infected patients by multilocus enzyme electrophoresis and DNA fingerprinting.

A total of 189 Candida albicans isolates have been typed by multilocus enzyme electrophoresis. The results obtained confirm the clonal mode of reproduction of C. albicans. The C. albicans populations found in the oropharynx of human immunodeficiency virus (HIV)-infected patients, in the oropharynx of healthy carriers, or in association with invasive candidiasis could not be distinguished. No clone or group of clones could be associated with the appearance of clinical disorders or with a reduced in vitro susceptibility to the antifungal agent fluconazole. Multiple and sequential oral isolates from 24 HIV-infected patients were also typed by restriction enzyme analysis with the enzymes EcoRI and HinfI and by use of the Ca3 repetitive probe. The results obtained by the combination of all three typing methods show that all but one patient each carried a unique major C. albicans clone in their oropharynx. The 21 patients with sequential isolates had the same C. albicans clones in their throats during recurrent oropharyngeal candidiasis episodes, independently of clinical status or of changes of in vitro susceptibility to fluconazole. Finally, several isolates of the same C. albicans clone found simultaneously in the oropharynx of a patient may present different levels of susceptibility to fluconazole.

Infrequent transmission of HIV-1 drug-resistant variants.

Transmission of drug-resistant variants is influenced by several factors, including the prevalence of drug resistance in the population of HIV-1-infected patients, HIV-1 RNA levels and transmission by recently infected patients. In order to evaluate the impact of these factors on the transmission of drug-resistant variants, we have defined the population of potential transmitters and compared their resistance profiles to those of newly infected patients. Sequencing of pol gene was performed in 220 recently infected patients and in 373 chronically infected patients with HIV-1 RNA >1000 copies/ml. Minimal and maximal drug-resistance profiles of potential transmitters were estimated by weighting resistance profiles of chronically infected patients with estimates of the Swiss HIV-1-infected population, the prevalence of exposure to antiviral drugs and the proportion of infections attributed to primary HIV infections. The drug-resistance prevalence in recently infected patients was 10.5% (one class drug resistance: 9.1%; two classes: 1.4%; three classes: 0%). Phylogenetic analysis revealed significant clustering for 30% of recent infections. The drug-resistance prevalence in chronically infected patients was 72.4% (one class: 29%; two classes: 27.6%; three classes: 15.8%). After adjustment, the risk of transmission relative to wild-type was reduced both for one class drug resistance (minimal and maximal estimates: odds ratio: 0.39, P<0.001; and odds ratio: 0.55, P=0.011, respectively), and for two to three class drug resistance (odds ratios: 0.05 and 0.07, respectively, P<0.001). Neither sexual behaviour nor HIV-1 RNA levels explained the low transmission of drug-resistant variants. These data suggest that drug-resistant variants and in particular multidrug-resistant variants have a substantially reduced transmission capacity.

Evolutionary trajectories of primate genes involved in HIV pathogenesis.

The current availability of five complete genomes of different primate species allows the analysis of genetic divergence over the last 40 million years of evolution. We hypothesized that the interspecies differences observed in susceptibility to HIV-1 would be influenced by the long-range selective pressures on host genes associated with HIV-1 pathogenesis. We established a list of human genes (n = 140) proposed to be involved in HIV-1 biology and pathogenesis and a control set of 100 random genes. We retrieved the orthologous genes from the genome of humans and of four nonhuman primates (Pan troglodytes, Pongo pygmaeus abeli, Macaca mulatta, and Callithrix jacchus) and analyzed the nucleotide substitution patterns of this data set using codon-based maximum likelihood procedures. In addition, we evaluated whether the candidate genes have been targets of recent positive selection in humans by analyzing HapMap Phase 2 single-nucleotide polymorphisms genotyped in a region centered on each candidate gene. A total of 1,064 sequences were used for the analyses. Similar median K(A)/K(S) values were estimated for the set of genes involved in HIV-1 pathogenesis and for control genes, 0.19 and 0.15, respectively. However, genes of the innate immunity had median values of 0.37 (P value = 0.0001, compared with control genes), and genes of intrinsic cellular defense had K(A)/K(S) values around or greater than 1.0 (P value = 0.0002). Detailed assessment allowed the identification of residues under positive selection in 13 proteins: AKT1, APOBEC3G, APOBEC3H, CD4, DEFB1, GML, IL4, IL8RA, L-SIGN/CLEC4M, PTPRC/CD45, Tetherin/BST2, TLR7, and TRIM5alpha. A number of those residues are relevant for HIV-1 biology. The set of 140 genes involved in HIV-1 pathogenesis did not show a significant enrichment in signals of recent positive selection in humans (intraspecies selection). However, we identified within or near these genes 24 polymorphisms showing strong signatures of recent positive selection. Interestingly, the DEFB1 gene presented signatures of both interspecies positive selection in primates and intraspecies recent positive selection in humans. The systematic assessment of long-acting selective pressures on primate genomes is a useful tool to extend our understanding of genetic variation influencing contemporary susceptibility to HIV-1.