A deep catalogue of protein-coding variation in 983,578 individuals.

Rare coding variants that significantly impact function provide insights into the biology of a gene1-3. However, ascertaining their frequency requires large sample sizes4-8. Here, we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. 23% of the Regeneron Genetics Center Million Exome data (RGC-ME) comes from non-European individuals of African, East Asian, Indigenous American, Middle Eastern, and South Asian ancestry. This catalogue includes over 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss-of-function, we identify 3,988 loss-of-function intolerant genes, including 86 that were previously assessed as tolerant and 1,153 lacking established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions depleted of missense variants despite being tolerant to pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this important resource of coding variation from the RGC-ME accessible via a public variant allele frequency browser.

Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Exome sequencing and analysis of 454,787 UK Biobank participants.

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing1 to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study2. We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

New Guinea has the world's richest island flora.

New Guinea is the world's largest tropical island and has fascinated naturalists for centuries1,2. Home to some of the best-preserved ecosystems on the planet3 and to intact ecological gradients-from mangroves to tropical alpine grasslands-that are unmatched in the Asia-Pacific region4,5, it is a globally recognized centre of biological and cultural diversity6,7. So far, however, there has been no attempt to critically catalogue the entire vascular plant diversity of New Guinea. Here we present the first, to our knowledge, expert-verified checklist of the vascular plants of mainland New Guinea and surrounding islands. Our publicly available checklist includes 13,634 species (68% endemic), 1,742 genera and 264 families-suggesting that New Guinea is the most floristically diverse island in the world. Expert knowledge is essential for building checklists in the digital era: reliance on online taxonomic resources alone would have inflated species counts by 22%. Species discovery shows no sign of levelling off, and we discuss steps to accelerate botanical research in the 'Last Unknown'8.

Activating STING1-dependent immune signaling in TP53 mutant and wild-type acute myeloid leukemia.

DNA methyltransferase inhibitors (DNMTis) reexpress hypermethylated genes in cancers and leukemias and also activate endogenous retroviruses (ERVs), leading to interferon (IFN) signaling, in a process known as viral mimicry. In the present study we show that in the subset of acute myeloid leukemias (AMLs) with mutations in TP53, associated with poor prognosis, DNMTis, important drugs for treatment of AML, enable expression of ERVs and IFN and inflammasome signaling in a STING-dependent manner. We previously reported that in solid tumors poly ADP ribose polymerase inhibitors (PARPis) combined with DNMTis to induce an IFN/inflammasome response that is dependent on STING1 and is mechanistically linked to generation of a homologous recombination defect (HRD). We now show that STING1 activity is actually increased in TP53 mutant compared with wild-type (WT) TP53 AML. Moreover, in TP53 mutant AML, STING1-dependent IFN/inflammatory signaling is increased by DNMTi treatment, whereas in AMLs with WT TP53, DNMTis alone have no effect. While combining DNMTis with PARPis increases IFN/inflammatory gene expression in WT TP53 AML cells, signaling induced in TP53 mutant AML is still several-fold higher. Notably, induction of HRD in both TP53 mutant and WT AMLs follows the pattern of STING1-dependent IFN and inflammatory signaling that we have observed with drug treatments. These findings increase our understanding of the mechanisms that underlie DNMTi + PARPi treatment, and also DNMTi combinations with immune therapies, suggesting a personalized approach that statifies by TP53 status, for use of such therapies, including potential immune activation of STING1 in AML and other cancers.

Facultative mycorrhization in a fern (Struthiopteris spicant L. Weiss) is bound to light intensity.

BACKGROUND: The establishment of mycorrhizal relationships between a fungus and a plant typically enhances nutrient and water uptake for the latter while securing a carbon source for the fungus. However, under a particular set of environmental conditions, such as low availability of light and abundant nutrients in the soil, the resources invested in the maintenance of the fungi surpass the benefits obtained by the host. In those cases, facultative mycorrhizal plants are capable of surviving without symbiosis. Facultative mycorrhization in ferns has been overlooked until now. The present study measured the response of Struthiopteris spicant L. Weiss, and its root-associated fungi to different levels of light and nutrient availability in terms of growth, mycorrhizal presence, and leaf nutrient content. This fern species exhibits a great tolerance to variable light, nutrient, and pH conditions, and it has been found with and without mycorrhizae. We conducted a greenhouse experiment with 80 specimens of S. spicant and three factors (Light, Phosphorus, and Nitrogen) resulting in eight treatments. RESULTS: We found a significant influence of the factor light on fungal community composition, plant biomass, and nutrient accumulation. Departing from a lack of colonization at the initial stage, plants showed a remarkable increment of more than 80% in the arbuscular mycorrhizal fungi (AMF) richness and abundance in their roots when grown under high light conditions, compared with the ones in low light. We also observed an upward trend of C:P and C:N ratios and the above- and belowground biomass production when AMF abundance increased. Furthermore, the compositional analysis of the whole fungal communities associated with S. spicant roots revealed clear differences among low-light and high-light treatments. CONCLUSIONS: This study is the first to investigate the importance of light and nutrient availability in determining fern-AMF relationships. We confirmed that Struthiopteris spicant is a facultative mycorrhizal plant. The composition and diversity of AMF found in the roots of this fern are strongly influenced by light and less by nutrient conditions. Our study shows that ferns respond very sensitively to changes in environmental factors, leading to shifts in the associated mycorrhizal communities.

