Genome Scan for Predisposing Loci for Distal Interphalangeal Joint Osteoarthritis: Evidence for a Locus on 2q

The genetic contribution to common forms of osteoarthritis (OA) is well established but poorly understood. We performed a genome scan, using 302 markers for loci predisposing to distal interphalangeal joint (DIP) OA. To minimize genetic heterogeneity in our study sample, we identified siblings with a severe, radiologically defined phenotype from the nationwide registers of Finland. In the initial genome scan, linkage analysis in 27 sibships gave a pairwise LOD score (Z) >1.00 with nine of the screening markers. In the second stage, additional markers and family members were genotyped in these chromosomal regions. On 2q12-q13, IL1R1 resulted in Z=2.34 at recombination fraction (theta) 0, allowing a dominant mode of inheritance. Association analysis of markers D2S2264, IL1R1, D2S373, and D2S1789 jointly provided some evidence for a shared haplotype among the affected individuals (P value of.012). Also, multipoint nonparametric linkage analysis yielded a P value of.0001 near the locus IL1R1 and P=.0007 approximately 20 cM telomeric near marker D2S1399, which, in two-point analysis, gave Z=1.48 (straight theta=. 02). This chromosomal region on 2q harbors the interleukin 1 gene cluster and, thus, represents a good candidate region for inflammatory and autoimmune disorders. Three additional chromosomal regions-4q26-q27, 7p15-p21, and Xcen-also provided some evidence for linkage, and further analyses would be justified to clarify their potential involvement in the genetic predisposition to DIP OA.     

A Genome Scan for Selection Signatures in Pigs

Identifying signatures of selection can provide a straightforward insight into the mechanism of artificial selection and further uncover the causal genes related to the phenotypic variation. Based on Illumina Porcine60KSNP chip data, four complementary methods, Long-Range Haplotype (LRH), Tajima’s D, Cross Population Extend Haplotype Homozygosity Test (XPEHH) and F_ST, were implemented in this study to detect the selection signatures in the whole genome of one typical Chinese indigenous breed, Rongchang, one Chinese cultivated breed, Songliao, and two western breeds, Landrace and Yorkshire. False Discovery Rate (FDR) was implemented to control the false positive rates. In our study, a total of 159, 127, 179 and 159 candidate selection regions with average length of 0.80 Mb, 0.73 Mb, 0.78 Mb and 0.73 Mb were identified in Landrace, Rongchang, Songliao and Yorkshire, respectively, that span approximately 128.00 Mb, 92.38 Mb, 130.30 Mb and 115.40 Mb and account for approximately 3.74–5.33% of genome across all autosomes. The selection regions of 11.52 Mb shared by Landrace and Yorkshire were the longest when chosen pairs from the pool of the four breeds were examined. The overlaps between Yorkshire and Songliao, approximately 9.20 Mb, were greater than those of Yorkshire and Rongchang. Meanwhile, the overlaps between Landrace and Songliao were greater than those of Landrace and Rongchang but less than those of Songliao and Ronchang. Bioinformatics analysis showed that the genes/QTLs relevant to fertility, coat color, and ear morphology were found in candidate selection regions. Some genes, such as LEMD3, MC1R, KIT, TRHR etc. that were reported under selection, were confirmed in our study, and this analysis also demonstrated the diversity of breeds.     

A Genome Scan for Positive Selection in Thoroughbred Horses

Thoroughbred horses have been selected for exceptional racing performance resulting in system-wide structural and functional adaptations contributing to elite athletic phenotypes. Because selection has been recent and intense in a closed population that stems from a small number of founder animals Thoroughbreds represent a unique population within which to identify genomic contributions to exercise-related traits. Employing a population genetics-based hitchhiking mapping approach we performed a genome scan using 394 autosomal and X chromosome microsatellite loci and identified positively selected loci in the extreme tail-ends of the empirical distributions for (1) deviations from expected heterozygosity (Ewens-Watterson test) in Thoroughbred (n = 112) and (2) global differentiation among four geographically diverse horse populations (F_ST). We found positively selected genomic regions in Thoroughbred enriched for phosphoinositide-mediated signalling (3.2-fold enrichment; P<0.01), insulin receptor signalling (5.0-fold enrichment; P<0.01) and lipid transport (2.2-fold enrichment; P<0.05) genes. We found a significant overrepresentation of sarcoglycan complex (11.1-fold enrichment; P<0.05) and focal adhesion pathway (1.9-fold enrichment; P<0.01) genes highlighting the role for muscle strength and integrity in the Thoroughbred athletic phenotype. We report for the first time candidate athletic-performance genes within regions targeted by selection in Thoroughbred horses that are principally responsible for fatty acid oxidation, increased insulin sensitivity and muscle strength: ACSS1 (acyl-CoA synthetase short-chain family member 1), ACTA1 (actin, alpha 1, skeletal muscle), ACTN2 (actinin, alpha 2), ADHFE1 (alcohol dehydrogenase, iron containing, 1), MTFR1 (mitochondrial fission regulator 1), PDK4 (pyruvate dehydrogenase kinase, isozyme 4) and TNC (tenascin C). Understanding the genetic basis for exercise adaptation will be crucial for the identification of genes within the complex molecular networks underlying obesity and its consequential pathologies, such as type 2 diabetes. Therefore, we propose Thoroughbred as a novel in vivo large animal model for understanding molecular protection against metabolic disease.     

