A step-by-step tutorial to use HierFstat to analyse populations hierarchically structured at multiple levels.

The populations of parasites and infectious agents are most of the time structured in complex hierarchy that lies beyond the classical nested design described by Wright's F-statistics (F(IS), F(ST) and F(IT)). In this note we propose a user-friendly step-by-step notice for using recent software (HierFstat) that computes and test fixation indices for any hierarchical structure. We add some tricks and tips for some special data kind (haploid, single locus), some other procedure (bootstrap over loci) and how to handle crossed factors.

Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.

Prominent use of distal 5’ transcription start sites and discovery of a large number of additional exons in ENCODE regions

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

A multicenter study to validate the reproducibility of MSI testing with a panel of 5 quasimonomorphic mononucleotide repeats.

Microsatellite instability (MSI) testing in clinics is becoming increasingly widespread; therefore, there is an urgent need for methodology standardization and the availability of quality control. This study is aimed to assess the interlaboratory reproducibility of MSI testing in archive samples by using a panel of 5 recently introduced, mononucleotide repeats (MNR). The quality control involved 8 European institutions. Participants were supplied with DNA extracted from 15 archive colon carcinoma samples and from the corresponding normal tissues. Every group was asked to assess the MSI status of the samples by using the BAT25, BAT26, NR21, NR24, and NR27 mononucleotide markers. Four institutions repeated the analysis using the NCI reference panel to confirm the results obtained with the MNR markers. The overall concordance among institutions for MSI analyses at single locus level was 97.7% when using the MNR panel and 95.0% with the NCI one. The laboratories obtained a full agreement in scoring the MSI status of each patient sample, both using the mononucleotide and the NCI marker sets. With the NCI marker set, however, concordance was lowered to 85.7% when considering the MSI-Low phenotype. Concordance between the 2 panels in scoring the MSI status of each sample was complete if no discrimination was made between MSI-Stable and MSI-L, whereas it dropped to 76.7% if MSI-L was considered. In conclusion, the use of the MNR panel seems to be a robust approach that yields a very high level of reproducibility. The results obtained with the 5 MNR are diagnostically consistent with those obtained by the use of the NCI markers, except for the MSI-Low phenotype.

Genetic causes of transitions from sexual reproduction to asexuality in plants and animals.

The persistence of sexual reproduction in the face of competition from asexual invaders is more likely if asexual lineages are produced infrequently or have low fitness. The generation rate and success of new asexual lineages will be influenced by the proximate mechanisms underlying transitions to asexuality. As such, characterization of these mechanisms can help explain the distribution of reproductive modes among natural populations. Here, we synthesize the literature addressing proximate causes of transitions from sexual to asexual reproduction in plants and animals. In cyclical and facultatively asexual taxa, individual mutations can cause obligate asexuality. The evolution of asexuality in obligately sexual groups is more complex, requiring the simultaneous acquisition of two traits generally controlled by different genetic factors: unreduced gamete formation and spontaneous development of unfertilized gametes. At least three 'pre-adaptations' could favour transitions to obligate asexuality in obligate sexuals. First, linkage among loci affecting separate key components of asexuality facilitates its spread, with evidence for these linkage blocks in plants. Second, asexuality should evolve more readily in haplodiploids; support for this hypothesis comes from two examples where a single locus causes transitions to asexuality. Third, standing genetic variation for the production of unreduced gametes could facilitate transitions to asexuality, but whether the ability to produce unreduced gametes contributes to the evolution of obligate asexuality remains unclear. We close by reviewing the associations between asexuality, hybridization and polyploidy, and argue that current data suggest that hybridization is more likely to play a causal role in transitions to asexuality than polyploidy.

Duplication and concerted evolution in a master sex determiner under balancing selection.

