Above and beyond state-of-the-art approaches to investigate sequence data: summary of methods and results from the population-based association group at the Genetic Analysis Workshop 19

This paper summarizes the contributions from the Population-Based Association group at the Genetic Analysis Workshop 19. It provides an overview of the new statistical approaches tried out by group members in order to take best advantage of population-based sequence data.Although contributions were highly heterogeneous regarding the applied quality control criteria and the number of investigated variants, several technical issues were identified, leading to practical recommendations. Preliminary analyses revealed that Hurdle-negative binomial regression is a promising approach to investigate the distribution of allele counts instead of called genotypes from sequence data. Convergence problems, however, limited the use of this approach, creating a technical challenge shared by environment-stratified models used to investigate rare variant-environment interactions, as well as by rare variant haplotype analyses using well-established public software. Estimates of relatedness and population structure strongly depended on the allele frequency of selected variants for inference. Another practical recommendation was that dissenting probability values from standard and small-sample tests of a particular hypothesis may reflect a lack of validity of large-sample approximations. Novel statistical approaches that integrate evolutionary information showed some advantage to detect weak genetic signals, and Bayesian adjustment for confounding was able to efficiently estimate causal genetic effects. Haplotype association methods may constitute a valuable complement of collapsing approaches for sequence data. This paper reports on the experience of members of the Population-Based Association group with several novel, promising approaches to preprocessing and analyzing sequence data, and to following up identified association signals.     

Longitudinal analytical approaches to genetic data

Background

Longitudinal phenotypic data provides a rich potential resource for genetic studies which may allow for greater understanding of variants and their covariates over time. Herein, we review 3 longitudinal analytical approaches from the Genetic Analysis Workshop 19 (GAW19). These contributions investigated both genome-wide association (GWA) and whole genome sequence (WGS) data from odd numbered chromosomes on up to 4 time points for blood pressure–related phenotypes. The statistical models used included generalized estimating equations (GEEs), latent class growth modeling (LCGM), linear mixed-effect (LME), and variance components (VC). The goal of these analyses was to test statistical approaches that use repeat measurements to increase genetic signal for variant identification.

Results

Two analytical methods were applied to the GAW19: GWA using real phenotypic data, and one approach to WGS using 200 simulated replicates. The first GWA approach applied a GEE-based model to identify gene-based associations with 4 derived hypertension phenotypes. This GEE model identified 1 significant locus, GRM7, which passed multiple test corrections for 2 hypertension-derived traits. The second GWA approach employed the LME to estimate genetic associations with systolic blood pressure (SBP) change trajectories identified using LCGM. This LCGM method identified 5 SBP trajectories and association analyses identified a genome-wide significant locus, near ATOX1 (p?=?1.0E^?8). Finally, a third VC-based model using WGS and simulated SBP phenotypes that constrained the β coefficient for a genetic variant across each time point was calculated and compared to an unconstrained approach. This constrained VC approach demonstrated increased power for WGS variants of moderate effect, but when larger genetic effects were present, averaging across time points was as effective.

Conclusion

In this paper, we summarize 3 GAW19 contributions applying novel statistical methods and testing previously proposed techniques under alternative conditions for longitudinal genetic association. We conclude that these approaches when appropriately applied have the potential to: (a) increase statistical power; (b) decrease trait heterogeneity and standard error; (c) decrease computational burden in WGS; and (d) have the potential to identify genetic variants influencing subphenotypes important for understanding disease progression.

     

Predicting existing targets for new drugs base on strategies for missing interactions

Background

There has been paid more and more attention to supervised classification models in the area of predicting drug-target interactions (DTIs). However, in terms of classification, unavoidable missing DTIs in data would cause three issues which have not yet been addressed appropriately by former approaches. Directly labeled as negatives (non-DTIs), missing DTIs increase the confusion of positives (DTIs) and negatives, aggravate the imbalance between few positives and many negatives, and are usually discriminated as highly-scored false positives, which influence the existing measures sharply.

Results

Under the framework of local classification model (LCM), this work focuses on the scenario of predicting how possibly a new drug interacts with known targets. To address the first two issues, two strategies, Spy and Super-target, are introduced accordingly and further integrated to form a two-layer LCM. In the bottom layer, Spy-based local classifiers for protein targets are built by positives, as well as reliable negatives identified among unlabeled drug-target pairs. In the top layer, regular local classifiers specific to super-targets are built with more positives generated by grouping similar targets and their interactions. Furthermore, to handle the third issue, an additional performance measure, Coverage, is presented for assessing DTI prediction. The experiments based on benchmark datasets are finally performed under five-fold cross validation of drugs to evaluate this approach. The main findings are concluded as follows. (1) Both two individual strategies and their combination are effective to missing DTIs, and the combination wins the best. (2) Having the advantages of less confusing decision boundary at the bottom layer and less biased decision boundary at the top layer, our two-layer LCM outperforms two former approaches. (3) Coverage is more robust to missing interactions than other measures and is able to evaluate how far one needs to go down the list of targets to cover all the proper targets of a drug.