Substantial Genetic Effects Involved in Determination of Circulating Levels of Calciotropic Hormones in Human Pedigrees

This paper reports the results of a series of univariate and bivariate statistical genetic analyses that were performed on a sample of nuclear and more complex pedigrees (N = 567 individuals) of an ethnically homogenous White population. Our major objectives were: (1) To quantitatively evaluate the extent and pattern of the putative genetic effects on plasma level variation and covariation of the intact parathyroid hormone (PTH) and 25-hydroxyvitamin D [25(OH)D]; (2) To evaluate the extent of the possible genetic covariation between each of the two calciotropic hormones and two important bone mass traits, namely radiographic hands bone mineral density (BMD) and cortical index (CI). Variance component analysis, as implemented in the statistical package FISHER unambigously demonstrated that in addition to age, genetic factors contribute significantly to interindividual variation of both calciotropic hormones (37.5% for PTH and 53.3% for 25(OH)D). Complex segregation analysis strongly suggested the involvement of major gene effects into the determination of 25(OH)D levels, but was not clear cut with respect to PTH. Significant correlations between circulating levels of study hormones were found (-0.146, P < 0.05 in men and -0.194, P < 0.01 in women). However, no genetic correlation was revealed between PTH and 25(OH)D plasma concentrations. Bivariate analyses showed that familial cross correlations between PTH and BMD and CI measured at the bones of the hand were consistently statistically significant, suggesting moderate, but detectable pleiotropic genetic effects. The corresponding genetic correlations were -0.461 +/- 0.153 and -0.223 +/- 0.113, respectively. Circulating levels of 25(OH)D showed neither phenotypic nor genetic correlation with any of the BMD or CI variation.     

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.     

A New Genetic Locus in Sinorhizobium meliloti Is Involved in Stachydrine Utilization

Stachydrine, a betaine released by germinating alfalfa seeds, functions as an inducer of nodulation genes, a catabolite, and an osmoprotectant in Sinorhizobium meliloti. Two stachydrine-inducible genes were found in S. meliloti 1021 by mutation with a Tn5-luxAB promoter probe. Both mutant strains (S10 and S11) formed effective alfalfa root nodules, but neither grew on stachydrine as the sole carbon and nitrogen source. When grown in the absence or presence of salt stress, S10 and S11 took up [¹⁴C]stachydrine as well as wild-type cells did, but neither used stachydrine effectively as an osmoprotectant. In the absence of salt stress, both S10 and S11 took up less [¹⁴C]proline than wild-type cells did. S10 and S11 appeared to colonize alfalfa roots normally in single-strain tests, but when mixed with the wild-type strain, their rhizosphere counts were reduced more than 50% (P ≤ 0.01) relative to the wild type. These results suggest that stachydrine catabolism contributes to root colonization. DNA sequence analysis identified the mutated locus in S11 as putA, and the luxAB fusion in that gene was induced by proline as well as stachydrine. DNA that restored the capacity of mutant S10 to catabolize stachydrine contained a new open reading frame, stcD. All data are consistent with the concept that stcD codes for an enzyme that produces proline by demethylation of N-methylproline, a degradation product of stachydrine.     

Genome Wide Association Identifies Novel Loci Involved in Fungal Communication