The transformer (tra) gene is a key regulator in the signalling hierarchy controlling all aspects of somatic sexual differentiation in Drosophila and other insects. Here, we show that six of the seven sequenced ants have two copies of tra. Surprisingly, the two paralogues are always more similar within species than among species. Comparative sequence analyses indicate that this pattern is owing to the ongoing concerted evolution after an ancestral duplication rather than independent duplications in each of the six species. In particular, there was strong support for inter-locus recombination between the paralogues of the ant Atta cephalotes. In the five species where the location of paralogues is known, they are adjacent to each other in four cases and separated by only few genes in the fifth case. Because there have been extensive genomic rearrangements in these lineages, this suggests selection acting to conserve their synteny. In three species, we also find a signature of positive selection in one of the paralogues. In three bee species where information is available, the tra gene is also duplicated, the copies are adjacent and in at least one species there was recombination between paralogues. These results suggest that concerted evolution plays an adaptive role in the evolution of this gene family.

The dermatophyte species Arthroderma benhamiae: intraspecies variability and mating behaviour.

Arthroderma benhamiae is a zoophilic dermatophyte belonging to the Trichophyton mentagrophytes species complex. Here, a population of A. benhamiae wild strains from the same geographical area (Switzerland) was studied by comparing their morphology, assessing their molecular variability using internal transcribed spacer (ITS) and 28S rRNA gene sequencing, and evaluating their interfertility. Sequencing of the ITS region and of part of the 28S rRNA gene revealed the existence of two infraspecific groups with markedly different colony phenotypes: white (group I) and yellow (group II), respectively. For all strains, the results of mating type identification by PCR, using HMG (high-mobility group) and α-box genes in the mating type locus as targets, were in total accordance with the results of mating type identification by strain confrontation experiments. White-phenotype strains were of mating type + (mt+) or mating type - (mt-), whilst yellow-phenotype strains were all mt-. White and yellow strains were found to produce fertile cleistothecia after mating with A. benhamiae reference tester strains, which belonged to a third group intermediate between groups I and II. However, no interfertility was observed between yellow strains and white strains of mt+. A significant result was that white strains of mt- were able to mate and produce fertile cleistothecia with the white A. benhamiae strain CBS 112371 (mt+), the genome of which has recently been sequenced and annotated. This finding should offer new tools for investigating the biology and genetics of dermatophytes using wild-type strains.

Effects of epigenetic regulatory elements on telomeric gene expression

