Strong Population Structure Research Articles

The Evolution of Wisteria Vein Mosaic Virus: A Case Study Approach to Track the Emergence of New Potyvirus Threats.

Wisteria vein mosaic virus (WVMV, Potyvirus wisteriae), a virus belonging to the genus Potyvirus, is responsible for Wisteria vein mosaic disease (WMD), a severe disease that affects Wisteria, a genus of garden plants acclaimed worldwide. Although probably originating in the Far East, WVMV infection was first reported in the US, and subsequently in numerous countries. Following the first molecular detection of an Italian isolate, WVMV Bari, its full-length genome was achieved using NGS barcoding technology. A PhyML phylogenetic analysis, supported by clustering algorithm validation, identified a clear separation between two phylogroups. One major clade comprised WVMV strains isolated from Wisteria spp. A second clade grouped three highly divergent strains, at the borderline species threshold, all found in non-wisteria hosts. Relying on a Relative Time Dated Tips (RTDT) molecular clock, the first emergence of WVMV clades has been traced back to around the 17th century. A network inference analysis confirmed the sharp separation between the two host-related phylogroups, also highlighting the presence of potential intermediate variants. Inter-population genetic parameters revealed a very high genetic differentiation in both populations, which was made reliable by statistically significant permutation tests. The migrant number (Nm) and fixation index (FST) evidenced a restricted gene flow and strong population structures. According to the dN/dS ratio and negative neutrality tests, it was derived that purifying selection at the expense of non-silent variants is underway within WVMV populations. Targeting WVMV evolutionary traits, the present effort raised interesting questions about the underestimated potential of this culpably neglected species to spread in economically relevant crops. The main intention of our study is, therefore, to propose an evolution-based analysis approach that serves as a case study to investigate how other potyviruses or newly emerging viruses may spread.

Population genetics of Macrognathus siamensis (Synbranchiformes: Mastacembelidae): Implications for non-migratory fishery resources in the Mekong River basin

The spotted spiny eel, Macrognathus siamensis is an economically important freshwater fish in the Mekong River basin, which is undergoing dramatic biodiversity changes due to anthropogenic impacts. The species is non-migratory, with a moderate larval duration that facilitates downstream genetic connectivity. Restriction site-associated DNA sequence datasets that includes 3736 and 1244 single nucleotide polymorphisms representing putatively neutral loci, from all geographic sites and without the Tachileik site, indicate strong population structure in this sedentary fish, with a significant isolation-by-distance signature. This structure reflects isolation of tributary populations and downstream dispersal of these distinct populations to mix in the Mekong mainstem and contributes to increased genetic diversity in the lower reaches of the Mekong basin. Genetic data indicates that the dispersal of eggs and larvae downstream is not impeded by Khone Falls and that the falls prevent upstream dispersal of eggs and larvae. Despite this downstream connectivity, there is a strong population structure among above- and below-Khone Falls groups. An analysis of outlier loci putatively under selection provides evidence that the difference between connected metapopulations above and below the falls is due to the fundamental differences in environmental regimes. We hypothesize that instead of a physical barrier, the falls simply represents a demarcation between previously recognized upstream limited-floodplain environments and downstream extensive-floodplain environments. These findings emphasize the need to monitor genetic diversity of key sedentary fishery species to assess whether proposed dams along the river could affect genetic diversity downstream. This genetic diversity is particularly important for the resilience of downstream populations because of the many environmental impacts caused by dams, land use, and climate change. Early detection of reduced downstream genetic diversity could trigger remediation programs to ensure resilience and continued fisheries productivity of important non-migratory fish resources.

Single nucleotide polymorphism (SNP) analysis reveals ancestry and genetic diversity of cultivated and wild grapevines in Croatia

