Dimensionality and modularity of adaptive variation: Divergence in threespine stickleback from diverse environments

In nature, populations are subjected to a wide variety of environmental conditions that affect fitness and induce adaptive or plastic responses in traits, resulting in phenotypic divergence between populations. The dimensionality of that divergence, however, remains contentious. At the extremes, some contend that populations diverge along a single axis of trait covariance with greatest availability of heritable variation, even if this does not lead a population directly to its fitness optimum. Those at the other extreme argue that selection can push populations towards their fitness optima along multiple phenotype axes simultaneously, resulting in divergence in numerous dimensions. Here, we address this debate using populations of threespine stickleback (Gasterosteus aculeatus) in the Cook Inlet region of southern Alaska from lakes with contrasting ecological conditions. We calculated effective dimensionality of divergence in several trait suites (defensive, swimming, and trophic) thought to be under correlated selection pressures, as well as across all traits. We also tested for integration among the trait suites and between each trait suite and the environment. We found that populations in the Cook Inlet radiation exhibit dimensionality of phenotype high enough to preclude a single axis of divergence.

The Human Phenotype Ontology in 2021

The Human Phenotype Ontology (HPO, https://hpo.jax.org) was launched in 2008 to provide a comprehensive logical standard to describe and computationally analyze phenotypic abnormalities found in human disease. The HPO is now a worldwide standard for phenotype exchange. The HPO has grown steadily since its inception due to considerable contributions from clinical experts and researchers from a diverse range of disciplines. Here, we present recent major extensions of the HPO for neurology, nephrology, immunology, pulmonology, newborn screening, and other areas. For example, the seizure subontology now reflects the International League Against Epilepsy (ILAE) guidelines and these enhancements have already shown clinical validity. We present new efforts to harmonize computational definitions of phenotypic abnormalities across the HPO and multiple phenotype ontologies used for animal models of disease. These efforts will benefit software such as Exomiser by improving the accuracy and scope of cross-species phenotype matching. The computational modeling strategy used by the HPO to define disease entities and phenotypic features and distinguish between them is explained in detail.We also report on recent efforts to translate the HPO into indigenous languages. Finally, we summarize recent advances in the use of HPO in electronic health record systems.

A multiple-phenotype imputation procedure as a method for prediction of cheese-making efficiency in Spanish Assaf sheep

Sheep milk is mainly intended to manufacture a wide variety of high-quality cheeses. The ovine cheese industry would benefit from an improvement, through genetic selection, of traits related to the milk coagulation properties (MCPs) and cheese yield-related traits, broadly denoted as “cheese-making traits.” Considering that routine measurements of these traits needed for genetic selection are expensive and time-consuming, this study aimed to evaluate the accuracy of a cheese-making phenotype imputation method based on the information from official milk control records combined with the pH of the milk. For this study, we analyzed records of milk production traits, milk composition traits, and measurements of cheese-making traits available from a total of 1,145 dairy ewes of the Spanish Assaf sheep breed. Cheese-making traits included five related to the MCPs and two cheese yield-related traits. The milk and cheese-making phenotypes were adjusted for significant effects based on a general linear model. The adjusted phenotypes were used to define a multiple-phenotype imputation procedure for the cheese-making traits based on multivariate normality and Markov chain Monte Carlo sampling. Five of the seven cheese-making traits considered in this study achieved a prediction accuracy of 0.60 computed as the correlation between the adjusted phenotypes and the imputed phenotypes. Particularly the logarithm of curd-firming time since rennet addition (logK20) (0.68), which has been previously suggested as a potential candidate trait to improve the cheese ability in this breed, and the logarithm of the ratio between the rennet clotting time and the curd firmness at 60 min (logRCT/A60) (0.65), which has been defined by other studies as an indicator trait of milk coagulation efficiency. This study represents a first step toward the possible use of the phenotype imputation of cheese-making traits to develop a practical methodology for the dairy sheep industry to impute cheese-making traits only based on the analysis of a milk sample without the need of pedigree information. This information could be also used in future planning of specific breeding programs considering the importance of the cheese-making efficiency in dairy sheep and highlights the potential of phenotype imputation to leverage sample size on expensive, hard-to-measure phenotypes.

A phenome-wide association study of four syndromic genes reveals pleiotropic effects of common and rare variants in the general population

The clinical evaluation of a genetic syndrome relies upon recognition of a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. However, individuals displaying a few phenotypes of a complex syndrome may not meet criteria for clinical diagnosis or genetic testing. Here, we present a phenome-wide association study (PheWAS) approach to systematically explore pleiotropy of common and rare alleles in genes associated with four well-described syndromic diseases (Alagille (AS), Marfan (MS), DiGeorge (DS), and Noonan (NS) syndromes) in the general population.Using human phenotype ontology (HPO) terms, we systematically mapped 60 phenotypes related to AS, MS, DS and NS in 337,198 unrelated white British from the UK Biobank (UKBB) based on their hospital admission records, self-administrated questionnaires, and physiological measurements. We performed logistic regression adjusting for age, sex, and the first 5 genetic principal components, for each phenotype and each variant in the target genes (JAG1, TBX1, FBN1, PTPN11, NOTCH2, and MAP2K1) and performed a gene burden testing.Overall, we observed multiple phenotype-genotype correlations, such as the association between variation in JAG1, FBN1, PTPN11 and SOS2 with diastolic and systolic blood pressure; and pleiotropy among multiple variants in syndromic genes. For example, rs11066309 in PTPN11 was significantly associated with a lower body mass index, an increased risk of hypothyroidism and a smaller size for gestational age, all in concordance with NS-related phenotypes. Similarly, rs589668 in FBN1 was associated with an increase in body height and blood pressure, and a reduced body fat percentage as observed in Marfan syndrome.Our findings suggest that the spectrum of associations of common and rare variants in genes involved in syndromic diseases can be extended to individual phenotypes within the general population.Author SummaryStandard medical evaluation of genetic syndromes relies upon recognizing a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. This may lead to missing diagnoses in patients with silent or a low expressed form of the syndrome. Here we take advantage of a rich electronic health record, various phenotypic measurements, and genetic information in 337,198 unrelated white British from the UKBB, to study the relation between single syndromic disease phenotypes and genes related to syndromic disease. We show multiple phenotype-genotypes associations in concordance with phenotypes variations found in syndromic diseases. For example, we show that mutation in FBN1 was associated with high standing/sitting height ratio and reduced body fat percentage as observed in individuals with Marfan syndrome. Our findings suggest that common and rare alleles in SD genes are causative of individual component phenotypes present in a general population; further research is needed to characterize the pleiotropic effect of alleles in syndromic genes in persons without the syndromic disease.

Meta-MultiSKAT: Multiple phenotype meta-analysis for region-based association test

The power of genetic association analyses can be increased by jointly meta-analyzing multiple correlated phenotypes. Here, we develop a meta-analysis framework, Meta-MultiSKAT, that uses summary statistics to test for association between multiple continuous phenotypes and variants in a region of interest. Our approach models the heterogeneity of effects between studies through a kernel matrix and performs a variance component test for association. Using a genotype kernel, our approach can test for rare-variants and the combined effects of both common and rare-variants. To achieve robust power, within Meta-MultiSKAT, we developed fast and accurate omnibus tests combining different models of genetic effects, functional genomic annotations, multiple correlated phenotypes and heterogeneity across studies. Additionally, Meta-MultiSKAT accommodates situations where studies do not share exactly the same set of phenotypes or have differing correlation patterns among the phenotypes. Simulation studies confirm that Meta-MultiSKAT can maintain type-I error rate at exome-wide level of 2.5×10−6. Further simulations under different models of association show that Meta-MultiSKAT can improve power of detection from 23% to 38% on average over single phenotype-based meta-analysis approaches. We demonstrate the utility and improved power of Meta-MultiSKAT in the meta-analyses of four white blood cell subtype traits from the Michigan Genomics Initiative (MGI) and SardiNIA studies.

Psychopathologic structure of bipolar disorders: exploring dimensional phenotypes, their relationships, and their associations with bipolar I and II disorders

BackgroundGiven its diverse disease courses and symptom presentations, multiple phenotype dimensions with different biological underpinnings are expected with bipolar disorders (BPs). In this study, we aimed to identify lifetime BP psychopathology dimensions. We also explored the differing associations with bipolar I (BP-I) and bipolar II (BP-II) disorders.MethodsWe included a total of 307 subjects with BPs in the analysis. For the factor analysis, we chose six variables related to clinical courses, 29 indicators covering lifetime symptoms of mood episodes, and 6 specific comorbid conditions. To determine the relationships among the identified phenotypic dimensions and their effects on differentiating BP subtypes, we applied structural equation modeling.ResultsWe selected a six-factor solution through scree plot, Velicer's minimum average partial test, and face validity evaluations; the six factors were cyclicity, depression, atypical vegetative symptoms, elation, psychotic/irritable mania, and comorbidity. In the path analysis, five factors excluding atypical vegetative symptoms were associated with one another. Cyclicity, depression, and comorbidity had positive associations, and they correlated negatively with psychotic/irritable mania; elation showed positive correlations with cyclicity and psychotic/irritable mania. Depression, cyclicity, and comorbidity were stronger in BP-II than in BP-I, and they contributed significantly to the distinction between the two disorders.ConclusionsWe identified six phenotype dimensions; in addition to symptom features of manic and depressive episodes, various comorbidities and high cyclicity constructed separate dimensions. Except for atypical vegetative symptoms, all factors showed a complex interdependency and played roles in discriminating BP-II from BP-I.