Copy number variation of the putative speciation genes in the European house mouse hybrid zone

Hybrid zones have long been described as "windows on the evolutionary process", and studying them has become even more important since the advance in the genome analysis tools. The hybrid zone between two subspecies of the house mouse (Mus musculus musculus and Mus m. domesticus) is a unique model speciation system to study fine scale interactions of recently diverged genomes. Here, we explore the role of gene Copy Number Variation in shaping the barrier to introgression in the hybrid zone within a previously established transect in Central Europe. The CNV of seven pre-selected candidate genes was determined via droplet-digital PCR and analyzed in the context of ~500k SNPs, with the ancestral population (i.e. musculus or domesticus) of every SNP allele previously inferred in the admixed individuals (Baird et al., in prep.). The copy numbers of five genes were clearly associated with the prevalence of either musculus or domesticus genomes across the hybrid zone. In three cases, the highest and/or outlying levels of association were observed at or very close to the annotated positions of the respective gene amplicons, demonstrating the power of our approach in confirming the reference locations of copy number variants. Notably, several other reference locations were recognized as positive outliers in the association with particular CNV genes, possibly representing the extra gene copies and/or their epistatic interaction sites.

Chloroplast and Nuclear Genetic Diversity Explain the Limited Distribution of Endangered and Endemic Thuja sutchuenensis in China

Narrow-ranged species face challenges from natural disasters and human activities, and to address why species distributes only in a limited region is of great significance. Here we investigated the genetic diversity, gene flow, and genetic differentiation in six wild and three cultivated populations of Thuja sutchuenensis, a species that survive only in the Daba mountain chain, using chloroplast simple sequence repeats (cpSSR) and nuclear restriction site-associated DNA sequencing (nRAD-seq). Wild T. sutchuenensis populations were from a common ancestral population at 203 ka, indicating they reached the Daba mountain chain before the start of population contraction at the Last Interglacial (LIG, ∼120–140 ka). T. sutchuenensis populations showed relatively high chloroplast but low nuclear genetic diversity. The genetic differentiation of nRAD-seq in any pairwise comparisons were low, while the cpSSR genetic differentiation values varied with pairwise comparisons of populations. High gene flow and low genetic differentiation resulted in a weak isolation-by-distance effect. The genetic diversity and differentiation of T. sutchuenensis explained its survival in the Daba mountain chain, while its narrow ecological niche from the relatively isolated and unique environment in the “refugia” limited its distribution.

Genetic Contributors of Incident Stroke in 10,700 African Americans With Hypertension: A Meta-Analysis From the Genetics of Hypertension Associated Treatments and Reasons for Geographic and Racial Differences in Stroke Studies

Background: African Americans (AAs) suffer a higher stroke burden due to hypertension. Identifying genetic contributors to stroke among AAs with hypertension is critical to understanding the genetic basis of the disease, as well as detecting at-risk individuals.Methods: In a population comprising over 10,700 AAs treated for hypertension from the Genetics of Hypertension Associated Treatments (GenHAT) and Reasons for Geographic and Racial Differences in Stroke (REGARDS) studies, we performed an inverse variance-weighted meta-analysis of incident stroke. Additionally, we tested the predictive accuracy of a polygenic risk score (PRS) derived from a European ancestral population in both GenHAT and REGARDS AAs aiming to evaluate cross-ethnic performance.Results: We identified 10 statistically significant (p < 5.00E-08) and 90 additional suggestive (p < 1.00E-06) variants associated with incident stroke in the meta-analysis. Six of the top 10 variants were located in an intergenic region on chromosome 18 (LINC01443-LOC644669). Additional variants of interest were located in or near the COL12A1, SNTG1, PCDH7, TMTC1, and NTM genes. Replication was conducted in the Warfarin Pharmacogenomics Cohort (WPC), and while none of the variants were directly validated, seven intronic variants of NTM proximal to our target variants, had a p-value <5.00E-04 in the WPC. The inclusion of the PRS did not improve the prediction accuracy compared to a reference model adjusting for age, sex, and genetic ancestry in either study and had lower predictive accuracy compared to models accounting for established stroke risk factors. These results demonstrate the necessity for PRS derivation in AAs, particularly for diseases that affect AAs disproportionately.Conclusion: This study highlights biologically plausible genetic determinants for incident stroke in hypertensive AAs. Ultimately, a better understanding of genetic risk factors for stroke in AAs may give new insight into stroke burden and potential clinical tools for those among the highest at risk.

Lower promoter activity of the ST8SIA2 gene has been favored in evolving human collective brains

