Census Size Research Articles

Amphibian decline in a Central European forest and the importance of woody debris for population persistence

Achieving an understanding of species’ temporal and spatial abundance patterns is a prerequisite for proper conservation and management decisions. Little is known of the ecological indicators of census size or biomass of terrestrial amphibian communities and historical estimates of population size are rare. Here we took advantage of abundance and biomass data collected in 1967/1968 for a leaf-litter amphibian community in Niepołomice Forest (NF) in south-central Poland. Using the same technique (removal sampling), we censused the amphibian community of this forest 50 years later (2016/2017), spatially matching a subset of current plots with the historical sampling sites. Averaged over all plots, we revealed at least a fourfold decline in amphibian abundance from 2817 ± 776 individuals/ha in 1967 to 674 ± 508 individuals/ha currently, paralleled by a decline in biomass from roughly 15.5 ± 1.7 to 3.8 ± 3.1 kg/ha. The decline was more striking (over 15-fold) for current plots placed in the same forest division as the historical plots (mean counts of 157 ± 103 individuals/ha corresponding in biomass to 1.0 ± 0.5 kg/ha). We also recorded a change in species composition between surveys, and propose that conversion of aquatic habitat and forest succession have decreased amphibian diversity and abundance in NF. Further, we estimated densities of anuran amphibians in forest plots differing in management regime and local habitat characteristics for contemporary plots (all historical plots were placed in managed forest). We found that unmanaged, old growth forest plots with coarse woody debris on the forest floor held on average twice as many anurans compared to plots in managed stands. Average body condition of the most common species, Rana temporaria, measured as a scaled mass index, was enhanced in plots situated in old growth forest. Our findings suggest that the preservation of old growth forests containing coarse woody debris may boost anuran abundance, biomass and body condition, and has the potential to buffer against long-term demographic decline. Moreover, the retention of deadwood could be a means of increasing the conservation potential of managed forests for terrestrial amphibian communities.

Genetic effects of demographic bottleneck and recovery in Kinkazan Island and mainland populations of Japanese macaques (Macaca fuscata).

Populations of Japanese macaques were significantly reduced in most areas from the 1900s to the 1960s and then recovered mainly in the northeastern part of Honshu. A drastic reduction in population size reduces genetic variability through a bottleneck effect. Demographic expansion after the reduction that accumulates new mutations can reduce the bottleneck effects or drive the recovery of genetic variability. We examined the genetic status of a small island population (Kinkazan Island) and a larger mainland population (southern Tohoku) of Japanese macaques that experienced recent demographic bottlenecks and recovery using eight microsatellite loci. The two populations were significantly genetically different from each other. The Kinkazan population exhibited lower genetic variability, remarkable evidence of bottleneck (i.e., significant heterozygosity excess and lower frequency of rare alleles), and a considerably smaller effective population size based on genetic data than based on the current census size. These results indicate that the genetic status has not completely recovered from the demographic bottleneck despite a full recovery in census size on Kinkazan Island. New mutations might rarely have accumulated because of the small carrying capacity of the island. Therefore, the genetic variability of the population would have been restrained by the severe bottleneck size, small carrying capacity, and long-term isolation. On the other hand, the bottleneck effect seems to be limited in the southern Tohoku population considering higher genetic variability, non-significant heterozygosity excess in many mutation conditions, and the highest frequency of rare alleles.

Changes in human effective population size overlap the beginning and end of a critical time in European medieval history, also characterized by the Black Death epidemic

The intersection of historiography and archaeology has long pondered over the impact of known historical events on census size. In recent times, genetic methods have successfully traced changes over time in the genetic size of a given population. Moreover, the correlation between genetic and census sizes of a population is contingent on several demographic assumptions that are relatively simple for our species. Our research endeavours to examine the changes in effective population size (Ne) in all human populations in the 1000 Genomes Project over the past two millennia. We compared our findings with estimates from historical censuses where available. Our investigation confirms what was already observed in France and reveals a common pattern found in most European populations, which manifests as a drastic population decrease beginning around the year 1300 and growth after the year 1600. This profile aligns well with known wars, famines, and epidemics that characterized these trying times in Europe. The most notable among them being the second plague epidemic, caused by Y. pestis, which in Europe commenced in 1347/8 and is also known as the "Black Death". Our findings demonstrate that changes in genetic population size through time can serve as a dependable proxy for census size, which is independent of potential biases in the written historical record. Consequently, we provide a robust estimate of the impact caused by the population crisis that followed the year 1300 on the European genomic landscape in light of previous results. Our study offers a new paradigm for interpreting the past and underscores the potential of genetic methods in reconstructing historical events.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples.

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Demographic history and conservation genomics of caribou (Rangifer tarandus) in Québec.

