Heterozygosity Estimates Research Articles

Effects of solute carrier family 27 member 1 (SLC27A1) genotypes on fat content and major fatty acids in colostrum and milk from Murrah and “Murrah × Carabao” crossbred buffaloes

This study analyzed the effect of Solute carrier family 27 member 1 (SLC27A1) genotypes on fat content and major fatty acids (lauric acid C12:0, myristic acid C14:0, palmitic acid C16:0, stearic acid C18:0, and oleic acid C18:1n-9) in milk and colostrum of 46 Murrah and “Murrah × Carabao” crossbred buffaloes at the Philippine Carabao Center - University of the Philippines Los Baños dairy herd. The SLC27A1 genotypes (CC and CT) were determined by the polymerase chain reaction-restricted fragment length polymorphism (PCR-RFLP) method using DNA extracted from hair follicles. This study found polymorphism in the SLC27A1 gene that were consistent with Hardy Weinberg’s law of equilibrium (HWE), with a polymorphic information content (PIC) and heterozygosity (H) estimate of 0.2915 and 0.3542, respectively, for Murrah; and 0.3219 and 0.4032, respectively for “Murrah × Carabao” crossbreeds. In Murrah buffaloes, the CT genotype was significantly associated (P<0.05) with higher colostrum yield and milk fat content, compared to CC. In “Murrah × Carabao” crossbreds, CT was associated with higher C12:0, C14:0, and C16:0, but lower C18:1n-9 in colostrum; and lower C16:0 in milk compared to CC. This study showed polymorphisms in SLC27A1 genotypes and their significant effects on colostrum yield and milk fat content in Murrah buffaloes and some major fatty acids in colostrum and milk from “Murrah × Carabao” crosses.

Low-coverage sequencing and Wahlund effect severely bias estimates of inbreeding, heterozygosity and effective population size in North American wolves.

vonHoldt etal. ((2024), Molecular Ecology, 33, e17231) (vH24) used low-coverage (average ~ 7X read depth) restriction site-associated DNA sequence data to estimate individual inbreeding and heterozygosity, and recent effective population size (Ne), in Great Lakes (GL) and Northern Rocky Mountain (RM) wolves. They concluded that RM heterozygosity rapidly declined between 1991 and 2020, and that Ne declined substantially in GL and RM over the last 50 generations. Here, we evaluate the effects of low sequence coverage and sampling strategy on vH24's findings and provide general recommendations for using sequence data to evaluate inbreeding, heterozygosity and Ne. Low-coverage sequencing resulted in downwardly biased estimates of individual inbreeding and heterozygosity, and an erroneous large temporal decline in RM heterozygosity due to declining read depth through time. Additionally, vH24's sampling strategy-which combined individuals from several genetically differentiated populations and across at least eight wolf generations-is expected to have resulted in severe downward bias in estimates of recent Ne for RM. We recommend using high-coverage sequence data ( 15-20X) to estimate inbreeding and heterozygosity. Carefully filtering individuals, loci and genotypes, and using genotype imputation or likelihoods can help to minimise bias when low-coverage sequence data must be used. For estimation of contemporary Ne, the marginal benefits of using more than 103-104 loci are small, so aggressive filtering of loci with low average read depth potentially can retain most individuals without sacrificing much precision. Individuals are relatively more valuable than loci because analyses of contemporary Ne should focus on roughly single-generation samples from local breeding populations.

Estimates of heterozygosity from single nucleotide polymorphism markers are context-dependent and oftenwrong.

Genetic diversity is frequently described using heterozygosity, particularly in a conservation context. Often, it is estimated using single nucleotide polymorphisms (SNPs); however, it has been shown that heterozygosity values calculated from SNPs can be biased by both study design and filtering parameters. Though solutions have been proposed to address these issues, our own work has found them to be inadequate in some circumstances. Here, we aimed to improve the reliability and comparability of heterozygosity estimates, specifically by investigating how sample size and missing data thresholds influenced the calculation of autosomal heterozygosity (heterozygosity calculated from across the genome, i.e. fixed and variable sites). We also explored how the standard practice of tri- and tetra-allelic site exclusion could bias heterozygosity estimates and influence eventual conclusions relating to genetic diversity. Across three distinct taxa (a frog, Litoria rubella; a tree, Eucalyptus microcarpa; and a grasshopper, Keyacris scurra), we found heterozygosity estimates to be meaningfully affected by sample size and missing data thresholds, partly due to the exclusion of tri- and tetra-allelic sites. These biases were inconsistent both between species and populations, with more diverse populations tending to have their estimates more severely affected, thus having potential to dramatically alter interpretations of genetic diversity. We propose a modified framework for calculating heterozygosity that reduces bias and improves the utility of heterozygosity as a measure of genetic diversity, whilst also highlighting the need for existing population genetic pipelines to be adjusted such that tri- and tetra-allelic sites be included in calculations.