BackgroundCroatia is a geographically small country with a remarkable diversity of cultivated and spontaneous grapevines. Local germplasm has been characterised by microsatellite markers, but a detailed analysis based on single nucleotide polymorphisms (SNPs) is still lacking. Here we characterize the genetic diversity of 149 accessions from three germplasm repositories and four natural sites using 516,101 SNPs to identify complete parent-offspring trios and their relations with spontaneous populations, offering a proof-of-concept for the use of reduced-representation genome sequencing in population genetics and genome-wide association studies (GWAS).ResultsPrincipal component analysis revealed a clear discontinuity between cultivated (V. vinifera subsp. sativa) and spontaneous grapevines, supporting the notion that the latter represent local populations of the wild progenitor (V. vinifera subsp. sylvestris). ADMIXTURE identified three ancestry components. Two sativa components are alternatively predominant in cultivars grown either in northern Adriatic Croatia and Continental Croatia or in Dalmatia (i.e. central and southern Adriatic Croatia). A sylvestris component, which is predominant in accessions from spontaneous populations, is a minor ancestry component in cultivated accessions. TREEMIX provided evidence of unidirectional migration from the vineyards to natural sites, suggesting that gene flow has gone preferentially from the introduced domesticated germplasm into local wild populations rather than vice versa. Identity-by-descent analysis indicated an extensive kinship network, including 14 complete parent-offspring trios, involving only cultivated accessions, six full-sibling relationships and invalidated a presumed pedigree of one of the most important varieties in Croatia, ‘Plavac Mali’. Despite this strong population structure, significant association was found between 143 SNPs and berry skin colour and between 2 SNPs and leaf hairiness, across two previously known genomic regions.ConclusionsThe clear genetic separation between Croatian cultivars and sylvestris ruled out the hypothesis that those cultivars originated from local domestication events. On the other hand, the evidence of a crop-to-wild gene flow signals the need for an urgent adoption of conservation strategies that preserve the residual genetic integrity of wild relatives. The use of this reduced-representation genome sequencing protocol in grapevine enables an accurate pedigree reconstruction and can be recommended for GWAS experiments.

Resequencing reveals population structure and genetic diversity in Tibetan sheep

BackgroundThe Tibetan sheep is one of the three major primitive sheep breeds in China, representing a unique and high-quality genetic resource in the Qinghai-Tibet Plateau and neighboring high-altitude regions, exhibiting exceptional adaptability to high-altitude climatic environments. However, research on the genetic relationships among different populations of Tibetan sheep at the whole-genome level remains insufficient. This study aims to explore the population structure and historical dynamics among 11 Tibetan sheep populations, accurately assess the genetic diversity within the populations, and providing a theoretical basis for the development of targeted genetic breeding strategies for Tibetan sheep.ResultsIn this study, a total of 10,884,454 high-quality SNPs were obtained. All Tibetan sheep populations exhibited varying degrees of linkage disequilibrium, with similar decay rates; among them, the WT population showed the fastest decay, while the TS population exhibited the slowest decay rate. Analyses using Tajima’s D and π indicated that the genetic diversity levels of the Tibetan sheep populations are generally low. Fst results revealed that most populations exhibited moderate to low levels of genetic differentiation. The effective population size among Tibetan sheep populations showed an increasing trend over time. The evolutionary relationships among Tibetan sheep populations reflect the correlation between their geographical locations and genomic genetic distances, while also indirectly confirming the impact of historical activities such as early human migration, admixture, fusion, and expansion on the population sizes and distributions of Tibetan sheep.ConclusionsThe results indicate that the genetic diversity levels and genetic differentiation among Tibetan sheep populations are relatively low. In this study, we identified the genetic characteristics of Tibetan sheep populations, which exhibit low levels of diversity, genetic differentiation, and a strong population structure. A deeper genomic exploration of the population structure and diversity status of Tibetan sheep populations will provide theoretical support for subsequent genetic breeding and diversity conservation efforts.

Examining the Effects of Environment, Geography, and Elevation on Patterns of DNA Methylation Across Populations of Two Widespread Bumble Bee Species.