Light, rather than circadian rhythm, regulates gas exchange in ferns and lycophytes.

Circadian regulation plays a vital role in optimizing plant responses to the environment. However, while circadian regulation has been extensively studied in angiosperms, very little is known for lycophytes and ferns, leaving a gap in our understanding of the evolution of circadian rhythms across the plant kingdom. Here, we investigated circadian regulation in gas exchange through stomatal conductance and photosynthetic efficiency in a phylogenetically broad panel of 21 species of lycophytes and ferns over a 46 h period under constant light and a selected few under more natural conditions with day-night cycles. No rhythm was detected under constant light for either lycophytes or ferns, except for two semi-aquatic species of the family Marsileaceae (Marsilea azorica and Regnellidium diphyllum), which showed rhythms in stomatal conductance. Furthermore, these results indicated the presence of a light-driven stomatal control for ferns and lycophytes, with a possible passive fine-tuning through leaf water status adjustments. These findings support previous evidence for the fundamentally different regulation of gas exchange in lycophytes and ferns compared to angiosperms, and they suggest the presence of alternative stomatal regulations in Marsileaceae, an aquatic family already well known for numerous other distinctive physiological traits. Overall, our study provides evidence for heterogeneous circadian regulation across plant lineages, highlighting the importance of broad taxonomic scope in comparative plant physiology studies.

Hyper-Cryptic radiation of a tropical montane plant lineage.

Species are seen as the fundamental unit of biotic diversity, and thus their delimitation is crucial for defining measures for diversity assessments and studying evolution. Differences between species have traditionally been associated with variation in morphology. And yet, the discovery of cryptic diversity suggests that the evolution of distinct lineages does not necessarily involve morphological differences. Here, we analyze 1,684,987 variant sites and over 4,000 genes for more than 400 samples to show how a tropical montane plant lineage (Geonoma undata species complex) is composed of numerous unrecognized genetic groups that are not morphologically distinct. We find that 11 to 14 clades do not correspond to the three currently recognized species. Most clades are genetically different and geographic distance and topography are the most important factors determining this genetic divergence. The genetic structure of this lineage does not match its morphological variation. Instead, this species complex constitutes the first example of a hyper-cryptic plant radiation in tropical mountains.

Spatial patterns and climatic drivers of phylogenetic structure of regional liverwort assemblages in China.

BACKGROUND AND AIMS: Latitudinal diversity gradients have been intimately linked to the tropical niche conservatism hypothesis, which posits that there has been a strong filter due to the challenges faced by ancestral tropical lineages to adapt to low temperatures and colonize extra-tropical regions. In liverworts, species richness is higher towards the tropics, but the centers of diversity of the basal lineages are distributed across extra-tropical regions, pointing to the colonization of tropical regions by phylogenetically clustered assemblages of species of temperate origin. Here, we test this hypothesis through analyses on the relationship between macroclimatic variation and phylogenetic diversity in Chinese liverworts. METHODS: Phylogenetic diversity metrics and their standardized effect sizes for liverworts in each of the 28 regional floras at the province level in China were related to latitude and six climate variables using regression analysis. We conducted variation partitioning analyses to determine the relative importance of each group of climatic variables. KEY RESULTS: We find that the number of species decreases with latitude, whereas phylogenetic diversity shows the reverse pattern, and that phylogenetic diversity is more strongly correlated with temperature-related variables, compared with precipitation-related variables. CONCLUSIONS: We interpret the opposite patterns observed in phylogenetic diversity and species richness in terms of a more recent origin of tropical diversity coupled with higher extinctions in temperate regions.

Global patterns and climatic determinants of phylogenetic structure of regional fern floras.