Giants among Cnidaria: Large Nuclear Genomes and Rearranged Mitochondrial Genomes in Siphonophores.

Siphonophores (Cnidaria: Hydrozoa) are abundant predators found throughout the ocean and are important constituents of the global zooplankton community. They range in length from a few centimeters to tens of meters. They are gelatinous, fragile, and difficult to collect, so many aspects of the biology of these roughly 200 species remain poorly understood. To survey siphonophore genome diversity, we performed Illumina sequencing of 32 species sampled broadly across the phylogeny. Sequencing depth was sufficient to estimate nuclear genome size from k-mer spectra in six specimens, ranging from 0.7 to 2.3 Gb, with heterozygosity estimates between 0.69% and 2.32%. Incremental k-mer counting indicates k-mer peaks can be absent with nearly 20× read coverage, suggesting minimum genome sizes range from 1.4 to 5.6 Gb in the 25 samples without peaks in the k-mer spectra. This work confirms most siphonophore nuclear genomes are large relative to the genomes of other cnidarians, but also identifies several with reduced size that are tractable targets for future siphonophore nuclear genome assembly projects. We also assembled complete mitochondrial genomes for 33 specimens from these new data, indicating a conserved gene order shared among nonsiphonophore hydrozoans, Cystonectae, and some Physonectae, revealing the ancestral mitochondrial gene order of siphonophores. Our results also suggest extensive rearrangement of mitochondrial genomes within other Physonectae and in Calycophorae. Though siphonophores comprise a small fraction of cnidarian species, this survey greatly expands our understanding of cnidarian genome diversity. This study further illustrates both the importance of deep phylogenetic sampling and the utility of k-mer-based genome skimming in understanding the genomic diversity of a clade.

Analysis of locus estimates of allelic diversity of STR markers in a sample of breeding bulls

Genotypes of the 11 DNA microsatellite loci of 84 bulls of seven breeds were used to evaluate 14 indicators of allelic diversity/differentiation. Traditional and multidimensional statistical methods were applied to the data matrices from the original and transformed estimates (11×14). Estimates of heterozygosity had coefficients of variability of 8-14 %, the number of alleles per locus and indicators of differentiation of breeds at the level of 20-26 %, fixation indices – 38-44 %. Statistically significant Kendall correlations (0.8-1.0) between indicators of allelic richness and heterozygosity, fixation indices, and differentiation indicators were established. The variability of the transformed estimates of diversity/differentiation indicators by loci was in the range of 6-32 %. Including by loci Eth3, Tgla122, Eth225, Bm2113 – 6-12 %, loci Inra23, Tgla126, Eth10 – 15-20 %, loci Tgla227, Sps115, Tgla53, Bm1824 – 28-32 %. The nonparametric Mann-Whitney-Wilcoxon test showed statistically significant differences in the medians of the Eth3 locus with the Bm2113 locus, the Tgla126 locus with the Eth3, Inra23, Tgla122, Eth225, Bm2113, Bm1824, Eth10 loci. The principal component analysis (PCA) identified two components with a total information content of 95,2 %. The first one took into account 59.4 % of the total variance, had the highest loads in intra-breed diversity data and was defined as an «alpha component». The second accounted for 35.8 % of the total variance, had the highest loads in inter-breed differentiation data and was defined as a «beta component». 2D-PCA-ordination showed that a characteristic grouping of loci took place for the analyzed breeds (samples), loci and measures of diversity. Loci Tgla227 and Tgla53 formed group A, group B – loci Tgla122, Eth225, Eth10, group C – loci Inra23, Bm2113 and Bm1824. The loci of the conditional group D (Eth3, Tgla126, Sps115) were defined as «untypical». Validation of ordination was confirmed by calculations on reduced data (dimension 11×7) and the method of non-metric multidimensional scaling (nMDS). The consistency of ordinations according to the Procrust test was 96 % (pperm <0.001). A similar classification of loci was obtained by cluster analysis (UPGMA) with butstrap probabilities of cluster: A – 73, B – 100, C – 73, D – 47 %. The distances and similarity indicators (S) between the profiles of loci and the «true» summary estimates for 11 loci were calculated. Loci Tgla126 and Sps115 had S ≈ 40 %, loci Tgla53 and Bm1824 – at the level of 60 %, loci Inra23, Tgla227 and Bm2113 – 70-75 %, loci Eth3, Tgla122, Eth225 and Eth10 – 84-88 %. The average absolute deviation of the estimates of diversity indicators for the four loci with S≥84 % from the «true» estimates was 3.4 %, for the four loci with S≤60 % – 12.4 %. According to component scores, a general diversity index, γLV, was calculated for each locus. Its correlation with the estimates of the Shannon/Sherwin′s γ-diversity with a 95 % probability value was in the range of 0.73-0.98, Kendall's rank correlation was 0.67 (pvalue = 0.005). The conducted research makes a certain contribution to the expansion of tools for processing molecular genetic data in the analysis of allelic diversity in subdivided populations.