Understanding the myriad avenues through which spatial and environmental factors shape evolution is a major focus in biological research. From a molecular perspective, much work has been focused on genomic sequence variation; however, recently there has been increased interest in how epigenetic variation may be shaped by different variables across the landscape. DNA methylation has been of particular interest given that it is dynamic and can alter gene expression, potentially offering a path for a rapid response to environmental change. We utilized whole genome enzymatic methyl sequencing to evaluate the distribution of CpG methylation across the genome and to analyze patterns of spatial and environmental association in the methylomes of two broadly distributed montane bumble bees (Bombus vancouverensis Cresson and Bombus vosnesenskii Radoszkowski) across elevational gradients in the western US. Methylation patterns in both species are similar at the genomic scale with ∼1% of CpGs being methylated and most methylation being found in exons. At the landscape scale, neither species exhibited strong spatial or population structuring in patterns of methylation, although some weak relationships between methylation and distance or environmental variables were detected. Differential methylation analysis suggests a stronger environment association in B. vancouverensis given the larger number of differentially methylated CpG's compared to B. vosnesenskii. We also observed only a handful of genes with both differentially methylated CpGs and previously detected environmentally associated outlier SNPs. Overall results reveal a weak but present pattern in variation in methylation over the landscape in both species.

Variation of the 3'RR1 HS1.2 Enhancer and Its Genomic Context.

In humans, the HS1.2 enhancer in the Ig heavy-chain locus is modular, with length polymorphism. Previous studies have shown the following features for this variation: (i) strong population structuring; (ii) association with autoimmune diseases; and (iii) association with developmental changes in Ig expression. The HS1.2 region could then be considered as a contributor to inter-individual diversity in humoral response in adaptive immunity. We experimentally determined the HS1.2-length class genotype in 72 of the 1000 Genomes CEU cell lines and assigned the HS1.2 alleles to haplotypes defined by 18 landmark SNPs. We also sequenced the variable portion and ~200 bp of the flanking DNA of 34 HS1.2 alleles. Furthermore, we computationally explored the ability of different allelic arrangements to bind transcription factors. Non-random association between HS1.2 and Gm allotypes in the European population clearly emerged. We show a wealth of variation in the modular composition of HS1.2, with five SNPs further contributing to diversity. Longer alleles offer more potential sites for binding but, for same-length alleles, SNP variation creates/destroys potential binding sites. Altogether, the arrangements of modules and SNP alleles both inside and outside HS1.2 denote an organization of diversity far from randomness. In the context of the strong divergence of human populations for this genomic region and the reported disease associations, our results suggest that selective forces shaped the pattern of its diversity.

Modeling the transmission of antibiotic-resistant Enterobacterales in the community: A systematic review

BackgroundAntibiotic-resistant Enterobacterales (ARE) are a public health threat worldwide. Dissemination of these opportunistic pathogens has been largely studied in hospitals. Despite high prevalence of asymptomatic colonization in the community in some regions of the world, less is known about ARE acquisition and spread in this setting. As explaining the community ARE dynamics has not been straightforward, mathematical models can be key to explore underlying phenomena and further evaluate the impact of interventions to curb ARE circulation outside of hospitals. MethodsWe conducted a systematic review of mathematical modeling studies focusing on the transmission of AR-E in the community, excluding models only specific to hospitals. We extracted model features (population, setting), formalism (compartmental, individual-based), biological hypotheses (transmission, infection, antibiotic impact, resistant strain specificities) and main findings. We discussed additional mechanisms to be considered, open scientific questions, and most pressing data needs. ResultsWe identified 18 modeling studies focusing on the human transmission of ARE in the community (n=11) or in both community and hospital (n=7). Models aimed at (i) understanding mechanisms driving resistance dynamics; (ii) identifying and quantifying transmission routes; or (iii) evaluating public health interventions to reduce resistance. To overcome the difficulty of reproducing observed ARE dynamics in the community using the classical two-strains competition model, studies proposed to include mechanisms such as within-host strain competition or a strong host population structure. Studies inferring model parameters from longitudinal carriage data were mostly based on models considering the ARE strain only. They showed differences in ARE carriage duration depending on the acquisition mode: returning travelers have a significantly shorter carriage duration than discharged hospitalized patient or healthy individuals. Interestingly, predictions across models regarding the success of public health interventions to reduce ARE rates depended on pathogens, settings, and antibiotic resistance mechanisms. For E. coli, reducing person-to-person transmission in the community had a stronger effect than reducing antibiotic use in the community. For Klebsiella pneumoniae, reducing antibiotic use in hospitals was more efficient than reducing community use. ConclusionsThis study raises the limited number of modeling studies specifically addressing the transmission of ARE in the community. It highlights the need for model development and community-based data collection especially in low- and middle-income countries to better understand acquisition routes and their relative contribution to observed ARE levels. Such modeling will be critical to correctly design and evaluate public health interventions to control ARE transmission in the community and further reduce the associated infection burden.

