Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.     

Rapid assessment of single-copy nuclear DNA variation in diverse species

We investigated the use of PCR primers designed to conserved exons within nuclear DNA to amplify potentially variable regions such as introns or hypervariable exons from a wide range of species. We then explored various approaches to assay population-level variation in these PCR products. Primers designed to amplify regions within the histone H2AF, myoglobin , MHC DQA , and aldolase (ALD) genes gave clean amplifications in diverse mammals (DQA) , and in birds, reptiles and mammals ( aldolase, H2AF, myoglobin ). The sequenced PCR products generally, but not always, confirmed that the correct locus had been amplified. Several primer sets produced smaller size fragments consistent with preferential amplification of intronless pseudogenes; this was confirmed by sequencing seal and reptile H2AF PCR products. Digestion with randomly selected four-base recognizing enzymes detected variation in some cases but not in others. In species/gene combinations with either low (e.g. seal H2AF, ALD-A ) or high (e.g. skink ALD-1 ) nucleotide diversity it was more efficient to sequence a small number of distantly related individuals (e.g. one per geographic population) and from these data to identify informative or potentially informative restriction enzymes for 'targeted' digestion. We conclude that for studies of population-level variation, the optimal approach is to use a battery of primers for initial PCR of both mtDNA and scnDNA loci, select those that give clean amplifications, and sequence one sample from each population to (i) confirm gene identity, (ii) estimate the amount of variation and, (iii) search for diagnostic restriction sites. This will allow determination of the most efficient approach for a large-scale study.     

Using molecular and ecological data to diagnose endangered populations of the puritan tiger beetle Cicindela puritana

Populations of the puritan tiger beetle Cicindela puritana in the eastern United States were found to be highly threatened at the Connecticut River, whereas several large populations on the western shore and newly discovered populations on the eastern shore of the Chesapeake Bay appeared to be less endangered. We assessed if the disjunct C. puritana subgroups are genetically distinct and therefore should be treated as separate units for conservation purposes. A total of 13 individuals from the Connecticut River and 27 individuals from the Chesapeake Bay were each analysed by sequencing of up to 837 base pairs of mitochondrial DNA per individual. Five different haplotypes could be distinguished. In a phylogenetic analysis of these DNA sequences that included four related Cicindela species as out-groups, haplotypes from the Chesapeake Bay represent a distinct clade. The conservation status of these populations was evaluated using a phylogenetic approach based on cladistic analysis and the framework of the phylogenetic species concept. According to this analysis, beetles from the Connecticut River and the Chesapeake Bay have to be considered as independent units. Populations from the eastern and western shore of Chesapeake Bay are not split in more than one unit using the same criteria, although they exhibited some degree of genetic subdivision. The results from the mtDNA analysis were corroborated by ecological parameters in that the Chesapeake Bay populations can be distinguished from all congeners by their different tat association.     

Does hunting affect the demography and genetic structure of the greywing francolin Francolinus africanus?

Hunted and unhunted populations of greywing francolin Francolinus africanus have been studied in the eastern Cape Province of South Africa in order to understand the effects of hunting on the demography and genetic structure of these populations. Greywing population density cycled annually for both hunted and unhunted populations. However, there was an apparent pulse of immigration of sub-dominant birds, and earlier reproduction, in the hunted populations immediately after the winter hunting season. Average levels of allozyme heterozygosity (H) for hunted and unhunted populations were both 0.076, and although the proportion of polymorphic loci per sample and the mean number of alleles per locus for each sample were lower for the hunted populations than for the unhunted populations, these differences were not significant. However, the hunted populations displayed higher levels of outbreeding (lower F _IS and F _IT values) than those for unhunted populations. Therefore, it is concluded that although greywing francolin populations contain relatively high levels of genetic heterogeneity, it is probably the increased levels of local immigration following hunting which reduces the effects of any reduction in genetic variation due to a decrease in local population size from hunting.     

Genetic inference of orchid population dynamics on different-aged lava flows in Costa Rica

Colonization is a fundamental ecological process that is important for the persistence of species, particularly when a changing environment necessitates range shifts. Vacant habitats available for colonization often arise from landscape disturbance. Colonization and population expansion processes can be inferred by examining the levels and spatial distribution of genetic variation of plant populations with known disturbance histories. Samples (N = 690) of the terrestrial orchid, Epidendrum radicans, were collected from five lava flow sites on the slopes of Volcán Arenal in Costa Rica that last experienced major eruptions in 1968 and 1992. Individuals were also sampled (N = 188) from four regional populations. Samples were characterized using 15 nuclear genetic markers and analyzed using population genetics statistics. Genetic diversity within sites was moderate (H_e = 0.092–0.192). Contrary to expectation, diversity tended to be lower on the older lava flows (0.131 vs. 0.172) which may reflect their more sheltered topography that restricted pollen/seed immigration, and/or greater intra- and interspecific competition. Genetic diversity measures indicate that the lava flows were colonized by numerous individuals that likely originated from multiple sources while spatial genetic structure (SGS) statistics indicate that most recruitment in the study sites subsequent to colonization resulted from in situ reproduction and localized seed deposition. Younger sites had significantly greater SGS over larger distances which reflects fewer reproductive events, and less spatial and temporal overlap of seed shadows relative to the older sites. Clones were also generally smaller on the younger sites (≤3 m vs. ≤8 m).     