Pangolin Genomes Offer Key Insights and Resources for the World's Most Trafficked Wild Mammals.

Pangolins form a group of scaly mammals that are trafficked at record numbers for their meat and purported medicinal properties. Despite their conservation concern, knowledge of their evolution is limited by a paucity of genomic data. We aim to produce exhaustive genomic resources that include 3,238 orthologous genes and whole-genome polymorphisms to assess the evolution of all eight extant pangolin species. Robust orthologous gene-based phylogenies recovered the monophyly of the three genera and highlighted the existence of an undescribed species closely related to Southeast Asian pangolins. Signatures of middle Miocene admixture between an extinct, possibly European, lineage and the ancestor of Southeast Asian pangolins, provide new insights into the early evolutionary history of the group. Demographic trajectories and genome-wide heterozygosity estimates revealed contrasts between continental versus island populations and species lineages, suggesting that conservation planning should consider intraspecific patterns. With the expected loss of genomic diversity from recent, extensive trafficking not yet realized in pangolins, we recommend that populations be genetically surveyed to anticipate any deleterious impact of the illegal trade. Finally, we produce a complete set of genomic resources that will be integral for future conservation management and forensic endeavors for pangolins, including tracing their illegal trade. These comprise the completion of whole-genomes for pangolins through the hybrid assembly of the first reference genome for the giant pangolin (Smutsia gigantea) and new draft genomes (∼43x-77x) for four additional species, as well as a database of orthologous genes with over 3.4 million polymorphic sites.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa.

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference.

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

Pathogenic Variants Associated with Rare Monogenic Diseases Established in Ancient Neanderthal and Denisovan Genome-Wide Data

Ancient anatomically modern humans (AMHs) encountered other archaic human species, most notably Neanderthals and Denisovans, when they left Africa and spread across Europe and Asia ~60,000 years ago. They interbred with them, and modern human genomes retain DNA inherited from these interbreeding events. High quality (high coverage) ancient human genomes have recently been sequenced allowing for a direct estimation of individual heterozygosity, which has shown that genetic diversity in these archaic human groups was very low, indicating low population sizes. In this study, we analyze ten ancient human genome-wide data, including four sequenced with high-coverage. We screened these ancient genome-wide data for pathogenic mutations associated with monogenic diseases, and established unusual aggregation of pathogenic mutations in individual subjects, including quadruple homozygous cases of pathogenic variants in the PAH gene associated with the condition phenylketonuria in a ~120,000 years old Neanderthal. Such aggregation of pathogenic mutations is extremely rare in contemporary populations, and their existence in ancient humans could be explained by less significant clinical manifestations coupled with small community sizes, leading to higher inbreeding levels. Our results suggest that pathogenic variants associated with rare diseases might be the result of introgression from other archaic human species, and archaic admixture thus could have influenced disease risk in modern humans.

Genetic diversity and population structure of a Peruvian cattle herd using SNP data.

