Identification of Medically Actionable Secondary Findings in the 1000 Genomes

The American College of Medical Genetics and Genomics (ACMG) recommends that clinical sequencing laboratories return secondary findings in 56 genes associated with medically actionable conditions. Our goal was to apply a systematic, stringent approach consistent with clinical standards to estimate the prevalence of pathogenic variants associated with such conditions using a diverse sequencing reference sample. Candidate variants in the 56 ACMG genes were selected from Phase 1 of the 1000 Genomes dataset, which contains sequencing information on 1,092 unrelated individuals from across the world. These variants were filtered using the Human Gene Mutation Database (HGMD) Professional version and defined parameters, appraised through literature review, and examined by a clinical laboratory specialist and expert physician. Over 70,000 genetic variants were extracted from the 56 genes, and filtering identified 237 variants annotated as disease causing by HGMD Professional. Literature review and expert evaluation determined that 7 of these variants were pathogenic or likely pathogenic. Furthermore, 5 additional truncating variants not listed as disease causing in HGMD Professional were identified as likely pathogenic. These 12 secondary findings are associated with diseases that could inform medical follow-up, including cancer predisposition syndromes, cardiac conditions, and familial hypercholesterolemia. The majority of the identified medically actionable findings were in individuals from the European (5/379) and Americas (4/181) ancestry groups, with fewer findings in Asian (2/286) and African (1/246) ancestry groups. Our results suggest that medically relevant secondary findings can be identified in approximately 1% (12/1092) of individuals in a diverse reference sample. As clinical sequencing laboratories continue to implement the ACMG recommendations, our results highlight that at least a small number of potentially important secondary findings can be selected for return. Our results also confirm that understudied populations will not reap proportionate benefits of genomic medicine, highlighting the need for continued research efforts on genetic diseases in these populations.     

Evolution of the 28S ribosomal RNA gene in anurans: regions of variability and their phylogenetic implications

Fifteen restriction sites were mapped to the 28S ribosomal RNA gene of individuals representing 54 species of frogs, two species of salamanders, a caecilian, and a lungfish. Eight of these sites were present in all species examined, and two were found in all but one species. Alignment of these conserved restriction sites revealed, among anuran 28S rRNA genes, five regions of major length variation that correspond to four of 12 previously identified divergent domains of this gene. One of the divergent domains (DD8) consists of two regions of length variation separated by a short segment that is conserved at least throughout tetrapods. Most of the insertions, deletions, and restriction-site variations identified in the 28S gene will require sequence-level analysis for a detailed reconstruction of their history. However, an insertion in DD9 that is coextensive with frogs in the suborder Neobatrachia, a BstEII site that is limited to representatives of two leptodactylid subfamilies, and a deletion in DD10 that is found only in three ranoid genera are probably synapomorphies.      

Functional and phenotypic changes associated with the in vitro development of human monocytes into activated macrophages

Freshly isolated human peripheral blood monocytes had minimal cytotoxic effect in vitro on the schistosomula of Schistosoma mansoni. However, stimulation of the cells with either interferon gamma (IFN) or specific anti-parasite antiserum caused an increase in cytotoxicity. Additionally, the normal development of monocytes into macrophages over 7 days was associated with a sharp increase in cytotoxicity. The non-cytotoxic monocytes were compared with activated macrophages to assess whether cytotoxicity was associated with changes in immunophenotype. As monocytes developed into macrophages there were marked increases in transferrin receptors (HB21), macrophage cellular integrin (3.9), and Fc receptors (KB61). A further three markers showed increased expression in 7-day-old macrophages stimulated by IFN, namely a high affinity Fc gamma receptor (10.1), MHC Class II (1B5) and tumour necrosis factor (TNF).     

Anthocyanin production in chl-rich and chl-poor seedlings

