POOR PERFORMANCE OF BOOTSTRAP CONFIDENCE INTERVALS FOR THE LOCATION OF A QUANTITATIVE TRAIT LOUCS

The aim of many genetic studies is to locate the genomic regions (called quantitative trait loci, QTLs) that contribute to variation in a quantitative trait (such as body weight). Confidence intervals for the locations of QTLs are particularly important for the design of further experiments to identify the gene or genes responsible for the effect. Likelihood support intervals are the most widely used method to obtain confidence intervals for QTL location, but the non-parametric bootstrap has also been recommended. Through extensive computer simulation, we show that bootstrap confidence intervals are poorly behaved and so should not be used in this context. The profile likelihood (or LOD curve) for QTL location has a tendency to peak at genetic markers, and so the distribution of the maximum likelihood estimate (MLE) of QTL location has the unusual feature of point masses at genetic markers; this contributes to the poor behavior of the bootstrap. Likelihood support intervals and approximate Bayes credible intervals, on the other hand, are shown to behave appropriately.

USE OF HIDDEN MARKOV MODELS FOR QTL MAPPING

An important aspect of the QTL mapping problem is the treatment of missing genotype data. If complete genotype data were available, QTL mapping would reduce to the problem of model selection in linear regression. However, in the consideration of loci in the intervals between the available genetic markers, genotype data is inherently missing. Even at the typed genetic markers, genotype data is seldom complete, as a result of failures in the genotyping assays or for the sake of economy (for example, in the case of selective genotyping, where only individuals with extreme phenotypes are genotyped). We discuss the use of algorithms developed for hidden Markov models (HMMs) to deal with the missing genotype data problem.

A polymorphism within the equine CRISP3 gene is associated with stallion fertility in Hanoverian warmblood horses

Fertility of stallions is of high economic importance, especially for large breeding organisations and studs. Breeding schemes with respect to fertility traits and selection of stallions at an early stage may be improved by including molecular genetic markers associated with traits. The genes coding for equine cysteine-rich secretory proteins (CRISPs) are promising candidate genes because previous studies have shown that CRISPs play a role in the fertilising ability of male animals. We have previously characterised the three equine CRISP genes and identified a non-synonymous polymorphism in the CRISP1 gene. In this study, we report one non-synonymous polymorphism in the CRISP2 gene and four non-synonymous polymorphisms in the CRISP3 gene. All six CRISP polymorphisms were genotyped in 107 Hanoverian breeding stallions. Insemination records of stallions were used to analyse the association between CRISP polymorphisms and fertility traits. Three statistical models were used to evaluate the influence of single mutations, genotypes and haplotypes of the polymorphisms. The CRISP3 AJ459965:c.+622G>A SNP leading to the amino acid substitution E208K was significantly associated with the fertility of stallions. Stallions heterozygous for the CRISP3 c.+622G>A SNP had lower fertility than homozygous stallions (P = 0.0234). The pregnancy rate per cycle in these stallions was estimated to be approximately 7% lower than in stallions homozygous at this position.

DEVELOPMENT AND APPLICATIONS OF GENOMIC RESOURCES FOR TWO NORTH AMERICAN HARDWOOD TAXA

