Genome‐scale metabolic modeling reveals key features of a minimal gene set

Mesoplasma florum, a fast‐growing near‐minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein‐coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species‐specific constraints to produce the functional iJL208 genome‐scale model (GEM) of metabolism. Genome‐wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI‐syn3.0 and its parent JCVI‐syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping‐stone toward model‐driven whole‐genome engineering.     

Low base‐substitution mutation rate and predominance of insertion‐deletion events in the acidophilic bacterium Acidobacterium capsulatum

Analyses of spontaneous mutation have shown that total genome‐wide mutation rates are quantitatively similar for most prokaryotic organisms. However, this view is mainly based on organisms that grow best around neutral pH values (6.0–8.0). In particular, the whole‐genome mutation rate has not been determined for an acidophilic organism. Here, we have determined the genome‐wide rate of spontaneous mutation in the acidophilic Acidobacterium capsulatum using a direct and unbiased method: a mutation‐accumulation experiment followed by whole‐genome sequencing. Evaluation of 69 mutation accumulation lines of A. capsulatum after an average of ~2900 cell divisions yielded a base‐substitution mutation rate of 1.22 × 10⁻¹⁰ per site per generation or 4 × 10⁻⁴ per genome per generation, which is significantly lower than the consensus value (2.5−4.6 × 10⁻³) of mesothermophilic (~15–40°C) and neutrophilic (pH 6–8) prokaryotic organisms. However, the insertion‐deletion rate (0.43 × 10⁻¹⁰ per site per generation) is high relative to the base‐substitution mutation rate. Organisms with a similar effective population size and a similar expected effect of genetic drift should have similar mutation rates. Because selection operates on the total mutation rate, it is suggested that the relatively high insertion‐deletion rate may be balanced by a low base‐substitution rate in A. capsulatum, with selection operating on the total mutation rate.     

Engineering a genome‐reduced bacterium to eliminate Staphylococcus aureus biofilms in vivo

Bacteria present a promising delivery system for treating human diseases. Here, we engineered the genome‐reduced human lung pathogen Mycoplasma pneumoniae as a live biotherapeutic to treat biofilm‐associated bacterial infections. This strain has a unique genetic code, which hinders gene transfer to most other bacterial genera, and it lacks a cell wall, which allows it to express proteins that target peptidoglycans of pathogenic bacteria. We first determined that removal of the pathogenic factors fully attenuated the chassis strain in vivo. We then designed synthetic promoters and identified an endogenous peptide signal sequence that, when fused to heterologous proteins, promotes efficient secretion. Based on this, we equipped the chassis strain with a genetic platform designed to secrete antibiofilm and bactericidal enzymes, resulting in a strain capable of dissolving Staphylococcus aureus biofilms preformed on catheters in vitro, ex vivo, and in vivo. To our knowledge, this is the first engineered genome‐reduced bacterium that can fight against clinically relevant biofilm‐associated bacterial infections.     

Genome‐wide analysis reveals demographic and life‐history patterns associated with habitat modification in landlocked,deep‐spawning sockeye salmon (Oncorhynchus nerka)

Human‐mediated habitat fragmentation in freshwater ecosystems can negatively impact genetic diversity, demography, and life history of native biota, while disrupting the behavior of species that are dependent on spatial connectivity to complete their life cycles. In the Alouette River system (British Columbia, Canada), dam construction in 1928 impacted passage of anadromous sockeye salmon (Oncorhynchus nerka), with the last records of migrants occurring in the 1930s. Since that time, O. nerka persisted as a resident population in Alouette Reservoir until experimental water releases beginning in 2005 created conditions for migration; two years later, returning migrants were observed for the first time in ~70 years, raising important basic and applied questions regarding life‐history variation and population structure in this system. Here, we investigated the genetic distinctiveness and population history of Alouette Reservoir O. nerka using genome‐wide SNP data (n = 7,709 loci) collected for resident and migrant individuals, as well as for neighboring anadromous sockeye salmon and resident kokanee populations within the Fraser River drainage (n = 312 individuals). Bayesian clustering and principal components analyses based on neutral loci revealed five distinct clusters, largely associated with geography, and clearly demonstrated that Alouette Reservoir resident and migrant individuals are genetically distinct from other O. nerka populations in the Fraser River drainage. At a finer level, there was no clear evidence for differentiation between Alouette Reservoir residents and migrants; although we detected eight high‐confidence outlier loci, they all mapped to sex chromosomes suggesting that differences were likely due to uneven sex ratios rather than life history. Taken together, these data suggest that contemporary Alouette Reservoir O. nerka represents a landlocked sockeye salmon population, constituting the first reported instance of deep‐water spawning behavior associated with this life‐history form. This finding punctuates the need for reassessment of conservation status and supports ongoing fisheries management activities in Alouette Reservoir.     