ST8SIA2 is an important molecule regulating expression of the phenotype involved in schizophrenia. Lowered promoter activity of the ST8SIA2 gene is considered to be protective against schizophrenia by conferring tolerance to psychosocial stress. Here, we examined the promoter-type composition of anatomically modern humans (AMHs) and archaic humans (AHs; Neanderthals and Denisovans), and compared the promoter activity at the population level (population promoter activity; PPA) between them. In AMHs, the TCT-type, showing the second lowest promoter activity, was most prevalent in the ancestral population of non-Africans. However, the detection of only the CGT-type from AH samples and recombination tracts in AH sequences showed that the CGT- and TGT-types, exhibiting the two highest promoter activities, were common in AH populations. Furthermore, interspecies gene flow occurred into AMHs from AHs and into Denisovans from Neanderthals, influencing promoter-type compositions independently in both AMHs and AHs. The difference of promoter-type composition makes PPA unique in each population. East and Southeast Asian populations show the lowest PPA. This results from the selective increase of the CGC-type, showing the lowest promoter activity, in these populations. Every non-African population shows significantly lower PPA than African populations, resulting from the TCT-type having the highest prevalence in the ancestral population of non-Africans. In addition, PPA reduction is also found among subpopulations within Africa via a slight increase of the TCT-type. These findings indicate a trend toward lower PPA in the spread of AMHs, interpreted as a continuous adaptation to psychosocial stress arising in migration. This trend is considered as genetic tuning for the evolution of collective brains. The inferred promoter-type composition of AHs differed markedly from that of AMHs, resulting in higher PPA in AHs than in AMHs. This suggests that the trend toward lower PPA is a unique feature in AMH spread.

Theoretical Analysis of Principal Components in an Umbrella Model of Intraspecific Evolution

Principal component analysis (PCA) is one of the most frequently-used approach to describe population structure from multilocus genotype data. Regarding geographic range expansions of modern humans, interpretations of PCA have, however, been questioned, as there is uncertainty about the wave-like patterns that have been observed in principal components. It has indeed been argued that wave-like patterns are mathematical artifacts that arise generally when PCA is applied to data in which genetic differentiation increases with geographic distance. Here, we present an alternative theory for the observation of wave-like patterns in PCA. We study a coalescent model -- the umbrella model -- for the diffusion of genetic variants. The model is based on a hierarchy of splits from an ancestral population without any particular geographical structure. In the umbrella model, splits occur almost continuously in time, giving birth to small daughter populations at a regular pace. Our results provide detailed mathematical descriptions of eigenvalues and eigenvectors for the PCA of sampled genomic sequences under the model. Removing variants uniquely represented in the sample, the PCA eigenvectors are defined as cosine functions of increasing periodicity, reproducing wave-like patterns observed in equilibrium isolation-by-distance models. Including rare variants in the analysis, the eigenvectors corresponding to the largest eigenvalues exhibit complex wave shapes. The accuracy of our predictions is further investigated with coalescent simulations. Our analysis supports the hypothesis that highly structured wave-like patterns could arise from genetic drift only, and may not always be artificial outcomes of spatially structured data. Genomic data related to the peopling of the Americas are reanalyzed in the light of our new theory.

Genomic insights into the recent population history of Mapuche Native Americans

The last few years have witnessed an explosive generation of genomic data from ancient and modern Native American populations. These data shed light on key demographic shifts that occurred in geographically diverse territories of South America, such as the Andean highlands, Southern Patagonia and the Amazon basin. We used genomic data to study the recent population history of the Mapuche, who are the major Native population from the Southern Cone (Chile and Argentina). We found evidence of specific shared genetic ancestry between the Mapuche and ancient populations from Southern Patagonia, Central Chile and the Argentine Pampas. Despite previous evidence of cultural influence of Inca and Tiwanaku polities over the Mapuche, we did not find evidence of specific shared ancestry between them, nor with Amazonian groups. We estimated the effective population size dynamics of the Mapuche ancestral population during the last millennia, identifying a population bottleneck around 1650 AD, coinciding with a period of Spaniards invasions into the territory inhabited by the Mapuche. Finally, we show that admixed Chileans underwent post-admixture adaptation in their Mapuche subancestry component in genes related with lipid metabolism, suggesting adaptation to scarce food availability.

Host Species is Linked to Pathogen Genotype for the Amphibian Chytrid Fungus (Batrachochytrium dendrobatidis)

Host-pathogen specificity can arise from certain selective environments mediated by both the host and pathogen. Therefore, understanding the degree to which host species identity is correlated with pathogen genotype can help reveal historical host-pathogen dynamics. One animal disease of particular concern is chytridiomycosis, typically caused by the global panzootic lineage of the amphibian chytrid fungus ( Batrachochytrium dendrobatidis , Bd), termed the Bd-GPL. This pathogen lineage has caused devastating declines in amphibian communities around the world. However, the origin of Bd-GPL and the fine-scale transmission dynamics of this lineage have remained a mystery. This is especially the case in North America where Bd-GPL is widespread, but disease outbreaks occur sporadically. Herein, we use Bd genetic data collected throughout the United States from amphibian skin swab and cultured isolate samples to investigate Bd genetic patterns. We highlight two case studies in Pennsylvania and Nevada where Bd-GPL genotypes are strongly correlated with host species identity. Specifically, in some localities bullfrogs ( Rana catesbeiana ) are infected with Bd-GPL lineages that are distinct from those infecting other sympatric amphibian species. Overall, we reveal a previously unknown association of Bd genotype with host species and identify the eastern United States as a Bd diversity hotspot and potential ancestral population for Bd-GPL.