New-generation sequencing technologies, among them SNP chips for massive genotyping, are useful for the effective management of genetic resources. To date, molecular studies in Peruvian cattle are still scarce. For the first time, the genetic diversity and population structure of a reproductive nucleus cattle herd of four commercial breeds from a Peruvian institution were determined. This nucleus comprises Brahman (N = 9), Braunvieh (N = 9), Gyr (N = 5), and Simmental (N = 15) breeds. Additionally, samples from a locally adapted creole cattle, the Arequipa Fighting Bull (AFB, N = 9), were incorporated. Female individuals were genotyped with the GGPBovine100K and males with the BovineHD. Quality control, and the proportion of polymorphic SNPs, minor allele frequency, expected heterozygosity, observed heterozygosity, and inbreeding coefficient were estimated for the five breeds. Admixture, principal component analysis (PCA), and discriminant analysis of principal components (DAPC) were performed. Also, a dendrogram was constructed using the Neighbor-Joining clustering algorithm. The genetic diversity indices in all breeds showed a high proportion of polymorphic SNPs, varying from 51.42% in Gyr to 97.58% in AFB. Also, AFB showed the highest expected heterozygosity estimate (0.41 ± 0.01), while Brahman the lowest (0.33 ± 0.01). Besides, Braunvieh possessed the highest observed heterozygosity (0.43 ± 0.01), while Brahman the lowest (0.37 ± 0.02), indicating that Brahman was less diverse. According to the molecular variance analysis, 75.71% of the variance occurs within individuals, whereas 24.29% occurs among populations. The pairwise genetic differentiation estimates (FST) between breeds showed values that ranged from 0.08 (Braunvieh vs. AFB) to 0.37 (Brahman vs. Braunvieh). Similarly, pairwise Reynold's distance ranged from 0.09 (Braunvieh vs. AFB) to 0.46 (Brahman vs. Braunvieh). The dendrogram, similar to the PCA, identified two groups, showing a clear separation between Bos indicus (Brahman and Gyr) and B. taurus breeds (Braunvieh, Simmental, and AFB). Simmental and Braunvieh grouped closely with the AFB cattle. Similar results were obtained for the population structure analysis with K = 2. The results from this study would contribute to the appropriate management, avoiding loss of genetic variability in these breeds and for future improvements in this nucleus. Additional work is needed to speed up the breeding process in the Peruvian cattle system.

Pipeline to Design Inbred Lines and F1 Hybrids of Leaf Chicory (Radicchio) Using Male Sterility and Genotyping-by-Sequencing.

Chicory, a horticultural crop cultivated worldwide, presents many botanical varieties and local biotypes. Among these, cultivars of the Italian radicchio group of the pure species Cichorium intybus L. and its interspecific hybrids with Cichorium endivia L.-as the "Red of Chioggia" biotype-includes several phenotypes. This study uses a pipeline to address the marker-assisted breeding of F1 hybrids: it presents the genotyping-by-sequencing results of four elite inbred lines using a RADseq approach and an original molecular assay based on CAPS markers for screening mutants with nuclear male sterility in the radicchio of Chioggia. A total of 2953 SNP-carrying RADtags were identified and used to compute the actual estimates of homozygosity and overall genetic similarity and uniformity of the populations, as well as to determine their genetic distinctiveness and differentiation. Molecular data were further used to investigate the genomic distribution of the RADtags among the two Cichorium species, allowing their mapping in 1131 and 1071 coding sequences in chicory and endive, respectively. Paralleling this, an assay to screen the genotype at the male sterility locus Cims-1 was developed to discriminate wild-type and mutant alleles of the causative gene myb80-like. Moreover, a RADtag mapped close to this genomic region proved the potential application of this method for future marker-assisted selection tools. Finally, after combining the genotype information of the core collection, the best 10 individuals from each inbred line were selected to compute the observed genetic similarity as a measure of uniformity as well as the expected homozygosity and heterozygosity estimates scorable by the putative progenies derived from selfing (pollen parent) and full-sibling (seed parent) or pair-wise crossing (F1 hybrids). This predictive approach was conducted as a pilot study to understand the potential application of RADseq in the fine tuning of molecular marker-assisted breeding strategies aimed at the development of inbred lines and F1 hybrids in leaf chicory.

An improved reference genome for Trifolium subterraneum L. provides insight into molecular diversity and intra-specific phylogeny.