Species-specific dynamics may cause deviations from general biogeographical predictions - evidence from a population genomics study of a New Guinean endemic passerine bird family (Melampittidae).

The family Melampittidae is endemic to New Guinea and consists of two monotypic genera: Melampitta lugubris (Lesser Melampitta) and Megalampitta gigantea (Greater Melampitta). Both Melampitta species have scattered and disconnected distributions across New Guinea in the central mountain range and in some of the outlying ranges. While M. lugubris is common and found in most montane regions of the island, M. gigantaea is elusive and known from only six localities in isolated pockets on New Guinea with very specific habitats of limestone and sinkholes. In this project, we apply museomics to determine the population structure and demographic history of these two species. We re-sequenced the genomes of all seven known M. gigantaea samples housed in museum collections as well as 24 M. lugubris samples from across its distribution. By comparing population structure between the two species, we investigate to what extent habitat dependence, such as in M. gigantaea, may affect population connectivity. Phylogenetic and population genomic analyses, as well as acoustic variation revealed that M. gigantaea consists of a single population in contrast to M. lugubris that shows much stronger population structure across the island. We suggest a recent collapse of M. gigantaea into its fragmented habitats as an explanation to its unexpected low diversity and lack of population structure. The deep genetic divergences between the M. lugubris populations on the Vogelkop region, in the western central range and the eastern central range, respectively, suggests that these three populations should be elevated to full species level. This work sheds new light on the mechanisms that have shaped the intriguing distribution of the two species within this family and is a prime example of the importance of museum collections for genomic studies of poorly known and rare species.

Detecting inbreeding depression in structured populations

Measuring inbreeding and its consequences on fitness is central for many areas in biology including human genetics and the conservation of endangered species. However, there is no consensus on the best method, neither for quantification of inbreeding itself nor for the model to estimate its effect on specific traits. We simulated traits based on simulated genomes from a large pedigree and empirical whole-genome sequences of human data from populations with various sizes and structures (from the 1,000 Genomes project). We compare the ability of various inbreeding coefficients ([Formula: see text]) to quantify the strength of inbreeding depression: allele-sharing, two versions of the correlation of uniting gametes which differ in the weight they attribute to each locus and two identical-by-descent segments-based estimators. We also compare two models: the standard linear model and a linear mixed model (LMM) including a genetic relatedness matrix (GRM) as random effect to account for the nonindependence of observations. We find LMMs give better results in scenarios with population or family structure. Within the LMM, we compare three different GRMs and show that in homogeneous populations, there is little difference among the different [Formula: see text] and GRM for inbreeding depression quantification. However, as soon as a strong population or family structure is present, the strength of inbreeding depression can be most efficiently estimated only if i) the phenotypes are regressed on [Formula: see text] based on a weighted version of the correlation of uniting gametes, giving more weight to common alleles and ii) with the GRM obtained from an allele-sharing relatedness estimator.

Genetic basis of local adaptation in the cold-tolerant mangrove Kandelia obovata.

Understanding the genetic basis of local adaption is crucial in the context of global climate change. Mangroves, as salt-tolerant trees and shrubs in the intertidal zone of tropical and subtropical coastlines, are particularly vulnerable to climate change. Kandelia obovata, the most cold-tolerant mangrove species, has undergone ecological speciation from its cold-intolerant counterpart, Kandelia candel, with geographic separation by the South China Sea. In this study, we conducted whole-genome re-sequencing of K. obovata populations along China's southeast coast, to elucidate the genetic basis responsible for mangrove local adaptation to climate. Our analysis revealed a strong population structure among the three K. obovata populations, with complex demographic histories involving population expansion, bottleneck, and gene flow. Genome-wide scans unveiled pronounced patterns of selective sweeps in highly differentiated regions among pairwise populations, with stronger signatures observed in the northern populations compared to the southern population. Additionally, significant genotype-environment associations for temperature-related variables were identified, while no associations were detected for precipitation. A set of 39 high-confidence candidate genes underlying local adaptation of K. obovata were identified, which are distinct from genes under selection detected by comparison between K. obovata and its cold-intolerant relative K. candel. These results significantly contribute to our understanding of the genetic underpinnings of local adaptation in K. obovata and provide valuable insights into the evolutionary processes shaping the genetic diversity of mangrove populations in response to climate change.