Hardwoods comprise about half of the biomass of forestlands in North America and present many uses including economic, ecological and aesthetic functions. Forest trees rely on the genetic variation within tree populations to overcome the many biotic, abiotic, anthropogenic factors which are further worsened by climate change, that threaten their continued survival and functionality. To harness these inherent genetic variations of tree populations, informed knowledge of the genomic resources and techniques, which are currently lacking or very limited, are imperative for forest managers. The current study therefore aimed to develop genomic microsatellite markers for the leguminous tree species, honey locust, Gleditsia triacanthos L. and test their applicability in assessing genetic variation, estimation of gene flow patterns and identification of a full-sib mapping population. We also aimed to test the usefulness of already developed nuclear and gene-based microsatellite markers in delineation of species and taxonomic relationships between four of the taxonomically difficult Section Lobatae species (Quercus coccinea, Q. ellipsoidalis, Q. rubra and Q. velutina. We recorded 100% amplification of G. triacanthos genomic microsatellites developed using Illumina sequencing techniques in a panel of seven unrelated individuals with 14 of these showing high polymorphism and reproducibility. When characterized in 36 natural population samples, we recorded 20 alleles per locus with no indication for null alleles at 13 of the 14 microsatellites. This is the first report of genomic microsatellites for this species. Honey locust trees occur in fragmented populations of abandoned farmlands and pastures and is described as essentially dioecious. Pollen dispersal if the main source of gene flow within and between populations with the ability to offset the effects of random genetic drift. Factors known to influence gene include fragmentation and degree of isolation, which make the patterns gene flow in fragmented populations of honey locust a necessity for their sustainable management. In this follow-up study, we used a subset of nine of the 14 developed gSSRs to estimate gene flow and identify a full-sib mapping population in two isolated fragments of honey locust. Our analyses indicated that the majority of the seedlings (65-100% - at both strict and relaxed assignment thresholds) were sired by pollen from outside the two fragment populations. Only one selfing event was recorded confirming the functional dioeciousness of honey locust and that the seed parents are almost completely outcrossed. From the Butternut Valley, TN population, pollen donor genotypes were reconstructed and used in paternity assignment analyses to identify a relatively large full-sib family comprised of 149 individuals, proving the usefulness of isolated forest fragments in identification of full-sib families. In the Ames Plantation stand, contemporary pollen dispersal followed a fat-tailed exponential-power distribution, an indication of effective gene flow. Our estimate of δ was 4,282.28 m, suggesting that insect pollinators of honey locust disperse pollen over very long distances. The high proportion of pollen influx into our sampled population implies that our fragment population forms part of a large effectively reproducing population. The high tendency of oak species to hybridize while still maintaining their species identity make it difficult to resolve their taxonomic relationships. Oaks of the section Lobatae are famous in this regard and remain unresolved at both morphological and genetic markers. We applied 28 microsatellite markers including outlier loci with potential roles in reproductive isolation and adaptive divergence between species to natural populations of four known interfertile red oaks, Q. coccinea, Q. ellpsoidalis, Q. rubra and Q. velutina. To better resolve the taxonomic relationships in this difficult clade, we assigned individual samples to species, identified hybrids and introgressive forms and reconstructed phylogenetic relationships among the four species after exclusion of genetically intermediate individuals. Genetic assignment analyses identified four distinct species clusters, with Q. rubra most differentiated from the three other species, but also with a comparatively large number of misclassified individuals (7.14%), hybrids (7.14%) and introgressive forms (18.83%) between Q. ellipsoidalis and Q. velutina. After the exclusion of genetically intermediate individuals, Q. ellipsoidalis grouped as sister species to the largely parapatric Q. coccinea with high bootstrap support (91 %). Genetically intermediate forms in a mixed species stand were located proximate to both potential parental species, which supports recent hybridization of Q. velutina with both Q. ellipsoidalis and Q. rubra. Analyses of genome-wide patterns of interspecific differentiation can provide a better understanding of speciation processes and taxonomic relationships in this taxonomically difficult group of red oak species.

Selection on floral traits through male fertility in a natural plant population

Most studies on selection in plants estimate female fitness components and neglect male mating success, although the latter might also be fundamental to understand adaptive evolution. Information from molecular genetic markers can be used to assess determinants of male mating success through parentage analyses. We estimated paternal selection gradients on floral traits in a large natural population of the herb Mimulus guttatus using a paternity probability model and maximum likelihood methods. This analysis revealed more significant selection gradients than a previous analysis based on regression of estimated male fertilities on floral traits. There were differences between results of univariate and multivariate analyses most likely due to the underlying covariance structure of the traits. Multivariate analysis, which corrects for the covariance structure of the traits, indicated that male mating success declined with distance from and depended on the direction to the mother plants. Moreover, there was directional selection for plants with fewer open flowers which have smaller corollas, a smaller anther-stigma separation, more red dots on the corolla and a larger fluctuating asymmetry therein. For most of these traits, however, there was also stabilizing selection indicating that there are intermediate optima for these traits. The large number of significant selection gradients in this study shows that even in relatively large natural populations where not all males can be sampled, it is possible to detect significant paternal selection gradients, and that such studies can give us valuable information required to better understand adaptive plant evolution.

Classification and pharmacological treatment of preschool wheezing: changes since 2008.

Since the publication of the European Respiratory Society Task Force report in 2008, significant new evidence has become available on the classification and management of preschool wheezing disorders. In this report, an international consensus group reviews this new evidence and proposes some modifications to the recommendations made in 2008. Specifically, the consensus group acknowledges that wheeze patterns in young children vary over time and with treatment, rendering the distinction between episodic viral wheeze and multiple-trigger wheeze unclear in many patients. Inhaled corticosteroids remain first-line treatment for multiple-trigger wheeze, but may also be considered in patients with episodic viral wheeze with frequent or severe episodes, or when the clinician suspects that interval symptoms are being under reported. Any controller therapy should be viewed as a treatment trial, with scheduled close follow-up to monitor treatment effect. The group recommends discontinuing treatment if there is no benefit and taking favourable natural history into account when making decisions about long-term therapy. Oral corticosteroids are not indicated in mild-to-moderate acute wheeze episodes and should be reserved for severe exacerbations in hospitalised patients. Future research should focus on better clinical and genetic markers, as well as biomarkers, of disease severity.

Major histocompatibility complex and other allergy-related candidate genes associated with insect bite hypersensitivity in Icelandic horses

Insect bite hypersensitivity (IBH) is an allergic dermatitis of horses caused by bites of insects. IBH is a multifactorial disease with contribution of genetic and environmental factors. Candidate gene association analysis of IBH was performed in a group of 89 Icelandic horses all born in Iceland and imported to Europe. Horses were classified in IBH-affected and non-affected based on clinical signs and history of recurrent dermatitis, and on the results of an in vitro sulfidoleukotriene (sLT)-release assay with Culicoides nubeculosus and Simulium vittatum extract. Different genetic markers were tested for association with IBH by the Fisher's exact test. The effect of the major histocompatibility complex (MHC) gene region was studied by genotyping five microsatellites spanning the MHC region (COR112, COR113, COR114, UM011 and UMN-JH34-2), and exon 2 polymorphisms of the class II Eqca-DRA gene. Associations with Eqca-DRA and COR113 were identified (p < 0.05). In addition, a panel of 20 single nucleotide polymorphisms (SNPs) in 17 candidate allergy-related genes was tested. During the initial screen, no marker from the panel was significantly (p < 0.05) associated with IBH. Five SNPs associated with IBH at p < 0.10 were therefore used for analysis of combined genotypes. Out of them, SNPs located in the genes coding for the CD14 receptor (CD14), interleukin 23 receptor (IL23R), thymic stromal lymphopoietin (TSLP) and transforming growth factor beta 3 (TGFB3) molecules were associated with IBH as parts of complex genotypes. These results are supported by similar associations and by expression data from different horse populations and from human studies.

Sex-specific clines support incipient speciation in a common European mammal

Hybrid zones provide excellent opportunities to study processes and mechanisms underlying reproductive isolation and speciation. Here we investigated sex-specific clines of molecular markers in hybrid zones of morphologically cryptic yet genetically highly-diverged evolutionary lineages of the European common vole (Microtus arvalis). We analyzed the position and width of four secondary contact zones along three independent transects in the region of the Alps using maternally (mitochondrial DNA) and paternally (Y-chromosome) inherited genetic markers. Given male-biased dispersal in the common vole, a selectively neutral secondary contact would show broader paternal marker clines than maternal ones. In a selective case, for example, involving a form of Haldane’s rule, Y-chromosomal clines would not be expected to be broader than maternal markers because they are transmitted by the heterogametic sex and thus gene flow would be restricted. Consistent with the selective case, paternal clines were significantly narrower or at most equal in width to maternal clines in all contact zones. In addition, analyses using maximum likelihood cline-fitting detected a shift of paternal relative to maternal clines in three of four contact zones. These patterns suggest that processes at the contact zones in the common vole are not selectively neutral, and that partial reproductive isolation is already established between these evolutionary lineages. We conclude that hybrid zone movement, sexual selection and/or genetic incompatibilities are likely associated with an unusual unidirectional manifestation of Haldane’s rule in this common European mammal.

Population structure, inbreeding and local adaptation within an endangered riverine specialist: the nase (Chondrostoma nasus)

Chondrostoma nasus is a cyprinid fish with highly specialized, ecologically and geographically distinct, ontogenetic trophic niches. Nase population numbers across their Swiss range have shown massive declines and many localized extinctions. Here we integrate data from different genetic markers with phenotypic and demographic data to survey patterns of neutral and adaptive genetic diversity in all extant (and one extinct) Swiss nase populations, with the aim to delineate intraspecific conservation units (CUs) and to inform future population management strategies. We discovered two major genetically and geographically distinct population groupings. The first population grouping comprises nase inhabiting rivers flowing into Lake Constance; the second comprises nase populations from Rhine drainages below Lake Constance. Within these clusters there is generally limited genetic differentiation among populations. Genomic outlier scans based on 256–377 polymorphic AFLP loci revealed little evidence of local adaptation both within and among population clusters, with the exception of one candidate locus identified in scans involving the inbred Schanzengraben population. However, significant phenotypic differentiation in body shape between certain populations suggests a need for more intensive future studies of local adaptation. Our data strongly suggests that the two major population groups should be treated as distinct CUs, with any supplemental stocking and reintroductions sourced only from within the range of the CU concerned.

Genetics of upper and lower airway diseases in the horse.

Genetic predispositions for guttural pouch tympany, recurrent laryngeal neuropathy and recurrent airway obstruction (RAO) are well documented. There is also evidence that exercise-induced pulmonary haemorrhage and infectious diseases of the respiratory tract in horses have a genetic component. The clinical expression of equine respiratory diseases with a genetic basis results from complex interactions between the environment and the genetic make-up of each individual horse. The genetic effects are likely to be due to variations in several genes, i.e. they are polygenic. It is therefore unlikely that single gene tests will be diagnostically useful in these disorders. Genetic profiling panels, combining several genetic factors with an assessment of environmental risk factors, may have greater value, but much work is still needed to uncover diagnostically useful genetic markers or even causative variants for equine respiratory diseases. Nonetheless, chromosomal regions associated with guttural pouch tympany, recurrent laryngeal neuropathy and RAO have been identified. The association of RAO with other hypersensitivities and with resistance to intestinal parasites requires further study. This review aims to provide an overview of the available data and current thoughts on the genetics of equine airway diseases.

Evaluation of genome-wide loci of iron metabolism in hereditary hemochromatosis identifies PCSK7 as a host risk factor of liver cirrhosis

Genome-wide association studies (GWAS) have revealed genetic determinants of iron metabolism, but correlation of these with clinical phenotypes is pending. Homozygosity for HFE C282Y is the predominant genetic risk factor for hereditary hemochromatosis (HH) and may cause liver cirrhosis. However, this genotype has a low penetrance. Thus, detection of yet unknown genetic markers that identify patients at risk of developing severe liver disease is necessary for better prevention. Genetic loci associated with iron metabolism (TF, TMPRSS6, PCSK7, TFR2 and Chr2p14) in recent GWAS and liver fibrosis (PNPLA3) in recent meta-analysis were analyzed for association with either liver cirrhosis or advanced fibrosis in 148 German HFE C282Y homozygotes. Replication of associations was sought in additional 499 Austrian/Swiss and 112 HFE C282Y homozygotes from Sweden. Only variant rs236918 in the PCSK7 gene (proprotein convertase subtilisin/kexin type 7) was associated with cirrhosis or advanced fibrosis (P = 1.02 × 10(-5)) in the German cohort with genotypic odds ratios of 3.56 (95% CI 1.29-9.77) for CG heterozygotes and 5.38 (95% CI 2.39-12.10) for C allele carriers. Association between rs236918 and cirrhosis was confirmed in Austrian/Swiss HFE C282Y homozygotes (P = 0.014; ORallelic = 1.82 (95% CI 1.12-2.95) but not in Swedish patients. Post hoc combined analyses of German/Swiss/Austrian patients with available liver histology (N = 244, P = 0.00014, ORallelic = 2.84) and of males only (N = 431, P = 2.17 × 10(-5), ORallelic = 2.54) were consistent with the premier finding. Association between rs236918 and cirrhosis was not confirmed in alcoholic cirrhotics, suggesting specificity of this genetic risk factor for HH. PCSK7 variant rs236918 is a risk factor for cirrhosis in HH patients homozygous for the HFE C282Y mutation.

Taking personalized medicine seriously: Biomarker approaches in phase IIb/III studies in major depression and schizophrenia

The success rate in the development of psychopharmacological compounds is insufficient. Two main reasons for failure have been frequently identified: 1) treating the wrong patients and 2) using the wrong dose. This is potentially based on the known heterogeneity among patients, both on a syndromal and a biological level. A focus on personalized medicine through better characterization with biomarkers has been successful in other therapeutic areas. Nevertheless, obstacles toward this goal that exist are 1) the perception of a lack of validation, 2) the perception of an expensive and complicated enterprise, and 3) the perception of regulatory hurdles. The authors tackle these concerns and focus on the utilization of biomarkers as predictive markers for treatment outcome. The authors primarily cover examples from the areas of major depression and schizophrenia. Methodologies covered include salivary and plasma collection of neuroendocrine, metabolic, and inflammatory markers, which identified subgroups of patients in the Netherlands Study of Depression and Anxiety. A battery of vegetative markers, including sleep-electroencephalography parameters, heart rate variability, and bedside functional tests, can be utilized to characterize the activity of a functional system that is related to treatment refractoriness in depression (e.g., the renin-angiotensin-aldosterone system). Actigraphy and skin conductance can be utilized to classify patients with schizophrenia and provide objective readouts for vegetative activation as a functional marker of target engagement. Genetic markers, related to folate metabolism, or folate itself, has prognostic value for the treatment response in patients with schizophrenia. Already, several biomarkers are routinely collected in standard clinical trials (e.g., blood pressure and plasma electrolytes), and appear to be differentiating factors for treatment outcome. Given the availability of a wide variety of markers, the further development and integration of such markers into clinical research is both required and feasible in order to meet the benefit of personalized medicine. This article is based on proceedings from the "Taking Personalized Medicine Seriously-Biomarker Approaches in Phase IIb/III Studies in Major Depression and Schizophrenia" session, which was held during the 10th Annual Scientific Meeting of the International Society for Clinical Trials Meeting (ISCTM) in Washington, DC, February 18 to 20, 2014.