Understanding how genomes encode complex cellular and organismal behaviors has become the outstanding challenge of modern genetics. Unlike classical screening methods, analysis of genetic variation that occurs naturally in wild populations can enable rapid, genome-scale mapping of genotype to phenotype with a medium-throughput experimental design. Here we describe the results of the first genome-wide association study (GWAS) used to identify novel loci underlying trait variation in a microbial eukaryote, harnessing wild isolates of the filamentous fungus Neurospora crassa. We genotyped each of a population of wild Louisiana strains at 1 million genetic loci genome-wide, and we used these genotypes to map genetic determinants of microbial communication. In N. crassa, germinated asexual spores (germlings) sense the presence of other germlings, grow toward them in a coordinated fashion, and fuse. We evaluated germlings of each strain for their ability to chemically sense, chemotropically seek, and undergo cell fusion, and we subjected these trait measurements to GWAS. This analysis identified one gene, NCU04379 (cse-1, encoding a homolog of a neuronal calcium sensor), at which inheritance was strongly associated with the efficiency of germling communication. Deletion of cse-1 significantly impaired germling communication and fusion, and two genes encoding predicted interaction partners of CSE1 were also required for the communication trait. Additionally, mining our association results for signaling and secretion genes with a potential role in germling communication, we validated six more previously unknown molecular players, including a secreted protease and two other genes whose deletion conferred a novel phenotype of increased communication and multi-germling fusion. Our results establish protein secretion as a linchpin of germling communication in N. crassa and shed light on the regulation of communication molecules in this fungus. Our study demonstrates the power of population-genetic analyses for the rapid identification of genes contributing to complex traits in microbial species.     

Identification of a Genetic Locus in Pseudomonas aureofaciens Involved in Fungal Inhibition

In iron-rich conditions, Pseudomonas aureofaciens PA147-2 produces an antibiotic-like compound that inhibits the growth of a plant fungal pathogen, Aphanomyces euteiches. To contribute to the potential use of PA147-2 as a biocontrol organism, we report the identification of a genetic locus important for antibiotic biosynthesis. Mutants defective for fungal inhibition (Af^-) were generated by Tn5 mutagenesis. Southern hybridization of total DNAs from three Af^- mutants indicated that loss of fungal inhibition was due to a single Tn5 insertion in each mutant. Restriction mapping of the mutation points showed that in two mutants the Tn5 insertions were in the same 16.0-kb EcoRI fragment and were separated by 2.1 kb. A genomic library of PA147-2 was constructed and screened by using a region of DNA flanking the Tn5 insertion in one mutant (PA109) as a probe to recover complementing cosmids. Three cosmids containing a 16.0-kb EcoRI fragment complementary to the two mutants were recovered. Allele replacement by homologous recombination with putative complementing cosmids restored one mutant to antifungal activity against A. euteiches. Southern analysis of the complemented mutants confirmed that allele replacement had occurred between cosmid DNA and Tn5. The wild-type 16.0-kb EcoRI fragment was cloned from the cosmid and complemented the two mutants to antifungal activity. An antifungal compound was isolated from PA147-2 grown on solid medium. Antifungal activity correlated to a peak on high-pressure liquid chromatography analysis. Under the same growth and extraction conditions, the antifungal activity seen in PA147-2 was absent in two Af^- mutants. Furthermore, absence of an antifungal compound in each mutant correlated to the absence of the wild-type “antifungal” peak on high-pressure liquid chromatography analysis.     

Genome Scan of M. tuberculosis Infection and Disease in Ugandans

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is an enduring public health problem globally, particularly in sub-Saharan Africa. Several studies have suggested a role for host genetic susceptibility in increased risk for TB but results across studies have been equivocal. As part of a household contact study of Mtb infection and disease in Kampala, Uganda, we have taken a unique approach to the study of genetic susceptibility to TB, by studying three phenotypes. First, we analyzed culture confirmed TB disease compared to latent Mtb infection (LTBI) or lack of Mtb infection. Second, we analyzed resistance to Mtb infection in the face of continuous exposure, defined by a persistently negative tuberculin skin test (PTST-); this outcome was contrasted to LTBI. Third, we analyzed an intermediate phenotype, tumor necrosis factor-alpha (TNFα) expression in response to soluble Mtb ligands enriched with molecules secreted from Mtb (culture filtrate). We conducted a full microsatellite genome scan, using genotypes generated by the Center for Medical Genetics at Marshfield. Multipoint model-free linkage analysis was conducted using an extension of the Haseman-Elston regression model that includes half sibling pairs, and HIV status was included as a covariate in the model. The analysis included 803 individuals from 193 pedigrees, comprising 258 full sibling pairs and 175 half sibling pairs. Suggestive linkage (p<10⁻³) was observed on chromosomes 2q21-2q24 and 5p13-5q22 for PTST-, and on chromosome 7p22-7p21 for TB; these findings for PTST- are novel and the chromosome 7 region contains the IL6 gene. In addition, we replicated recent linkage findings on chromosome 20q13 for TB (p = 0.002). We also observed linkage at the nominal α = 0.05 threshold to a number of promising candidate genes, SLC11A1 (PTST- p = 0.02), IL-1 complex (TB p = 0.01), IL12BR2 (TNFα p = 0.006), IL12A (TB p = 0.02) and IFNGR2 (TNFα p = 0.002). These results confirm not only that genetic factors influence the interaction between humans and Mtb but more importantly that they differ according to the outcome of that interaction: exposure but no infection, infection without progression to disease, or progression of infection to disease. Many of the genetic factors for each of these stages are part of the innate immune system.     