The loss of genetic diversity is a challenge many species are facing, with genomics being a potential tool to inform and prioritize decision-making. Most caribou (Rangifer tarandus) populations have experienced significant recent declines throughout Québec, Canada, and are considered of concern, threatened or endangered. Here, we calculated the ancestral and contemporary patterns of genomic diversity of five representative caribou populations and applied a comparative population genomics framework to assess the interplay between demographic events and genomic diversity. We first calculated a caribou specific mutation rate, μ, by extracting orthologous genes from related ungulates and estimating the rate of synonymous mutations. Whole genome re-sequencing was then completed on 67 caribou: from these data we calculated nucleotide diversity, θ π and estimated the coalescent or ancestral effective population size (N e), which ranged from 12,030 to 15,513. When compared to the census size, N C, the endangered Gaspésie Mountain caribou population had the highest ancestral N e:N C ratio which is consistent with recent work suggesting high ancestral N e:N C is of conservation concern. In contrast, values of contemporary N e, estimated from linkage-disequilibrium, ranged from 11 to 162, with Gaspésie having among the highest contemporary N e:N C ratio. Importantly, classic conservation genetics theory would predict this population to be of less concern based on this ratio. Interestingly, F varied only slightly between populations, and despite evidence of bottlenecks across the province, runs of homozygosity were not abundant in the genome. Tajima's D estimates mirrored the demographic models and current conservation status. Our study highlights how genomic patterns are nuanced and potentially misleading if viewed only through a contemporary lens; we argue a holistic conservation genomics view should integrate ancestral N e and Tajima's D into management decisions.

On the feasibility of estimating contemporary effective population size (Ne) for genetic conservation and monitoring of forest trees

Estimates of contemporary effective population size (Ne) can provide valuable information for genetic conservation and monitoring, pinpointing populations at higher risk of genetic erosion, decreased fitness, maladaptation and, ultimately, demographic decline. There are however potential limitations in the application of commonly employed genetic estimators of contemporary Ne to widespread forest tree populations. Genetic isolation by distance within populations, small and spatially restricted samples, among-population gene flow and overlapping generations are factors that can potentially affect the accuracy of marker-based estimates of contemporary Ne, depending on the demographic scenario. In particular, we illustrate the uncertainty faced by forest researchers and managers when interpreting contemporary Ne estimates obtained for continuously distributed tree populations with large census size N. To that end, we use previously published genotypic data of 21 Pinus pinaster populations, including distinct sampling schemes, together with a widely used method based on linkage disequilibrium patterns observed in a single (one-time) population sample. We hypothesize that spatially restricted sampling might be the main putative factor behind the apparently low Ne/N estimates obtained for the large and continuously-distributed populations studied here. Because of its statistically-inherent assessment difficulty, we call for caution when interpreting marker-based estimates of contemporary Ne for monitoring widely-distributed tree populations or small genetic conservation units embedded within large continuous tree populations.

Genetic Analysis and Status of Brown Bear Sub-Populations in Three National Parks of Greece Functioning as Strongholds for the Species' Conservation.