Subterranean clover (Trifolium subterraneum L., Ts) is a geocarpic, self-fertile annual forage legume with a compact diploid genome (n = x = 8, 544 Mb/1C). Its resilience and climate adaptivity have made it an economically important species in Mediterranean and temperate zones. Using the cultivar Daliak, we generated higher resolution sequence data, created a new genome assembly TSUd_3.0, and conducted molecular diversity analysis for copy number variant (CNV) and single-nucleotide polymorphism (SNP) among 36 cultivars. TSUd_3.0 substantively improves prior genome assemblies with new Hi-C and long-read sequence data, covering 531 Mb, containing 41,979 annotated genes and generating a 94.4% BUSCO score. Comparative genomic analysis among select members of the tribe Trifolieae indicated TSUd 3.0 corrects six assembly-error inversion/duplications and confirmed phylogenetic relationships. Its synteny with T. pratense, T. repens, Medicago truncatula and Lotus japonicus genomes were assessed, with the more distantly related T. repens and M. truncatula showing higher levels of co-linearity with Ts than between Ts and its close relative T. pratense. Resequencing of 36 cultivars discovered 7,789,537 SNPs subsequently used for genomic diversity assessment and sequence-based clustering. Heterozygosity estimates ranged from 1% to 21% within the 36 cultivars and may be influenced by admixture. Phylogenetic analysis supported subspecific genetic structure, although it indicates four or five groups, rather than the three recognized subspecies. Furthermore, there were incidences where cultivars characterized as belonging to a particular subspecies clustered with another subspecies when using genomic data. These outcomes suggest that further investigation of Ts sub-specific classification using molecular and morpho-physiological data is needed to clarify these relationships. This upgraded reference genome, complemented with comprehensive sequence diversity analysis of 36 cultivars, provides a platform for future gene functional analysis of key traits, and genome-based breeding strategies for climate adaptation and agronomic performance. Pangenome analysis, more in-depth intra-specific phylogenomic analysis using the Ts core collection, and functional genetic and genomic studies are needed to further augment knowledge of Trifolium genomes.

Evaluation of Inbred Lines and Their Tester Hybrids and Estimation of Heterozygosity for Some Quantitative Traits in Maize

Nine inbred lines of maize were sowing :(OH, ZP- 301, unn440, Inbreed12, Zp- 607, SH, Ik-58, IK-8 and R-153). The first five was used as a tester (males) and the reminder used as line (females) during at 2014 season (spring) according to method (Line X Tester) to obtain 20 crosses, while the evaluation season was in Shiwan province / Kirkuk governorate of parents and its crosses and the commercial variety (Moring) at spring of 2015 using Random Complete Block Design (RCBD) design with three replicates. The studied traits were: (No. of days to 50% male flowering and female flowering, plant height, ear leaf area, No. of the rows ear-1, No. of the grains row-1, No. of the grains ear-1, 100 grains weight and grain yield plant-1). The most results were summarized as follows: The genotypes (parents, crosses and commercial variety) showed significant differences at 1% level probe for all studied except days to 50% female flowering for parents. Crosses gave highly (Inbred12 X IK-58) for ear leaf area (681.06 cm) and the cross (R-153 X Inbred 12) for both No. of the grains row-1 and grain yield plant-1 (760.23 grain) and (175.457gm) respectively, the cross (Ik-58 x OH) for 100 grains weight (30.387gm). The crosses (IK-58 X OH), (SH X Zp-301), (Zp-301 X Ik-8), (Ik-8 X Un-44052), (IK-8 X Inbred12), (R-153 X Inbred12) and (SH X Zp-607) gave desired significant heterosis over the better parents as well as the commercial variety for most traits included the grain yield plant-1. The parents SH, Zp-607 and Zp-301 to be used in single cross production and synthetic varieties, also the crosses (SH X OH), (SH X Un-44052), and (R-153 X Un-44052) in double cross production and compound crosses because these hybrids put in different main and secondary groups through cluster analysis.

Transitioning from microsatellites to SNP-based microhaplotypes in genetic monitoring programmes: Lessons from paired data spanning 20 years.