Phenotypic and Molecular Diversity of Wild Populations of Acca sellowiana (Berg.) Burret in the Southern Area of Natural Distribution

Acca sellowiana is a subtropical tree in the myrtle family (Myrtaceae) native to southern Brazil and northeastern Uruguay. It is recognized for its value as a fruit-bearing, ornamental, and medicinal species. Based on distinctive characteristics of fruits, seeds, and leaves, as well as its geographical distribution pattern, two variants of the species are distinguished: the “Brazilian type” and the “Uruguayan type”. The objective of this study was to characterize, for the first time, the diversity of 202 individuals from four wild populations in Uruguay, representative of the species’ most southern natural distribution. Twenty-three morphological descriptors (leaf, flower, and fruit) and 204 RAPD (Random Amplified Polymorphic DNA) markers were used. The morphological data collected validated the main criteria that distinguish “Uruguayan type” populations from “Brazilian type” populations, such as lower seed weight and fruit size, thin and slightly rough skin, high pulp percentage, and hairy white abaxial leaf surfaces. Analyses of both morphological and molecular data indicated wide diversity and strong population structuring, which is relevant information for designing conservation plans, sustainable utilization, and genetic improvement of the plant genetic resources of this species.

Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant.

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Genetic and geographic population structure in the malaria vector, Anopheles farauti, provides a candidate system for pioneering confinable gene-drive releases

Indoor insecticide applications are the primary tool for reducing malaria transmission in the Solomon Archipelago, a region where Anopheles farauti is the only common malaria vector. Due to the evolution of behavioural resistance in some An. farauti populations, these applications have become less effective. New malaria control interventions are therefore needed in this region, and gene-drives provide a promising new technology. In considering developing a population-specific (local) gene-drive in An. farauti, we detail the species’ population genetic structure using microsatellites and whole mitogenomes, finding many spatially confined populations both within and between landmasses. This strong population structure suggests that An. farauti would be a useful system for developing a population-specific, confinable gene-drive for field release, where private alleles can be used as Cas9 targets. Previous work on Anopheles gambiae has used the Cardinal gene for the development of a global population replacement gene-drive. We therefore also analyse the Cardinal gene to assess whether it may be a suitable target to engineer a gene-drive for the modification of local An. farauti populations. Despite the extensive population structure observed in An. farauti for microsatellites, only one remote island population from Vanuatu contained fixed and private alleles at the Cardinal locus. Nonetheless, this study provides an initial framework for further population genomic investigations to discover high-frequency private allele targets in localized An. farauti populations. This would enable the development of gene-drive strains for modifying localised populations with minimal chance of escape and may provide a low-risk route to field trial evaluations.

Continent-island boundary and environment-shaped evolution in the marine amphipod Ampithoe marcuzzii complex (Crustacea: Eumalacostraca: Ampithoidae)

Marine amphipods are crustaceans that lack a larval phase and consequently have low dispersion rates. Despite that, these crustaceans present a remarkable ability to be transported by rafting on natural floating substrata, especially macroalgae, where they find shelter, food and a mating ground. The species Ampithoe marcuzzii is widely distributed throughout the western Atlantic Ocean. Here, it was used as a model to study seascape genomics and phylogeography in invertebrates with low dispersion capacities. We anticipated that the lineages would present isolation-by-distance patterns. However, surface currents and other abiotic variables could facilitate connectivity among distant sites. Based on mitochondrial and nuclear genes, SNPs, and environmental associations, we observed the presence of a species complex within A. marcuzzii, separating mainland and insular populations. Each species showed an independent evolutionary history, with a strong latitudinal population structure and evidence of isolation-by-distance and isolation-by-environment, characterizing the 'continent' species. Historical expansion and environmental variables were observed associated with the southeastern population, and ecological niche modeling corroborated the region as a paleorefuge. Conversely, populations from 'islands' presented complicated evolutionary histories, with closer localities genetically isolated and distant localities connected. These findings indicate that insular populations with low dispersion capacity might be more susceptible to spatial connectivity by floating substrata and to changes in surface currents. In contrast, mainland populations might be more vulnerable to local climate changes due to lack of gene flow.