Genetics of recurrent airway obstruction (RAO)

Recurrent airway obstruction (RAO) is a multifactorial and polygenic disease. Affected horses are typically 7 years of age or older and show exercise intolerance, increased breathing effort, coughing, airway neutrophilia, mucus accumulation and hyperreactivity as well as cholinergic bronchospasm. The environmental factors responsible are predominantly allergens and irritants in haydust, but the immunological mechanisms underlying RAO are still unclear. Several studies have demonstrated a familiar predisposition for RAO and it is now proven that the disease has a genetic basis. In offspring, the risk of developing RAO is 3-fold increased when one parent is affected and increases to almost 5-fold when both parents have RAO. Segregation analysis in two high-prevalence families demonstrated a high heritability and a complex inheritance with several major genes. A whole genomescan showed chromosome-wide significant linkage of seven chromosomal regions with RAO. Of the microsatellites, which were located near atopy candidate genes, those in a region of chromosome 13 harboring the IL4R gene were strongly associated with the RAO phenotype in the offspring of one RAO-affected stallion. Furthermore, IgE-levels are influenced by hereditary factors in the horse, and we have evidence that RAO-affected offspring of the same stallion have increased levels of specific IgE against moldspore allergens. The identification of genetic markers and ultimately of the responsible genes will not only allow for an improved prophylaxis, i.e. early identification of susceptible individuals and avoidance of high-risk matings, but also improve our ability to find new therapeutic targets and to optimize existing treatments.

