Liuyang Wang

History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae)

How Quaternary climatic oscillations affected range distributions and intraspecific divergence of alpine plants on the Qinghai‐Tibetan Plateau (QTP) remains largely unknown. Here, we report a survey of chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variation aimed at exploring the phylogeographical history of the QTP alpine endemic Aconitum gymnandrum. We sequenced three cpDNA fragments (rpl20–rps12 intergenic spacer, the trnV intron and psbA‐trnH spacer) and also the nuclear (ITS) region in 245 individuals from 23 populations sampled throughout the species’ range. Two distinct lineages, with eastern and western geographical distributions respectively, were identified from a phylogenetic analysis of ITS sequence variation. Based on a fast substitution rate, these were estimated to have diverged from each other in the early Pleistocene approximately 1.45 Ma. The analysis of cpDNA variation identified nine chlorotypes that clustered into two major clades that were broadly congruent in geographical distribution with the two ITS lineages. The east–west split of cpDNA divergence was supported by an amova which partitioned approximately half of the total variance between these two groups of populations. Analysis of the spatial distribution of chlorotypes showed that each clade was subdivided into two groups of populations such that a total of four population groups existed in the species. It is suggested that these different groups derive from four independent glacial refugia that existed during the Last Glacial Maximum (LGM), and that three of these refugia were located at high altitude on the QTP platform itself at that time. Coalescent simulation of chlorotype genealogies supported both an early Pleistocene origin of the two main cpDNA clades and also the ‘four‐refugia’ hypothesis during the LGM. Two previous phylogeographical studies of QTP alpine plants indicated that such plants retreated to refugia at the eastern/south‐eastern plateau edge during the LGM and/or previous glacial maxima. However, the results for A. gymnandrum suggest that at least some of these cold‐tolerant species may have also survived centrally on the QTP platform throughout the Quaternary.

Genome-wide architecture of reproductive isolation in a naturally occurring hybrid zone between Mus musculus musculus and M. m. domesticus

Studies of a hybrid zone between two house mouse subspecies (Mus musculus musculus and M. m. domesticus) along with studies using laboratory crosses reveal a large role for the X chromosome and multiple autosomal regions in reproductive isolation as a consequence of disrupted epistasis in hybrids. One limitation of previous work has been that most of the identified genomic regions have been large. The goal here is to detect and characterize precise genomic regions underlying reproductive isolation. We surveyed 1401 markers evenly spaced across the genome in 679 mice collected from two different transects. Comparisons between transects provide a means for identifying common patterns that likely reflect intrinsic incompatibilities. We used a genomic cline approach to identify patterns that correspond to epistasis. From both transects, we identified contiguous regions on the X chromosome in which markers were inferred to be involved in epistatic interactions. We then searched for autosomal regions showing the same patterns and found they constitute about 5% of autosomal markers. We discovered substantial overlap between these candidate regions underlying reproductive isolation and QTL for hybrid sterility identified in laboratory crosses. Analysis of gene content in these regions suggests a key role for several mechanisms, including the regulation of transcription, sexual conflict and sexual selection operating at both the postmating prezygotic and postzygotic stages of reproductive isolation. Taken together, these results indicate that speciation in two recently diverged (c. 0.5 Ma) house mouse subspecies is complex, involving many genes dispersed throughout the genome and associated with distinct functions.

HIV-1 gp120 Vaccine Induces Affinity Maturation in both New and Persistent Antibody Clonal Lineages

ABSTRACT Most antibodies that broadly neutralize HIV-1 are highly somatically mutated in antibody clonal lineages that persist over time. Here, we describe the analysis of human antibodies induced during an HIV-1 vaccine trial (GSK PRO HIV-002) that used the clade B envelope (Env) gp120 of clone W6.1D (gp120W6.1D). Using dual-color antigen-specific sorting, we isolated Env-specific human monoclonal antibodies (MAbs) and studied the clonal persistence of antibodies in the setting of HIV-1 Env vaccination. We found evidence of VH somatic mutation induced by the vaccine but only to a modest level (3.8% ± 0.5%; range 0 to 8.2%). Analysis of 34 HIV-1-reactive MAbs recovered over four immunizations revealed evidence of both sequential recruitment of naïve B cells and restimulation of previously recruited memory B cells. These recombinant antibodies recapitulated the anti-HIV-1 activity of participant serum including pseudovirus neutralization and antibody-dependent cell-mediated cytotoxicity (ADCC). One antibody (3491) demonstrated a change in specificity following somatic mutation with binding of the inferred unmutated ancestor to a linear C2 peptide while the mutated antibody reacted only with a conformational epitope in gp120 Env. Thus, gp120W6.1D was strongly immunogenic but over four immunizations induced levels of affinity maturation below that of broadly neutralizing MAbs. Improved vaccination strategies will be needed to drive persistent stimulation of antibody clonal lineages to induce affinity maturation that results in highly mutated HIV-1 Env-reactive antibodies.

Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone

