Demonstration of loss of heterozygosity by single-nucleotide polymorphism microarray analysis and alterations in strain morphology in Candida albicans strains during infection

Candida albicans is a diploid yeast with a predominantly clonal mode of reproduction, and no complete sexual cycle is known. As a commensal organism, it inhabits a variety of niches in humans. It becomes an opportunistic pathogen in immunocompromised patients and can cause both superficial and disseminated infections. It has been demonstrated that genome rearrangement and genetic variation in isolates of C. albicans are quite common. One possible mechanism for generating genome-level variation among individuals of this primarily clonal fungus is mutation and mitotic recombination leading to loss of heterozygosity (LOH). Taking advantage of a recently published genome-wide single-nucleotide polymorphism (SNP) map (A. Forche, P. T. Magee, B. B. Magee, and G. May, Eukaryot. Cell 3:705-714, 2004), an SNP microarray was developed for 23 SNP loci residing on chromosomes 5, 6, and 7. It was used to examine 21 strains previously shown to have undergone mitotic recombination at the GAL1 locus on chromosome 1 during infection in mice. In addition, karyotypes and morphological properties of these strains were evaluated. Our results show that during in vivo passaging, LOH events occur at observable frequencies, that such mitotic recombination events occur independently in different loci across the genome, and that changes in karyotypes and alterations of phenotypic characteristics can be observed alone, in combination, or together with LOH.     

Development of mapped simple sequence repeat markers from common bean (Phaseolus vulgaris L.) based on genome sequences of a Chinese landrace and diversity evaluation

Microsatellite or single sequence repeat (SSR) markers have been commonly used in genetic research in many crop species, including common bean (Phaseolus vulgaris L.). A limited number of existing SSR markers have been designed from high-throughput sequencing of the genome, warranting the exploitation of new SSR markers from genomic regions. In this paper, we sequenced total DNA from the genotype Hong Yundou with a 454-FLX pyrosequencer and found numerous SSR loci. Based on these, a large number of SSR markers were developed and 90 genomic-SSR markers with clear bands were tested for mapping and diversity detection. The new SSR markers proved to be highly polymorphic for molecular polymorphism, with an average polymorphism information content value of 0.44 in 131 Chinese genotypes and breeding lines, effective for distinguishing Andean and Mesoamerican genotypes. In addition, we integrated 85 primers of the 90 polymorphism markers into the bean map using an F2 segregating population derived from Hong Yundou crossed with Jingdou. The distribution of SSR markers among 11 chromosomes was not random and tended to cluster on the linkage map, with 14 new markers mapped on chromosome Pv01, whereas only four loci were located on chromosome Pv04. Overall, these new markers have potential for genetic mapping, genetic diversity studies and map-based cloning in common bean.     

Generalized T2 test for genome association studies

Recent progress in the development of single-nucleotide polymorphism (SNP) maps within genes and across the genome provides a valuable tool for fine-mapping and has led to the suggestion of genomewide association studies to search for susceptibility loci for complex traits. Test statistics for genome association studies that consider a single marker at a time, ignoring the linkage disequilibrium between markers, are inefficient. In this study, we present a generalized T2 statistic for association studies of complex traits, which can utilize multiple SNP markers simultaneously and considers the effects of multiple disease-susceptibility loci. This generalized T2 statistic is a corollary to that originally developed for multivariate analysis and has a close relationship to discriminant analysis and common measure of genetic distance. We evaluate the power of the generalized T2 statistic and show that power to be greater than or equal to those of the traditional chi2 test of association and a similar haplotype-test statistic. Finally, examples are given to evaluate the performance of the proposed T2 statistic for association studies using simulated and real data.     

Comparative mapping of a region on chromosome 10 containing QTL for reproduction in swine