Many long-term genetic monitoring programmes began before next-generation sequencing became widely available. Older programmes can now transition to new marker systems usually consisting of 1000s of SNP loci, but there are still important questions about comparability, precision, and accuracy of key metrics estimated using SNPs. Ideally, transitioned programmes should capitalize on new information without sacrificing continuity of inference across the time series. We combined existing microsatellite-based genetic monitoring information with SNP-based microhaplotypes obtained from archived samples of Rio Grande silvery minnow (Hybognathus amarus) across a 20-year time series to evaluate point estimates and trajectories of key genetic metrics. Demographic and genetic monitoring bracketed multiple collapses of the wild population and included cases where captive-born repatriates comprised the majority of spawners in the wild. Even with smaller sample sizes, microhaplotypes yielded comparable and in some cases more precise estimates of variance genetic effective population size, multilocus heterozygosity and inbreeding compared to microsatellites because many more microhaplotype loci were available. Microhaplotypes also recorded shifts in allele frequencies associated with population bottlenecks. Trends in microhaplotype-based inbreeding metrics were associated with the fraction of hatchery-reared repatriates to the wild and should be incorporated into future genomic monitoring. Although differences in accuracy and precision of some metrics were observed between marker types, biological inferences and management recommendations were consistent.

Diploid-dominant life cycles characterize the early evolution of Fungi

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

Mating behavior and reproductive morphology predict macroevolution of sex allocation in hermaphroditic flatworms

BackgroundSex allocation is the distribution of resources to male or female reproduction. In hermaphrodites, this concerns an individual’s resource allocation to, for example, the production of male or female gametes. Macroevolutionary studies across hermaphroditic plants have revealed that the self-pollination rate and the pollination mode are strong predictors of sex allocation. Consequently, we expect similar factors such as the selfing rate and aspects of the reproductive biology, like the mating behaviour and the intensity of postcopulatory sexual selection, to predict sex allocation in hermaphroditic animals. However, comparative work on hermaphroditic animals is limited. Here, we study sex allocation in 120 species of the hermaphroditic free-living flatworm genus Macrostomum. We ask how hypodermic insemination, a convergently evolved mating behaviour where sperm are traumatically injected through the partner’s epidermis, affects the evolution of sex allocation. We also test the commonly-made assumption that investment into male and female reproduction should trade-off. Finally, we ask if morphological indicators of the intensity of postcopulatory sexual selection (female genital complexity, male copulatory organ length, and sperm length) can predict sex allocation.ResultsWe find that the repeated evolution of hypodermic insemination predicts a more female-biased sex allocation (i.e., a relative shift towards female allocation). Moreover, transcriptome-based estimates of heterozygosity reveal reduced heterozygosity in hypodermically mating species, indicating that this mating behavior is linked to increased selfing or biparental inbreeding. Therefore, hypodermic insemination could represent a selfing syndrome. Furthermore, across the genus, allocation to male and female gametes is negatively related, and larger species have a more female-biased sex allocation. Finally, increased female genital complexity, longer sperm, and a longer male copulatory organ predict a more male-biased sex allocation.ConclusionsSelfing syndromes have repeatedly originated in plants. Remarkably, this macroevolutionary pattern is replicated in Macrostomum flatworms and linked to repeated shifts in reproductive behavior. We also find a trade-off between male and female reproduction, a fundamental assumption of most theories of sex allocation. Beyond that, no theory predicts a more female-biased allocation in larger species, suggesting avenues for future work. Finally, morphological indicators of more intense postcopulatory sexual selection appear to predict more intense sperm competition.

Correlations between microsatellite DNA heterozygosity and reproductive traits in Large White sows

The main aim of this paper was to determine whether heterozygosity (assessed using microsatellite genotypes) was correlated with the reproductive traits in sows. The study was conducted on two herds of sows of the Large White sows breed at the Limited Liability Company “Tavriys’ki svyni” (Kherson region, Ukraine) and the Agricultural Private Enterprise “Techmet-Yug” (Mykolayiv region, Ukraine). During the study, we used eleven microsatellite loci recommended by International Society for Animal Genetics (ISAG) – S0101, S0155, S0228, S0355, S0386, Sw24, Sw72, Sw240, Sw857, Sw936 and Sw951. The litter records included information on the total number of piglets born (TNB), number of piglets born alive (NBA), number of stillborn piglets (NSB), frequency of stillborn piglets (FSB), litter size at weaning (NW) in the first five parities. Individual heterozygosity estimates (for each microsatellite loci separately) and microsatellite multilocus heterozygosity (for all used loci) estimates (MLH) were used in our analysis. ANOVA was used to determine the relationship of the dependent effects (reproductive traits) to single locus heterozygosity using two classes: 0 (for homozygous individuals) and 1 (for heterozygous individuals). In addition, the following indicators were calculated for each genotype: the squared distances (d2) between alleles within an individual for each microsatellite loci and mean squared distances (mean d2) between alleles within an individual for 11 microsatellite loci. Spearman’s rank correlation coefficients were used to measure the association between d2 (for each microsatellite loci) and reproductive traits in sows. ANOVA on reproductive traits of sows belonging to different MLH and mean d2 classes was also undertaken. For sows from the Agricultural Private Enterprise “‘Techmet-Yug” were observed negative associations between heterozygosity and reproductive traits. We conclude that care should be taken when crossing between different breeders (English and Hungarian selection) to avoid outbreeding depression.

