Reference Intervals Revisited: A Novel Model for Population-Based Reference Intervals, Using a Small Sample Size and Biological Variation Data.

BACKGROUND: Conventional population-based reference intervals (popRIs) are established on the ranking of single measurement results from at least 120 reference individuals. In this study, we aimed to explore a new model for popRIs, utilizing biological variation (BV) data to define the reference interval (RI) limits and compared BV-based popRI from different sample sizes with previously published conventional popRIs from the same population. METHODS: The model is based on defining the population set point (PSP) from single-measurement results of a group of reference individuals and using the total variation around the PSP, derived from the combination of BV and analytical variation, to define the RI limits. Using data from 143 reference individuals for 48 clinical chemistry and hematology measurands, BV-based popRIs were calculated for different sample sizes (n = 16, n = 30, and n = 120) and considered acceptable if they covered 90% of the population. In addition, simulation studies were performed to estimate the minimum number of required reference individuals. RESULTS: The median ratio of the BV-based to conventional RI ranges was 0.98. The BV-based popRIs calculated from the different samples were similar, and most met the coverage criterion. For 25 measurands ≤16 reference individuals and for 23 measurands >16 reference individuals were required to estimate the PSP. CONCLUSIONS: The BV-based popRI model delivered robust RIs for most of the included measurands. This new model requires a smaller group of reference individuals than the conventional popRI model and can be implemented if reliable BV data are available.

Revisiting the steatosis-associated fibrosis estimator score in young Asian subjects with steatotic liver disease and consideration for population variability.

BACKGROUND AND AIM: The steatosis-associated fibrosis estimator (SAFE) score has been developed to distinguish clinically significant fibrosis in patients with steatotic liver disease (SLD). However, validation of its performance in Asian subjects is limited. This study aimed to evaluate the performance of the SAFE score in Asian subjects with biopsy-proven SLD and in different subgroups according to age, sex, and body mass index. METHODS: We retrospectively analyzed 6383 living liver donors who underwent a liver biopsy between 2005 and 2023. Of these, 1551 subjects with biopsy-proven SLD were included. The performance of the SAFE score was evaluated using areas under the curve and compared with those of the nonalcoholic fatty liver disease fibrosis score (NFS) and fibrosis-4 index (FIB-4). RESULTS: The prevalence of clinically significant fibrosis in the cohort was 2.2%. The proportion of subjects with a "low-risk" SAFE score was the highest (91.0%), followed by those with "intermediate-risk" (7.8%) and "high-risk" (1.2%) scores. The prevalence of fibrosis in subjects with low-risk, intermediate-risk, and high-risk scores was 1.6%, 6.6%, and 21.1%, respectively. The SAFE outperformed FIB-4 and NFS (area under the curve: 0.70 vs 0.64 for both NFS and FIB-4). However, it showed low diagnostic accuracy and sensitivity (27%) at the low cutoff (SAFE < 0) in subjects aged 30-39 years (fibrosis: 1.2%), despite having a high negative predictive value (0.99). CONCLUSION: While the SAFE score demonstrates superior performance compared with other noninvasive tests in Asian subjects with SLD, its performance varies across age groups. In younger subjects, particularly, its performance may be more limited.

Among-population variation in drought responses is consistent across life stages but not between native and non-native ranges.

Understanding how widespread species adapt to variation in abiotic conditions across their ranges is fundamental to ecology. Insight may come from studying how among-population variation (APV) in the common garden corresponds with the environmental conditions of source populations. However, there are no such studies comparing native vs non-native populations across multiple life stages. We examined APV in the performance and functional traits of 59 Conyza canadensis populations, in response to drought, across large aridity gradients in the native (North America) and non-native (Eurasia) ranges in three experiments. Our treatment (dry vs wet) was applied at the recruitment, juvenile, and adult life stages. We found contrasting patterns of APV in drought responses between the two ranges. In the native range, plant performance was less reduced by drought in populations from xeric than mesic habitats, but such relationship was not apparent for non-native populations. These range-specific patterns were consistent across the life stages. The weak adaptive responses of non-native populations indicate that they can become highly abundant even without complete local adaptation to abiotic environments and suggest that long-established invaders may still be evolving to the abiotic environment. These findings may explain lag times in invasions and raise concern about future expansions.

Variation and plasticity in life-history traits and fitness of wild Arabidopsis thaliana populations are not related to their genotypic and ecological diversity.