Conclusions

Proposing two strategies and one performance measure, this work has addressed the issues derived from missing interactions, which cause confusing and biased decision boundaries in classifiers, as well as the inappropriate measure of predicting performance, in the scenario of predicting interactions between new drugs and known targets.

     

Expression of an extracellular ribonuclease gene increases resistance to <Emphasis Type="Italic">Cucumber mosaic virus</Emphasis> in tobacco

Background

The apoplast plays an important role in plant defense against pathogens. Some extracellular PR-4 proteins possess ribonuclease activity and may directly inhibit the growth of pathogenic fungi. It is likely that extracellular RNases can also protect plants against some viruses with RNA genomes. However, many plant RNases are multifunctional and the direct link between their ribonucleolytic activity and antiviral defense still needs to be clarified. In this study, we evaluated the resistance of Nicotiana tabacum plants expressing a non-plant single-strand-specific extracellular RNase against Cucumber mosaic virus.

Results

Severe mosaic symptoms and shrinkage were observed in the control non-transgenic plants 10 days after inoculation with Cucumber mosaic virus (CMV), whereas such disease symptoms were suppressed in the transgenic plants expressing the RNase gene. In a Western blot analysis, viral proliferation was observed in the uninoculated upper leaves of control plants, whereas virus levels were very low in those of transgenic plants. These results suggest that resistance against CMV was increased by the expression of the heterologous RNase gene.

Conclusion

We have previously shown that tobacco plants expressing heterologous RNases are characterized by high resistance to Tobacco mosaic virus. In this study, we demonstrated that elevated levels of extracellular RNase activity resulted in increased resistance to a virus with a different genome organization and life cycle. Thus, we conclude that the pathogen-induced expression of plant apoplastic RNases may increase non-specific resistance against viruses with RNA genomes.

     

Comparison of human and mouse E-selectin binding to Sialyl-Lewis<Superscript>x</Superscript>

Background

During inflammation, leukocytes are captured by the selectin family of adhesion receptors lining blood vessels to facilitate exit from the bloodstream. E-selectin is upregulated on stimulated endothelial cells and binds to several ligands on the surface of leukocytes. Selectin:ligand interactions are mediated in part by the interaction between the lectin domain and Sialyl-Lewis x (sLe^x), a tetrasaccharide common to selectin ligands. There is a high degree of homology between selectins of various species: about 72 and 60 % in the lectin and EGF domains, respectively. In this study, molecular dynamics, docking, and steered molecular dynamics simulations were used to compare the binding and dissociation mechanisms of sLe^x with mouse and human E-selectin. First, a mouse E-selectin homology model was generated using the human E-selectin crystal structure as a template.

Results

Mouse E-selectin was found to have a greater interdomain angle, which has been previously shown to correlate with stronger binding among selectins. sLe^x was docked onto human and mouse E-selectin, and the mouse complex was found to have a higher free energy of binding and a lower dissociation constant, suggesting stronger binding. The mouse complex had higher flexibility in a few key residues. Finally, steered molecular dynamics was used to dissociate the complexes at force loading rates of 2000–5000 pm/ps². The mouse complex took longer to dissociate at every force loading rate and the difference was statistically significant at 3000 pm/ps². When sLe^x-coated microspheres were perfused through microtubes coated with human or mouse E-selectin, the particles rolled more slowly on mouse E-selectin.

Conclusions

Both molecular dynamics simulations and microsphere adhesion experiments show that mouse E-selectin protein binds more strongly to sialyl Lewis x ligand than human E-selectin. This difference was explained by a greater interdomain angle for mouse E-selectin, and greater flexibility in key residues. Future work could introduce similar amino acid substitutions into the human E-selectin sequence to further modulate adhesion behavior.

     

Case report of whole genome sequencing in the XY female: identification of a novel <Emphasis Type="Italic">SRY</Emphasis> mutation and revision of a misdiagnosis of androgen insensitivity syndrome

Background

The 46,XY female is characterised by a male karyotype and female phenotype arising due to any interruption in the sexual development pathways in utero. The cause is usually genetic and various genes are implicated.

Case presentation

Herein we describe a 46,XY woman who was first diagnosed with androgen insensitivity syndrome (testicular feminisation) at 18 years; however, this was later questioned due to the presence of intact Müllerian structures. The clinical phenotype suggested several susceptibility genes including SRY, DHH, NR5A1, NR0B1, AR, AMH, and AMHR2. To study candidate genes simultaneously, we performed whole genome sequencing. This revealed a novel and likely pathogenic missense variant (p.Arg130Pro, c.389G>C) in SRY, one of the major genes implicated in complete gonadal dysgenesis, hence securing this condition over androgen insensitivity syndrome as the cause of the patient’s disorder of sexual development.

Conclusion

This case highlights the emerging clinical utility of whole genome sequencing as a tool in differentiating disorders of sexual development.