Transgene expression in eukaryotic cells strongly depends on the locus of integration in the host genome and often results in limited transcription level because of unfavorable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which are believed to act on the chromatin structure and may protect transgenes from this so-called position effect. Despite being extensively used to increase transgene expression, the mechanism of action of many of these elements remains largely unknown. Here we evaluated different epigenetic regulatory DNA elements for their ability to protect transgene transcription at telomeres, a defined chromatin environment associated to low or inconsistent expression caused by the Telomere Position Effect (TPE). For the assessment of the effects of epigenetic regulators at telomeres, a novel dual reporter system had to be designed. Telomeric integration of the newly-developed dual reporter system carrying different epigenetic regulators showed that MARs (Matrix Attachment Regions), a UCOE (Ubiquitous Chromatin-Opening Element) or the chicken cHS4 insulator have strong barrier activity which prevented TPE from spreading toward the centromere, resulting in stable and in some cases increased expression of a telomeric-distal reporter gene. In addition, MARs and STAR element 40 resulted in an increase of cells expressing the telomeric-proximal reporter gene, suggesting also an anti-silencing effect. Chromatin immunoprecipitation assays revealed that at telomeres MARs actively promote the deposition of euchromatic histone marks, especially acetylation of both histone H3 and H4, which might be involved in MARs' barrier and transcriptional activator activities. Differently, the chromatin in proximity of the UCOE element was depleted of several repressive chromatin marks, such as trimethylated lysine 9 and lysine 27 on histone H3 and trimethylated lysine 20 of histone H4, possibly favoring the preservation of an open chromatin structure at the integration site. We conclude that epigenetic regulatory elements that may be used to enhance and sustain transgene expression have all a specific epigenetic signature which might be at the basis of their mechanism of action, and that a combination of different classes of epigenetic regulators might be advantageous when high levels of protein expression are required. - L'expression des transgènes dans les cellules eucaryotes est fortement influencée par leur site d'intégration dans le génome. En effet, une structure chromatinienne défavorable au niveau du site d'intégration peut fortement limiter le degré d'expression d'un transgène. Il existe toutefois des séquences d'ADN qui, en agissant sur la structure de la chromatine, permettent de limiter cet effet de position et, par conséquent, de promouvoir l'expression soutenue d'un transgène. Ces éléments génomiques, connus comme régulateurs épigénétiques, sont largement utilisés dans plusieurs domaines où une expression élevée et soutenue est requise, malgré un mode de fonctionnement parfois méconnu. Dans cette étude, j'ai évalué la capacité de différents régulateurs épigénétiques à protéger la transcription de transgènes au niveau des télomères, régions chromatiniennes bien définies qui ont été associées à un fort effet de silençage, connu comme «effet de position télomérique». Pour cela, un nouveau système à deux gènes rapporteurs a été développé. Lorsque des MARs (Matrix Attachment Regions, séquences d'ADN pouvant s'associer à la matrice nucléaire), un UCOE (Ubiquitous Chromatin-Opening Element, élément pouvant ouvrir la chromatine) ou l'isolateur génétique cHS4 (dérivé du locus de la β-globine de poulet) sont placés entre les deux gènes rapporteurs, une forte activité barrière bloquant la propagation de la chromatine répressive télomérique est observée, résultant en un plus grand nombre de cellules exprimant le gène télomérique-distal. D'autre part, une augmentation du nombre de cellules exprimant le gène télomérique-proximal, observée en présence des éléments MAR et STAR 40 (Stabilizing Anti-Repressor element 40, un élément pouvant prévenir la répression génique), suggère aussi un faible effet anti-silençage pour ces éléments. Des expériences d'immunoprécipitation de la chromatine démontrent qu'au télomère, les MARs favorisent l'assemblage de marqueurs de la chromatine active, surtout l'acétylation des histones H3 et H4, qui pourraient être à la base de l'activité barrière et de celle d'activateur transcriptionel. Différemment, la chromatine à proximité de l'élément UCOE est particulièrement pauvre en marqueurs de la chromatine silencieuse, comme la trimethylation des lysines 9 et 27 de l'histone H3, ainsi que la trimethylation de la lysine 20 de l'histone H4. Cela suggère que UCOE pourrait préserver une structure chromatinienne ouverte au site d'intégration, favorisant l'expression des gènes à sa proximité. En conclusion, les régulateurs épigénétiques analysés lors de cette étude ont tous montré une signature épigénétique spécifique qui pourrait être à la base de leurs mécanismes de fonctionnement, suggérant aussi qu'une utilisation d'éléments épigénétiques de classe différente dans un même vecteur d'expression pourrait être avantageuse lorsque de hauts et soutenus niveaux d'expression sont nécessaires.

Heterozygote advantage and the maintenance of polymorphism for multilocus traits.

Explaining how polymorphism is maintained in the face of selection remains a puzzle since selection tends to erode genetic variation. Provided an infinitely large unsubdivided population and no frequency-dependance of selective values, heterozygote advantage is the text book explanation for the maintenance of polymorphism when selection acts at a diallelic locus. Here, we investigate whether this remains true when selection acts at multiple diallelic loci. We use five different definitions of heterozygote advantage that largely cover this concept for multiple loci. Using extensive numerical simulations, we found no clear associations between the presence of any of the five definitions of heterozygote advantage and the maintenance of polymorphism at all loci. The strength of the association decreases as the number of loci increases or as recombination decreases. We conclude that heterozygote advantage cannot be a general mechanism for the maintenance of genetic polymorphism at multiple loci. These findings suggest that a correlation between the number of heterozygote loci and fitness is not warranted on theoretical ground.