In order to optimize the appropriate conservation actions for the brown bear (Ursus arctos L.) population in Greece, we estimated the census (Nc) and effective (Ne) population size as well as the genetic status of brown bear sub-populations in three National Parks (NP): Prespa (MBPNP), Pindos (PINDNP), and Rhodopi (RMNP). The Prespa and Pindos sub-populations are located in western Greece and the Rhodopi population is located in eastern Greece. We extracted DNA from 472 hair samples and amplified through PCR 10 microsatellite loci. In total, 257 of 472 samples (54.5%) were genotyped for 6–10 microsatellite loci. Genetic analysis revealed that the Ne was 35, 118, and 61 individuals in MBPNP, PINDNP, and RMNP, respectively, while high levels of inbreeding were found in Prespa and Rhodopi but not in Pindos. Moreover, analysis of genetic structure showed that the Pindos population is genetically distinct, whereas Prespa and Rhodopi show mutual overlaps. Finally, we found a notable gene flow from Prespa to Rhodopi (10.19%) and from Rhodopi to Prespa (14.96%). Therefore, targeted actions for the conservation of the bears that live in the abovementioned areas must be undertaken, in order to ensure the species’ viability and to preserve the corridors that allow connectivity between the bear sub-populations in Greece.

Phylogenomics of the world’s otters

Comparative whole-genome analyses hold great power to illuminate commonalities and differences in the evolution of related species that share similar ecologies. The mustelid subfamily Lutrinae includes 13 currently recognized extant species of otters,1-5 a semiaquatic group whose evolutionary history is incompletely understood. We assembled a dataset comprising 24 genomes from all living otter species, 14 of which were newly sequenced. We used this dataset to infer phylogenetic relationships and divergence times, to characterize patterns of genome-wide genealogical discordance, and to investigate demographic history and current genomic diversity. We found that genera Lutra, Aonyx, Amblonyx, and Lutrogale form a coherent clade that should be synonymized under Lutra, simplifying the taxonomic structure of the subfamily. The poorly known tropical African Aonyx congicus and the more widespread Aonyx capensis were found to be reciprocally monophyletic (having diverged 440,000 years ago), supporting the validity of the former as a distinct species. We observed variable changes in effective population sizes over time among otters within and among continents, although several species showed similar trends of expansions and declines during the last 100,000 years. This has led to different levels of genomic diversity assessed by overall heterozygosity, genome-wide SNV density, and run of homozygosity burden. Interestingly, there were cases in which diversity metrics were consistent with the current threat status (mostly based on census size), highlighting the potential of genomic data for conservation assessment. Overall, our results shed light on otter evolutionary history and provide a framework for further in-depth comparative genomic studies targeting this group.

Estimation of the census (Nc) and effective (Ne) population size of a wild mandrill (Mandrillus sphinx) horde in the Lopé National Park, Gabon using a non-invasive genetic approach

Estimation of the census (Nc) and effective (Ne) population size of a wild mandrill (Mandrillus sphinx) horde in the Lopé National Park, Gabon using a non-invasive genetic approach

How Can We Resolve Lewontin's Paradox?

We discuss the genetic, demographic, and selective forces that are likely to be at play in restricting observed levels of DNA sequence variation in natural populations to a much smaller range of values than would be expected from the distribution of census population sizes alone—Lewontin’s Paradox. While several processes that have previously been strongly emphasized must be involved, including the effects of direct selection and genetic hitchhiking, it seems unlikely that they are sufficient to explain this observation without contributions from other factors. We highlight a potentially important role for the less-appreciated contribution of population size change; specifically, the likelihood that many species and populations may be quite far from reaching the relatively high equilibrium diversity values that would be expected given their current census sizes.

Sibship and parentage reconstruction as genetic tool for designing and monitoring plant translocations

Estimating contemporary effective population size (Ne), sibship and parentage can be useful for designing and monitoring plant translocations. Contemporary Ne is an important indicator of population viability, as it considers several demographic factors, e.g. sex ratio, mate compatibility, inbreeding, nonrandom mating, and family size. Sibship and parentage analyses allow not only estimating genetic relatedness of transplants and post‐translocation generations, but also selfing rates, parent contributions to offspring and contemporary seed and pollen dispersal. Sibship assignment method was used to reconstruct offspring pedigree and parentage and to estimate contemporary Ne, selfing rate, and contemporary seed and pollen dispersal in wild (used as seed sources for translocation) and translocated populations of the self‐compatible, insect‐pollinated Dianthus deltoides (Caryophyllaceae). Source populations show very small Ne estimates (5–21) and high selfing rates (28.6–68.9%), despite high census sizes. The first post‐translocation generation (recruits and transplant seed progeny) showed small contemporary Ne estimates (19–53) and substantial selfing rates (8.4–22.3%), despite high flowering sizes, overall genetic diversity and multiple paternity. Parentage assignments revealed high levels of admixture between sources, restricted seed dispersal and successful—often restricted—pollen dispersal in the post‐translocation generation. Pollinator service, crosses between selfed or related transplants and floral display, may have contributed to selfing and pollen dispersal patterns in translocated populations. The present findings support the interest of sibship and parentage reconstruction as genetic tool for designing and monitoring plant translocations, emphasizing the importance of assessing the genetic quality of the transplants, and contributing to a more precise evaluation of population extinction risk.