Privacy-preserving genomic testing in the clinic: a model using HIV treatment

PURPOSE The implementation of genomic-based medicine is hindered by unresolved questions regarding data privacy and delivery of interpreted results to health-care practitioners. We used DNA-based prediction of HIV-related outcomes as a model to explore critical issues in clinical genomics. METHODS We genotyped 4,149 markers in HIV-positive individuals. Variants allowed for prediction of 17 traits relevant to HIV medical care, inference of patient ancestry, and imputation of human leukocyte antigen (HLA) types. Genetic data were processed under a privacy-preserving framework using homomorphic encryption, and clinical reports describing potentially actionable results were delivered to health-care providers. RESULTS A total of 230 patients were included in the study. We demonstrated the feasibility of encrypting a large number of genetic markers, inferring patient ancestry, computing monogenic and polygenic trait risks, and reporting results under privacy-preserving conditions. The average execution time of a multimarker test on encrypted data was 865 ms on a standard computer. The proportion of tests returning potentially actionable genetic results ranged from 0 to 54%. CONCLUSIONS The model of implementation presented herein informs on strategies to deliver genomic test results for clinical care. Data encryption to ensure privacy helps to build patient trust, a key requirement on the road to genomic-based medicine.Genet Med advance online publication 14 January 2016Genetics in Medicine (2016); doi:10.1038/gim.2015.167.

