Pheno2Geno - High-throughput generation of genetic markers and maps from molecular phenotypes for crosses between inbred strains

Background

Genetic markers and maps are instrumental in quantitative trait locus (QTL) mapping in segregating populations. The resolution of QTL localization depends on the number of informative recombinations in the population and how well they are tagged by markers. Larger populations and denser marker maps are better for detecting and locating QTLs. Marker maps that are initially too sparse can be saturated or derived de novo from high-throughput omics data, (e.g. gene expression, protein or metabolite abundance). If these molecular phenotypes are affected by genetic variation due to a major QTL they will show a clear multimodal distribution. Using this information, phenotypes can be converted into genetic markers.

Results

The Pheno2Geno tool uses mixture modeling to select phenotypes and transform them into genetic markers suitable for construction and/or saturation of a genetic map. Pheno2Geno excludes candidate genetic markers that show evidence for multiple possibly epistatically interacting QTL and/or interaction with the environment, in order to provide a set of robust markers for follow-up QTL mapping.We demonstrate the use of Pheno2Geno on gene expression data of 370,000 probes in 148 A. thaliana recombinant inbred lines. Pheno2Geno is able to saturate the existing genetic map, decreasing the average distance between markers from 7.1 cM to 0.89 cM, close to the theoretical limit of 0.68 cM (with 148 individuals we expect a recombination every 100/148=0.68 cM); this pinpointed almost all of the informative recombinations in the population.

Conclusion

The Pheno2Geno package makes use of genome-wide molecular profiling and provides a tool for high-throughput de novo map construction and saturation of existing genetic maps. Processing of the showcase dataset takes less than 30 minutes on an average desktop PC. Pheno2Geno improves QTL mapping results at no additional laboratory cost and with minimum computational effort. Its results are formatted for direct use in R/qtl, the leading R package for QTL studies. Pheno2Geno is freely available on CRAN under “GNU GPL v3”. The Pheno2Geno package as well as the tutorial can also be found at: http://pheno2geno.nl.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0475-6) contains supplementary material, which is available to authorized users.     

Relationship of aerobic and anaerobic parameters with 400 m front crawl swimming performance

The aims of the present study were to investigate the relationship of aerobic and anaerobic parameters with 400 m performance, and establish which variable better explains long distance performance in swimming. Twenty-two swimmers (19.1±1.5 years, height 173.9±10.0 cm, body mass 71.2±10.2 kg; 76.6±5.3% of 400 m world record) underwent a lactate minimum test to determine lactate minimum speed (LMS) (i.e., aerobic capacity index). Moreover, the swimmers performed a 400 m maximal effort to determine mean speed (S400m), peak oxygen uptake (V.O2PEAK) and total anaerobic contribution (C_ANA). The C_ANA was assumed as the sum of alactic and lactic contributions. Physiological parameters of 400 m were determined using the backward extrapolation technique (V.O2PEAK and alactic contributions of C_ANA) and blood lactate concentration analysis (lactic anaerobic contributions of C_ANA). The Pearson correlation test and backward multiple regression analysis were used to verify the possible correlations between the physiological indices (predictor factors) and S400m (independent variable) (p < 0.05). Values are presented as mean ± standard deviation. Significant correlations were observed between S400m (1.4±0.1 m·s^-1) and LMS (1.3±0.1 m·s^-1; r = 0.80), V.O2PEAK (4.5±3.9 L·min^-1; r = 0.72) and C_ANA (4.7±1.5 L·O₂; r= 0.44). The best model constructed using multiple regression analysis demonstrated that LMS and V.O2PEAK explained 85% of the 400 m performance variance. When backward multiple regression analysis was performed, C_ANA lost significance. Thus, the results demonstrated that both aerobic parameters (capacity and power) can be used to predict 400 m swimming performance.     

Somatic mosaic activating mutations in PIK3CA cause CLOVES syndrome

Congenital lipomatous overgrowth with vascular, epidermal, and skeletal anomalies (CLOVES) is a sporadically occurring, nonhereditary disorder characterized by asymmetric somatic hypertrophy and anomalies in multiple organs. We hypothesized that CLOVES syndrome would be caused by a somatic mutation arising during early embryonic development. Therefore, we employed massively parallel sequencing to search for somatic mosaic mutations in fresh, frozen, or fixed archival tissue from six affected individuals. We identified mutations in PIK3CA in all six individuals, and mutant allele frequencies ranged from 3% to 30% in affected tissue from multiple embryonic lineages. Interestingly, these same mutations have been identified in cancer cells, in which they increase phosphoinositide-3-kinase activity. We conclude that CLOVES is caused by postzygotic activating mutations in PIK3CA. The application of similar sequencing strategies will probably identify additional genetic causes for sporadically occurring, nonheritable malformations.     

Use of polyethyleneimine polymer in cell culture as attachment factor and lipofection enhancer

Background  

Several cell lines and primary cultures benefit from the use of positively charged extracellular matrix proteins or polymers that enhance their ability to attach to culture plates. Polyethyleneimine is a positively charged polymer that has gained recent attention as a transfection reagent. A less known use of this cationic polymer as an attachment factor was explored with several cell lines.     

An AFLP-based genetic linkage map of Plasmodium chabaudi chabaudi

Background

Plasmodium chabaudi chabaudi can be considered as a rodent model of human malaria parasites in the genetic analysis of important characters such as drug resistance and immunity. Despite the availability of some genome sequence data, an extensive genetic linkage map is needed for mapping the genes involved in certain traits.

Methods

The inheritance of 672 Amplified Fragment Length Polymorphism (AFLP) markers from two parental clones (AS and AJ) of P. c. chabaudi was determined in 28 independent recombinant progeny clones. These, AFLP markers and 42 previously mapped Restriction Fragment Length Polymorphism (RFLP) markers (used as chromosomal anchors) were organized into linkage groups using Map Manager software.