Knowledge of relationships between phylogenetic structure of a biological assemblage and ecological factors that drive the variation of phylogenetic structure among regions is crucial for understanding the causes of variation in taxonomic composition and richness among regions, but this knowledge is lacking for the global flora of ferns. Here, we fill this critical knowledge gap. We divided the globe into 392 geographic units on land, collated species lists of ferns for each geographic unit, and used different phylogenetic metrics (tip- vs basal-weighted) reflecting different evolutionary depths to quantify phylogenetic structure. We then related taxonomic and phylogenetic structure metrics to six climatic variables for ferns as a whole and for two groups of ferns (old clades vs polypods) reflecting different evolutionary histories across the globe and within each continental region. We found that when old clades and polypods were considered separately, temperature-related variables explained more variation in these metrics than did precipitation-related variables in both groups. When analyses were conducted for continental regions separately, this pattern holds in most cases. Climate extremes have a stronger relationship with phylogenetic structure of ferns than does climate seasonality. Climatic variables explained more variation in phylogenetic structure at deeper evolutionary depths.

Phylogenomic evolutionary insights in the fern family Gleicheniaceae.

The pantropical fern family Gleicheniaceae comprises approximately 157 species. Seven genera are currently recognized in the family, although their monophyly is still uncertain due to low sampling in phylogenetic studies. We examined the monophyly of the genera through extended sampling, using the first phylogenomic inference of the family including data from both nuclear and plastid genomes. Seventy-six samples were sequenced (70 Gleicheniaceae species and six outgroups) using high throughput sequencing, including all seven currently recognized genera. Plastid and nuclear data were recovered and assembled; the nuclear data was phased to reduce paralogy as well as hybrid noise in the final recovered topology. Maximum likelihood trees were built for each locus, and a concatenated dataset was built for both datasets. A species tree based on a multispecies coalescent model was generated, and divergence time analyses performed. We here present the first genomic phylogenetic inferences concerning Gleicheniaceae, confirming the monophyly of most genera except Sticherus, which we recovered as paraphyletic. Although most of the extant genera of Gleicheniaceae originated during the Mesozoic, several genera show Neogene and even Quaternary diversifications, and our results suggest that reticulation and polyploidy may have played significant roles during this diversification. However, some genera, such as Rouxopteris and Stromatopteris, appear to represent evolutionary relicts.

Splitting one species into 22: an unusual tripling of molecular, morphological, and geographical differentiation in the fern family Didymochlaenaceae (Polypodiales).

The pantropical fern genus Didymochlaena (Didymochlaenaceae) has long been considered to contain one species only. Recent studies have resolved this genus/family as either sister to the rest of eupolypods I or as the second branching lineage of eupolypods I, and have shown that this genus is not monospecific, but the exact species diversity is unknown. In this study, a new phylogeny is reconstructed based on an expanded taxon sampling and six molecular markers. Our major results include: (i) Didymochlaena is moderately or weakly supported as sister to the rest of eupolypods I, highlighting the difficulty in resolving the relationships of this important fern lineage in the polypods; (ii) species in Didymochlaena are resolved into a New World clade and an Old World clade, and the latter further into an African clade and an Asian-Pacific clade; (iii) an unusual tripling of molecular, morphological and geographical differentiation in Didymochlaena is detected, suggesting single vicariance or dispersal events in individual regions and no evidence for reversals at all, followed by allopatric speciation at more or less homogeneous rates; (iv) evolution of 18 morphological characters is inferred and two morphological synapomorphies defining the family are recognized-the elliptical sori and fewer than 10 sori per pinnule, the latter never having been suggested before; (v) based on morphological and molecular variation, 22 species in the genus are recognized contrasting with earlier estimates of between one and a few; and (vi) our biogeographical analysis suggests an origin for Didymochlaena in the latest Jurassic-earliest Cretaceous and the initial diversification of the extant lineages in the Miocene-all but one species diverged from their sisters within the last 27 Myr, in most cases associated with allopatric speciation owing to geologic and climatic events, or dispersal.

Influence of the Degree of Swelling on the Stiffness and Toughness of Microgel-Reinforced Hydrogels.

The stiffness and toughness of conventional hydrogels decrease with increasing degree of swelling. This behavior makes the stiffness-toughness compromise inherent to hydrogels even more limiting for fully swollen ones, especially for load-bearing applications. The stiffness-toughness compromise of hydrogels can be addressed by reinforcing them with hydrogel microparticles, microgels, which introduce the double network (DN) toughening effect into hydrogels. However, to what extent this toughening effect is maintained in fully swollen microgel-reinforced hydrogels (MRHs) is unknown. Herein, it is demonstrated that the initial volume fraction of microgels contained in MRHs determines their connectivity, which is closely yet nonlinearly related to the stiffness of fully swollen MRHs. Remarkably, if MRHs are reinforced with a high volume fraction of microgels, they stiffen upon swelling. By contrast, the fracture toughness linearly increases with the effective volume fraction of microgels present in the MRHs regardless of their degree of swelling. These findings provide a universal design rule for the fabrication of tough granular hydrogels that stiffen upon swelling and hence, open up new fields of use of these hydrogels.