Genome Assembly and Population Resequencing Reveal the Geographical Divergence of 'Shanmei'(Rubus corchorifolius)

Rubus corchorifolius (Shanmei or mountain berry, 2n =14) is widely distributed in China, and its fruit has high nutritional and medicinal values. Here, we report a high-quality chromosome-scale genome assembly of Shanmei, with a size of 215.69 Mb and encompassing 26696 genes. Genome comparisons among Rosaceae species show that Shanmei and Fupenzi(Rubus chingii Hu) are most closely related, and then is blackberry (Rubus occidentalis). Further resequencing of 101 samples of Shanmei collected from four regions in provinces of Yunnan, Hunan, Jiangxi and Sichuan in South China reveals that the Hunan population of Shanmei possesses the highest diversity and may represent the relatively more ancestral population. Moreover, the Yunnan population undergoes strong selection based on nucleotide diversity, linkage disequilibrium and the historical effective population size analyses. Furthermore, genes from candidate genomic regions that show strong divergence are significantly enriched in flavonoid biosynthesis and plant hormone signal transduction, indicating the genetic basis of adaptation of Shanmei to the local environments. The high-quality genome sequences and the variome dataset of Shanmei provide valuable resources for breeding applications and for elucidating the genome evolution and ecological adaptation of Rubus species.

Phylogeography of the desert scorpion Mesobuthus mongolicus illuminates a route out of Central Asia and Junggar Basin as a corridor for faunal exchange between Central and East Asia

A comprehensive understanding of the ecological, evolutionary and genetic consequences of climate changes requires integration of information from different geographic regions and ecosystems. However, a clear knowledge gap exists in arid biota of Asian drylands, where climate change and biological evolution demonstrate high regional specificity. We report here our phylogeographic study on the desert scorpion Mesobuthus mongolicus, attempting to narrow the gap. Range wide samples for one mitochondrial gene and three nuclear loci were subjected to phylogenetic inferences, molecular dating, Bayesian phylogeographic and demographic reconstructions, and population genetic analyses. Ecological niche modeling was employed to generate models of potential distribution range at present and during the last glacial maximum (LGM). Our results revealed that M. mongolicus represents a single coherent lineage that diverged with its most closely related lineage from Central Asia about 1.36 Ma and underwent radiation ever since. Bayesian phylogeographic reconstruction suggested the ancestral population of M. mongolicus dispersed gradually eastward through the Junngar Basin to the Gobi region during the late Pleistocene, strengthening the view that the Junggar Basin has constituted an important corridor for faunal exchange between Central Asian and East Asia. We found that the desert scorpion underwent demographic expansion while experiencing range contraction during the LGM. Development of sandy deserts in northwest China incurred by Eurasian aridification might have opened up substantial empty niches that sustained population expansion during the LGM. Our results suggested that organisms' ecological adaptations have acted as an important determinant for their phylogeographic and demographic responses to past climate changes.

Variation in nuclear genome size within the Eisenia nordenskioldi complex (Lumbricidae, Annelida)

The size of the nuclear genome in eukaryotes is mostly determined by mobile elements and noncoding sequences and may vary within wide limits. It can differ significantly both among higher-order taxa and closely related species within a genus; genome size is known to be uncorrelated with organism complexity (the so-called C-paradox). Less is known about intraspecific variation of this parameter. Typically, genome size is stable within a species, and the known exceptions turn out be cryptic taxa. The Eisenia nordenskioldi complex encompasses several closely related earthworm species. They are widely distributed in the Urals, Siberia, and the Russian Far East, as well as adjacent regions. This complex is characterized by significant morphological, chromosomal, ecological, and genetic variation. The aim of our study was to estimate the nuclear genome size in several genetic lineages of the E. nordenskioldi complex using flow cytometry. The genome size in different genetic lineages differed strongly, which supports the hypothesis that they are separate species. We found two groups of lineages, with small (250–500 Mbp) and large (2300–3500 Mbp) genomes. Moreover, different populations within one lineage also demonstrated variation in genome size (15–25 %). We compared the obtained data to phylogenetic trees based on transcriptome data. Genome size in ancestral population was more likely to be big. It increased or decreased independently in different lineages, and these processes could be associated with changes in genome size and/or transition to endogeic lifestyle.