A covering method for detecting genetic associations between rare variants and common phenotypes

Genome wide association (GWA) studies, which test for association between common genetic markers and a disease phenotype, have shown varying degrees of success. While many factors could potentially confound GWA studies, we focus on the possibility that multiple, rare variants (RVs) may act in concert to influence disease etiology. Here, we describe an algorithm for RV analysis, RareCover. The algorithm combines a disparate collection of RVs with low effect and modest penetrance. Further, it does not require the rare variants be adjacent in location. Extensive simulations over a range of assumed penetrance and population attributable risk (PAR) values illustrate the power of our approach over other published methods, including the collapsing and weighted-collapsing strategies. To showcase the method, we apply RareCover to re-sequencing data from a cohort of 289 individuals at the extremes of Body Mass Index distribution (NCT00263042). Individual samples were re-sequenced at two genes, FAAH and MGLL, known to be involved in endocannabinoid metabolism (187Kbp for 148 obese and 150 controls). The RareCover analysis identifies exactly one significantly associated region in each gene, each about 5 Kbp in the upstream regulatory regions. The data suggests that the RVs help disrupt the expression of the two genes, leading to lowered metabolism of the corresponding cannabinoids. Overall, our results point to the power of including RVs in measuring genetic associations.     

Integration of Brassica A genome genetic linkage map between Brassica napus and B. rapa.

An integrated linkage map between B. napus and B. rapa was constructed based on a total of 44 common markers comprising 41 SSR (33 BRMS, 6 Saskatoon, and 2 BBSRC) and 3 SNP/indel markers. Between 3 and 7 common markers were mapped onto each of the linkage groups A1 to A10. The position and order of most common markers revealed a high level of colinearity between species, although two small regions on A4, A5, and A10 revealed apparent local inversions between them. These results indicate that the A genome of Brassica has retained a high degree of colinearity between species, despite each species having evolved independently after the integration of the A and C genomes in the amphidiploid state. Our results provide a genetic integration of the Brassica A genome between B. napus and B. rapa. As the analysis employed sequence-based molecular markers, the information will accelerate the exploitation of the B. rapa genome sequence for the improvement of oilseed rape.     

A versatile computational pipeline for bacterial genome annotation improvement and comparative analysis, with Brucella as a use case

We present a bacterial genome computational analysis pipeline, called GenVar. The pipeline, based on the program GeneWise, is designed to analyze an annotated genome and automatically identify missed gene calls and sequence variants such as genes with disrupted reading frames (split genes) and those with insertions and deletions (indels). For a given genome to be analyzed, GenVar relies on a database containing closely related genomes (such as other species or strains) as well as a few additional reference genomes. GenVar also helps identify gene disruptions probably caused by sequencing errors. We exemplify GenVar's capabilities by presenting results from the analysis of four Brucella genomes. Brucella is an important human pathogen and zoonotic agent. The analysis revealed hundreds of missed gene calls, new split genes and indels, several of which are species specific and hence provide valuable clues to the understanding of the genome basis of Brucella pathogenicity and host specificity.     

Use of parents, sibs, and unrelated controls for detection of associations between genetic markers and disease.

Detecting the association between genetic markers and complex diseases can be a critical first step toward identification of the genetic basis of disease. Misleading associations can be avoided by choosing as controls the parents of diseased cases, but the availability of parents often limits this design to early-onset disease. Alternatively, sib controls offer a valid design. A general multivariate score statistic is presented, to detect the association between a multiallelic genetic marker locus and affection status; this general approach is applicable to designs that use parents as controls, sibs as controls, or even unrelated controls whose genotypes do not fit Hardy-Weinberg proportions or that pool any combination of these different designs. The benefit of this multivariate score statistic is that it will tend to be the most powerful method when multiple marker alleles are associated with affection status. To plan these types of studies, we present methods to compute sample size and power, allowing for varying sibship sizes, ascertainment criteria, and genetic models of risk. The results indicate that sib controls have less power than parental controls and that the power of sib controls can be increased by increasing either the number of affected sibs per sibship or the number of unaffected control sibs. The sample-size results indicate that the use of sib controls to test for associations, by use of either a single-marker locus or a genomewide screen, will be feasible for markers that have a dominant effect and for common alleles having a recessive effect. The results presented will be useful for investigators planning studies using sibs as controls.     

The Arabidopsis genome sequence as a tool for genome analysis in Brassicaceae. A comparison of the Arabidopsis and Capsella rubella genomes