Evolutionary history of modern Samoans.

Archaeological studies estimate the initial settlement of Samoa at 2,750 to 2,880 y ago and identify only limited settlement and human modification to the landscape until about 1,000 to 1,500 y ago. At this point, a complex history of migration is thought to have begun with the arrival of people sharing ancestry with Near Oceanic groups (i.e., Austronesian-speaking and Papuan-speaking groups), and was then followed by the arrival of non-Oceanic groups during European colonialism. However, the specifics of this peopling are not entirely clear from the archaeological and anthropological records, and is therefore a focus of continued debate. To shed additional light on the Samoan population history that this peopling reflects, we employ a population genetic approach to analyze 1,197 Samoan high-coverage whole genomes. We identify population splits between the major Samoan islands and detect asymmetrical gene flow to the capital city. We also find an extreme bottleneck until about 1,000 y ago, which is followed by distinct expansions across the islands and subsequent bottlenecks consistent with European colonization. These results provide for an increased understanding of Samoan population history and the dynamics that inform it, and also demonstrate how rapid demographic processes can shape modern genomes.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population.

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.

Hydroelastomers: soft, tough, highly swelling composites.

Inspired by the cellular design of plant tissue, we present an approach to make versatile, tough, highly water-swelling composites. We embed highly swelling hydrogel particles inside tough, water-permeable, elastomeric matrices. The resulting composites, which we call hydroelastomers, combine the properties of their parent phases. From their hydrogel component, the composites inherit the ability to highly swell in water. From the elastomeric component, the composites inherit excellent stretchability and fracture toughness, while showing little softening as they swell. Indeed, the fracture properties of the composite match those of the best-performing, tough hydrogels, exhibiting fracture energies of up to 10 kJ m-2. Our composites are straightforward to fabricate, based on widely-available materials, and can easily be molded or extruded to form shapes with complex swelling geometries. Furthermore, there is a large design space available for making hydroelastomers, since one can use any hydrogel as the dispersed phase in the composite, including hydrogels with stimuli-responsiveness. These features make hydroelastomers excellent candidates for use in soft robotics and swelling-based actuation, or as shape-morphing materials, while also being useful as hydrogel replacements in other fields.

The relationship between chlorophyllous spores and mycorrhizal associations in ferns: evidence from an evolutionary approach.

PREMISE: Approximately 14% of all fern species have physiologically active chlorophyllous spores that are much more short-lived than the more common and dormant achlorophyllous spores. Most chlorophyllous-spored species (70%) are epiphytes and account for almost 37% of all epiphytic ferns. Chlorophyllous-spored ferns are also overrepresented among fern species in habitats with waterlogged soils, of which nearly 60% have chlorophyllous spores. Ferns in these disparate habitat types also have a low incidence of mycorrhizal associations. We therefore hypothesized that autotrophic chlorophyllous spores represent an adaptation of ferns to habitats with scarce mycorrhizal associations. METHODS: We evaluated the coevolution of chlorophyllous spores and mycorrhizal associations in ferns and their relation to habitat type using phylogenetic comparative methods. RESULTS: Although we did not find support for the coevolution of spore type and mycorrhizal associations, we did find that chlorophyllous spores and the absence of mycorrhizal associations have coevolved with epiphytic and waterlogged habitats. Transition rates to epiphytic and waterlogged habitats were significantly higher in species with chlorophyllous spores compared to achlorophyllous lineages. CONCLUSIONS: Spore type and mycorrhizal associations appear to play important roles in the radiation of ferns into different habitat types. Future work should focus on clarifying the functional significance of these associations.

Does the Size of Microgels Influence the Toughness of Microgel-Reinforced Hydrogels?

Rapid advances in the biomedical field increasingly often demand soft materials that can be processed into complex 3D shapes while being able to reliably bear significant loads. Granular hydrogels have the potential to serve as artificial tissues because they can be 3D printed into complex shapes and their composition can be tuned over short length scales. Unfortunately, granular hydrogels are typically soft such that they cannot be used for load-bearing applications. To address this shortcoming, individual microgels can be connected through a percolating network, such that they introduce the double network toughening mechanism into granular hydrogels. However, the influence of the microgel size and concentration on the processing and toughness of microgel-reinforced hydrogels (MRHs) remains to be elucidated. Here, it is demonstrated that processing and toughness depend on the inter-microgel connectivity, while the stress at break is solely dependent on the microgel size. These findings offer an in-depth understanding of how liquid- and paste-like precursors containing soft, deformable microgels can be processed into bulk microstructured soft materials and how the size and concentration of these microgels influence the mechanical properties of microgel-reinforced hydrogels.