POCUS: mining genomic sequence annotation to predict disease genes

Here we present POCUS (prioritization of candidate genes using statistics), a novel computational approach to prioritize candidate disease genes that is based on over-representation of functional annotation between loci for the same disease. We show that POCUS can provide high (up to 81-fold) enrichment of real disease genes in the candidate-gene shortlists it produces compared with the original large sets of positional candidates. In contrast to existing methods, POCUS can also suggest counterintuitive candidates.     

Microsatellite DNA markers for the rainbow darter,Etheostoma caeruleum (Percidae), and their potential utility for other darter species

The rainbow darter Etheostoma caeruleum is a small fish in the perch family (Percidae) that is adapted to fast‐flowing streams in eastern North America. It is relatively sensitive to habitat degradation and is widely used as a sentinel of stream condition. To provide a complementary tool for assessing the integrity of stream ecosystems, 16 highly polymorphic tetranucleotide microsatellite markers were identified for these darters. Between four and 16 loci were found to be useful in five congeneric species. These markers will be useful for characterizing population genetic structure and diversity of rainbow darters and related fishes.     

Biodiversity of 52 chicken populations assessed by microsatellite typing of DNA pools

In a project on the biodiversity of chickens funded by the European Commission (EC), eight laboratories collaborated to assess the genetic variation within and between 52 populations from a wide range of chicken types. Twenty-two di-nucleotide microsatellite markers were used to genotype DNA pools of 50 birds from each population. The polymorphism measures for the average, the least polymorphic population (inbred C line) and the most polymorphic population (Gallus gallus spadiceus) were, respectively, as follows: number of alleles per locus, per population: 3.5, 1.3 and 5.2; average gene diversity across markers: 0.47, 0.05 and 0.64; and proportion of polymorphic markers: 0.91, 0.25 and 1.0. These were in good agreement with the breeding history of the populations. For instance, unselected populations were found to be more polymorphic than selected breeds such as layers. Thus DNA pools are effective in the preliminary assessment of genetic variation of populations and markers. Mean genetic distance indicates the extent to which a given population shares its genetic diversity with that of the whole tested gene pool and is a useful criterion for conservation of diversity. The distribution of population-specific (private) alleles and the amount of genetic variation shared among populations supports the hypothesis that the red jungle fowl is the main progenitor of the domesticated chicken.     

Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken

Background

Toll like receptors (TLR) play the central role in the recognition of pathogen associated molecular patterns (PAMPs). Mutations in the TLR1, TLR2 and TLR4 genes may change the ability to recognize PAMPs and cause altered responsiveness to the bacterial pathogens.

Results

The study presents association between TLR gene mutations and increased susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection. Novel mutations in TLR genes (TLR1- Ser150Gly and Val220Met; TLR2 – Phe670Leu) were statistically correlated with the hindrance in recognition of MAP legends. This correlation was confirmed subsequently by measuring the expression levels of cytokines (IL-4, IL-8, IL-10, IL-12 and IFN-γ) in the mutant and wild type moDCs (mocyte derived dendritic cells) after challenge with MAP cell lysate or LPS. Further in silico analysis of the TLR1 and TLR4 ectodomains (ECD) revealed the polymorphic nature of the central ECD and irregularities in the central LRR (leucine rich repeat) motifs.

Conclusion

The most critical positions that may alter the pathogen recognition ability of TLR were: the 9^th amino acid position in LRR motif (TLR1–LRR10) and 4^th residue downstream to LRR domain (exta-LRR region of TLR4). The study describes novel mutations in the TLRs and presents their association with the MAP infection.     

A computational investigation of kinetoplastid trans-splicing

Trans-splicing is an unusual process in which two separate RNA strands are spliced together to yield a mature mRNA. We present a novel computational approach which has an overall accuracy of 82% and can predict 92% of known trans-splicing sites. We have applied our method to chromosomes 1 and 3 of Leishmania major, with high-confidence predictions for 85% and 88% of annotated genes respectively. We suggest some extensions of our method to other systems.     

Dieback and mortality of South African fynbos shrubs is likely driven by a novel pathogen and pathogen‐induced hydraulic failure