Limited genetic differentiation of Mycetomoellerius mikromelanos in Parque National Soberanía,Panama: Implications for queen dispersal

The coevolutionary relationship between fungus-growing ants (Formicidae: Attini: Attina) and their symbionts has been well studied in the Panamanian rain forests. To further understand the ecological context of these evolutionary relationships, we have examined the population-genetic structure of the fungus-growing ant species Mycetomoellerius mikromelanos Cardenas, Schultz, Adams 2021 in the Panama Canal Zone. We specifically investigated the presence of population structure, the significance of geographic features (i.e., creeks) limiting gene flow, and relatedness between ant colonies. To accomplish this, we genotyped 85 ant colonies from nine creeks across an approximately 30 km transect in Parque National Soberanía, Panama, using double-digest restriction-site-associated DNA sequencing. We did not find distinct population structure using two genetic clustering methods; however, we did detect an effect of isolation by distance. Furthermore, related colonies were frequently detected on the same creek or neighboring creeks, and some at further geographic distances. Collectively, these findings demonstrate that new colonies tend to establish on natal creeks and occasionally on distant creeks following long-distance dispersal events. We discuss how population-genetic patterns reveal the natural history of M. mikromelanos in Parque National Soberanía and how these results fit into the context of fungus-growing ant mutualisms. Abstract in Spanish is available with online material.     

59个Y-STR基因座的法医学多态性应用效能评估

为了评估法医学DNA数据库建设中涉及的59个Y-STR基因座的遗传多态性和法医学应用效能,通过AGCU Y SUPP PLUS试剂盒和AGCU Y37试剂盒检测374个广东汉族无关男性个体,将检测出的59个Y-STR基因座按照突变率的高低进行分类组合和统计分析。结果显示,59个Y-STR基因座联合运用在374个无关男性个体中检出了374个单倍型,其中44个中低突变Y-STR组合、15个高快突变Y-STR组合分别检出373和372个单倍型。59个Y-STR的基因多态性数值分布在0.055 1 (DYS645)~0.958 0 (DYF387S1 a/b)之间。结果表明,这59个Y-STR在广东汉族群体中均具有良好的多态性,按中低突变Y-STR组合和高快突变Y-STR组合研发新的检测体系可更好地满足法医实践的不同需求。     

snpR: User friendly population genomics for SNP data sets with categorical metadata

The analysis of genomic data can be an intimidating process, particularly for researchers who are not experienced programmers. Commonly used analyses are spread across many programs, each requiring their own specific input formats, and so data must often be repeatedly reorganized and transformed into new formats. Analyses often require splitting data according to metadata variables such as population or family, which can be challenging to manage in large data sets. Here, we introduce snpR, a user-friendly data analysis package in R for processing SNP genomic data. snpR is designed to automate data subsetting and analyses across categorical metadata while also streamlining repeated analyses by integrating approaches contained in many different packages in a single ecosystem. snpR facilitates iterative and efficient analyses centred on a single R object for an entire analysis pipeline.     

Inferring landscape resistance to gene flow when genetic drift is spatially heterogeneous

In connectivity models, land cover types are assigned cost values characterizing their resistance to species movements. Landscape genetic methods infer these values from the relationship between genetic differentiation and cost distances. The spatial heterogeneity of population sizes, and consequently genetic drift, is rarely included in this inference although it influences genetic differentiation. Similarly, migration rates and population spatial distributions potentially influence this inference. Here, we assessed the reliability of cost value inference under several migration rates, population spatial patterns and degrees of population size heterogeneity. Additionally, we assessed whether considering intra-population variables, here using gravity models, improved the inference when drift is spatially heterogeneous. We simulated several gene flow intensities between populations with varying local sizes and spatial distributions. We then fit gravity models of genetic distances as a function of (i) the ‘true’ cost distances driving simulations or alternative cost distances, and (ii) intra-population variables (population sizes, patch areas). We determined the conditions making the identification of the ‘true’ costs possible and assessed the contribution of intra-population variables to this objective. Overall, the inference ranked cost scenarios reliably in terms of similarity with the ‘true’ scenario (cost distance Mantel correlations), but this ‘true’ scenario rarely provided the best model goodness of fit. Ranking inaccuracies and failures to identify the ‘true’ scenario were more pronounced when migration was very restricted (<4 dispersal events/generation), population sizes were most heterogeneous and some populations were spatially aggregated. In these situations, considering intra-population variables helps identify cost scenarios reliably, thereby improving cost value inference from genetic data.