BACKGROUND: Despite its implications for population dynamics and evolution, the relationship between genetic and phenotypic variation in wild populations remains unclear. Here, we estimated variation and plasticity in life-history traits and fitness of the annual plant Arabidopsis thaliana in two common garden experiments that differed in environmental conditions. We used up to 306 maternal inbred lines from six Iberian populations characterized by low and high genotypic (based on whole-genome sequences) and ecological (vegetation type) diversity. RESULTS: Low and high genotypic and ecological diversity was found in edge and core Iberian environments, respectively. Given that selection is expected to be stronger in edge environments and that ecological diversity may enhance both phenotypic variation and plasticity, we expected genotypic diversity to be positively associated with phenotypic variation and plasticity. However, maternal lines, irrespective of the genotypic and ecological diversity of their population of origin, exhibited a substantial amount of phenotypic variation and plasticity for all traits. Furthermore, all populations harbored maternal lines with canalization (robustness) or sensitivity in response to harsher environmental conditions in one of the two experiments. CONCLUSIONS: Overall, we conclude that the environmental attributes of each population probably determine their genotypic diversity, but all populations maintain substantial phenotypic variation and plasticity for all traits, which represents an asset to endure in changing environments.

Multilocus pathogenic variants contribute to intrafamilial clinical heterogeneity: a retrospective study of sibling pairs with neurodevelopmental disorders.

BACKGROUND: Multilocus pathogenic variants (MPVs) are genetic changes that affect multiple gene loci or regions of the genome, collectively leading to multiple molecular diagnoses. MPVs may also contribute to intrafamilial phenotypic variability between affected individuals within a nuclear family. In this study, we aim to gain further insights into the influence of MPVs on a disease manifestation in individual research subjects and explore the complexities of the human genome within a familial context. METHODS: We conducted a systematic reanalysis of exome sequencing data and runs of homozygosity (ROH) regions of 47 sibling pairs previously diagnosed with various neurodevelopmental disorders (NDD). RESULTS: We found siblings with MPVs driven by long ROH regions in 8.5% of families (4/47). The patients with MPVs exhibited significantly higher FROH values (p-value = 1.4e-2) and larger total ROH length (p-value = 1.8e-2). Long ROH regions mainly contribute to this pattern; the siblings with MPVs have a larger total size of long ROH regions than their siblings in all families (p-value = 6.9e-3). Whereas the short ROH regions in the siblings with MPVs are lower in total size compared to their sibling pairs with single locus pathogenic variants (p-value = 0.029), and there are no statistically significant differences in medium ROH regions between sibling pairs (p-value = 0.52). CONCLUSION: This study sheds light on the significance of considering MPVs in families with affected sibling pairs and the role of ROH as an adjuvant tool in explaining clinical variability within families. Identifying individuals carrying MPVs may have implications for disease management, identification of possible disease risks to different family members, genetic counseling and exploring personalized treatment approaches.

How do stochastic processes and genetic threshold effects explain incomplete penetrance and inform causal disease mechanisms?

Incomplete penetrance is the rule rather than the exception in Mendelian disease. In syndromic monogenic disorders, phenotypic variability can be viewed as the combination of incomplete penetrance for each of multiple independent clinical features. Within genetically identical individuals, such as isogenic model organisms, stochastic variation at molecular and cellular levels is the primary cause of incomplete penetrance according to a genetic threshold model. By defining specific probability distributions of causal biological readouts and genetic liability values, stochasticity and incomplete penetrance provide information about threshold values in biological systems. Ascertainment of threshold values has been achieved by simultaneous scoring of relatively simple phenotypes and quantitation of molecular readouts at the level of single cells. However, this is much more challenging for complex morphological phenotypes using experimental and reductionist approaches alone, where cause and effect are separated temporally and across multiple biological modes and scales. Here I consider how causal inference, which integrates observational data with high confidence causal models, might be used to quantify the relative contribution of different sources of stochastic variation to phenotypic diversity. Collectively, these approaches could inform disease mechanisms, improve predictions of clinical outcomes and prioritize gene therapy targets across modes and scales of gene function. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Analytical performance specifications based on biological variation data - considerations, strengths and limitations.