Results

614 AFLP markers formed linkage groups assigned to 10 of 14 chromosomes, and 12 other linkage groups not assigned to known chromosomes. The genetic length of the genome was estimated to be about 1676 centiMorgans (cM). The mean map unit size was estimated to be 13.7 kb/cM. This was slightly less then previous estimates for the human malaria parasite, Plasmodium falciparum

Conclusion

The P. c. chabaudi genetic linkage map presented here is the most extensive and highly resolved so far available for this species. It can be used in conjunction with the genome databases of P. c chabaudi, P. falciparum and Plasmodium yoelii to identify genes underlying important phenotypes such as drug resistance and strain-specific immunity.     

Phylogenetic relationships within the Alcidae (Charadriiformes: Aves) inferred from total molecular evidence

The Alcidae is a unique assemblage of Northern Hemisphere seabirds that forage by "flying" underwater. Despite obvious affinities among the species, their evolutionary relationships are unclear. We analyzed nucleotide sequences of 1,045 base pairs of the mitochondrial cytochrome b gene and allelic profiles for 37 allozyme loci in all 22 extant species. Trees were constructed on independent and combined data sets using maximum parsimony and distance methods that correct for superimposed changes. Alternative methods of analysis produced only minor differences in relationships that were supported strongly by bootstrapping or standard error tests. Combining sequence and allozyme data into a single analysis provided the greatest number of relationships receiving strong support. Addition of published morphological and ecological data did not improve support for any additional relationship. All analyses grouped species into six distinct lineages: (1) the dovekie (Alle alle) and auks, (2) guillemots, (3) brachyramphine murrelets, (4) synthliboramphine murrelets, (5) true auklets, and (6) the rhinoceros auklet (Cerorhinca monocerata) and puffins. The two murres (genus Uria) were sister taxa, and the black guillemot (Cepphus grylle) was basal to the other guillemots. The Asian subspecies of the marbled murrelet (Brachyramphus marmoratus perdix) was the most divergent brachyramphine murrelet, and two distinct lineages occurred within the synthliboramphine murrelets. Cassin's auklet (Ptychoramphus aleuticus) and the rhinoceros auklet were basal to the other auklets and puffins, respectively, and the Atlantic (Fratercula arctica) and horned (Fratercula corniculata) puffins were sister taxa. Several relationships among tribes, among the dovekie and auks, and among the auklets could not be resolved but resembled "star" phylogenies indicative of adaptive radiations at different depths within the trees.      

Stress protein of cyanobacteria CP36: interaction with photoactive complexes and formation of supramolecular structures

Cyanobacterium Anacystis nidulans R2, Synechocystis sp. PCC 6803 (wild-type strain and mutants Delta2 and Delta3 lacking PSII and PSI, respectively), and Synechocystis sp. BO 9201 synthesize the pigment--protein complex CP36 (CPIV-4, CP43') under iron deficiency in the medium. Accumulation of CP36 is accompanied by structural reorganizations in the photosynthetic membranes. Integrating mean times of excitation relaxation (quenching) are 2.2 nsec (CP36), 1 nsec (PSI), and 420 psec (PSII in Fm state). The energy migration between CP36 and the photosystems can be described by a model of a one-layer ring of CP36 around core-complexes. The excitation from CP36 to PSI is transferred within <10 psec. The energy transfer from CP36 to PSII occurs during 170 psec. Cells with low content of CP36 probably contain only a latent fraction of unbound to phycobilisomes PSII which is the analog of PSIIbeta of higher plants. In PSI there are four binding sites for CP36 monomers per RC. PSII can bind up to 32 molecules of CP36 per RC. Cells with a large amount of CP36 contain monomer form of PSII core-complex which can bind eight tetramers of CP36 (8 binding sites). In conditions of iron deficiency only one monomer of a dimer PSII core-complex is destroyed and released chlorophyll is accumulated in CP36. Accumulation of CP36 in A. nidulans cells can be accompanied by membrane stacking which is similar to the stacking in chlorophyll b-containing organisms. The stacking can occur in the region of localization of PSII latent fraction bound to CP36. The membrane stacking shields PSII stromal surfaces from the aqueous phase for activation of electron transfer on the acceptor side of PSII.     

Relationship of In Vivo MR Parameters to Histopathological and Molecular Characteristics of Newly Diagnosed,Nonenhancing Lower-Grade Gliomas

The goal of this research was to elucidate the relationship between WHO 2016 molecular classifications of newly diagnosed, nonenhancing lower grade gliomas (LrGG), tissue sample histopathology, and magnetic resonance (MR) parameters derived from diffusion, perfusion, and ¹H spectroscopic imaging from the tissue sample locations and the entire tumor. A total of 135 patients were scanned prior to initial surgery, with tumor cellularity scores obtained from 88 image-guided tissue samples. MR parameters were obtained from corresponding sample locations, and histograms of normalized MR parameters within the T2 fluid-attenuated inversion recovery lesion were analyzed in order to evaluate differences between subgroups. For tissue samples, higher tumor scores were related to increased normalized apparent diffusion coefficient (nADC), lower fractional anisotropy (nFA), lower cerebral blood volume (nCBV), higher choline (nCho), and lower N-acetylaspartate (nNAA). Within the T2 lesion, higher tumor grade was associated with higher nADC, lower nFA, and higher Cho to NAA index. Pathological analysis confirmed that diffusion and metabolic parameters increased and perfusion decreased with tumor cellularity. This information can be used to select targets for tissue sampling and to aid in making decisions about treating residual disease.