Genomic erosion in a demographically recovered bird species during conservation rescue.

The pink pigeon (Nesoenas mayeri) is an endemic species of Mauritius that has made a remarkable recovery after a severe population bottleneck in the 1970s to early 1990s. Prior to this bottleneck, an ex situ population was established from which captive‐bred individuals were released into free‐living subpopulations to increase population size and genetic variation. This conservation rescue led to rapid population recovery to 400–480 individuals, and the species was twice downlisted on the International Union for the Conservation of Nature (IUCN) Red List. We analyzed the impacts of the bottleneck and genetic rescue on neutral genetic variation during and after population recovery (1993–2008) with restriction site‐associated sequencing, microsatellite analyses, and quantitative genetic analysis of studbook data of 1112 birds from zoos in Europe and the United States. We used computer simulations to study the predicted changes in genetic variation and population viability from the past into the future. Genetic variation declined rapidly, despite the population rebound, and the effective population size was approximately an order of magnitude smaller than census size. The species carried a high genetic load of circa 15 lethal equivalents for longevity. Our computer simulations predicted continued inbreeding will likely result in increased expression of deleterious mutations (i.e., a high realized load) and severe inbreeding depression. Without continued conservation actions, it is likely that the pink pigeon will go extinct in the wild within 100 years. Conservation rescue of the pink pigeon has been instrumental in the recovery of the free‐living population. However, further genetic rescue with captive‐bred birds from zoos is required to recover lost variation, reduce expression of harmful deleterious variation, and prevent extinction. The use of genomics and modeling data can inform IUCN assessments of the viability and extinction risk of species, and it helps in assessments of the conservation dependency of populations.

Optimum contribution selection (OCS) analyses prompted successful conservation actions for Faroese Horse population

The Faroese horse, an endangered indigenous horse breed, is a part of the cultural and societal heritage of the Faroe Islands. Population history describes a severe bottleneck, prompting for quantification of the genetic diversity (level of inbreeding, probability of gene origin, effective population size) and assessment of sustainable conservation potential (Optimum Contribution Selection, OCS) of the Faroese horse population. The pedigree completeness (PCI) of the Faroese horse is adequate for a realistic estimation of the level of inbreeding (PCI5 = 0.96). In concordance with the known population history, the average inbreeding is exceptionally high; in the last cohort, it was equal to 26.8%. An estimate of the effective population size, based on individual increase in inbreeding and coancestry, accounting for the whole population history, was eight. OCS offers a tool to understand and control the increase in the average relationships in the population. Within a fixed number of matings, the repetitive use of stallions resulted in the lowest level of average relationships. Successful follow-up of mating schemes planned together with a holistic assessment of the suitability of an individual as a breeding candidate, will minimize the increase in inbreeding in future generations and maximize the possibility to increase the census size of the Faroese horse population.

Divided by the range: evidence for geographic isolation of the highly mobile Emu (Dromaius novaehollandiae)