Screening by chlorophyll (Chl) affects photoconversion rates and photoequilibrium ratios of phytochrome in vivo and may cause distortion of the action spectra of photomorphogenesis (N Kazarinova-Fukshansky, M Seyfried, E Schäfer 1985 Photochem Photobiol 41: 689-702). Inhibitors that reduce the Chl content of seedlings are sometimes used in photomorphogenesis research to decrease the effects of Chl screening on the state of phytochrome in vivo. Streptomycin is one of the inhibitors that can be used for this purpose. The effects of streptomycin on phytochrome-mediated anthocyanin accumulation in young seedlings are significantly different in closely related systems. The use of `Chl-bleachers' in photomorphogenesis studies may produce undesirable side effects. At the level of the expression of a photoregulated response, the effects of differences in the state of phytochrome between water-grown Chl-rich and inhibitor-treated Chl-poor seedlings may be difficult to evaluate because they may be masked by the effects of the inhibitor on the response.     

Characterization of cross-bridge elasticity and kinetics of cross-bridge cycling during force development in single smooth muscle cells

Force development in smooth muscle, as in skeletal muscle, is believed to reflect recruitment of force-generating myosin cross-bridges. However, little is known about the events underlying cross-bridge recruitment as the muscle cell approaches peak isometric force and then enters a period of tension maintenance. In the present studies on single smooth muscle cells isolated from the toad (Bufo marinus) stomach muscularis, active muscle stiffness, calculated from the force response to small sinusoidal length changes (0.5% cell length, 250 Hz), was utilized to estimate the relative number of attached cross-bridges. By comparing stiffness during initial force development to stiffness during force redevelopment immediately after a quick release imposed at peak force, we propose that the instantaneous active stiffness of the cell reflects both a linearly elastic cross-bridge element having 1.5 times the compliance of the cross-bridge in frog skeletal muscle and a series elastic component having an exponential length-force relationship. At the onset of force development, the ratio of stiffness to force was 2.5 times greater than at peak isometric force. These data suggest that, upon activation, cross-bridges attach in at least two states (i.e., low-force-producing and high-force-producing) and redistribute to a steady state distribution at peak isometric force. The possibility that the cross-bridge cycling rate was modulated with time was also investigated by analyzing the time course of tension recovery to small, rapid step length changes (0.5% cell length in 2.5 ms) imposed during initial force development, at peak force, and after 15 s of tension maintenance. The rate of tension recovery slowed continuously throughout force development following activation and slowed further as force was maintained. Our results suggest that the kinetics of force production in smooth muscle may involve a redistribution of cross-bridge populations between two attached states and that the average cycling rate of these cross-bridges becomes slower with time during contraction.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

15N resonance assignments of oxidized and reducedChromatium vinosum high-potential iron protein

The¹⁵N resonances in reduced and oxidizedChromatium vinosum high-potential iron protein have been assigned by use of¹H-¹H COSY spectra and¹H-¹⁵N HMQC, HMQC-COSY, and HMQC-NOESY spectra. Unambiguous assignment of 70 of 85 backbone¹⁵N resonances in the reduced protein and 62 of 85 resonances in the oxidized protein are made, as are 12 of 21 side-chain¹⁵N resonances.     

Prediction of near-term climate change impacts on UK wheat quality and the potential for adaptation through plant breeding

Wheat is a major crop worldwide, mainly cultivated for human consumption and animal feed. Grain quality is paramount in determining its value and downstream use. While we know that climate change threatens global crop yields, a better understanding of impacts on wheat end-use quality is also critical. Combining quantitative genetics with climate model outputs, we investigated UK-wide trends in genotypic adaptation for wheat quality traits. In our approach, we augmented genomic prediction models with environmental characterisation of field trials to predict trait values and climate effects in historical field trial data between 2001 and 2020. Addition of environmental covariates, such as temperature and rainfall, successfully enabled prediction of genotype by environment interactions (G × E), and increased prediction accuracy of most traits for new genotypes in new year cross validation. We then extended predictions from these models to much larger numbers of simulated environments using climate scenarios projected under Representative Concentration Pathways 8.5 for 2050–2069. We found geographically varying climate change impacts on wheat quality due to contrasting associations between specific weather covariables and quality traits across the UK. Notably, negative impacts on quality traits were predicted in the East of the UK due to increased summer temperatures while the climate in the North and South-west may become more favourable with increased summer temperatures. Furthermore, by projecting 167,040 simulated future genotype–environment combinations, we found only limited potential for breeding to exploit predictable G × E to mitigate year-to-year environmental variability for most traits except Hagberg falling number. This suggests low adaptability of current UK wheat germplasm across future UK climates. More generally, approaches demonstrated here will be critical to enable adaptation of global crops to near-term climate change.     

Use of Plant Strategy Ordination,DCA and ANOVA to elucidate relationships among habitats of Salix planifolia and Salix monticola

Abstract. Detrended Correspondence Analysis (DCA) was used on a data set of 112 species from 14 stands dominated by Salix planifolia and/or S. monticola. Environmental variables were determined which might be important to explain differences in stand vegetation composition, and willow stature (height), in willow carrs in the Rocky Mountain National Park. Correlation of environmental variables with DCA stand scores indicate that peat depth and soil redox potential are highly related to the primary DCA axis. Soils in stands of S. planifolia are more reduced in the early summer and have greater peat depths, than soils in stands of S. monticola. The second DCA axis is related to a water chemistry gradient of hydrogen ion, sulphate, and ammonium concentration. Four environmental variables selected by correlation analysis were regressed on stand scores from the first two DCA axes. Regression coefficients of August soil redox potential and soil water pH were significantly non-zero on the primary DCA axis. In addition to overall vegetation relationships elucidated by DCA, the relationships between willow stature and environmental variables were tested by simultaneous ANOVAs. ANOVA results and DCA ordinations indicate that S. planifolia and S. monticola plants achieve the greatest stature in the least reduced soils, and in shallow peat. Competitive, stress tolerant, and ruderal plant strategy theory applied to the vegetation data produced an arrangement of stands highly correlated to the DCA Axis 1 stand arrangement. These correlations indicate that the primary axis of vegetation structure can be extracted either by species composition, or species morphology/physiology relationships. Concurrent use of DCA and plant strategy theory is suggested as a tool for the prediction of community composition.