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.

Rapid antagonistic coevolution between strains of the social amoeba Dictyostelium discoideum.

Social groups face a fundamental problem of overcoming selfish individuals capable of destroying cooperation. In the social amoeba Dictyostelium discoideum, there is evidence that some clones ('cheaters') contribute disproportionately to the viable spores in a fruiting body while avoiding the dead stalk cell fate. It remains unclear, however, whether this cheating is actually the product of selection. Here, I report the results of an experimental evolution study designed to test whether clones of D. discoideum will evolve resistance to cheating in the laboratory with genetic variation created only through spontaneous mutation. Two strains, one green fluorescent protein (GFP)-labelled and one wild-type, were allowed to grow and develop together before the wild-type strain was removed and replaced with a naïve strain evolving in parallel. Over the course of 10 social generations, the GFP-labelled strain reliably increased its representation in the spores relative to control populations that had never experienced the competitor. This competitive advantage extended to the non-social, vegetative growth portion of the life cycle, but not to pairwise competition with two other strains. These results indicate strong antagonism between strains, mediated by ample mutational variation for cheating and also suggest that arms races between strains in the wild may be common.

Melanin-based coloration is a nondirectionally selected sex-specific signal of offspring development in the Alpine swift

Two mutually exclusive hypotheses have been put forward to explain the evolution and adaptive function of melanin-based color traits. According to sexual selection theory melanism is a directionally selected signal of individual quality, whereas theory on the maintenance of genetic polymorphism proposes that alternative melanin-based variants achieve equal fitness. Alpine swift (Apus melba) males and females have a conspicuous patch of white feathers on the breast with their rachis varying continuously from white to black, and hence the breast varies from white to striated. If this trait is a sexually selected signal of quality, its expression should be condition dependent and the degree of melanism directionally selected. If variation in melanism is a polymorphism, its expression should be genetically determined and fitness of melanin-based variants equal. We experimentally tested these predictions by exchanging eggs or hatchlings between randomly chosen nests and by estimating survival and reproduction in relation to melanism. We found that breast melanism is heritable and that the environment and body condition do not significantly influence its expression. Between 5 and 50 days of age nestlings were heavier and their wings longer when breast feathers of their biological father were blacker, and they also fledged at a younger age. This shows that aspects of offspring quality covary positively with the degree of melanism. However, this did not result in directional selection because nestling survival and recruitment in the local breeding population were not associated with father breast melanism. Furthermore, adult survival, age at first reproduction and probability of skipping reproduction did not covary with the degree of melanism. Genetic variation in breast melanism is therefore maintained either because nonmelanic males achieve fitness similar to melanic males via a different route than producing fast-growing offspring, or because the advantage of producing fast-growing offspring is not sufficiently pronounced to result in directional selection.

Common variants in the ATP2B1 gene are associated with susceptibility to hypertension: the Japanese Millennium Genome Project.

Hypertension is one of the most common complex genetic disorders. We have described previously 38 single nucleotide polymorphisms (SNPs) with suggestive association with hypertension in Japanese individuals. In this study we extend our previous findings by analyzing a large sample of Japanese individuals (n=14 105) for the most associated SNPs. We also conducted replication analyses in Japanese of susceptibility loci for hypertension identified recently from genome-wide association studies of European ancestries. Association analysis revealed significant association of the ATP2B1 rs2070759 polymorphism with hypertension (P=5.3×10(-5); allelic odds ratio: 1.17 [95% CI: 1.09 to 1.26]). Additional SNPs in ATP2B1 were subsequently genotyped, and the most significant association was with rs11105378 (odds ratio: 1.31 [95% CI: 1.21 to 1.42]; P=4.1×10(-11)). Association of rs11105378 with hypertension was cross-validated by replication analysis with the Global Blood Pressure Genetics consortium data set (odds ratio: 1.13 [95% CI: 1.05 to 1.21]; P=5.9×10(-4)). Mean adjusted systolic blood pressure was highly significantly associated with the same SNP in a meta-analysis with individuals of European descent (P=1.4×10(-18)). ATP2B1 mRNA expression levels in umbilical artery smooth muscle cells were found to be significantly different among rs11105378 genotypes. Seven SNPs discovered in published genome-wide association studies were also genotyped in the Japanese population. In the combined analysis with replicated 3 genes, FGF5 rs1458038, CYP17A1, rs1004467, and CSK rs1378942, odds ratio of the highest risk group was 2.27 (95% CI: 1.65 to 3.12; P=4.6×10(-7)) compared with the lower risk group. In summary, this study confirmed common genetic variation in ATP2B1, as well as FGF5, CYP17A1, and CSK, to be associated with blood pressure levels and risk of hypertension.