Several quantitative trait loci (QTL) for important reproductive traits (age of puberty, ovulation rate, nipple number and plasma FSH) have been identified on the long arm of porcine chromosome 10. Bi-directional chromosome painting has shown that this region is homologous to human chromosome 10p. Because few microsatellite or type I markers have been placed on SSC10, we wanted to increase the density of known ESTs mapped in this region of the porcine genome. Genes were chosen for their position on human chromosome 10, sequence availability from the TIGR pig gene indices, and their potential as a candidate gene. The PCR primers were designed to amplify across introns or 3'-UTR to maximize single nucleotide polymorphism (SNP) discovery. Parents of the mapping population (one sire and seven dams) were amplified and sequenced to find informative markers. The SNPs were genotyped using primer extension and mass spectrometry. These amplification products were also used to probe a BAC library (RPCI-44, Roswell Park Cancer Institute) for positive clones and screened for microsatellites. Six genes from human chromosome 10p (AKR1C2, PRKCQ, ITIH2, ATP5C1, PIP5K2A and GAD2) were mapped in the MARC swine mapping population. Gene order was conserved within these markers from centromere to telomere of porcine chromosome 10q, as compared with human chromosome 10p. Four of these genes (PIP5K2A, ITIH2, GAD2 and AKR1C2), which map under QTL, are potential candidate genes. Identification of porcine homologues near important QTL and development of a comparative map for this chromosome will allow further fine- mapping and positional cloning of candidate genes affecting reproductive traits.     

Characterization of a divergent chromosome region in the willow warbler Phylloscopus trochilus using avian genomic resources

Genome scans have made it possible to find outlier markers thought to have been influenced by divergent selection in almost any wild population. However, the lack of genomic information in nonmodel species often makes it difficult to associate these markers with certain genes or chromosome regions. Furthermore, the extent of linkage disequilibrium (LD) in the genome will determine the density of markers required to identify the genes under selection. In this study, we investigated a chromosome region in the willow warbler Phylloscopus trochilus surrounding a single marker previously identified in a genome scan. We first located the marker in the assembled genome of another species, the zebra finch Taeniopygia guttata, and amplified surrounding sequences in Fennoscandian willow warblers. Within an investigated chromosome region of 7.3 Mb as mapped to the zebra finch genome, we observed elevated genetic differentiation between a southern and a northern population across a 2.5-Mb interval comprising numerous coding genes. Within the southern and northern populations, higher values of LD were mostly found between SNPs within the same locus, but extended across distantly situated loci when the analyses were restricted to sampling sites showing intermediate allele frequencies of southern and northern alleles. Our study shows that cross-species genome information is a useful resource to obtain candidate sequences adjacent to outlier markers in nonmodel species.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

SELECTION ON OVERDOMINANT GENES MAINTAINS HETEROZYGOSITY ALONG MULTIPLE CHROMOSOMES IN A CLONAL LINEAGE OF HONEY BEE

Correlations between fitness and genome‐wide heterozygosity (heterozygosity‐fitness correlations, HFCs) have been reported across a wide range of taxa. The genetic basis of these correlations is controversial: do they arise from genome‐wide inbreeding (“general effects”) or the “local effects” of overdominant loci acting in linkage disequilibrium with neutral loci? In an asexual thelytokous lineage of the Cape honey bee (Apis mellifera capensis), the effects of inbreeding have been homogenized across the population, making this an ideal system in which to detect overdominant loci, and to make inferences about the importance of overdominance on HFCs in general. Here we investigate the pattern of zygosity along two chromosomes in 42 workers from the clonal Cape honey bee population. On chromosome III (which contains the sex‐locus, a gene that is homozygous‐lethal) and chromosome IV we show that the pattern of zygosity is characterized by loss of heterozygosity in short regions followed by the telomeric restoration of heterozygosity. We infer that at least four selectively overdominant genes maintain heterozygosity on chromosome III and three on chromosome IV via local effects acting on neutral markers in linkage disequilibrium. We conclude that heterozygote advantage and local effects may be more common and evolutionarily significant than is generally appreciated.     

Variants of the anti-Müllerian hormone gene in a compound heterozygote with the persistent Müllerian duct syndrome and his family