The annotated Arabidopsis genome sequence was exploited as a tool for carrying out comparative analyses of the Arabidopsis and Capsella rubella genomes. Comparison of a set of random, short C. rubella sequences with the corresponding sequences in Arabidopsis revealed that aligned protein-coding exon sequences differ from aligned intron or intergenic sequences in respect to the degree of sequence identity and the frequency of small insertions/deletions. Molecular-mapped markers and expressed sequence tags derived from Arabidopsis were used for genetic mapping in a population derived from an interspecific cross between Capsella grandiflora and C. rubella. The resulting eight Capsella linkage groups were compared to the sequence maps of the five Arabidopsis chromosomes. Fourteen colinear segments spanning approximately 85% of the Arabidopsis chromosome sequence maps and 92% of the Capsella genetic linkage map were detected. Several fusions and fissions of chromosomal segments as well as large inversions account for the observed arrangement of the 14 colinear blocks in the analyzed genomes. In addition, evidence for small-scale deviations from genome colinearity was found. Colinearity between the Arabidopsis and Capsella genomes is more pronounced than has been previously reported for comparisons between Arabidopsis and different Brassica species.     

A race between tumor immunoescape and genome maintenance selects for optimum levels of (epi)genetic instability

The human immune system functions to provide continuous body-wide surveillance to detect and eliminate foreign agents such as bacteria and viruses as well as the body's own cells that undergo malignant transformation. To counteract this surveillance, tumor cells evolve mechanisms to evade elimination by the immune system; this tumor immunoescape leads to continuous tumor expansion, albeit potentially with a different composition of the tumor cell population ("immunoediting"). Tumor immunoescape and immunoediting are products of an evolutionary process and are hence driven by mutation and selection. Higher mutation rates allow cells to more rapidly acquire new phenotypes that help evade the immune system, but also harbor the risk of an inability to maintain essential genome structure and functions, thereby leading to an error catastrophe. In this paper, we designed a novel mathematical framework, based upon the quasispecies model, to study the effects of tumor immunoediting and the evolution of (epi)genetic instability on the abundance of tumor and immune system cells. We found that there exists an optimum number of tumor variants and an optimum magnitude of mutation rates that maximize tumor progression despite an active immune response. Our findings provide insights into the dynamics of tumorigenesis during immune system attacks and help guide the choice of treatment strategies that best inhibit diverse tumor cell populations.     

RFLP- and RAPD-based genetic relationships of seven diploid species of Avena with the A genome.

Relatively few molecular analyses are available for diploid oat species, which constitute the majority of the wild species of Avena and, therefore, the principal natural reservoir of variability. The present work reports an RAPD-(random amplified polymorphic DNA) and RFLP-(restriction fragment length polymorphism) based study of the intra- and interspecific variability of seven diploid A-genome oat species. Both types of markers resulted in valid tools for identifying polymorphisms both within and between species. The two statistical analyses, UPGMA (unweighted pair group method, arithmetic mean) and PCoA (principal coordinate analysis), computed on the basis of genetic similarities estimated from RAPDs and RFLPs, showed that the different accessions grouped according to species, but the similarity coefficients were consistently higher in the RFLP analysis. Furthermore, slight differences were observed in the intra- and interspecific relationships found with the two types of markers. This may support the hypothesis that the polymorphisms revealed by the two types of markers may associate with regions of the genome having different evolutionary rates. The relationships among species are not identical to those deduced from previous karyotypic and morphological studies, thus suggesting a partially different evolutionary pathway in oat speciation.     

RetroPred: A tool for prediction, classification and extraction of non-LTR retrotransposons (LINEs & SINEs) from the genome by integrating PALS, PILER, MEME and ANN

The problem of predicting non-long terminal repeats (LTR) like long interspersed nuclear elements (LINEs) and short interspersed nuclear elements (SINEs) from the DNA sequence is still an open problem in bioinformatics. To elevate the quality of annotations of LINES and SINEs an automated tool "RetroPred" was developed. The pipeline allowed rapid and thorough annotation of non-LTR retrotransposons. The non-LTR retrotransposable elements were initially predicted by Pairwise Aligner for Long Sequences (PALS) and Parsimonious Inference of a Library of Elementary Repeats (PILER). Predicted non-LTR elements were automatically classified into LINEs and SINEs using ANN based on the position specific probability matrix (PSPM) generated by Multiple EM for Motif Elicitation (MEME). The ANN model revealed a superior model (accuracy = 78.79 +/- 6.86 %, Q(pred) = 74.734 +/- 17.08 %, sensitivity = 84.48 +/- 6.73 %, specificity = 77.13 +/- 13.39 %) using four-fold cross validation. As proof of principle, we have thoroughly annotated the location of LINEs and SINEs in rice and Arabidopsis genome using the tool and is proved to be very useful with good accuracy. Our tool is accessible at http://www.juit.ac.in/RepeatPred/home.html.     