Evolution of assortative mating following selective introgression of pigmentation genes between two Drosophila species

Adaptive introgression is ubiquitous in animals, but experimental support for its role in driving speciation remains scarce. In the absence of conscious selection, admixed laboratory strains of Drosophila asymmetrically and progressively lose alleles from one parental species and reproductive isolation against the predominant parent ceases after 10 generations. Here, we selectively introgressed during 1 year light pigmentation genes of D. santomea into the genome of its dark sibling D. yakuba, and vice versa. We found that the pace of phenotypic change differed between the species and the sexes and identified through genome sequencing common as well as distinct introgressed loci in each species. Mating assays showed that assortative mating between introgressed flies and both parental species persisted even after 4 years (~60 generations) from the end of the selection. Those results indicate that selective introgression of as low as 0.5% of the genome can beget morphologically distinct and reproductively isolated strains, two prerequisites for the delimitation of new species. Our findings hence represent a significant step toward understanding the genome‐wide dynamics of speciation‐through‐introgression.     

Roosting ecology of endangered plant‐roosting bats on Okinawa Island: Implications for bat‐friendly forestry practices

Roosting information is crucial to guiding bat conservation and bat‐friendly forestry practices. The Ryukyu tube‐nosed bat Murina ryukyuana (Endangered) and Yanbaru whiskered bat Myotis yanbarensis (Critically Endangered) are forest‐dwelling bats endemic to the central Ryukyu Archipelago, Japan. Despite their threatened status, little is known about the roosting ecology of these species and the characteristics of natural maternity roosts are unknown. To inform sustainable forestry practices and conservation management, we radio‐tracked day roosts of both species in the subtropical forests of Okinawa''s Kunigami Village District. We compared roost and roost site characteristics statistically between M. ryukyuana nonmaternity roosts (males or nonreproductive females), maternity roosts, and all M. yanbarensis roosts. Generalized linear models were used to investigate roost site selection by M. ryukyuana irrespective of sex and age class. Lastly, we compiled data on phenology from this and prior studies. Nonreproductive M. ryukyuana roosted alone and primarily in understory foliage. Murina ryukyuana maternity roosts were limited to stands >50 years old, and ~60% were in foliage. Myotis yanbarensis roosted almost entirely in cavities along gulch bottoms and only in stands >70 years old (~1/3 of Kunigami''s total forest area). Murina ryukyuana maternity roosts were higher (4.3 ± 0.6 m) than conspecific nonmaternity roosts (2.3 ± 0.5 m; p < .001) and M. yanbarensis roosts (2.7 ± 0.5 m; not significant). Model results were inconclusive. Both species appear to be obligate plant roosters throughout their life cycle, but the less flexible roosting preferences of M. yanbarensis may explain its striking rarity. To conserve these threatened bats, we recommend the following forestry practices: (a) reduce clearing of understory vegetation, (b) refrain from removing trees along streams, (c) promote greater tree cavity densities by protecting old‐growth forests and retaining snags, and (d) refrain from removing trees or understory between April and July, while bats are pupping.     

Using fractal self‐similarity to increase precision of shrub biomass estimates

We show that aerial tips are self‐similar fractals of whole shrubs and present a field method that applies this fact to improves accuracy and precision of biomass estimates of tall‐shrubs, defined here as those with diameter at root collar (DRC) ≥ 2.5 cm. Power function allometry of biomass to stem diameter generates a disproportionate prediction error that increases rapidly with diameter. Thus, biomass should be modeled as a single measure of stem diameter only if stem diameter is less than a threshold D_max. When stem diameter exceeds D_max, then the stem internode should be treated as a conic frustrum requiring two additional measures: a second, node‐adjacent diameter and a length. If the second diameter is less than D_max, then the power function allometry can be applied to the aerial tip; otherwise an additional internode is measured. This “two‐component” allometry—internodes as frustra and aerial tips as shrubs—can reduce estimated biomass error propagated to the plot‐level by as much as 50% or more where very large shrubs are present D_max is any diameter such that the ratio of single‐component to two‐component uncertainty exceeds the ratio of two‐component to single‐component measurement time. Guidelines for estimating D_max based on pilot field data are provided. Tall shrubs are increasing in abundance and distribution across Arctic, alpine, boreal, and dryland ecosystems. Estimating their biomass is important for both ecological studies and carbon accounting. Reducing field‐sample prediction error increases precision in multi‐stage modeling because additional measures efficiently improve plot‐level biomass precision, reducing uncertainty for shrub biomass estimates.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