Theory predicts that naturally occurring hybrid zones between genetically distinct taxa can move over space and time as a result of selection and/or demographic processes, with certain types of hybrid zones being more or less likely to move. Determining whether a hybrid zone is stationary or moving has important implications for understanding evolutionary processes affecting interactions in hybrid populations. However, direct observations of hybrid zone movement are difficult to make unless the zone is moving rapidly. Here, evidence for movement in the house mouse Mus musculus domesticus × Mus musculus musculus hybrid zone is provided using measures of LD and haplotype structure among neighbouring SNP markers from across the genome. Local populations of mice across two transects in Germany and the Czech Republic were sampled, and a total of 1301 mice were genotyped at 1401 markers from the nuclear genome. Empirical measures of LD provide evidence for extinction and (re)colonization in single populations and, together with simulations, suggest hybrid zone movement because of either geography‐dependent asymmetrical dispersal or selection favouring one subspecies over the other.

Functional Organization of the Genome May Shape the Species Boundary in the House Mouse

Genomic features such as rate of recombination and differentiation have been suggested to play a role in species divergence. However, the relationship of these phenomena to functional organization of the genome in the context of reproductive isolation remains unexplored. Here, we examine genomic characteristics of the species boundaries between two house mouse subspecies (Mus musculus musculus/M. m. domesticus). These taxa form a narrow semipermeable zone of secondary contact across Central Europe. Due to the incomplete nature of reproductive isolation, gene flow in the zone varies across the genome. We present an analysis of genomic differentiation, rate of recombination, and functional composition of genes relative to varying amounts of introgression. We assessed introgression using 1,316 autosomal single nucleotide polymorphism markers, previously genotyped in hybrid populations from three transects. We found a significant relationship between amounts of introgression and both genomic differentiation and rate of recombination with genomic regions of reduced introgression associated with higher genomic differentiation and lower rates of recombination, and the opposite for genomic regions of extensive introgression. We also found a striking functional polarization of genes based on where they are expressed in the cell. Regions of elevated introgression exhibit a disproportionate number of genes involved in signal transduction functioning at the cell periphery, among which olfactory receptor genes were found to be the most prominent group. Conversely, genes expressed intracellularly and involved in DNA binding were the most prevalent in regions of reduced introgression. We hypothesize that functional organization of the genome is an important driver of species divergence.

Comparative analyses of clinical and environmental populations of Cryptococcus neoformans in Botswana

Cryptococcus neoformans var. grubii (Cng) is the most common cause of fungal meningitis, and its prevalence is highest in sub‐Saharan Africa. Patients become infected by inhaling airborne spores or desiccated yeast cells from the environment, where the fungus thrives in avian droppings, trees and soil. To investigate the prevalence and population structure of Cng in southern Africa, we analysed isolates from 77 environmental samples and 64 patients. We detected significant genetic diversity among isolates and strong evidence of geographic structure at the local level. High proportions of isolates with the rare MATa allele were observed in both clinical and environmental isolates; however, the mating‐type alleles were unevenly distributed among different subpopulations. Nearly equal proportions of the MATa and MATα mating types were observed among all clinical isolates and in one environmental subpopulation from the eastern part of Botswana. As previously reported, there was evidence of both clonality and recombination in different geographic areas. These results provide a foundation for subsequent genomewide association studies to identify genes and genotypes linked to pathogenicity in humans.

Human genetic and metabolite variation reveals that methylthioadenosine is a prognostic biomarker and an inflammatory regulator in sepsis.