Summary The human genome contains a large number of interspersed simple repeat sequences that are variable in length and can therefore serve as highly informative, polymorphic markers. Typing procedures include conventional multilocus and single locus probing, and polymerase chain reaction aided analysis. We have identified simple sequences in a cosmid clone stemming from the human Y chromosome and consisting of (gata)_n repeats. We have compared these with two equivalent simple repeat loci from chromosome 12. After amplifying the tandemly repeated motifs, we detected between four and eight different alleles at each of the three loci. Codominant inheritance of the alleles was established in family studies and the informativity of the simple repeat loci was determined by typing unrelated individuals. The polymorphisms are suitable for application in linkage studies, practical forensic case work, deficiency cases in paternity determination, and for studying ethnological questions. The mutational mechanisms that bring about changes in simple repeats located both on the autosomes and on the sex chromosomes, are discussed.Professor Dr. Otto Prokop (Humboldt-Universität Berlin) on the occasion of his 70th birthday     

Comparative mapping in grasses. Wheat relationships

Conventionally, the genetics of species of the family Gramineae have been studied separately. Comparative mapping using DNA markers offers a method of combining the research efforts in each species. In this study, we developed consensus maps for members of the Triticeae tribe (Triticum aestivum, T. tauschii, andHordeum spp.) and compared them to rice, maize and oat. The aneuploid stocks available in wheat are invaluable for comparative mapping because almost every DNA fragment can be allocated to a chromosome arm, thus preventing erroneous conclusions about probes that could not be mapped due to a lack of polymorphism between mapping parents. The orders of the markers detected by probes mapped in rice, maize and oat were conserved for 93, 92 and 94% of the length of Triticeae consensus maps, respectively. The chromosome segments duplicated within the maize genome by ancient polyploidization events were identified by homoeology of segments from two maize chromosomes to regions of one Triticeae chromosome. Homoeologous segments conserved across Triticeae species, rice, maize, and oat can be identified for each Triticeae chromosome. Putative orthologous loci for several simply inherited and quantitatively inherited traits in Gramineae species were identified.Communicated by H. Saedler     

Characterization of AT-rich microsatellites in common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.)

Polymorphism of microsatellite markers is often associated with the simple sequence repeat motif targeted. AT-rich microsatellites tend to be highly variable and this appears to be notable, especially in legume genomes. To analyze the value of AT-rich microsatellites for common bean (Phaseolus vulgaris L.), we developed a total of 85 new microsatellite markers, 74 of which targeted ATA or other AT-rich motif loci and 11 of which were made for GA, CA or CAC motif loci. We evaluated the loci for the level of allelic diversity in comparison to previously characterized microsatellites using a panel of 18 standard genotypes and genetically mapped any loci polymorphic in the DOR364 × G19833 population. The majority of the microsatellites produced single bands and detected single loci, however, 15 of the AT-rich microsatellites produced multiple or double banding patterns; while only one of the GA or CA-rich microsatellites did. The polymorphism information content (PIC) values averaged 0.892 and 0.600 for the AT and ATA motif microsatellites, respectively, but only 0.140 for the CA-rich microsatellites. GA microsatellites, which had a large average number of repeats, had high to intermediate PIC, averaging 0.706. A total of 45 loci could be genetically mapped and distribution of the loci across the genome was skewed towards non-distal locations with a greater prevalence of loci on linkage groups b02, b09 and b11. AT-rich microsatellites were found to be a useful source of polymorphic markers for mapping and diversity assessment in common bean that appears to uncover higher diversity than other types of simple sequence repeat markers.     

Standard linkage and association methods identify the mechanism of four susceptibility genes for a simulated complex disease