Evaluation of adding a second marker to overcome Staphylococcus aureus spa typing homoplasies.

The utility of sequencing a second highly variable locus in addition to the spa gene (e.g., double-locus sequence typing [DLST]) was investigated to overcome limitations of a Staphylococcus aureus single-locus typing method. Although adding a second locus seemed to increase discriminatory power, it was not sufficient to definitively infer evolutionary relationships within a single multilocus sequence type (ST-5).

Diversity of staphylococcal cassette chromosome mec elements in predominant methicillin-resistant Staphylococcus aureus clones in a small geographic area.

Recent population genetic studies suggest that staphylococcal cassette chromosome mec (SCCmec) was acquired much more frequently than previously thought. In the present study, we aimed to investigate the diversity of SCCmec elements in a local methicillin-resistant Staphylococcus aureus (MRSA) population. Each MRSA isolate (one per patient) recovered in the Vaud canton of Switzerland from January 2005 to December 2008 was analyzed by the double-locus sequence typing (DLST) method and SCCmec typing. DLST analysis indicated that 1,884/2,036 isolates (92.5%) belong to four predominant clones. As expected from the local spread of a clone, most isolates within clones harbored an identical SCCmec type. However, three to seven SCCmec types have been recovered in every predominant DLST clone, suggesting that some of these elements might have been acquired locally. This pattern could also be explained by distinct importations of related isolates into the study region. The addition of a third highly variable locus to further increase the discriminatory power of typing as well as epidemiological data suggested that most ambiguous situations were explained by the second hypothesis. In conclusion, our study showed that even if the acquisition of new SCCmec elements at a local level likely occurs, it does not explain all the diversity observed in the study region.

Recessive multiple epiphyseal dysplasia (rMED) with homozygosity for C653S mutation in the DTDST gene - phenotype, molecular diagnosis and surgical treatment of habitual dislocation of multilayered patella: case report.

BACKGROUND: Multiple epiphyseal dysplasia (MED) is one of the more common generalised skeletal dysplasias. Due to its clinical heterogeneity diagnosis may be difficult. Mutations of at least six separate genes can cause MED. Joint deformities, joint pain and gait disorders are common symptoms. CASE PRESENTATION: We report on a 27-year-old male patient suffering from clinical symptoms of autosomal recessive MED with habitual dislocation of a multilayered patella on both sides, on the surgical treatment and on short-term clinical outcome. Clinical findings were: bilateral hip and knee pain, instability of femorotibial and patellofemoral joints with habitual patella dislocation on both sides, contractures of hip, elbow and second metacarpophalangeal joints. Main radiographic findings were: bilateral dislocated multilayered patella, dysplastic medial tibial plateaus, deformity of both femoral heads and osteoarthritis of the hip joints, and deformity of both radial heads. In the molecular genetic analysis, the DTDST mutation g.1984T > A (p.C653S) was found at the homozygote state. Carrier status was confirmed in the DNA of the patient's parents. The mutation could be considered to be the reason for the patient's disease. Surgical treatment of habitual patella dislocation with medialisation of the tibial tuberosity led to an excellent clinical outcome. CONCLUSIONS: The knowledge of different phenotypes of skeletal dysplasias helps to select genes for genetic analysis. Compared to other DTDST mutations, this is a rather mild phenotype. Molecular diagnosis is important for genetic counselling and for an accurate prognosis. Even in case of a multilayered patella in MED, habitual patella dislocation could be managed successfully by medialisation of the tibial tuberosity.

A limited role for beta-selection during gamma delta T cell development.