ABSTRACT Conserving evolutionary processes is becoming increasingly important in conservation management as environmental changes continue to threaten wild populations. Characterising genetically distinct populations and assessing connectivity across the landscape enables wildlife managers to prioritise conservation efforts with limited resources. In the NSW North Coast bioregion of Australia, one of the last remaining coastal populations of the Emu (Dromaius novaehollandiae) has been state-listed as an Endangered Population, owing to its geographic isolation and small census size. Using mitochondrial and nuclear genetic markers, we examined the spatial genetic structure and diversity of Emus across south-eastern Australia. For the NSW North Coast population, we estimated the effective population size () and carried out simulations to predict future levels of genetic variation. We show significant genetic divergence between the NSW North Coast Emu and other localities based on thousands of highly resolving nuclear markers. Among NSW North Coast Emus, we found less genetic diversity and a critically low-effective population size ( = 14.84 and 22.49 based on independent methods). Our simulations predict that the of the NSW North Coast Emu population is insufficient to maintain genetic diversity and the population may be at risk of inbreeding depression. Incorporating genetic data into the design of captive-release and translocation projects would refine management plans for this locally important population and monitor risks to its long-term survival.

Limited gene flow and pronounced population genetic structure of Eastern Massasauga (Sistrurus catenatus) in a Midwestern prairie remnant.

As anthropogenic changes continue to ecologically stress wildlife, obtaining measures of gene flow and genetic diversity are crucial for evaluating population trends and considering management and conservation strategies for small, imperiled populations. In our study, we conducted a molecular assessment to expand on previous work to elucidate patterns of diversity and connectivity in the remaining disjunct Eastern Massasauga Rattlesnake (Sistrurus catenatus) hibernacula in Illinois. We assayed genetic data for 327 samples collected during 1999-2015 from the Carlyle Lake study area across 21 microsatellite loci. We found hibernacula formed distinct genetic clusters corresponding to the three main study areas (Dam Recreation Areas, Eldon Hazlet State Park, and South Shore State Park). Genetic structuring and low estimates of dispersal indicated that connectivity among these study areas is limited and each is demographically independent. Hibernacula exhibited moderate levels of heterozygosity (0.60-0.73), but estimates of effective population size (5.2-41.0) were low and track census sizes generated via long-term mark-recapture data. Hibernacula at Carlyle Lake, which represent the only Eastern Massasauga remaining in Illinois, are vulnerable to future loss of genetic diversity through lack of gene flow as well as demographic and environmental stochastic processes. Our work highlights the need to include population-level genetic data in recovery planning and suggests that recovery efforts should focus on managing the three major study areas as separate conservation units in order to preserve and maintain long-term adaptive potential of these populations. Specific management goals should include improving connectivity among hibernacula, maintaining existing wet grassland habitat, and minimizing anthropogenic sources of mortality caused by habitat management (e.g., mowing, prescribed fire) and recreational activities. Our molecular study provides additional details about demographic parameters and connectivity at Carlyle Lake that can be used to guide recovery of Eastern Massasauga in Illinois and throughout its range.

Origin and Diversification of the Saguaro Cactus (Carnegiea gigantea): A Within-Species Phylogenomic Analysis.

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these, we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects postglacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects. [Phylogenomics; phylogeography; rooting; Sonoran Desert.].

Pedigree analysis and estimates of effective breeding size characterize sea lamprey reproductive biology.

The sea lamprey (Petromyzon marinus) is an invasive species in the Great Lakes and the focus of a large control and assessment program. Current assessment methods provide information on the census size of spawning adult sea lamprey in a small number of streams, but information characterizing reproductive success of spawning adults is rarely available. We used RAD‐capture sequencing to genotype single nucleotide polymorphism (SNP) loci for ~1600 sea lamprey larvae collected from three streams in northern Michigan (Black Mallard, Pigeon, and Ocqueoc Rivers). Larval genotypes were used to reconstruct family pedigrees, which were combined with Gaussian mixture analyses to identify larval age classes for estimation of spawning population size. Two complementary estimates of effective breeding size (N b), as well as the extrapolated minimum number of spawners (N s), were also generated for each cohort. Reconstructed pedigrees highlighted inaccuracies of cohort assignments from traditionally used mixture analyses. However, combining genotype‐based pedigree information with length‐at‐age assignment of cohort membership greatly improved cohort identification accuracy. Population estimates across all three streams sampled in this study indicate a small number of successfully spawning adults when barriers were in operation, implying that barriers limited adult spawning numbers but were not completely effective at blocking access to spawning habitats. Thus, the large numbers of larvae present in sampled systems were a poor indicator of spawning adult abundance. Overall, pedigree‐based N b and N s estimates provide a promising and rapid assessment tool for sea lamprey and other species.