Taxonomy, phylogeny and molecular epidemiology of Echinococcus multilocularis: From fundamental knowledge to health ecology.

Alveolar echinococcosis, caused by the tapeworm Echinococcus multilocularis, is one of the most severe parasitic diseases in humans and represents one of the 17 neglected diseases prioritised by the World Health Organisation (WHO) in 2012. Considering the major medical and veterinary importance of this parasite, the phylogeny of the genus Echinococcus is of considerable importance; yet, despite numerous efforts with both mitochondrial and nuclear data, it has remained unresolved. The genus is clearly complex, and this is one of the reasons for the incomplete understanding of its taxonomy. Although taxonomic studies have recognised E. multilocularis as a separate entity from the Echinococcus granulosus complex and other members of the genus, it would be premature to draw firm conclusions about the taxonomy of the genus before the phylogeny of the whole genus is fully resolved. The recent sequencing of E. multilocularis and E. granulosus genomes opens new possibilities for performing in-depth phylogenetic analyses. In addition, whole genome data provide the possibility of inferring phylogenies based on a large number of functional genes, i.e. genes that trace the evolutionary history of adaptation in E. multilocularis and other members of the genus. Moreover, genomic data open new avenues for studying the molecular epidemiology of E. multilocularis: genotyping studies with larger panels of genetic markers allow the genetic diversity and spatial dynamics of parasites to be evaluated with greater precision. There is an urgent need for international coordination of genotyping of E. multilocularis isolates from animals and human patients. This could be fundamental for a better understanding of the transmission of alveolar echinococcosis and for designing efficient healthcare strategies.