T cells belong to two distinct lineages expressing either alpha beta or gamma delta TCR. During alpha beta T cell development, it is clearly established that productive rearrangement at the TCR beta locus in immature precursor cells leads to the expression of a pre-TCR complex. Signaling through the pre-TCR results in the selective proliferation and maturation of TCR beta+ cells, a process that is known as beta-selection. However, the potential role of beta-selection during gamma delta T cell development is controversial. Whereas PCR-RFLP and sequencing techniques have provided evidence for a bias toward in-frame VDJ beta rearrangements in gamma delta cells (consistent with beta-selection), gamma delta cells apparently develop normally in mice that are unable to assemble a pre-TCR complex due to a deficiency in TCR beta or pT alpha genes. In this report, we have directly addressed the physiologic significance of beta-selection during gamma delta cell development in normal mice by quantitating intracellular TCR beta protein in gamma delta cells and correlating its presence with cell cycle status. Our results indicate that beta-selection plays a significant (although limited) role in gamma delta cell development by selectively amplifying a minor subset of gamma delta precursor cells with productively rearranged TCR beta genes.

A large deletion in the adRP gene PRPF31: evidence that haploinsufficiency is the cause of disease.

PURPOSE: To report a large deletion that encompasses more than 90% of PRPF31 gene and two other neighboring genes in their entirety in an adRP pedigree that appears to show only the typical clinical features of retinitis pigmentosa. METHODS: To identify PRPF31 mutation in a dominant RP family (ADRP2) previously linked to the RP11 locus, the 14 exons of PRPF31 were screened for mutations by direct sequencing. To investigate the possibility of a large deletion, microsatellite markers near PRPF31 gene were analyzed by non-denaturing PAGE. RESULTS: Initial screening of PRPF31 gene in the ADRP2 family did not reveal an obvious mutation. A large deletion was however suspected due to lack of heterozygosity for nearly all PRPF31 intragenic single nucleotide polymorphysm (SNPs). In order to estimate the size of the deletion, SNPs and microsatellite markers spanning and flanking PRPF31 were analyzed in the entire ADRP2 family. Haplotype analysis with the above markers suggested a deletion of approximately 30 kb that included the putative promoter region of a novel gene OSCAR, the entire genomic content of genes NDUFA3, TFPT and more than 90% of PRPF31 gene. Sequence analysis of the region flanking the potential deletion showed a high presence of Alu elements implicating Alu mediated recombination as the mechanism responsible for this event. CONCLUSIONS: This mutation provides evidence that haploinsufficiency rather than aberrant function of mutated proteins is the cause of disease in these adRP patients with mutations in PRPF31 gene.

Androgenetic alopecia: identification of four genetic risk loci and evidence for the contribution of WNT signaling to its etiology.

The pathogenesis of androgenetic alopecia (AGA, male-pattern baldness) is driven by androgens, and genetic predisposition is the major prerequisite. Candidate gene and genome-wide association studies have reported that single-nucleotide polymorphisms (SNPs) at eight different genomic loci are associated with AGA development. However, a significant fraction of the overall heritable risk still awaits identification. Furthermore, the understanding of the pathophysiology of AGA is incomplete, and each newly associated locus may provide novel insights into contributing biological pathways. The aim of this study was to identify unknown AGA risk loci by replicating SNPs at the 12 genomic loci that showed suggestive association (5 × 10(-8)<P<10(-5)) with AGA in a recent meta-analysis. We analyzed a replication set comprising 2,759 cases and 2,661 controls of European descent to confirm the association with AGA at these loci. Combined analysis of the replication and the meta-analysis data identified four genome-wide significant risk loci for AGA on chromosomes 2q35, 3q25.1, 5q33.3, and 12p12.1. The strongest association signal was obtained for rs7349332 (P=3.55 × 10(-15)) on chr2q35, which is located intronically in WNT10A. Expression studies in human hair follicle tissue suggest that WNT10A has a functional role in AGA etiology. Thus, our study provides genetic evidence supporting an involvement of WNT signaling in AGA development.