Anthropogenic fragmentation increases risk of genetic decline in the threatened orchid Platanthera leucophaea.

Protecting biodiversity requires an understanding of how anthropogenic changes impact the genetic processes associated with extinction risk. Studies of the genetic changes due to anthropogenic fragmentation have revealed conflicting results. This is likely due to the difficulty in isolating habitat loss and fragmentation, which can have opposing impacts on genetic parameters. The well‐studied orchid, Platanthera leucophaea, provides a rich dataset to address this issue, allowing us to examine range‐wide genetic changes. Midwestern and Northeastern United States. We sampled 35 populations of P. leucophaea that spanned the species’ range and varied in patch composition, degree of patch isolation, and population size. From these populations we measured genetic parameters associated with increased extinction risk. Using this combined dataset, we modeled landscape variables and population metrics against genetic parameters to determine the best predictors of increased extinction risk. All genetic parameters were strongly associated with population size, while development and patch isolation showed an association with genetic diversity and genetic structure. Genetic diversity was lowest in populations with small census sizes, greater urbanization pressures (habitat loss), and small patch area. All populations showed moderate levels of inbreeding, regardless of size. Contrary to expectation, we found that critically small populations had negative inbreeding values, indicating non‐random mating not typically observed in wild populations, which we attribute to selection for less inbred individuals. The once widespread orchid, Platanthera leucophaea, has suffered drastic declines and extant populations show changes in the genetic parameters associated with increased extinction risk, especially smaller populations. Due to the important correlation with risk and habitat loss, we advocate continued monitoring of population sizes by resource managers, while the critically small populations may need additional management to reverse genetic declines.

Codon Usage Bias in Phytoplankton

Non-random usage of synonymous codons, known as “codon bias”, has been described in many organisms, from bacteria to Drosophila, but little is known about it in phytoplankton. This phenomenon is thought to be driven by selection for translational efficiency. As the efficacy of selection is proportional to the effective population size, species with large population sizes, such as phytoplankton, are expected to have strong codon bias. To test this, we measured codon bias in 215 strains from Haptophyta, Chlorophyta, Ochrophyta (except diatoms that were studied previously), Dinophyta, Cryptophyta, Ciliophora, unicellular Rhodophyta and Chlorarachniophyta. Codon bias is modest in most groups, despite the astronomically large population sizes of marine phytoplankton. The strength of the codon bias, measured with the effective number of codons, is the strongest in Haptophyta and the weakest in Chlorarachniophyta. The optimal codons are GC-ending in most cases, but several shifts to AT-ending codons were observed (mainly in Ochrophyta and Ciliophora). As it takes a long time to reach a new equilibrium after such shifts, species having AT-ending codons show a lower frequency of optimal codons compared to other species. Genetic diversity, calculated for species with more than three strains sequenced, is modest, indicating that the effective population sizes are many orders of magnitude lower than the astronomically large census population sizes, which helps to explain the modest codon bias in marine phytoplankton. This study represents the first comparative analysis of codon bias across multiple major phytoplankton groups.

Forensic and genetic characterizations of diverse southern Thai populations based on 15 autosomal STRs

Southern Thailand is home to various populations; the Moklen, Moken and Urak Lawoi’ sea nomads and Maniq negrito are the minority, while the southern Thai groups (Buddhist and Muslim) are the majority. Although previous studies have generated forensic STR dataset for major groups, such data of the southern Thai minority have not been included; here we generated a regional forensic database of southern Thailand. We newly genotyped common 15 autosomal STRs in 184 unrelated southern Thais, including all minorities and majorities. When combined with previously published data of major southern Thais, this provides a total of 334 southern Thai samples. The forensic parameter results show appropriate values for personal identification and paternity testing; the probability of excluding paternity is 0.99999622, and the combined discrimination power is 0.999999999999999. Probably driven by genetic drift and/or isolation with small census size, we found genetic distinction of the Maniq and sea nomads from the major groups, which were closer to the Malay and central Thais than the other Thai groups. The allelic frequency results can strength the regional forensic database in southern Thailand and also provide useful information for anthropological perspective.