Identification of novel susceptibility loci for non‐syndromic cleft lip with or without cleft palate

Although several genome‐wide association studies (GWAS) of non‐syndromic cleft lip with or without cleft palate (NSCL/P) have been reported, more novel association signals are remained to be exploited. Here, we performed an in‐depth analysis of our previously published Chinese GWAS cohort study with replication in an extra dbGaP case‐parent trios and another in‐house Nanjing cohort, and finally identified five novel significant association signals (rs11119445: 3’ of SERTAD4, P = 6.44 × 10⁻¹⁴; rs227227 and rs12561877: intron of SYT14, P = 5.02 × 10⁻¹³ and 2.80 × 10⁻¹¹, respectively; rs643118: intron of TRAF3IP3, P = 4.45 × 10⁻⁶; rs2095293: intron of NR6A1, P = 2.98 × 10⁻⁵). The mean (standard deviation) of the weighted genetic risk score (wGRS) from these SNPs was 1.83 (0.65) for NSCL/P cases and 1.58 (0.68) for controls, respectively (P = 2.67 × 10⁻¹⁶). Rs643118 was identified as a shared susceptible factor of NSCL/P among Asians and Europeans, while rs227227 may contribute to the risk of NSCL/P as well as NSCPO. In addition, sertad4 knockdown zebrafish models resulted in down‐regulation of sox2 and caused oedema around the heart and mandibular deficiency, compared with control embryos. Taken together, this study has improved our understanding of the genetic susceptibility to NSCL/P and provided further clues to its aetiology in the Chinese population.     

Low genome‐wide divergence between two lizard populations with high adaptive phenotypic differentiation

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome‐wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with F_ST statistics significantly higher than those expected on the basis of overall, genome‐wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome‐wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with F_ST‐based outlier detection tests. Instead, our analyses confirmed the lack of genome‐wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.     

PIP2 depletion and altered endocytosis caused by expression of Alzheimer's disease‐protective variant PLCγ2 R522

Variants identified in genome‐wide association studies have implicated immune pathways in the development of Alzheimer’s disease (AD). Here, we investigated the mechanistic basis for protection from AD associated with PLCγ2 R522, a rare coding variant of the PLCG2 gene. We studied the variant''s role in macrophages and microglia of newly generated PLCG2‐R522‐expressing human induced pluripotent cell lines (hiPSC) and knockin mice, which exhibit normal endogenous PLCG2 expression. In all models, cells expressing the R522 mutation show a consistent non‐redundant hyperfunctionality in the context of normal expression of other PLC isoforms. This manifests as enhanced release of cellular calcium ion stores in response to physiologically relevant stimuli like Fc‐receptor ligation or exposure to Aβ oligomers. Expression of the PLCγ2‐R522 variant resulted in increased stimulus‐dependent PIP₂ depletion and reduced basal PIP₂ levels in vivo. Furthermore, it was associated with impaired phagocytosis and enhanced endocytosis. PLCγ2 acts downstream of other AD‐related factors, such as TREM2 and CSF1R, and alterations in its activity directly impact cell function. The inherent druggability of enzymes such as PLCγ2 raises the prospect of PLCγ2 manipulation as a future therapeutic approach in AD.     

Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

Low‐coverage reduced representation sequencing reveals subtle within‐island genetic structure in Aldabra giant tortoises

Aldabrachelys gigantea (Aldabra giant tortoise) is one of only two giant tortoise species left in the world and survives as a single wild population of over 100,000 individuals on Aldabra Atoll, Seychelles. Despite this large current population size, the species faces an uncertain future because of its extremely restricted distribution range and high vulnerability to the projected consequences of climate change. Captive‐bred A. gigantea are increasingly used in rewilding programs across the region, where they are introduced to replace extinct giant tortoises in an attempt to functionally resurrect degraded island ecosystems. However, there has been little consideration of the current levels of genetic variation and differentiation within and among the islands on Aldabra. As previous microsatellite studies were inconclusive, we combined low‐coverage and double‐digest restriction‐associated DNA (ddRAD) sequencing to analyze samples from 33 tortoises (11 from each main island). Using 5426 variant sites within the tortoise genome, we detected patterns of within‐island population structure, but no differentiation between the islands. These unexpected results highlight the importance of using genome‐wide genetic markers to capture higher‐resolution genetic structure to inform future management plans, even in a seemingly panmictic population. We show that low‐coverage ddRAD sequencing provides an affordable alternative approach to conservation genomic projects of non‐model species with large genomes.