A comparison between biochemical tests and the 'CAMP' reaction for identification of Aeromonas spp.

Aeromonas bacteria (110 strains) from a variety of clinical, food and environmental sources, were identified using routine biochemical tests. Concurrently they were tested aerobically and anaerobically for their ability to perform synergistic haemolysis with Staphylococcus aureus (the 'CAMP' reaction). Results did not support a reported observation that the 'CAMP' reaction can he used to facilitate speciation of Aeromonas bacteria.     

A strategy for constructing high-resolution genetic maps of the human genome: a genetic map of chromosome 17p, ordered with meiotic breakpoint-mapping panels.

Genetic linkage analyses with genotypic data obtained from four CEPH reference families initially assigned 24 new PCR-based markers to chromosome 17 and located the markers at specific intervals of an existing genetic map of chromosome 17p. Each marker was additionally genotyped with an ordered set of obligate, phase-known recombinant chromosomes. The breakpoint-mapping panels for each family consisted of two parents, one sib with a nonrecombinant chromosome, and one or more sibs with obligate recombinant chromosomes. The relative order of markers was determined by sorting segregation patterns of new markers and ordered anchor markers and by minimizing double-recombination events. Consistency of segregation patterns with multiple flanking loci constituted support for order. A genetic map of chromosome 17p was completed with 39 markers in 23 clusters, with an average space of 3 cM between clusters. The collection of informative genotypes was highly efficient, requiring fivefold fewer genotypes than would be collected with all the CEPH families. Given the availability of large numbers of highly informative PCR-based markers, meiotic breakpoint mapping should facilitate construction of a human genomic map with 1-cM resolution.     

A multinational Andean genetic and health program. VIII. Lung function changes with migration between altitudes.

Studies of lung function in high altitude populations have suggested the influence of hypoxic environment on the development of this characteristic independent of confounding variables such as ethnicity and habitual exercise. However, often the effect of altitude on vital capacity is greater in children than adults, suggesting that more than developmental adaptation is operative. Also selective migration could account for the similarity of migrants and permanent residents at a destination altitude. To explore these problems we studied the lung function (FVC, FEV1, PFR) of 377 individuals who had migrated between altitudes in northern Chile. Migrant measurements were adjusted to those of permanent residents of appropriate age, sex and height at the altitudes of origin and destination. The measurements were then related to ethnicity (Spanish-Aymara ancestry), occupation and permanence, the latter combining information on both age at migration to and length of stay at a destination altitude. Upward migration was associated with increased chest depth, FVC and FEV1, but not height or other chest measurements. Downward migration had no significant effect. The flow-dependent test PFR was so sensitive to observer variability and occupation that it was difficult to establish its relationship to permanence. Unlike the body measurements, lung function measurements (especially PFR) tended to deviate from permanent controls at the origin altitude in a direction suggestive of selective migration, nor was permanence itself independent of ethnicity and occupation. Because of these difficulties the question of developmental adaptation in lung function may not be answerable in cross-sectional studies like the present and previous efforts, but rather in longitudinal investigations in which the control is the individual him/herself.     

Microsatellite markers for genome analysis in Brassica. II. Assignment of rapeseed microsatellites to the A and C genomes and genetic mapping in Brassica oleracea L.