A Genomewide Scan Identifies Two Novel Loci Involved in Specific Language Impairment

Approximately 4% of English-speaking children are affected by specific language impairment (SLI), a disorder in the development of language skills despite adequate opportunity and normal intelligence. Several studies have indicated the importance of genetic factors in SLI; a positive family history confers an increased risk of development, and concordance in monozygotic twins consistently exceeds that in dizygotic twins. However, like many behavioral traits, SLI is assumed to be genetically complex, with several loci contributing to the overall risk. We have compiled 98 families drawn from epidemiological and clinical populations, all with probands whose standard language scores fall 1.5 SD below the mean for their age. Systematic genomewide quantitative-trait–locus analysis of three language-related measures (i.e., the Clinical Evaluation of Language Fundamentals–Revised [CELF-R] receptive and expressive scales and the nonword repetition [NWR] test) yielded two regions, one on chromosome 16 and one on 19, that both had maximum LOD scores of 3.55. Simulations suggest that, of these two multipoint results, the NWR linkage to chromosome 16q is the most significant, with empirical P values reaching 10⁻⁵, under both Haseman-Elston (HE) analysis (LOD score 3.55; P=.00003) and variance-components (VC) analysis (LOD score 2.57; P=.00008). Single-point analyses provided further support for involvement of this locus, with three markers, under the peak of linkage, yielding LOD scores >1.9. The 19q locus was linked to the CELF-R expressive-language score and exceeds the threshold for suggestive linkage under all types of analysis performed—multipoint HE analysis (LOD score 3.55; empirical P=.00004) and VC (LOD score 2.84; empirical P=.00027) and single-point HE analysis (LOD score 2.49) and VC (LOD score 2.22). Furthermore, both the clinical and epidemiological samples showed independent evidence of linkage on both chromosome 16q and chromosome 19q, indicating that these may represent universally important loci in SLI and, thus, general risk factors for language impairment.     

Genetic Evidence That the Sans Fille Locus Is Involved in Drosophila Sex Determination

Females homozygous for sans fille1621 (= fs(1)1621) have an abnormal germ line. Instead of producing eggs, the germ-line cells proliferate forming ovarian tumors or excessive numbers of nurse cells. The Sex-lethal gene product(s) regulate the branch point of the dosage compensation and sex determination pathways in the soma. The role of Sex-lethal in the germ line is not clear but the germ line of females homozygous for female sterile Sex-lethal alleles or germ-line clones of loss-of-function alleles are characterized by ovarian tumors. Females heterozygous for sans fille1621 or Sex-lethal are phenotypically wild type with respect to viability and fertility but females trans-heterozygous for sans fille1621 and Sex-lethal show ovarian tumors, somatic sexual transformations, and greatly reduced viability.     

Draft Genome Sequences of Two Streptococcus pyogenes Strains Involved in Abnormal Sharp Raised Scarlet Fever in China, 2011

A scarlet fever outbreak caused by Streptococcus pyogenes occurred in China in 2011. To determine the genomic features of the outbreak strains, we deciphered genomes of two strains isolated from the regions with the highest incidence rates. The sequences will provide valuable information for comprehensive study of mechanisms related to this outbreak.     

Cloning of Rhizobium meliloti Locus sxf C Involved in Symbiosis and Phage Sensitivity

A symbiotic mutant of Rhizobium meliloti Rmd438 (sxf C:: Tn5) which was phage resistant against RMP64, failed to utilize galactose as carbon source as reported earlier (21). The Bg/ll gene bank of wild type R. meliloti was mobilized into Rmd438 and a clone pSP676 which complemented for phage sensitivity was isolated. In order to characterize this clone, a Bg/ll and EcoRI map was constructed. The insert of 13.2 kb had three Bg/ll fragments of 4.0, 3.6 and 5.6 kb in this order. All three fragments were subcloned on the vector pRK290 and mobilized into Sxf-mutants. The complementation for phage sensitivity, symbiosis and galactose utilization properties are discussed.     