Analytical performance specifications (APS) are typically established through one of three models: (i) outcome studies, (ii) biological variation (BV), or (iii) state-of-the-art. Presently, The APS can, for most measurands that have a stable concentration, be based on BV. BV based APS, defined for imprecision, bias, total allowable error and allowable measurement uncertainty, are applied to many different processes in the laboratory. When calculating APS, it is important to consider the different APS formulae, for what setting they are to be applied and if they are suitable for the intended purpose. In this opinion paper, we elucidate the background, limitations, strengths, and potential intended applications of the different BV based APS formulas. When using BV data to set APS, it is important to consider that all formulae are contingent on accurate and relevant BV estimates. During the last decade, efficient procedures have been established to obtain reliable BV estimates that are presented in the EFLM biological variation database. The database publishes detailed BV data for numerous measurands, global BV estimates derived from meta-analysis of quality-assured studies of similar study design and automatic calculation of BV based APS.

Temporal phenotypic variation of spinach root traits and its relation to shoot performance.

The root system is important for the growth and development of spinach. To reveal the temporal variability of the spinach root system, root traits of 40 spinach accessions were measured at three imaging times (20, 30, and 43 days after transplanting) in this study using a non-destructive and non-invasive root analysis system. Results showed that five root traits were reliably measured by this system (RootViz FS), and two of which were highly correlated with manually measured traits. Root traits had higher variations than shoot traits among spinach accessions, and the trait of mean growth rate of total root length had the largest coefficients of variation across the three imaging times. During the early stage, only tap root length was weakly correlated with shoot traits (plant height, leaf width, and object area (equivalent to plant surface area)), whereas in the third imaging, root fresh weight, total root length, and root area were strongly correlated with shoot biomass-related traits. Five root traits (total root length, tap root length, total root area, root tissue density, and maximal root width) showed high variations with coefficients of variation values (CV ≥ 0.3, except maximal root width) and high heritability (H2 > 0.6) among the three stages. The 40 spinach accessions were classified into five subgroups with different growth dynamics of the primary and lateral roots by cluster analysis. Our results demonstrated the potential of in-situ phenotyping to assess dynamic root growth in spinach and provide new perspectives for biomass breeding based on root system ideotypes.

Predator selection on phenotypic variability of cryptic and aposematic moths.

Natural selection generally favours phenotypic variability in camouflaged organisms, whereas aposematic organisms are expected to evolve a more uniform warning coloration. However, no comprehensive analysis of the phenotypic consequences of predator selection in aposematic and cryptic species exists. Using state-of-the-art image analysis, we examine 2800 wing images of 82 moth species accessed via three online museum databases. We test whether anti-predator strategy (i.e., camouflage or aposematism) explains intraspecific variation in wing colour and pattern across northern hemisphere moths. In addition, we test two mutually non-exclusive, ecological hypotheses to explain variation in colour pattern: diel-activity or dietary-niche. In this work, taking into account phylogenetic relationships, moth phenotypic variability is best explained by anti-predator strategy with camouflaged moths being more variable in wing patterning than aposematic species.

Phylogeography and phenotypic wing shape variation in a damselfly across populations in Europe.

BACKGROUND: Describing geographical variation in morphology of organisms in combination with data on genetic differentiation and biogeography can provide important information on how natural selection shapes such variation. Here we study genetic structure using ddRAD seq and wing shape variation using geometric morphometrics in 14 populations of the damselfly Lestes sponsa along its latitudinal range in Europe. RESULTS: The genetic analysis showed a significant, yet relatively weak population structure with high genetic heterozygosity and low inbreeding coefficients, indicating that neutral processes contributed very little to the observed wing shape differences. The genetic analysis also showed that some regions of the genome (about 10%) are putatively shaped by selection. The phylogenetic analysis showed that the Spanish and French populations were the ancestral ones with northern Swedish and Finnish populations being the most derived ones. We found that wing shape differed significantly among populations and showed a significant quadratic (but weak) relationship with latitude. This latitudinal relationship was largely attributed to allometric effects of wing size, but non-allometric variation also explained a portion of this relationship. However, wing shape showed no phylogenetic signal suggesting that lineage-specific variation did not contribute to the variation along the latitudinal gradient. In contrast, wing size, which is correlated with body size in L. sponsa, had a strong negative correlation with latitude. CONCLUSION: Our results suggest a relatively weak population structure among the sampled populations across Europe, but a clear differentiation between south and north populations. The observed geographic phenotypic variation in wing shape may have been affected by different local selection pressures or environmental effects.

Phenotypic variation in biting behavior associated with differences in expression of olfactory genes in the vector mosquito, Aedes albopictus (Diptera: Culicidae).