Microsatellites are highly polymorphic and efficient markers for the analysis of plant genomes. Primer specificity, however, may restrict the applicability of these markers even between closely related species for comparative mapping studies. We have demonstrated that the majority of microsatellites identified in oilseed rape (Brassica napus L; AC genome) correspond to loci which can be easily assigned to the A and C progenitor genomes. A study with 63 primer pairs has shown that 54% detect two loci, one from each genome, while 25% and 21%, respectively, are either A or C genome-specific. The distribution of rapeseed microsatellites in the C genome was investigated by genetic mapping in Brassica oleracea L. Ninety two dinucleotide microsatellites were screened for polymorphism in an F₂ population derived from a cross between collard and cauliflower, for which an RFLP map has been constructed previously. Thirty three primer pairs (35.7%) have yielded either unspecific or no PCR products whereas the remaining primer pairs amplified one or more distinct loci. The level of polymorphism found in the mapping population was 49.2%. A total of 29 primer pairs disclosed 34 loci of which 31 are evenly distributed on 8 of the 9 B. oleracea linkage groups. For the remaining three markers linkage could not be established. Our results showed that microsatellite markers from the composite genome of B. napus can serve as a useful marker system in genetic studies and for plant-breeding objectives in B. oleracea. Received: 14 April 2000 / Accepted: 3 July 2000     

Thinopyrum bessarabicum: biochemical and cytological markers for the detection of genetic introgression in its hybrid derivatives with Triticum aestivum L.

Summary Thinopyrum bessarabicum (2n = 2x = 14, JJ) with its unique property of salt tolerance provides a potential means for the transfer of this important and complex trait into cultivated wheat through intergeneric hybridization. To accomplish this, diagnostic markers for detecting the presence of Th. bessarabicum chromosomes in a wheat background have to be established. The C-banded karyotype of Th. bessarabicum distinctly identifies individual Th. bessarabicum chromosomes and separates them from those of Triticum aestivum. Also, seven protein/isozymes, i.e., malate dehydrogenase, high-molecular-weight glutenin, Superoxide dismutase, grain esterase, glutamate oxaloacetate transaminase, -amylase and -amylase, were identified as being positive markers specific to Th. bessarabicum; these were also expressed in the T. aestivum/Th. bessarabicum amphiploid. These diagnostic biochemical markers could be useful in detecting and establishing homoeology of Th. bessarabicum chromosomes in T. aestivum/Th. bessarabicum intergeneric hybrid derivatives.     

A convenient and rapid method for genetic transformation of E. coli with plasmids

A convenient and rapid method for the genetic transformation of Escherichia coli with plasmids is proposed. By mixing the recipient cells and plasmid DNA and spreading them directly on selective medium plates containing Ca2+, the so-called 'plate transformation' could achieve almost the same transformation efficiency as the classical transformation method with calcium. The whole protocol takes only about 2 min, its simplicity compared favorably, not only to the usual protocol, but also to all other documented modifications.     

Development of a set of near isogenic and backcross recombinant inbred lines containing most of the Lycopersicon hirsutum genome in a L. esculentum genetic background: a tool for gene mapping and gene discovery.

A novel population consisted of a set of 99 near isogenic lines (NILs) and backcross recombinant inbred lines (BCRILs) derived from a cross between the cultivated tomato Lycopersicon esculentum cv. E6206 and L. hirsutum accession LA1777 is presented. Most of the lines contain a single defined introgression from L. hirsutum in the L. esculentum genetic background and together, the lines provide a coverage of more than the 85% of the L. hirsutum genome. These lines represent a new tool to uncover the genetic resources hidden in L. hirsutum as well as to study the genes responsible of its unique biology. Furthermore, the study of the allelic frequency and heterozygosity among BCRILs showed that specific genomic regions were likely subjected to unintentional selection pressures during the stock development. Genes involved in the reproductive behavior and (or) pollen viability are hypothesized to be responsible for these alterations.     

InterPro as a new tool for whole genome analysis. A comparative analysis of Mycobacterium tuberculosis, Bacillus subtilis and Escherichia coli as a case study

InterPro was developed as a new integrated documentation resource for protein families, domains and functional sites to rationalize the complementary efforts of the PROSITE, PRINTS, Pfam and ProDom database projects and has applications in computational functional classification of newly determined sequences lacking biochemical characterization and in comparative genome analysis. InterPro contains over 3500 entries, with more than 1000000 hits in SWISS-PROT and TrEMBL. The database is accessible for text- and sequence-based searches at http://www.ebi.ac.uk/interpro/. InterPro was used for whole proteome analysis of the pathogenic microorganism, Mycobacterium tuberculosis, and comparison with the predicted protein coding sequences of the complete genomes of Bacillus subtilis and Escherichia coli. 64.8% of the M. tuberculosis proteins in the proteome matched InterPro entries, and these could be classified according to function. The comparison with B. subtilis and E. coli provided information on the most common protein families and domains, and the most highly represented families in each organism. InterPro thus provides a useful tool for global views of whole proteomes and their compositions.