We examined whether extensive dry season dieback and mortality in a South African fynbos community were due to drought or pathogen attack. Plant dieback and mortality have been reported elsewhere in similar plant communities suggesting potential for a widespread climatic or biotic threat to this community. We collected tissue samples from Brunia noduliflora, the dominant plant in the community, and cultured them for potential plant pathogens. We also measured dry season predawn and midday water potentials of healthy and stressed plants and constructed pressure‐volume curves to assess turgor loss point. Plant stress and mortality were monitored over a 2‐year study period. Both healthy plants and plants that displayed moderate signs of stress had dry season predawn water potentials well above their turgor loss point suggesting plants were not water stressed. However, plants displaying >60% crown dieback had much lower water potentials (as low as ?12 MPa). A previously undescribed fungus (Pythium sp.) was isolated from the root vascular tissue of all stressed plants and was not present in healthy plants. The proximate cause of plant stress was likely pathogen‐induced, while the ultimate cause of plant death appears to be extreme water stress. The present study suggests that Brunia (Bruniceae), Leucadendron (Proteaceae) and Erica (Ericaceae), all emblematic and dominant genera within the diverse fynbos community, may be susceptible to Pythium infection. This may pose a serious threat to communities already threatened by climate change.     

A novel disorder caused by defective biosynthesis of N-linked oligosaccharides due to glucosidase I deficiency

Glucosidase I is an important enzyme in N-linked glycoprotein processing, removing specifically distal alpha-1,2-linked glucose from the Glc3Man9GlcNAc2 precursor after its en bloc transfer from dolichyl diphosphate to a nascent polypeptide chain in the endoplasmic reticulum. We have identified a glucosidase I defect in a neonate with severe generalized hypotonia and dysmorphic features. The clinical course was progressive and was characterized by the occurrence of hepatomegaly, hypoventilation, feeding problems, seizures, and fatal outcome at age 74 d. The accumulation of the tetrasaccharide Glc(alpha1-2)Glc(alpha1-3)Glc(alpha1-3)Man in the patient's urine indicated a glycosylation disorder. Enzymological studies on liver tissue and cultured skin fibroblasts revealed a severe glucosidase I deficiency. The residual activity was <3% of that of controls. Glucosidase I activities in cultured skin fibroblasts from both parents were found to be 50% of those of controls. Tissues from the patient subjected to SDS-PAGE followed by immunoblotting revealed strongly decreased amounts of glucosidase I protein in the homogenate of the liver, and a less-severe decrease in cultured skin fibroblasts. Molecular studies showed that the patient was a compound heterozygote for two missense mutations in the glucosidase I gene: (1) one allele harbored a G-->C transition at nucleotide (nt) 1587, resulting in the substitution of Arg at position 486 by Thr (R486T), and (2) on the other allele a T-->C transition at nt 2085 resulted in the substitution of Phe at position 652 by Leu (F652L). The mother was heterozygous for the G-->C transition, whereas the father was heterozygous for the T-->C transition. These base changes were not seen in 100 control DNA samples. A causal relationship between the alpha-glucosidase I deficiency and the disease is postulated.     

Response to genomic selection: The Bulmer effect and the potential of genomic selection when the number of phenotypic records is limiting

Background

Over the last ten years, genomic selection has developed enormously. Simulations and results on real data suggest that breeding values can be predicted with high accuracy using genetic markers alone. However, to reach high accuracies, large reference populations are needed. In many livestock populations or even species, such populations cannot be established when traits are difficult or expensive to record, or when the population size is small. The value of genomic selection is then questionable.

Methods

In this study, we compare traditional breeding schemes based on own performance or progeny information to genomic selection schemes, for which the number of phenotypic records is limiting. Deterministic simulations were performed using selection index theory. Our focus was on the equilibrium response obtained after a few generations of selection. Therefore, we first investigated the magnitude of the Bulmer effect with genomic selection.

Results

Results showed that the reduction in response due to the Bulmer effect is the same for genomic selection as for selection based on traditional BLUP estimated breeding values, and is independent of the accuracy of selection. The reduction in response with genomic selection is greater than with selection based directly on phenotypes without the use of pedigree information, such as mass selection. To maximize the accuracy of genomic estimated breeding values when the number of phenotypic records is limiting, the same individuals should be phenotyped and genotyped, rather than genotyping parents and phenotyping their progeny. When the generation interval cannot be reduced with genomic selection, large reference populations are required to obtain a similar response to that with selection based on BLUP estimated breeding values based on own performance or progeny information. However, when a genomic selection scheme has a moderate decrease in generation interval, relatively small reference population sizes are needed to obtain a similar response to that with selection on traditional BLUP estimated breeding values.

Conclusions

When the trait of interest cannot be recorded on the selection candidate, genomic selection schemes are very attractive even when the number of phenotypic records is limited, because traditional breeding requires progeny testing schemes with long generation intervals in those cases.     

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.