We evaluated miRNA and mRNA expression differences in head tissues between avid-biting vs. reluctant-biting Aedes albopictus (Skuse) females from a single population over a 20-min timescale. We found no differences in miRNA expression between avid vs. reluctant biters, indicating that translational modulation of blood-feeding behavior occurs on a longer timescale than mRNA transcription. In contrast, we detected 19 differentially expressed mRNAs. Of the 19 differentially expressed genes at the mRNA level between avid-biting vs. reluctant-biting A. albopictus, 9 are implicated in olfaction, consistent with the well-documented role of olfaction in mosquito host-seeking. Additionally, several of the genes that we identified as differentially expressed in association with phenotypic variation in biting behavior share similar functions with or are inferred orthologues of, genes associated with evolutionary variation in biting behaviors of Wyeomyia smithii (Coq.) and Culex pipiens (Lin.). A future goal is to determine whether these genes are involved in the evolutionary transition from a biting to a non-biting life history.

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison.

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Genome dosage alteration caused by chromosome pyramiding and shuffling effects on karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids.

KEY MESSAGE: We developed an array of Zea-Tripsacum tri-hybrid allopolyploids with multiple ploidies. We unveiled that changes in genome dosage due to the chromosomes pyramiding and shuffling of three species effects karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids. Polyploidy, or whole genome duplication, has played a major role in evolution and speciation. The genomic consequences of polyploidy have been extensively studied in many plants; however, the extent of chromosomal variation, genome dosage, phenotypic diversity, and heterosis in allopolyploids derived from multiple species remains largely unknown. To address this question, we synthesized an allohexaploid involving Zea mays, Tripsacum dactyloides, and Z. perennis by chromosomal pyramiding. Subsequently, an allooctoploid and an allopentaploid were obtained by hybridization of the allohexaploid with Z. perennis. Moreover, we constructed three populations with different ploidy by chromosomal shuffling (allopentaploid × Z. perennis, allohexaploid × Z. perennis, and allooctoploid × Z. perennis). We have observed 3 types of sexual reproductive modes and 2 types of asexual reproduction modes in the tri-species hybrids, including 2n gamete fusion (2n + n), haploid gamete fusion (n + n), polyspermy fertilization (n + n + n) or 2n gamete fusion (n + 2n), haploid gametophyte apomixis, and asexual reproduction. The tri-hybrids library presents extremely rich karyotype heterogeneity. Chromosomal compensation appears to exist between maize and Z. perennis. A rise in the ploidy of the trihybrids was linked to a higher frequency of chromosomal translocation. Variation in the degree of phenotypic diversity observed in different segregating populations suggested that genome dosage effects phenotypic manifestation. These findings not only broaden our understanding of the mechanisms of polyploid formation and reproductive diversity but also provide a novel insight into genome pyramiding and shuffling driven genome dosage effects and phenotypic diversity.

Characterization of Phenotypic Variability in Becker Muscular Dystrophy for Clinical Practice and Towards Trial Readiness: A Two-Years Follow up Study.

Background: Becker muscular dystrophy (BMD) is a dystrophinopathy due to in-frame mutations in the dystrophin gene (DMD) which determines a reduction of dystrophin at muscle level. BMD has a wide spectrum of clinical variability with different degrees of disability. Studies of natural history are needed also in view of up-coming clinical trials. Objectives: From an initial cohort of 32 BMD adult subjects, we present a detailed phenotypic characterization of 28 patients, then providing a description of their clinical natural history over the course of 12 months for 18 and 24 months for 13 of them. Methods: Each patient has been genetically characterized. Baseline, and 1-year and 2 years assessments included North Star Ambulatory Assessment (NSAA), timed function tests (time to climb and descend four stairs), 6-minute walk test (6MWT), Walton and Gardner-Medwin Scale and Medical Research Council (MRC) scale. Muscle magnetic resonance imaging (MRI) was acquired at baseline and in a subgroup of 9 patients after 24 months. Data on cardiac function (electrocardiogram, echocardiogram, and cardiac MRI) were also collected. Results and conclusions: Among the clinical heterogeneity, a more severe involvement is often observed in patients with 45-X del, with a disease progression over two years. The 6MWT appears sensitive to detect modification from baseline during follow up while no variation was observed by MRC testing. Muscle MRI of the lower limbs correlates with clinical parameters.Our study further highlights how the phenotypic variability of BMD adult patients makes it difficult to describe an uniform course and substantiates the need to identify predictive parameters and biomarkers to stratify patients.

Non-genetic phenotypic variability affects populations and communities in protist microcosms.

