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'.

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.

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.

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.

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.

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.

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.

Phenotypic variability in Joubert syndrome is partially explained by ciliary pathophysiology.

INTRODUCTION: Joubert syndrome (JS) arises from defects of primary cilia resulting in potential malformations of the brain, kidneys, eyes, liver, and limbs. Several of the 35+ genes associated with JS have recognized genotype/phenotype correlations, but most genes have not had enough reported individuals to draw meaningful conclusions. METHODS: A PubMed literature review identified 688 individuals with JS across 32 genes and 112 publications to bolster known genotype/phenotype relationships and identify new correlations. All included patients had the "molar tooth sign" and a confirmed genetic diagnosis. Individuals were categorized by age, ethnicity, sex and the presence of developmental disability/intellectual disability, hypotonia, abnormal eye movements, ataxia, visual impairment, renal impairment, polydactyly, and liver abnormalities. RESULTS: Most genes demonstrated unique phenotypic profiles. Grouping proteins based on physiologic interactions established stronger phenotypic relationships that reflect known ciliary pathophysiology. Age-stratified data demonstrated that end-organ disease is progressive in JS. Most genes demonstrated a significant skew towards having variants with either residual protein function or no residual protein function. CONCLUSION: This cohort demonstrates that clinically meaningful genotype/phenotype relationships exist within most JS-related genes and can be referenced to allow for more personalized clinical care.

Mutational spectrum and phenotypic variability of Duchenne muscular dystrophy and related disorders in a Bangladeshi population.

Duchenne muscular dystrophy (DMD) is a severe rare neuromuscular disorder caused by mutations in the X-linked dystrophin gene. Several mutations have been identified, yet the full mutational spectrum, and their phenotypic consequences, will require genotyping across different populations. To this end, we undertook the first detailed genotype and phenotype characterization of DMD in the Bangladeshi population. We investigated the rare mutational and phenotypic spectrum of the DMD gene in 36 DMD-suspected Bangladeshi participants using an economically affordable diagnostic strategy involving initial screening for exonic deletions in the DMD gene via multiplex PCR, followed by testing PCR-negative patients for mutations using whole exome sequencing. The deletion mapping identified two critical DMD gene hotspot regions (near proximal and distal ends, spanning exons 8-17 and exons 45-53, respectively) that comprised 95% (21/22) of the deletions for this population cohort. From our exome analysis, we detected two novel pathogenic hemizygous mutations in exons 21 and 42 of the DMD gene, and novel pathogenic recessive and loss of function variants in four additional genes: SGCD, DYSF, COL6A3, and DOK7. Our phenotypic analysis showed that DMD suspected participants presented diverse phenotypes according to the location of the mutation and which gene was impacted. Our study provides ethnicity specific new insights into both clinical and genetic aspects of DMD.

Thurston syndrome with thalassaemia: a rare case devising a novel molecular and phenotypic variation.

A male infant presented with progressive paleness of the body since 3 months of age. On examination, the child had pallor, microcephaly with dysmorphic facies (depressed nasal bridge, low set ears, retrognathia, high arched palate and tongue hamartoma). Postaxial polydactyly in bilateral hands and feet, broad great toes, with syndactyly of left fourth and fifth toes were present. The haemogram showed severe anaemia with a microcytic hypochromic picture. High-performance liquid chromatography (HPLC) was normal. However, the parents' HPLC was suggestive of beta thalassaemia trait. Whole-exome sequencing revealed Thurston syndrome with beta-thalassaemia in homozygous pattern with a novel mutation. It is a rare genetic syndrome exclusively found in the South Asian population. Due to the rarity, identification of this syndrome is often difficult and requires awareness among clinicians. However, it is important to diagnose the disorder accurately in order to provide appropriate genetic counselling and prognostication to the parents.

Phenotypic variations induced emergence of orientation order and morphology in Bacillus subtilis biofilm growth.

During the Bacillus subtilis biofilm growth on the solid MSgg substrate, the biofilm exhibits highly ordered structures such as matrix-producing-cell chains and Van Gogh bundles due to bacterial orientation order. These structures make the biofilm have strong mobility and environmental adaptability, thus making bacteria easier to survive and thrive in biofilms comparing to planktonic bacteria. We tested the behaviors of different phenotypes as well as their impacts on bacterial clusters: motile cells arrange disorderly, the biofilm made up of motile cells tends to be circular and isotropic; matrix-producing cells form cellular chains that guide motile cells along the chain to form a locally nematic phase, the morphology of the biofilm made up of both motile cells and matrix-producing cells is rendered irregular. Combining the results of a coarse-grained and individual-based model, we can control the biofilm growth through regulating environmental friction, bacterial growth rate and adhesion between cells.

Combined reference-free and multi-reference based GWAS uncover cryptic variation underlying rapid adaptation in a fungal plant pathogen.

Microbial pathogens often harbor substantial functional diversity driven by structural genetic variation. Rapid adaptation from such standing variation threatens global food security and human health. Genome-wide association studies (GWAS) provide a powerful approach to identify genetic variants underlying recent pathogen adaptation. However, the reliance on single reference genomes and single nucleotide polymorphisms (SNPs) obscures the true extent of adaptive genetic variation. Here, we show quantitatively how a combination of multiple reference genomes and reference-free approaches captures substantially more relevant genetic variation compared to single reference mapping. We performed reference-genome based association mapping across 19 reference-quality genomes covering the diversity of the species. We contrasted the results with a reference-free (i.e., k-mer) approach using raw whole-genome sequencing data in a panel of 145 strains collected across the global distribution range of the fungal wheat pathogen Zymoseptoria tritici. We mapped the genetic architecture of 49 life history traits including virulence, reproduction and growth in multiple stressful environments. The inclusion of additional reference genome SNP datasets provides a nearly linear increase in additional loci mapped through GWAS. Variants detected through the k-mer approach explained a higher proportion of phenotypic variation than a reference genome-based approach and revealed functionally confirmed loci that classic GWAS approaches failed to map. The power of GWAS in microbial pathogens can be significantly enhanced by comprehensively capturing structural genetic variation. Our approach is generalizable to a large number of species and will uncover novel mechanisms driving rapid adaptation of pathogens.

Reproductive individuality of clonal fish raised in near-identical environments and its link to early-life behavioral individuality.

Recent studies have documented among-individual phenotypic variation that emerges in the absence of apparent genetic and environmental differences, but it remains an open question whether such seemingly stochastic variation has fitness consequences. We perform a life-history experiment with naturally clonal fish, separated directly after birth into near-identical (i.e., highly standardized) environments, quantifying 2522 offspring from 152 broods over 280 days. We find that (i) individuals differ consistently in the size of offspring and broods produced over consecutive broods, (ii) these differences are observed even when controlling for trade-offs between brood size, offspring size and reproductive onset, indicating individual differences in life-history productivity and (iii) early-life behavioral individuality in activity and feeding patterns, with among-individual differences in feeding being predictive of growth, and consequently offspring size. Thus, our study provides experimental evidence that even when minimizing genetic and environmental differences, systematic individual differences in life-history measures and ultimately fitness can emerge.