Digest: Islands promote population differentiation in dispersive swallows.

How can highly dispersive species give rise to genetically distinct populations? This seemingly paradoxical pattern is common among insular birds, but not in those with continental distributions. Broyles and Myers sequence the genomes of almost 150 individuals from the island-dwelling Pacific swallow (Hirundo tahitica) and its continental counterpart, the welcome swallow (H. neoxena). They find strong population structure only among island populations and attribute this to a behavioral reduction in dispersal propensity following island colonization. However, wing shape remains consistent across populations, suggesting it might not accurately reflect dispersal propensity in this group. This study illustrates the interplay between dispersal, isolation, and divergence, offering insights into how geographic factors affect speciation and population differentiation on islands.

The association between phylogenetic lineage and the subclinical phenotype of pulmonary tuberculosis: A retrospective 2-cohort study

Subclinical pulmonary tuberculosis (PTB) is an asymptomatic disease state between established TB infection and symptomatic (clinical) TB disease. It is present in 20-25% of PTB patients in high-income countries. Mycobacterium tuberculosis complex (MTBC) genetic heterogeneity, and differential host immunological responses, have been implicated in its pathogenesis. To determine the association between MTBC lineage and PTB disease phenotype, we used two retrospective cohorts of PTB patients in Canada and two independent lineage attribution methods (DNA fingerprinting and genome sequencing). The first cohort, Cohort 1, consisted of consecutively diagnosed PTB patients between 2014 and 2020. The second, Cohort 2, consisted of newly-arrived foreign-born PTB patients who either were or were not referred for post-landing medical surveillance between 2004 and 2017. Univariable and multivariable logistic regression models were sequentially fitted to both cohorts, adjusting for age, sex, disease type, drug resistance and HIV. Evolution of radiographic features was correlated to lineage in Cohort 2. Cohort 1 and 2 included 874 (209 subclinical) and 111 (44 subclinical) patients, respectively. In both cohorts, subclinical patients were more likely than clinical patients to have relapse/retreatment disease, be smear-negative, have longer times-to-culture positivity and to harbor an ancestral MTBC lineage (Indo-Oceanic or Mycobacterium africanum). Relapse/retreatment disease and ancestral MTBC lineage were independent predictors of subclinical disease (ORs and 95% CIs in Cohort 1, 1.85 [1.07,3.28], p<0.029 and 2.30 [1.66,3.18], p<0.001, respectively, and Cohort 2, 5.74 [1.37-24.06], p<0.017 and 3.21 (1.29,7.97], p<0.012, respectively). The geographic distribution of Indo-Oceanic strains causing subclinical disease was uneven. Non-progressive lung disease was more common in patients infected with ancestral than modern lineages in Cohort 2, 56.0% vs 25.4%, p<0.005. MTBC lineage is a strong predictor of PTB disease phenotype. The genetic drivers of this association, and the relative contribution of other explanatory variables, are unknown.

Indigenous Australian genomes show deep structure and rich novel variation