Intraspecific trait variation (ITV), potentially driven by genetic and non-genetic mechanisms, can underlie variability in resource acquisition, individual fitness and ecological interactions. Impacts of ITV at higher levels of biological organizations are hence likely, but up-scaling our knowledge about ITV importance to communities and comparing its relative effects at population and community levels has rarely been investigated. Here, we tested the effects of genetic and non-genetic ITV on morphological traits in microcosms of protist communities by contrasting the effects of strains showing different ITV levels (i.e. trait averages and variance) on population growth, community composition and biomass production. We found that genetic and non-genetic ITV can lead to different effects on populations and communities across several generations. Furthermore, the effects of ITV declined across levels of biological organization: ITV directly altered population performance, with cascading but indirect consequences for community composition and biomass productivity. Overall, these results show that the drivers of ITV can have distinct effects on populations and communities, with cascading impacts on higher levels of biological organization that might mediate biodiversity-ecosystem functioning relationships.

Megaproject will chart human immune diversity.

Human Immunome Project's survey aims to improve drugs and vaccines.

Variable Phenotype of Congenital Corneal Opacities in Biallelic CYP1B1 Pathogenic Variants.

PURPOSE: The aim of this study is to describe the variable phenotype of congenital corneal opacities occurring in patients with biallelic CYP1B1 pathogenic variants. METHODS: A retrospective chart review was conducted to identify patients with congenital corneal opacities and CYP1B1 pathogenic variants seen at UPMC Children's Hospital of Pittsburgh. Ophthalmic examination, high-frequency ultrasound, anterior segment optical coherence tomography, histopathologic images, and details of genetic testing were reviewed. RESULTS: Three children were identified. All presented with raised intraocular pressure. Two patients showed bilateral limbus-to-limbus avascular corneal opacification that did not resolve with intraocular pressure control; 1 showed unilateral avascular corneal opacity with a crescent of clear cornea, iridocorneal adhesions, iridolenticular adhesions, and classical features of congenital glaucoma in the fellow eye (enlarged corneal diameter, Haab striae, and clearing of the corneal clouding with appropriate intraocular pressure control). The first 2 patients were visually rehabilitated with penetrating keratoplasty. Histopathology revealed distinct features: a variably keratinized epithelium; a thick but discontinuous Bowman-like layer with areas of disruption and abnormal cellularity; Descemet membrane, when observed, showed reduced endothelial cells; and no pathological changes of Haab striae were identified. Two patients had compound heterozygous pathogenic variants in CYP1B1 causing premature stop codons, whereas 1 was homozygous for a pathogenic missense variant. CONCLUSIONS: Congenital corneal opacities seen in biallelic CYP1B1 pathogenic variants have a variable phenotype. One is that commonly termed as Peters anomaly type 1 (with iridocorneal adhesions, with or without iridolenticular adhesions) and the other is a limbus-to-limbus opacity, termed CYP1B1 cytopathy. Clinicians should be aware of this phenotypic variability.

Digenic inheritance accounts for phenotypic variability in amelogenesis imperfecta.

Amelogenesis imperfecta (AI) represents a group of clinically and genetically heterogeneous disorders that affect enamel formation and mineralization. Although AI is commonly considered a monogenic disorder, digenic inheritance is rarely reported. In this study, we recruited two nonconsanguineous Chinese families exhibiting diverse phenotypes of enamel defects among affected family members. Digenic variants were discovered in both probands. In family 1, the proband inherited a paternal frameshift variant in LAMA3 (NM_198129.4:c.3712dup) and a maternal deletion encompassing the entire AMELX gene. This resulted in a combined hypoplastic and hypomineralized AI phenotype, which was distinct from the parents' manifestations. In family 2, whole-exome sequencing analysis revealed the proband carried a maternal heterozygous splicing variant in COL17A1 (NC_000010.11 (NM_000494.3): c.4156 + 2dup) and compound heterozygous variants in RELT (paternal: NM_032871.4:c.260A > T; maternal: NM_032871.4:c.521 T > G). These genetic changes caused the abundant irregular enamel defects observed in the proband, whereas other affected family members carrying heterozygous variants in both COL17A1 and RELT displayed only horizontal grooves as their phenotype. The pathogenicity of the novel COL17A1 splice site variant was confirmed through RT-PCR and minigene assay. This study enhances our understanding by highlighting the potential association between the co-occurrence of variants in two genes and variable phenotypes observed in AI patients.