Genome-Wide Investigation of the Multiple Origins Hypothesis for Deep-Spawning Kokanee Salmon (Oncorhynchus nerka) across its Pan-Pacific Distribution

Salmonids have emerged as important study systems for investigating molecular processes underlying parallel evolution given their tremendous life history variation. Kokanee, the resident form of anadromous sockeye salmon (Oncorhynchus nerka), have evolved multiple times across the species' pan-Pacific distribution, exhibiting multiple reproductive ecotypes including those that spawn in streams, on lake-shores, and at lake depths >50 m. The latter has only been detected in 5 locations in Japan and British Columbia, Canada. Here, we investigated the multiple origins hypothesis for deep-spawning kokanee, using 9721 single nucleotide polymorphisms distributed across the genome analyzed for the vast majority of known populations in Japan (Saiko Lake) and Canada (Anderson, Seton, East Barrière Lakes) relative to stream-spawning populations in both regions. We detected 397 outlier loci, none of which were robustly identified in paired-ecotype comparisons in Japan and Canada independently. Bayesian clustering and principal components analyses based on neutral loci revealed 6 distinct clusters, largely associated with geography or translocation history, rather than ecotype. Moreover, a high level of divergence between Canadian and Japanese populations, and between deep- and stream-spawning populations regionally, suggests the deep-spawning ecotype independently evolved on the 2 continents. On a finer level, Japanese kokanee populations exhibited low estimates of heterozygosity, significant levels of inbreeding, and reduced effective population sizes relative to Canadian populations, likely associated with transplantation history. Along with preliminary evidence for hybridization between deep- and stream-spawning ecotypes in Saiko Lake, these findings should be considered within the context of on-going kokanee fisheries management in Japan.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō.

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Unraveling the Genetic Diversity and Population Structure of Bangladeshi Indigenous Cattle Populations Using 50K SNP Markers.

Simple SummaryIndigenous cattle have extraordinary adaptation capability to diverse environments under low input production system. However, the population size is declining rapidly in Bangladesh due to massive imports of high yielding dairy breeds. The genetic diversity measures are important for assessing population architecture as well as for development of conservation strategies. The aim of this study was to investigate genetic variability and population structure of indigenous cattle genetic resources of Bangladesh using Illumina Bovine SNP50K BeadChip genotyped data. Similar to other zebu populations, low genetic diversity measures were found in Bangladeshi cattle populations. Our findings revealed their distinct genetic structure but showed low levels of genetic differentiation among the six indigenous cattle populations. Moreover, admixture and phylogenetic analysis highlighted historical gene flow among the studied populations. Altogether, our findings provide a comprehensive genomic information on indigenous cattle populations of Bangladesh that could be utilized in their future conservation and breeding research.Understanding the genetic basis of locally adapted indigenous cattle populations is essential to design appropriate strategies and programs for their genetic improvement and conservation. Here, we report genetic diversity measures, population differentiation, and structure of 218 animals sampled from six indicine cattle populations of Bangladesh. Animals were genotyped with Illumina Bovine SNP50K BeadChip along with genotyped data of 505 individuals included from 19 zebu and taurine breeds worldwide. The principal component analysis (PCA) showed clear geographic separation between taurine and indicine lineages where Bangladeshi indigenous cattle clustered with South Asian zebu populations. However, overlapped clusters in PCA, heterozygosity estimates, and Neighbor-Joining phylogenetic tree analysis revealed weak genetic differentiation among the indigenous cattle populations of Bangladesh. The admixture analysis at K = 5 and 9 suggests distinct genetic structure of the studied populations along with 1 to 4% of taurine ancestry. The effective population size suggested a limited pool of ancestors particularly for Sahiwal and North Bengal Grey cattle. In conclusion, these findings shed insights into the genetic architecture of six indigenous cattle populations of Bangladesh for the first time and suggested as distinct gene pools without potential admixture with zebu or taurine populations.