Genetic Evidence That the Ovo Locus Is Involved in Drosophila Germ Line Sex Determination

Zygotically contributed ovo gene product is required for the survival of female germ cells in Drosophila melanogaster. Trans-allelic combinations of weak and dominant ovo mutations (ovoD) result in viable germ cells that appear to be partially transformed from female to male sexual identity. The ovoD2 mutation is partially suppressed by many Sex-lethal alleles that affect the soma, while those that affect only the germ line fail to interact with ovoD2. One of two loss-of-function ovo alleles is suppressed by a loss-of-function Sex-lethal allele. Because ovo mutations are germ line dependent, it is likely that ovo is suppressed by way of communication between the somatic and germ lines. A loss-of-function allele of ovo is epistatic to germ line dependent mutations in Sex-lethal. The germ line dependent sex determination mutation, sans fille, and ovoD mutations show a dominant synergistic interaction resulting in partial transformation of germ line sexual identity. The ovo locus appears to be involved in germ line sex determination and is linked in some manner to sex determination in the soma.     

Quantitative Trait Locus Based Virulence Determinant Mapping of the HSV-1 Genome in Murine Ocular Infection: Genes Involved in Viral Regulatory and Innate Immune Networks Contribute to Virulence

Herpes simplex virus type 1 causes mucocutaneous lesions, and is the leading cause of infectious blindness in the United States. Animal studies have shown that the severity of HSV-1 ocular disease is influenced by three main factors; innate immunity, host immune response and viral strain. We previously showed that mixed infection with two avirulent HSV-1 strains (OD4 and CJ994) resulted in recombinants that exhibit a range of disease phenotypes from severe to avirulent, suggesting epistatic interactions were involved. The goal of this study was to develop a quantitative trait locus (QTL) analysis of HSV-1 ocular virulence determinants and to identify virulence associated SNPs. Blepharitis and stromal keratitis quantitative scores were characterized for 40 OD4:CJ994 recombinants. Viral titers in the eye were also measured. Virulence quantitative trait locus mapping (vQTLmap) was performed using the Lasso, Random Forest, and Ridge regression methods to identify significant phenotypically meaningful regions for each ocular disease parameter. The most predictive Ridge regression model identified several phenotypically meaningful SNPs for blepharitis and stromal keratitis. Notably, phenotypically meaningful nonsynonymous variations were detected in the UL24, UL29 (ICP8), UL41 (VHS), UL53 (gK), UL54 (ICP27), UL56, ICP4, US1 (ICP22), US3 and gG genes. Network analysis revealed that many of these variations were in HSV-1 regulatory networks and viral genes that affect innate immunity. Several genes previously implicated in virulence were identified, validating this approach, while other genes were novel. Several novel polymorphisms were also identified in these genes. This approach provides a framework that will be useful for identifying virulence genes in other pathogenic viruses, as well as epistatic effects that affect HSV-1 ocular virulence.     

A Genomic Scan for Selection Reveals Candidates for Genes Involved in the Evolution of Cultivated Sunflower (Helianthus annuus)

Genomic scans for selection are a useful tool for identifying genes underlying phenotypic transitions. In this article, we describe the results of a genome scan designed to identify candidates for genes targeted by selection during the evolution of cultivated sunflower. This work involved screening 492 loci derived from ESTs on a large panel of wild, primitive (i.e., landrace), and improved sunflower (Helianthus annuus) lines. This sampling strategy allowed us to identify candidates for selectively important genes and investigate the likely timing of selection. Thirty-six genes showed evidence of selection during either domestication or improvement based on multiple criteria, and a sequence-based test of selection on a subset of these loci confirmed this result. In view of what is known about the structure of linkage disequilibrium across the sunflower genome, these genes are themselves likely to have been targeted by selection, rather than being merely linked to the actual targets. While the selection candidates showed a broad range of putative functions, they were enriched for genes involved in amino acid synthesis and protein catabolism. Given that a similar pattern has been detected in maize (Zea mays), this finding suggests that selection on amino acid composition may be a general feature of the evolution of crop plants. In terms of genomic locations, the selection candidates were significantly clustered near quantitative trait loci (QTL) that contribute to phenotypic differences between wild and cultivated sunflower, and specific instances of QTL colocalization provide some clues as to the roles that these genes may have played during sunflower evolution.     