The simulated dataset of the Genetic Analysis Workshop 14 provided affection status and the presence or absence of 12 traits. It was determined that all affected individuals must have traits E, F and H (EFH phenotype) and they must also have either trait B (B subtype) or traits C, D, and G (CDG subtype). A genome screen was performed, and linkage peaks were identified on chromosomes 1, 3, 5, and 9 using microsatellite markers. Dense panels of single-nucleotide polymorphism (SNP) markers were ordered for each of the four linkage peaks. In each case, association analyses identified a single SNP that accounted for the linkage evidence. The SNP on chromosome 1 appeared to primarily influence the B subtype, while the SNPs on chromosomes 5 and 9 primarily influenced the CDG subtype. The chromosome 3 SNP had the strongest effect and influenced both subtypes, as well as the requisite EFH phenotype. Recognizing the two subtypes prior to linkage analysis was key to identifying these loci using only a single replicate. This highlights the need in real life situations for careful examination of the phenotypic data prior to genetic analysis.     

Differential divergence in autosomes and sex chromosomes is associated with intra‐island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome‐wide markers to investigate the diversification of the Reunion grey white‐eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping‐by‐sequencing and pooled RAD‐seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z‐linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.     

Allele-specific amplification in cancer revealed by SNP array analysis

Amplification, deletion, and loss of heterozygosity of genomic DNA are hallmarks of cancer. In recent years a variety of studies have emerged measuring total chromosomal copy number at increasingly high resolution. Similarly, loss-of-heterozygosity events have been finely mapped using high-throughput genotyping technologies. We have developed a probe-level allele-specific quantitation procedure that extracts both copy number and allelotype information from single nucleotide polymorphism (SNP) array data to arrive at allele-specific copy number across the genome. Our approach applies an expectation-maximization algorithm to a model derived from a novel classification of SNP array probes. This method is the first to our knowledge that is able to (a) determine the generalized genotype of aberrant samples at each SNP site (e.g., CCCCT at an amplified site), and (b) infer the copy number of each parental chromosome across the genome. With this method, we are able to determine not just where amplifications and deletions occur, but also the haplotype of the region being amplified or deleted. The merit of our model and general approach is demonstrated by very precise genotyping of normal samples, and our allele-specific copy number inferences are validated using PCR experiments. Applying our method to a collection of lung cancer samples, we are able to conclude that amplification is essentially monoallelic, as would be expected under the mechanisms currently believed responsible for gene amplification. This suggests that a specific parental chromosome may be targeted for amplification, whether because of germ line or somatic variation. An R software package containing the methods described in this paper is freely available at http://genome.dfci.harvard.edu/~tlaframb/PLASQ.     

Completion of a parasexual cycle in Candida albicans by induced chromosome loss in tetraploid strains

The human pathogenic fungus Candida albicans has traditionally been classified as a diploid, asexual organism. However, mating-competent forms of the organism were recently described that produced tetraploid mating products. In principle, the C.albicans life cycle could be completed via a sexual process, via a parasexual mechanism, or by both mechanisms. Here we describe conditions in which growth of a tetraploid strain of C.albicans on Saccharomyces cerevisiae 'pre-sporulation' medium induced efficient, random chromosome loss in the tetraploid. The products of chromosome loss were often strains that were diploid, or very close to diploid, in DNA content. If they inherited the appropriate MTL (mating-type like) loci, these diploid products were themselves mating competent. Thus, an efficient parasexual cycle can be performed in C.albicans, one that leads to the reassortment of genetic material in this organism. We show that this parasexual cycle-consisting of mating followed by chromosome loss-can be used in the laboratory for simple genetic manipulations in C.albicans.     

High Potential of a Transposon mPing as a Marker System in japonica �� japonica Cross in Rice

Although quantitative traits loci (QTL) analysis has been widely performed to isolate agronomically important genes, it has been difficult to obtain molecular markers between individuals with similar phenotypes (assortative mating). Recently, the miniature inverted-repeat transposable element mPing was shown to be active in the japonica strain Gimbozu EG4 where it had accumulated more than 1000 copies. In contrast, most other japonicas, including Nipponbare, have 50 or fewer mPing insertions in their genome. In this study we have exploited the polymorphism of mPing insertion sites to generate 150 PCR markers in a cross between the closely related japonicas, Nipponbare × Gimbozu (EG4). These new markers were distributed in genic regions of the whole genome and showed significantly higher polymorphism (150 of 183) than all other molecular markers tested including short sequence repeat markers (46 of 661). In addition, we performed QTL analysis with these markers using recombinant inbred lines derived from Nipponbare × Gimbozu EG4, and successfully mapped a locus involved in heading date on the short arm of chromosome 6. Moreover, we could easily map two novel loci involved in the culm length on the short arms of chromosomes 3 and 10.Key words: Linkage mapping, Transposon, japonica, Oryza sativa L., QTL analysis     