The substrate of the methionine salvage pathway, methylthioadenosine, is a predictor of death in sepsis patients. Sepsis is a deleterious inflammatory response to infection with high mortality. Reliable sepsis biomarkers could improve diagnosis, prognosis, and treatment. Integration of human genetics, patient metabolite and cytokine measurements, and testing in a mouse model demonstrate that the methionine salvage pathway is a regulator of sepsis that can accurately predict prognosis in patients. Pathway-based genome-wide association analysis of nontyphoidal Salmonella bacteremia showed a strong enrichment for single-nucleotide polymorphisms near the components of the methionine salvage pathway. Measurement of the pathway’s substrate, methylthioadenosine (MTA), in two cohorts of sepsis patients demonstrated increased plasma MTA in nonsurvivors. Plasma MTA was correlated with levels of inflammatory cytokines, indicating that elevated MTA marks a subset of patients with excessive inflammation. A machine-learning model combining MTA and other variables yielded approximately 80% accuracy (area under the curve) in predicting death. Furthermore, mice infected with Salmonella had prolonged survival when MTA was administered before infection, suggesting that manipulating MTA levels could regulate the severity of the inflammatory response. Our results demonstrate how combining genetic data, biomolecule measurements, and animal models can shape our understanding of disease and lead to new biomarkers for patient stratification and potential therapeutic targeting.

Refugial isolation and range expansions drive the genetic structure of Oxyria sinensis (Polygonaceae) in the Himalaya-Hengduan Mountains.

The formation of the Mekong-Salween Divide and climatic oscillations in Pleistocene were the main drivers for the contemporary diversity and genetic structure of plants in the Himalaya-Hengduan Mountains (HHM). To identify the relative roles of the two historical events in shaping population history of plants in HHM, we investigated the phylogeographic pattern of Oxyria sinensis, a perennial plant endemic to the HHM. Sixteen chloroplast haplotypes were identified and were clustered into three phylogenetic clades. The age of the major clades was estimated to be in the Pleistocene, falling into several Pleistocene glacial stages and postdating the formation of the Mekong-Salween Divide. Range expansions occurred at least twice in the early and middle Pleistocene, but the spatial genetic distribution rarely changed since the Last Glacial Maximum. Our results suggest that temporary mountain glaciers may act as barriers in promoting the lineage divergence in O. sinensis and that subsequential range expansions and secondary contacts might reshape the genetic distribution in geography and blur the boundary of population differentiation created in the earlier glacial stages. This study demonstrates that Pleistocene climatic change and mountain glaciers, rather than the Mekong-Salween Divide, play the primary role in shaping the spatial genetic structure of O. sinensis.

Confirmation of natural hybrids between Gentiana straminea and G. siphonantha (Gentianaceae) based on molecular evidence

A few individuals with intermediate morphology always appeared in the sympatric distributions of Gentiana straminea and G. siphonantha. These intermediate individuals were hypothesized to be the hybrids of two species after a careful evaluation of their morphological characteristics. To test this hypothesis, sequence comparison of the internal transcribed spacer (ITS) regions of the nuclear ribosomal and trnS (GCU)-trnG (UCC) intergenic spacer region of the chloroplast DNA from Gentiana straminea, G. siphonantha and the putative hybrids was performed. The results suggest that most intermediate individuals were the natural hybrids between G. straminea and G. siphonantha. In addition, we examined the sequence variation among the individuals of both parent species and analyzed the possibility leading to the incongruent identification in some individuals based on morphologic and molecular evidences, respectively. The intraspecific diversification of DNA fragments within both parent species and their high variability in hybrid swarms probably resulted from chloroplast genome recombination and incomplete lineage sorting during the early stages of speciation origin of the parent species.

minotaur: A platform for the analysis and visualization of multivariate results from genome scans with R Shiny.

Genome scans are widely used to identify ‘outliers’ in genomic data: loci with different patterns compared with the rest of the genome due to the action of selection or other nonadaptive forces of evolution. These genomic data sets are often high dimensional, with complex correlation structures among variables, making it a challenge to identify outliers in a robust way. The Mahalanobis distance has been widely used, but has the major limitation of assuming that data follow a simple parametric distribution. Here, we develop three new metrics that can be used to identify outliers in multivariate space, while making no strong assumptions about the distribution of the data. These metrics are implemented in the R package minotaur, which also includes an interactive web‐based application for visualizing outliers in high‐dimensional data sets. We illustrate how these metrics can be used to identify outliers from simulated genetic data and discuss some of the limitations they may face in application.