Complete Genome Sequences of Yersinia pestis from Natural Foci in China

Yersinia pestis, the causative agent of plague, is a deadly bacterium that affects humans. Strain D106004 was isolated from a new plague focus in Yulong County, China, in 2006. To gain insights into the epidemic origin, we have sequenced the genomes of D106004 and strains Z176003 and D182038, isolated from neighboring regions.This article describes genomic comparisons between three respective Yersinia pestis strains isolated from new natural plague foci in China. Y. pestis strain D106004 was isolated from Apodemus chevrieri in Yulong County in 2006, and its genome was compared to those of strain D182038 (isolated from A. chevrieri in 1982 from Jianchuan County) and strain Z176003 (isolated from Marmota himalayana in 1976 in Naqu [Tibet] County).Between 25 October 2005 and 2 November 2005, there was an outbreak of pneumonic plague in Yulong, which was identified as a new natural plague focus (13). The primary Y. pestis reservoirs associated with this outbreak were A. chevrieri, Eothenmys miletus, and Apodemus latronum, and the primary vectors associated with plague transmission were also identified as similar to what was observed in neighboring Jianchuan County (7). However, the Y. pestis strain identified metabolized maltose significantly differently than the previously described strains (6).Whole-genome shotgun and solexa methods were used, as previously described (3), to compare the Y. pestis D106004, D182038, and Z176003 sequences, which consisted of 475, 385, and 413 contigs, respectively, resulting in an average 9-fold coverage across the genomes. All isolates examined possessed a single circular chromosome with the three virulence plasmids (pMT, pCD, and pPCP) associated with classical Y. pestis strains. Automated gene modeling was carried out using the Glimmer3 software program (11) in addition to comparing the respective gene products using the Nt, Nr, KEGG, Swissprot, and COG databases using the basic local alignment search tool for proteins (BLASTP). Open reading frames (ORFs) in the respective 4,626,944-bp, 4,640,720-bp, and 4,553,586-bp genomes of strains D182038, D106004, and Z176003 were predicted to be of 3,642, 3,636, 3,543, and more than 300 bp in length. Strains D182038, D106004, and Z176003 each had six rRNA (16S-23S-5S) genes and 73 (D182038), 70 (D106004), or 68 (Z176003) tRNA genes predicted by the tRNAScan-SE server (9).Comparison of Y. pestis strains 91001 and KIM to Y. pestis strain CO92 identified genetic rearrangements (5, 10, 12) resulting from insertion sequences (2), and pulsed-field gel electrophoresis (PFGE) profile comparisons between D182038 and D106004 suggested that genomic variability of the Y. pestis strains from different foci was caused by genome rearrangement (16). According to our analyses, the Y. pestis strains isolated from the two foci have very different syntenic structures due to rearrangement, but they share high similarity between plates (8). In addition, a unique multiple-locus variable-number tandem repeat analysis (MLVA) type was defined for the strains isolated from Yulong, indicating a new clonal group. These results also suggested that the Yulong strains were closely related to the strains from the Qinghai-Tibet Plateau plague foci (15). Analysis of Y. pestis microevolution has been made possible by comparing single- nucleotide polymorphism (SNP) profiles as previously described (1, 4, 14).The availability of high-quality sequences is crucial in order to resolve the origins of the new strains isolated from natural plague foci.     

Filtering for Compound Heterozygous Sequence Variants in Non-Consanguineous Pedigrees

The identification of disease-causing mutations in next-generation sequencing (NGS) data requires efficient filtering techniques. In patients with rare recessive diseases, compound heterozygosity of pathogenic mutations is the most likely inheritance model if the parents are non-consanguineous. We developed a web-based compound heterozygous filter that is suited for data from NGS projects and that is easy to use for non-bioinformaticians. We analyzed the power of compound heterozygous mutation filtering by deriving background distributions for healthy individuals from different ethnicities and studied the effectiveness in trios as well as more complex pedigree structures. While usually more then 30 genes harbor potential compound heterozygotes in single exomes, this number can be markedly reduced with every additional member of the pedigree that is included in the analysis. In a real data set with exomes of four family members, two sisters affected by Mabry syndrome and their healthy parents, the disease-causing gene PIGO, which harbors the pathogenic compound heterozygous variants, could be readily identified. Compound heterozygous filtering is an efficient means to reduce the number of candidate mutations in studies aiming at identifying recessive disease genes in non-consanguineous families. A web-server is provided to make this filtering strategy available at www.gene-talk.de.