Identification of 11 polymorphic simple sequence repeat loci in the phytopathogenic fungus Fusarium pseudograminearum as a tool for genetic studies

Simple sequence repeat (SSR) markers for Fusarium pseudograminearum with 2 to 3 bp repeat motifs were identified by screening the genome database of the related species Fusarium graminearum. Twelve SSRs amplified single loci in both F. graminearum and F. pseudograminearum. Forty F. pseudograminearum and six F. graminearum individual isolates were screened to determine levels of polymorphism, with all SSRs displaying three to 14 alleles across all isolates. Eleven SSRs were polymorphic across F. pseudograminearum isolates tested proving the usefulness of genome databases of closely related species in identifying genetic markers.     

EST-SSR DNA polymorphism in durum wheat (Triticum durum L.) collections

SSRs derived from EST were molecular markers belonging to the transcribed region of the genome. Therefore, any polymorphism detected using EST-SSRs might reflect the better relationship among species or varieties. Using wheat EST-SSR markers, 60 durum wheat (Triticum durum L.) accessions from seven countries were investigated. Twenty-five primer pairs could amplify successfully in the 60 durum wheat accessions, of which tri-nucleotide repeats were the dominant type, and revealed 26 loci on all seven wheat homologous chromosome groups. A total of 87 eSSR alleles were detected, and the number of alleles detected by a single pair of primers ranged from 1 to 11, with an average of 3.3 alleles per locus. Higher numbers of alleles and PIC were identified on the B genome than those on the A genome.     

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.     

Levels of linkage disequilibrium in a wild bird population

Population-based mapping approaches are attractive for tracing the genetic background to phenotypic traits in wild species, given that it is often difficult to gather extensive and well-defined pedigrees needed for quantitative trait locus analysis. However, the feasibility of association or hitch-hiking mapping is dependent on the degree of linkage disequilibrium (LD) in the population, on which there is yet limited information for wild species. Here we use single nucleotide polymorphism (SNP) markers from 23 genes in a recently established linkage map of the Z chromosome of the collared flycatcher, to study the extent of LD in a natural bird population. In most but not all cases we find SNPs within the same intron (less than 500 bp) to be in perfect LD. However, LD then decays to background level at a distance 1cM or 400-500 kb. Although LD seems more extensive than in other species, if the observed pattern is representative for other regions of the genome and turns out to be a general feature of natural bird populations, dense marker maps might be needed for genome scans aimed at identifying association between marker and trait loci.     

Primary genetic linkage maps of the ascidian, Ciona intestinalis

For whole-genome analysis in a basal chordate (protochordate), we used F1 pseudo-testcross mapping strategy and amplified fragment length polymorphism (AFLP) markers to construct primary linkage maps of the ascidian tunicate Ciona intestinalis. Two genetic maps consisted of 14 linkage groups, in agreement with the haploid chromosome number, and contained 276 and 125 AFLP loci derived from crosses between British and Neapolitan individuals. The two maps covered 4218.9 and 2086.9 cM, respectively, with an average marker interval of 16.1 and 18.9 cM. We observed a high recombinant ratio, ranging from 25 to 49 kb/cM, which can explain the high degree of polymorphism in this species. Some AFLP markers were converted to sequence tagged sites (STSs) by sequence determination, in order to create anchor markers for the fragmental physical map. Our recombination tools provide basic knowledge of genetic status and whole genome organization, and genetic markers to assist positional cloning in C. intestinalis.