The Indigenous peoples of Australia have a rich linguistic and cultural history. How this relates to genetic diversity remains largely unknown because of their limited engagement with genomic studies. Here we analyse the genomes of 159 individuals from four remote Indigenous communities, including people who speak a language (Tiwi) not from the most widespread family (Pama–Nyungan). This large collection of Indigenous Australian genomes was made possible by careful community engagement and consultation. We observe exceptionally strong population structure across Australia, driven by divergence times between communities of 26,000–35,000 years ago and long-term low but stable effective population sizes. This demographic history, including early divergence from Papua New Guinean (47,000 years ago) and Eurasian groups1, has generated the highest proportion of previously undescribed genetic variation seen outside Africa and the most extended homozygosity compared with global samples. A substantial proportion of this variation is not observed in global reference panels or clinical datasets, and variation with predicted functional consequence is more likely to be homozygous than in other populations, with consequent implications for medical genomics2. Our results show that Indigenous Australians are not a single homogeneous genetic group and their genetic relationship with the peoples of New Guinea is not uniform. These patterns imply that the full breadth of Indigenous Australian genetic diversity remains uncharacterized, potentially limiting genomic medicine and equitable healthcare for Indigenous Australians.

Sample contamination explains evidence of close kin proximity in yellowfin tuna (Thunnus albacares) in the Western and Central Pacific Ocean

Anderson et al.’s 2019 paper ‘Close Kin Proximity in Yellowfin Tuna (Thunnus albacares) as a Driver of Population Genetic Structure in the Tropical Western and Central Pacific’ provided observations of genetically related individuals among sampled yellowfin tuna. Resampling of some individuals from the original study produced very different sequencing results when compared with the original dataset, one possible explanation of which is cross-contamination among the original samples. Re-analyses produced no indication of strong population structure, including a very slightly and consistent heterozygosity deficit with only one statistically significant pairwise FST value of 0.002. We identified only one pair of genetically similar individuals that could still be an artefact of lingering contamination. These new results therefore support a highly connected population of yellowfin tuna in the study area, although we do not rule out the potential for local adaptation driven by non-observed loci or genetic sub-structure operating at more contemporary scales.

Population genetics of Camellia granthamiana, an endangered plant species with extremely small populations in China.

Introduction: Camellia, the largest genus of Theaceae, is well-known for having high economic values. Camellia granthamiana demonstrates large beautiful flowers with some primitive characters, such as multiple large and persistent bracteoles and sepals, was listed as Vulnerable species on the IUCN Red List. Methods: In this study, we investigated all possible records of the species, and sampled four natural populations and five cultivated individuals. By applying shallow-genome sequencing for nine individuals and RAD-seq sequencing for all the sampled 77 individuals, we investigated population genetic diversity and population structure of the species. Results and discussion: The results showed that the population sampled from Fengkai, previously identified as C. albogigias, possessed different plastid genome from other species possibly due to plastid capture; the species possesses strong population structure possibly due to the effect of isolation by distance, habitat fragmentation, and self-crossing tendency of the species, whose effective population size declined quickly in the past 4,000 years. Nevertheless, C. granthamiana maintains a medium level of genetic diversity within population, and significant differentiation was observed among the four investigated populations, it is anticipated that more populations are expected to be found and all these extant populations should be taken into instant protection.

Eelgrass Zostera capensis populations in KwaZulu-Natal, South Africa, harbour distinct genomic signals despite limited geographical distance

Seagrasses are threatened by anthropogenic stressors and climate change, with numerous population declines reported. In South Africa, the eelgrass Zostera capensis is restricted to estuarine environments and has a disjunct distribution and declining status. With the majority of the distribution of Z. capensis on the west and southwest coasts of South Africa, the isolated easternmost populations in KwaZulu-Natal Province (KZN) are of particular interest. Following the extirpation of Z. capensis at Durban and St Lucia, only five populations remain, of which three (in the Amatikulu, Mlalazi and Mhlatuze estuaries) are situated <50 km apart. Previous molecular analyses showed strong population structure between Z. capensis populations, but the geographically close populations on the east coast were not included. In this study, using single nucleotide polymorphisms (SNPs) isolated from 31 individual sampled plants collected from three northern KZN estuaries, we provide evidence for distinct population clusters, with unique evolutionary signals. Zostera capensis in the Mlalazi Estuary has a low level of genomic diversity, likely as a result of a small, dynamic population unable to withstand prolonged freshwater exposure. Our results suggest that conservation efforts need to consider unique population signals even among geographically close populations, in particular within the context of restoration, where genomic compatibility may determine the persistence of restored populations.