Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet. 2011;12:499-510. [PMID: 21681211 DOI: 10.1038/nrg3012]

Plant Neighbourhood as Predictor of Allopatric Speciation in Sedges

Adaptation to different environments in geographically separated populations is key for allopatric speciation. Most research has focused on the effects of geographical isolation and abiotic factors, but disjunct populations frequently co-occur with different pools of species, favouring divergent adaptation and speciation. We show the importance of plant neighbourhood, compared to geographic and environmental factors, in the allopatric speciation of two closely related plants, Carex elata and Carex reuteriana. Both species share similar ecological requirements and inhabit river shores at medium to low altitudes across the Iberian Peninsula. We employed a multidisciplinary approach integrating abiotic, biotic and geographical factors, and genomic data (genotyping-by-sequencing) to infer the relative role of different evolutionary drivers. Abiotic factors were assessed based on 38 bioclimatic variables, biotic factors using the community of co-occurring plant species (1536 vegetation inventories), and geographical factors with a distance matrix based on geographic coordinates. Using regularised generalised linear models, we identified the key variables explaining distribution patterns. We also examined the relationships between inter-population genetic distances, and biotic, abiotic and geographic factors to understand the drivers of lineage splitting, revealing varying degrees of influence. Plant neighbourhood emerged as a stronger predictor of allopatric distributions than abiotic or geographic factors, with the largest effect observed in Carex elata, which exhibited the greatest population differentiation. These findings suggest that the biotic and microenvironmental factors influencing divergent plant neighbourhoods have significantly contributed to the differentiation of these taxa, providing new insights into the evolutionary processes shaping the origin and distribution of species.

Phylogenomics resolves the puzzling phylogeny of banded newts (genus Ommatotriton)

Resolving the order of speciation events that occurred in rapid succession is inherently hard and typically requires a phylogenomic approach. A case in point concerns the previously unresolved phylogeny of the three species of banded newt (genus Ommatotriton). We obtain c. 7k nuclear DNA markers using target enrichment by sequence capture and analyze the dataset using maximum likelihood inference of concatenated data with RAxML, summary multi-species coalescent analysis with ASTRAL and Bayesian species tree inference using a diffusion model with SNAPPER, and use TreeMix and PhyloNet to test for interspecific gene flow. All analyses recover three distinct species with no evidence of interspecific gene flow. All analyses retrieved the topology (O. nesterovi, (O. ophryticus, O. vittatus)), with high support. SNAPPER did show the tendency to get stuck in a local optimum, resulting in a different but still highly supported topology. Furthermore, we notice that fewer SNAPPER runs get stuck in a local optimum when we include an outgroup. Therefore, we recommend the exploration of multiple independent runs and the use of an outgroup with this approach. The banded newt radiation illustrates the use of genome-wide data to tackle formerly unresolved phylogenies.

Clinical exome sequencing for carrier screening in assisted reproductive technology and sperm donation

PURPOSE

To assess the efficacy of clinical exome sequencing (CES) in individuals involved in assisted reproductive technology (ART) or sperm donor programs, with a specific focus on its impact on clinical decision-making.

METHODS

A total of 3991 individuals without a family history of genetic disorders underwent CES targeting 5595 genes at a reproductive center between December 2022 and April 2024. The cohort comprised 217 sperm donors, 232 female recipients, and 1771 couples (3542 patients) undergoing ART with their own gametes. At-risk couples (ARCs) were identified when both partners had a pathogenic or likely pathogenic variant (P/LP) in the same autosomal recessive gene or X-linked variants in females. The analysis primarily examined carrier frequencies, reproductive choices, and outcomes of ARCs.

RESULTS

Among the 3991 individuals screened, 3895 (97.6%) were found to carry at least one P/LP variant, with an average carrier burden of 3.8 variants per individual, showing no significant disparity in carrier status between individuals with infertility and sperm donors/recipients. Within the screened couples, 9.3% were identified as ARCs and 2.3% opted for preimplantation genetic testing for monogenic diseases (PGT-M). As of now, 31 ARCs proceeded with the transfer of euploid and unaffected blastocysts, resulting in 8 healthy live births and 13 ongoing pregnancies.

CONCLUSION

The findings reveal a significant prevalence of carrier status for autosomal recessive and X-linked diseases, irrespective of fertility status. This highlights the critical importance of integrating genetic risk counseling and informed reproductive decision-making into fertility clinics and sperm banks.

Combining genotyping approaches improves resolution for association mapping: a case study in tropical maize under water stress conditions

Genome-wide Association Studies (GWAS) identify genome variations related to specific phenotypes using Single Nucleotide Polymorphism (SNP) markers. Genotyping platforms like SNP-Array or sequencing-based techniques (GBS) can genotype samples with many SNPs. These approaches may bias tropical maize analyses due to reliance on the temperate line B73 as the reference genome. An alternative is a simulated genome called "Mock," adapted to the population using bioinformatics. Recent studies show SNP-Array, GBS, and Mock yield similar results for population structure, heterotic groups definition, tester selection, and genomic hybrid prediction. However, no studies have examined the results generated by these different genotyping approaches for GWAS. This study aims to test the equivalence among the three genotyping scenarios in identifying significant effect genes in GWAS. To achieve this, maize was used as the model species, where SNP-Array genotyped 360 inbred lines from a public panel via the Affymetrix platform and GBS. The GBS data were used to perform SNP calling using the temperate inbred line B73 as the reference genome (GBS-B73) and a simulated genome "Mock" obtained in-silico (GBS-Mock). The study encompassed four above-ground traits with plants grown under two levels of water supply: well-watered (WW) and water-stressed (WS). In total, 46, 34, and 31 SNP were identified in the SNP-Array, GBS-B73, and GBS-Mock scenarios, respectively, across the two water levels, associated with the evaluated traits following the comparative analysis of each genotyping method individually. Overall, the identified candidate genes varied along the various scenarios but had the same functionality. Regarding SNP-Array and GBS-B73, genes with functional similarity were identified even without coincidence in the physical position of the SNPs. These genes and regions are involved in various processes and responses with applications in plant breeding. In terms of accuracy, the combination of genotyping scenarios compared to those isolated is feasible and recommended, as it increased all traits under both water conditions. In this sense, it is worth highlighting the combination of GBS-B73 and GBS-Mock scenarios, not only due to the increase in the resolution of GWAS results but also the reduction of costs associated with genotyping and the possibility of conducting genomic breeding methods.

Genomic insights into Thermomonas hydrothermalis: potential applications in industrial biotechnology

Thermomonas hydrothermalis, a thermophilic bacterium isolated from hot springs, exhibits unique genomic features that underpin its adaptability to extreme environments and its potential in industrial biotechnology. In this study, we present a comparative genomic analysis of two strains, DSM 14834 and HOT.CON.106, revealing distinct metabolic pathways and stress response mechanisms. The genome annotation highlighted strain-specific variations, such as enhanced motility and chemotaxis capabilities in HOT.CON.106 and a stronger genomic stability emphasis in DSM 14834. Comparative analysis with other Thermomonas species demonstrated that T. hydrothermalis possesses a unique genomic architecture, including genes for thermostable enzymes (e.g., amylases and pullulanases) and secondary metabolite biosynthesis. These enzymes and metabolites have significant industrial potential in high-temperature processes such as bioenergy production, bioplastics synthesis, and bioremediation. The findings underscore the relative differentiation between the strains and their broader implications for sustainable biotechnology, offering a basis for further exploration of thermophilic microorganisms in industrial applications.

Advances in Whole Genome Sequencing: Methods, Tools, and Applications in Population Genomics

With the rapid advancement of high-throughput sequencing technologies, whole genome sequencing (WGS) has emerged as a crucial tool for studying genetic variation and population structure. Utilizing population genomics tools to analyze resequencing data allows for the effective integration of selection signals with population history, precise estimation of effective population size, historical population trends, and structural insights, along with the identification of specific genetic loci and variations. This paper reviews current whole genome sequencing technologies, detailing primary research methods, relevant software, and their advantages and limitations within population genomics. The goal is to examine the application and progress of resequencing technologies in this field and to consider future developments, including deep learning models and machine learning algorithms, which promise to enhance analytical methodologies and drive further advancements in population genomics.

Genomic selection for resistance to one pathogenic strain of Vibrio splendidus in blue mussel Mytilus edulis

Introduction

The blue mussel is one of the major aquaculture species worldwide. In France, this species faces a significant threat from infectious disease outbreaks in both mussel farms and the natural environment over the past decade. Diseases caused by various pathogens, particularly Vibrio spp., have posed a significant challenge to the mussel industry. Genetic improvement of disease resistance can be an effective approach to overcoming this issue.

Methods

In this work, we tested genomic selection in the blue mussel (Mytilus edulis) to understand the genetic basis of resistance to one pathogenic strain of Vibrio splendidus (strain 14/053 2T1) and to predict the accuracy of selection using both pedigree and genomic information. Additionally, we performed a genome-wide association study (GWAS) to identify putative QTLs underlying disease resistance. We conducted an experimental infection involving 2,280 mussels sampled from 24 half-sib families containing each two full-sib families which were injected with V. splendidus. Dead and survivor mussels were all sampled, and among them, 348 dead and 348 surviving mussels were genotyped using a recently published multi-species medium-density 60K SNP array.

Results

From potentially 23.5K SNPs for M. edulis present on the array, we identified 3,406 high-quality SNPs, out of which 2,204 SNPs were successfully mapped onto the recently published reference genome. Heritability for resistance to V. splendidus was moderate ranging from 0.22 to 0.31 for a pedigree-based model and from 0.28 to 0.36 for a genomic-based model.

Discussion

GWAS revealed the polygenic architecture of the resistance trait in the blue mussel. The genomic selection models studied showed overall better performance than the pedigree-based model in terms of accuracy of breeding values prediction. This work provides insights into the genetic basis of resistance to V. splendidus and exemplifies the potential of genomic selection in family-based breeding programs in M. edulis.

Unveiling the Genetic Diversity and Demographic History of Coffea stenophylla in Sierra Leone Using Genotyping-By-Sequencing

Coffea stenophylla is a rare Coffea species boasting a flavor profile comparable to Arabica coffee (Coffea arabica) and has a good adaptability to lowland tropical climates. This species faces increasing threats from climate change, deforestation, and habitat fragmentation in its West African homeland. Using 1037 novel SNP markers derived from Genotyping-by-Sequencing (GBS), we revealed the presence of three distinct natural populations (mean Fst = 0.176) in Sierra Leone. Evidence of recent bottlenecks and small effective population size (118-140) was found across all three populations, reflecting the impact of recent anthropogenic disturbances on this species. Using a model-flexible inference approach, we unveiled a strong ancient bottleneck approximately 23,000 years ago, coinciding with the last glacial maximum (LGM), followed by post-glacial expansion and divergence into distinct genetic clusters. A comparative analysis between ex situ genebanks and natural populations detected a significant gap in genetic diversity, with two out of three natural populations missing from the ex situ genebank collection. These findings highlight the urgent need to improve conservation practices for C. stenophylla in Sierra Leone. The novel SNP markers developed in this study provided valuable tools to support future efforts in conservation and utilization of C. stenophylla genetic resources in West Africa.

Promises and challenges of crop translational genomics

Crop translational genomics applies breeding techniques based on genomic datasets to improve crops. Technological breakthroughs in the past ten years have made it possible to sequence the genomes of increasing numbers of crop varieties and have assisted in the genetic dissection of crop performance. However, translating research findings to breeding applications remains challenging. Here we review recent progress and future prospects for crop translational genomics in bringing results from the laboratory to the field. Genetic mapping, genomic selection and sequence-assisted characterization and deployment of plant genetic resources utilize rapid genotyping of large populations. These approaches have all had an impact on breeding for qualitative traits, where single genes with large phenotypic effects exert their influence. Characterization of the complex genetic architectures that underlie quantitative traits such as yield and flowering time, especially in newly domesticated crops, will require further basic research, including research into regulation and interactions of genes and the integration of genomic approaches and high-throughput phenotyping, before targeted interventions can be designed. Future priorities for translation include supporting genomics-assisted breeding in low-income countries and adaptation of crops to changing environments.

Development of Aegilops comosa and Aegilops caudata-specific molecular markers and fluorescence in situ hybridization probes based on specific-locus amplified fragment sequencing

As tertiary gene pools of wheat, Aegilops comosa and Ae. caudata contain many excellent genes/traits and gradually become important and noteworthy wild resources for wheat improvement worldwide. However, the lack of molecular markers and cytological probes with good specificity and high sensitivity limits the development and utilization of Triticum aestivum-Ae. comosa (Ta. Aeco)/Ae. caudata (Ta. Aeca) introgression lines. Using specific-locus amplified fragment sequencing, two Ae. comosa and one Ae. caudata accessions, Chinese Spring, and three Ta. Aeco and Ta. Aeca introgression lines each were sequenced to develop new molecular markers and cytological probes. After strict sequence comparison and verification in different materials, a total of 39 molecular markers specific to three chromosomes in Ae. comosa (nine, seven, and 10 markers for 1M, 2M, and 7M, respectively) and Ae. caudata (two, six, and five markers for 3C, 4C, and 5C, respectively) and 21 fluorescence in situ hybridization (FISH) probes (one centromeric probe with signals specific to the M chromosomes, two centromeric probes with signals in all the tested genomes, and six, eight, and four FISH probes specific to the M, C, and M, C, and U chromosomes, respectively) were successfully exploited. The newly developed molecular markers and cytological probes could be used in karyotype studies, centromere evolutionary analyses of Aegilops, and had the ability to detect the fusion centromeres and small-fragment translocations in introgression lines.

Persistent Gene Flow Suggests an Absence of Reproductive Isolation in an African Antelope Speciation Model

African antelope diversity is a globally unique vestige of a much richer world-wide Pleistocene megafauna. Despite this, the evolutionary processes leading to the prolific radiation of African antelopes are not well understood. Here, we sequenced 145 whole genomes from both subspecies of the waterbuck (Kobus ellipsiprymnus), an African antelope believed to be in the process of speciation. We investigated genetic structure and population divergence and found evidence of a mid-Pleistocene separation on either side of the eastern Great Rift Valley, consistent with vicariance caused by a rain shadow along the so-called "Kingdon's Line." However, we also found pervasive evidence of both recent and widespread historical gene flow across the Rift Valley barrier. By inferring the genome-wide landscape of variation among subspecies, we found 14 genomic regions of elevated differentiation, including a locus that may be related to each subspecies' distinctive coat pigmentation pattern. We investigated these regions as candidate speciation islands. However, we observed no significant reduction in gene flow in these regions, nor any indications of selection against hybrids. Altogether, these results suggest a pattern whereby climatically driven vicariance is the most important process driving the African antelope radiation and suggest that reproductive isolation may not set in until very late in the divergence process. This has a significant impact on taxonomic inference, as many taxa will be in a gray area of ambiguous systematic status, possibly explaining why it has been hard to achieve consensus regarding the species status of many African antelopes. Our analyses demonstrate how population genetics based on low-depth whole genome sequencing can provide new insights that can help resolve how far lineages have gone along the path to speciation.

High Crimean-Congo hemorrhagic fever incidence linked to greater genetic diversity and differentiation in Hyalomma marginatum populations in Türkiye

BACKGROUND

Ticks are crucial vectors of a wide range of pathogens, posing significant threats to human and animal health globally. Understanding the genetic basis of tick biology and host-parasite interactions is essential for developing effective control programs. This study investigates the fine-scale genetic structure of Hyalomma marginatum Koch, 1844, the primary vector of Crimean-Congo hemorrhagic fever (CCHF) in Türkiye. Despite its significant public health importance, information regarding its population structure and genetic diversity is quite limited.

METHODS

We used restriction site-associated DNA sequencing (RAD-Seq) to obtain genome-wide sequence data from 10 tick populations in Türkiye, collected from regions with low, moderate, and high incidence rates of CCHF. Based on these data, we determined population structure and diversity of populations using principal component analysis (PCA) and admixture analysis. Furthermore, we calculated pairwise FST and utilized discriminant analysis of principal components (DAPC) to understand genetic differentiation between populations.

RESULTS

PCA and admixture analysis indicated minimal genetic structure between populations, but we detected notable genetic differentiation and high genetic diversity from regions with high CCHF rates. Furthermore, our DAPC identified 31 significant single-nucleotide polymorphisms (SNPs) associated with regions with high CCHF incidence, with 25 SNPs located near genes involved in critical biological functions such as nucleic acid binding, transmembrane transport, and proteolysis. These findings suggest that genetic variations in these regions may confer adaptive advantages in environments with high pathogen loads.

CONCLUSIONS

This study provides the first comprehensive analysis of H. marginatum genetic diversity in Türkiye, revealing significant differentiation in populations from CCHF-endemic regions. These results underscore the importance of considering fine-scale genetic diversity to fully understand the drivers of genetic variation in ticks and their implications for vectorial capacity.

Genetic diversity assessment of cucumber landraces using molecular signatures

Genetic profiling of the biodiversity in cultivated crop plants is necessary to preserve important genes and utilize them in a breeding program. Cucumber is used as a model plant to study various characteristics of Cucurbitaceae. Its adaptation to a wide range of climatic conditions suggested analyzing the landraces. The present study was conducted to evaluate the differences, at the genetic level, among landraces spanning five continents. DNA extracted from fifty-six landraces selected from USDA germplasm bank to cover a global representative sample of world cucumber landraces was used for polymerase chain reaction using twenty-eight polymorphic expressed sequence tags simple sequence repeat (EST-SSR) markers. Twenty-eight EST-SSR markers covering all seven chromosomes yielded 98 bands with an average of 3.42 bands per marker. Polymorphic information content ranged from 0.00 (EC35) to 0.74 (EC17) with an average of 0.34. Six clusters provided an appropriate summary of the variation among the landraces, with the two largest groups including 32 (Asiatic) and 17 (European and American) landraces, respectively. Four small groups, three with two members, and one with one member (PI 525155-Egypt) were dissimilar to the two main groups. Landraces from the same region were often clustered together. Genetic similarity of the landraces was revealed by marker banding patterns. The locations of genetic diversity for cucumber landraces can be identified from this study.

Induced and natural variation affect traits independently in hybrid Populus

The genetic control of many plant traits can be highly complex. Both allelic variation (sequence change) and dosage variation (copy number change) contribute to a plant's phenotype. While numerous studies have investigated the effect of allelic or dosage variation, very few have documented both within the same system, leaving their relative contribution to phenotypic effects unclear. The Populus genome is highly polymorphic, and poplars are fairly tolerant of gene dosage variation. Here, using a previously established Populus hybrid F1 population, we assessed and compared the effect of natural allelic variation and induced dosage variation on biomass, phenology, and leaf morphology traits. We identified QTLs for many of these traits, but our results indicate limited overlap between the QTLs associated with natural allelic variation and induced dosage variation. Additionally, the integration of data from both allelic and dosage variation identifies a larger set of QTLs that together explain a larger percentage of the phenotypic variance. Finally, our results suggest that the effect of the large indels might mask that of allelic QTLs. Our study helps clarify the relationship between allelic and dosage variation and their effects on quantitative traits.

Genotyping by sequencing; a strategy for identification and mapping of induced mutation in newly developed wheat mutant lines

Exposing genetic material with physical mutagens can create novel genetic resources capable of combating different stresses. High throughput GBS-DArTseq™ assay was deployed to estimate genetic diversity of 33 newly developed stable wheat mutants as compared to the wild type. The identified 1,57,608 PAVs markers were randomly distributed across wheat chromosomes and sub-genomes with the highest number detected on Chr-7D (2877) and Chr-7B (2711). The B sub-genome contained the most PAVs followed by D and A-sub genome. Among mutant lines, Pb-M-2061 and Pb-M-59 had the highest PAV count, while Pb-M-605 and Pb-M-196 had the lowest. A total of 7,910 PAVs were consistently present over all replicates, with 3,252 specifically present in mutants and absent in wild type. The maximum PAVs (1480) were found in Pb-M-1027 and Pb-M-1323 (656). Functional characterization revealed that out of 3,252, 1,238 were found in wheat transcriptome database that contained 152 characterized and 1,196 uncharacterized genes. COGs and GO-terms analysis linked many PAVs with pathways involving signaling, metabolism and defense. Maximum number of gene-containing PAVs were identified in Pb-M-1027, Pb-M-2302 and Pb-M-1323 which were involved in tolerance to diseases and abiotic stresses, improved photosynthetic efficiency, larger grain size, increased grain yield and harvest index pathways. This study provides valuable insights into the genetic diversity and potential agronomic benefits of PAVs in wheat mutant lines. These findings can help molecular geneticist and breeders for exploiting the induced genetic diversity for unravelling the genetic circuits as well as exploiting in wheat breeding for developing resilient cultivars.

Internet-Based Abnormal Chromosomal Diagnosis During Pregnancy Using a Noninvasive Innovative Approach to Detecting Chromosomal Abnormalities in the Fetus: Scoping Review

BACKGROUND

Chromosomal abnormalities are genetic disorders caused by chromosome errors, leading to developmental delays, birth defects, and miscarriages. Currently, invasive procedures such as amniocentesis or chorionic villus sampling are mostly used, which carry a risk of miscarriage. This has led to the need for a noninvasive and innovative approach to detect and prevent chromosomal abnormalities during pregnancy.

OBJECTIVE

This review aims to describe and appraise the potential of internet-based abnormal chromosomal preventive measures as a noninvasive approach to detecting and preventing chromosomal abnormalities during pregnancy.

METHODS

A thorough review of existing literature and research on chromosomal abnormalities and noninvasive approaches to prenatal diagnosis and therapy was conducted. Electronic databases such as PubMed, Google Scholar, ScienceDirect, CENTRAL, CINAHL, Embase, OVID MEDLINE, OVID PsycINFO, Scopus, ACM, and IEEE Xplore were searched for relevant studies and articles published in the last 5 years. The keywords used included chromosomal abnormalities, prenatal diagnosis, noninvasive, and internet-based, and diagnosis.

RESULTS

The review of literature revealed that internet-based abnormal chromosomal diagnosis is a potential noninvasive approach to detecting and preventing chromosomal abnormalities during pregnancy. This innovative approach involves the use of advanced technology, including high-resolution ultrasound, cell-free DNA testing, and bioinformatics, to analyze fetal DNA from maternal blood samples. It allows early detection of chromosomal abnormalities, enabling timely interventions and treatment to prevent adverse outcomes. Furthermore, with the advancement of technology, internet-based abnormal chromosomal diagnosis has emerged as a safe alternative with benefits including its cost-effectiveness, increased accessibility and convenience, potential for earlier detection and intervention, and ethical considerations.

CONCLUSIONS

Internet-based abnormal chromosomal diagnosis has the potential to revolutionize prenatal care by offering a safe and noninvasive alternative to invasive procedures. It has the potential to improve the detection of chromosomal abnormalities, leading to better pregnancy outcomes and reduced risk of miscarriage. Further research and development in this field is needed to make this approach more accessible and affordable for pregnant women.

Genomic Analysis of Aedes aegypti in the Northern Chihuahuan Desert of Texas and Mexico

Background: Aedes aegypti, is the primary vector of dengue, Chikungunya, Zika, and yellow fever viruses. Both natural and human-impacted landscapes have selective pressures on Ae. aegypti, resulting in strong genomic structure even within close geographical distances. Materials and Methods: We assess the genetic structure of this medically important mosquito species at the northern leading edge of their distribution in Southwestern USA. Ae. aegypti were collected during 2017 in the urban communities of El Paso and Sparks, Texas (USA) and in the city of Ciudad Juárez, Mexico. Results: Thousands of nuclear loci were sequenced across 260 captured Ae. aegypti. First, we recovered the genetic structure of Ae. aegypti following geography, with all four major collection communities being genetically distinct. Importantly, we found population structure and genetic diversity that suggest rapid expansion through active-short distance dispersals, with Anapra being the likely source for the others. Next, tests of selection recovered eight functional genes across six outliers: calmodulin with olfactory receptor function; the protein superfamily C-type lectin with function in mosquito immune system and development; and TATA box binding protein with function in gene regulation. Conclusion: Despite these populations being documented in the early 2000s, we find that selective pressures on specific genes have already occurred and likely facilitate Ae. aegypti range expansion.

Streamlined whole-genome genotyping through NGS-enhanced thermal asymmetric interlaced (TAIL)-PCR

Whole-genome genotyping (WGG) stands as a pivotal element in genomic-assisted plant breeding. Nevertheless, sequencing-based approaches for WGG continue to be costly, primarily owing to the high expenses associated with library preparation and the laborious protocol. During prior development of foreground and background integrated genotyping by sequencing (FBI-seq), we discovered that any sequence-specific primer (SP) inherently possesses the capability to amplify a massive array of stable and reproducible non-specific PCR products across the genome. Here, we further improved FBI-seq by replacing the adapter ligated by Tn5 transposase with an arbitrary degenerate (AD) primer. The protocol for the enhanced FBI-seq unexpectedly mirrors a simplified thermal asymmetric interlaced (TAIL)-PCR, a technique that is widely used for isolation of flanking sequences. However, the improved TAIL-PCR maximizes the primer-template mismatched annealing capabilities of both SP and AD primers. In addition, leveraging of next-generation sequencing enhances the ability of this technique to assay tens of thousands of genome-wide loci for any species. This cost-effective, user-friendly, and powerful WGG tool, which we have named TAIL-PCR by sequencing (TAIL-peq), holds great potential for widespread application in breeding programs, thereby facilitating genome-assisted crop improvement.

Identification of Insertion and Deletion (InDel) Markers for Chickpea (Cicer arietinum L.) Based on Double-Digest Restriction Site-Associated DNA Sequencing

Enhancing the marker repository and the development of breeder-friendly markers in chickpeas is important in relation to chickpea genomics-assisted breeding applications. Insertion-deletion (InDel) markers are widely distributed across genomes and easily observed with specifically designed primers, leading to less time, cost, and labor requirements. In light of this, the present study focused on the identification and development of InDel markers through the use of double-digest restriction site-associated DNA sequencing (ddRADSeq) data from 20 chickpea accessions. Bioinformatic analysis identified 20,700 InDel sites, including 15,031 (72.61%) deletions and 5669 (27.39%) insertions, among the chickpea accessions. The InDel markers ranged from 1 to 25 bp in length, while single-nucleotide-length InDel markers were found to represent the majority of the InDel sites and account for 79% of the total InDel markers. However, we focused on InDel markers wherein the length was greater than a single nucleotide to avoid any read or alignment errors. Among all of the InDel markers, 96.1% were less than 10 bp, 3.6% were between 10 and 20 bp, and 0.3% were more than 20 bp in length. We examined the InDel markers that were 10 bp and longer for the development of InDel markers based on a consideration of the genomic distribution and low-cost genotyping with agarose gels. A total of 29 InDel regions were selected, and primers were successfully designed to evaluate their efficiency. Annotation analysis of the InDel markers revealed them to be found with the highest frequency in the intergenic regions (82.76%), followed by the introns (6.90%), coding sequences (6.90%), and exons (3.45%). Genetic diversity analysis demonstrated that the polymorphic information content of the markers varied from 0.09 to 0.37, with an average of 0.20. Taken together, these results showed the efficiency of InDel marker development for chickpea genetic and genomic studies using the ddRADSeq method. The identified markers might prove valuable for chickpea breeders.

Conservation genomic study of Hopea hainanensis (Dipterocarpaceae), an endangered tree with extremely small populations on Hainan Island, China

Introduction

Hopea hainanensis Merrill & Chun is considered a keystone and indicator species in the tropical lowland rainforests of Hainan Island. Owing to its high-quality timber, H. hainanensis has been heavily exploited, leading to its classification as a first-class national protected plant in China and a plant species with extremely small populations (PSESPs).

Methods

This study analyzed genome-wide single nucleotide polymorphisms obtained through restriction site-associated DNA sequencing from 78 adult trees across 10 H. hainanensis populations on Hainan Island.

Results and discussion

The nucleotide diversity of the sampled populations ranged from 0.00096 to 0.00138, which is lower than that observed in several other PSESPs and endangered tree species. Bayesian unsupervised clustering, principal component analysis, and neighbor-joining tree reconstruction identified three to five genetic clusters in H. hainanensis, most of which were geographically widespread and shared by multiple populations. Demographic history analysis based on pooled samples indicated that the decline in the H. hainanensis population began approximately 20,000 years ago, starting from an ancestral population size of approximately 10,000 individuals. The reduction in population size accelerated approximately 4,000 years ago and has continued to the present, resulting in a severely reduced population on Hainan Island. Intensified genetic drift in small and isolated H. hainanensis populations may contribute to moderate differentiation between some of them, as revealed by pairwise F st. In conclusion, our conservation genomic study confirms a severe population decline and an extremely low level of nucleotide variation in H. hainanensis on Hainan Island. These findings provide critical insights for the sustainable management and genetic restoration of H. hainanensis on Hainan Island.

Genome Sequencing of Lentinula edodes Revealed a Genomic Variant Block Associated with a Thermo-Tolerant Trait in Fruit Body Formation

The formation of multicellular fruiting bodies in basidiomycete mushrooms is a crucial developmental process for sexual reproduction and subsequent spore development. Temperature is one of the most critical factors influencing the phase transition for mushroom reproduction. During the domestication of mushrooms, traits related to fruiting bodies have significantly impacted agricultural adaptation and human preferences. Recent research has demonstrated that chromosomal variations, such as structural variants (SVs) and variant blocks (VBs), play crucial roles in agronomic traits and evolutionary processes. However, the lack of high-quality genomic information and important trait data have hindered comprehensive identification and characterization in Lentinula edodes breeding processes. In this study, the genomes of two monokaryotic L. edodes strains, characterized by thermo-tolerance and thermo-sensitivity during fruiting body formation, were reassembled at the chromosomal level. Comparative genomic studies of four thermo-tolerant and thermo-sensitive monokaryotic L. edodes strains identified a 0.56 Mbp variant block on chromosome 9. Genes associated with DNA repair or cellular response to DNA damage stimulus were enriched in this variant block. Finally, we developed eight CAPS markers from the variant block to discriminate the thermo-tolerant traits in L. edodes cultivars. Our findings show that the identified variant block is highly correlated with the thermo-tolerant trait for fruiting body formation and that alleles present in this block may have been artificially selected during L. edodes domestication.

Sequencing vs. amplification for the estimation of allele dosages in sugarcane (Saccharum spp.)

Premise

Detecting single-nucleotide polymorphisms (SNPs) in a cost-effective way is fundamental in any plant breeding pipeline. Here, we compare three genotyping techniques for their ability to reproduce the allele dosage of SNPs of interest in sugarcane (Saccharum spp.).

Methods

To identify a reproducible technique to estimate allele dosage for the validation of SNP markers, the correlation between Flex-Seq, kompetitive allele-specific PCR (KASP), and genotyping-by-sequencing and restriction site-associated DNA sequencing (GBS+RADseq) was determined for a set of 76 SNPs. To find alternative methodologies for allele dosage estimation, the KASP and Flex-Seq techniques were compared for the same set of SNPs. For the three techniques, a population of 53 genotypes from the diverse sugarcane panel of the Centro de Investigación de la Caña de Azúcar (Cenicaña), Colombia, was selected.

Results

The average Pearson correlation coefficients between GBS+RADseq and Flex-Seq, GBS+RADseq and KASP, and Flex-Seq and KASP were 0.62 ± 0.27, 0.38 ± 0.27, and 0.38 ± 0.30, respectively.

Discussion

Flex-Seq reproduced the allele dosages determined using GBS+RADseq with good levels of precision because of its depth of sequencing and ability to target specific positions in the genome. Additionally, Flex-Seq outperformed KASP by allowing the conversion of a higher number of SNPs and a more accurate estimation of the allele dosage. Flex-Seq has therefore become the genotyping methodology of choice for marker validation at Cenicaña.

Congruent patterns of cryptic cladogenesis revealed using RADseq and Sanger sequencing in a velvet worm species complex (Onychophora: Peripatopsidae: Peripatopsis sedgwicki)

In the present study, first generation DNA sequencing (mitochondrial cytochrome c oxidase subunit one, COI) and reduced-representative genomic RADseq data were used to understand the patterns and processes of diversification of the velvet worm, Peripatopsis sedgwicki species complex across its distribution range in South Africa. For the RADseq data, three datasets (two primary and one supplementary) were generated corresponding to 1,259-11,468 SNPs, in order to assess the diversity and phylogeography of the species complex. Tree topologies for the two primary datasets were inferred using maximum likelihood and Bayesian inferences methods. Phylogenetic analyses using the COI datasets retrieved four distinct, well-supported clades within the species complex. Five species delimitation methods applied to the COI data (ASAP, bPTP, bGMYC, STACEY and iBPP) all showed support for the distinction of the Fort Fordyce Nature Reserve specimens. In the main P. sedgwicki species complex, the species delimitation methods revealed a variable number of operational taxonomic units and overestimated the number of putative taxa. Divergence time estimates coupled with the geographic exclusivity of species and phylogeographic results suggest recent cladogenesis during the Plio/Pleistocene. The RADseq data were subjected to a principal components analysis and a discriminant analysis of principal components, under a maximum-likelihood framework. The latter results corroborate the four main clades observed using the COI data, however, applying additional filtering revealed additional diversity. The high overall congruence observed between the RADseq data and COI data suggest that first generation sequence data remain a cheap and effective method for evolutionary studies, although RADseq does provide a far greater resolution of contemporary temporo-spatial patterns.

Genomic divergence and mutation load in the Begonia masoniana complex from limestone karsts

Understanding genome-wide diversity, inbreeding, and the burden of accumulated deleterious mutations in small and isolated populations is essential for predicting and enhancing population persistence and resilience. However, these effects are rarely studied in limestone karst plants. Here, we re-sequenced the nuclear genomes of 62 individuals of the Begonia masoniana complex (B. liuyanii, B. longgangensis, B. masoniana and B. variegata) and investigated genomic divergence and genetic load for these four species. Our analyses revealed four distinct clusters corresponding to each species within the complex. Notably, there was only limited admixture between B. liuyanii and B. longgangensis occurring in overlapping geographic regions. All species experienced historical bottlenecks during the Pleistocene, which were likely caused by glacial climate fluctuations. We detected an asymmetric historical gene flow between group pairs within this timeframe, highlighting a distinctive pattern of interspecific divergence attributable to karst geographic isolation. We found that isolated populations of B. masoniana have limited gene flow, the smallest recent population size, the highest inbreeding coefficients, and the greatest accumulation of recessive deleterious mutations. These findings underscore the urgency to prioritize conservation efforts for these isolated population. This study is among the first to disentangle the genetic differentiation and specific demographic history of karst Begonia plants at the whole-genome level, shedding light on the potential risks associated with the accumulation of deleterious mutations over generations of inbreeding. Moreover, our findings may facilitate conservation planning by providing critical baseline genetic data and a better understanding of the historical events that have shaped current population structure of rare and endangered karst plants.

New insights into mycotoxin risk management through fungal population genetics and genomics

Mycotoxin contamination of food and feed is a major global concern. Chronic or acute dietary exposure to contaminated food and feed can negatively affect both human and animal health. Contamination occurs through plant infection by toxigenic fungi, primarily Aspergillus and Fusarium spp., either before or after harvest. Despite the application of various management strategies, controlling these pathogens remains a major challenge primarily because of their ability to adapt to environmental changes and selection pressures. Understanding the genetic structure of plant pathogen populations is pivotal for gaining new insights into their biology and epidemiology, as well as for understanding the mechanisms behind their adaptability. Such deeper understanding is crucial for developing effective and preemptive management strategies tailored to the evolving nature of pathogenic populations. This review focuses on the population-level variations within the two most economically significant toxigenic fungal genera according to space, host, and pathogenicity. Outcomes in terms of migration patterns, gene flow within populations, mating abilities, and the potential for host jumps are examined. We also discuss effective yet often underutilized applications of population genetics and genomics to address practical challenges in the epidemiology and disease control of toxigenic fungi.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

A chromosome-level genome assembly provides insights into the local adaptation of Tamarix austromongolica in the Yellow River Basin, China

Tamarix austromongolica is endemic to the Yellow River Basin and has adapted to diverse ecological settings in the region, including the arid areas of northwestern China and the saline soil regions of the Yellow River Delta. However, the genetic basis of its local adaptation remains unclear. We report a chromosome-level assembly of the T. austromongolica genome based on PacBio high-fidelity sequencing and Hi-C technology. The 12 pseudochromosomes cover 98.44% of the 1.32 Gb assembly, with a contig N50 of 52.57 Mb and a BUSCO score of 98.2%. The genome comprises 913.6 Mb (68.83%) of repetitive sequences and 22,374 protein-coding genes. Genome evolution analyses suggest that genes under positive selection and significantly expanded gene families have facilitated T. austromongolica's adaptability to diverse environmental factors and high resistance to diseases. Using genotyping-by-sequencing, we conducted population structure and selection analyses of 114 samples from 15 sites. Two genetic groups were identified, and 114 and 289 candidate genes were assigned to the populations of the northwestern and eastern parts of the Yellow River, respectively. Furthermore, we discovered numerous candidate genes associated with high-altitude adaptability and salt tolerance. This research provides valuable genomic resources for the evolutionary study and genetic breeding of tamarisk.

Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers

Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.

Quantitative Trait Loci Mapping and Comparative Transcriptome Analysis of Fruit Weight (FW) in Watermelon (Citrullus lanatus L.)

The watermelon (Citrullus lanatus L.) holds substantial economic value as a globally cultivated horticultural crop. However, the genetic architecture of watermelon fruit weight (FW) remains poorly understood. In this study, we used sh14-11 with small fruit and N14 with big fruit to construct 100 recombinant inbred lines (RILs). Based on whole-genome resequencing (WGR), 218,127 single nucleotide polymorphisms (SNPs) were detected to construct a high-quality genetic map. After quantitative trait loci (QTL) mapping, a candidate interval of 31-38 Mb on chromosome 2 was identified for FW. Simultaneously, the bulked segregant analysis (BSA) in the F2 population corroborated the identification of the same interval, encompassing the homologous gene linked to the known FW-related gene fas. Additionally, RNA-seq was carried out across 11 tissues from sh14-11 and N14, revealing expression profiles that identified 1695 new genes and corrected the annotation of 2941 genes. Subsequent differential expression analysis unveiled 8969 differentially expressed genes (DEGs), with 354 of these genes exhibiting significant differences across four key developmental stages. The integration of QTL mapping and differential expression analysis facilitated the identification of 14 FW-related genes, including annotated TGA and NAC transcription factors implicated in fruit development. This combined approach offers valuable insights into the genetic basis of FW, providing crucial resources for enhancing watermelon cultivation.

Phylogenetic relationships and genetic diversity of the Korean endemic Phedimus latiovalifolius (Crassulaceae) and its close relatives

Phedimus latiovalifolius (Y.N.Lee) D.C.Son & H.J.Kim is exclusively distributed in the high mountains in the Korean Peninsula, mainly along the Baekdudaegan mountain range. Despite its morphological and distributional distinction from other Phedimus Raf. species, its taxonomic identity and phylogenetic relationship with congeneric species remain unclear. This study employs genotyping-by-sequencing-derived genome-wide single nucleotide polymorphisms to establish the monophyly of P. latiovalifolius and its relationship with closely related species. Genetic diversity and population differentiation of P. latiovalifolius are also assessed to provide baseline genetic information for future conservation and management strategies. Our phylogenetic analyses robustly demonstrate the monophyletic nature of P. latiovalifolius, with P. aizoon (L.) 't Hart identified as its closest sister lineage. There is no genetic evidence supporting a hybrid origin of P. latiovalifolius from P. aizoon involving either P. ellacombeanus (Praeger) 't Hart or P. kamtschaticus (Fisch.) 't Hart. Population genetic analyses reveal two major groups within P. latiovalifolius. A higher genetic variation is observed in P. ellacombeanus than in the congeneric species. Notably, most of the genetic variation exists within P. latiovalifolius populations. Given its distribution and the potential role of Baekdudaegan as an East Asian Pleistocene refugia, P. latiovalifolius could be considered rare and endemic, persisting in the refugium across glacial/interglacial cycles.

Polyploidy and hybridization in the Mediterranean: unravelling the evolutionary history of Centaurium (Gentianaceae)

BACKGROUND AND AIMS

Polyploidy is considered one of the main mechanisms of plant evolution and speciation. In the Mediterranean Basin, polyploidy has contributed to making this region a biodiversity hotspot, along with its geological and climatic history and other ecological and biogeographical factors. The Mediterranean genus Centaurium (Gentianaceae) comprises ~25 species, of which 60 % are polyploids, including tetraploids and hexaploids. To date, the evolutionary history of centauries has been studied using Sanger sequencing phylogenies, which have been insufficient to fully understand the phylogenetic relationships in this lineage. The goal of this study is to gain a better understanding of the evolutionary history of Centaurium by exploring the mechanisms that have driven its diversification, specifically hybridization and polyploidy. We aim to identify the parentage of hybrid species, at the species or clade level, as well as assessing whether morphological traits are associated with particular ploidy levels.

METHODS

We sequenced RADseq markers from 42 samples of 28 Centaurium taxa, and performed phylogenomic analyses using maximum likelihood, summary coalescent SVDquartets and Neighbor-Net approaches. To identify hybrid taxa, we used PhyloNetworks and the fastSTRUCTURE algorithm. To infer the putative parental species of the allopolyploids, we employed genomic analyses (SNIPloid). The association between different traits and particular ploidy levels was explored with non-metric multidimensional scaling.

KEY RESULTS

Our phylogenetic analyses confirmed the long-suspected occurrence of recurrent hybridization. The allopolyploid origin of the tetraploid C. serpentinicola and the hexaploids C. mairei, C. malzacianum and C. centaurioides was also confirmed, unlike that of C. discolor. We inferred additional signatures of hybridization events within the genus and identified morphological traits differentially distributed in different ploidy levels.

CONCLUSIONS

This study highlights the important role that hybridization has played in the evolution of a Mediterranean genus such as Centaurium, leading to a polyploid complex, which facilitated its diversification and may exemplify that of other Mediterranean groups.

Different reference genomes determine different results: Comparing SNP calling in RAD-seq of Engelhardia roxburghiana using different reference genomes

Advances in next-generation sequencing (NGS) have significantly reduced the cost and improved the efficiency of obtaining single nucleotide polymorphism (SNP) markers, particularly through restriction site-associated DNA sequencing (RAD-seq). Meanwhile, the progression in whole genome sequencing has led to the utilization of an increasing number of reference genomes in SNP calling processes. This study utilized RAD-seq data from 242 individuals of Engelhardia roxburghiana, a tropical tree of the walnut family (Juglandaceae), with SNP calling conducted using the STACKS pipeline. We aimed to compare both reference-based approaches, namely, employing a closely related species as the reference genome versus the species itself as the reference genome, to evaluate their respective merits and limitations. Our findings indicate a substantial discrepancy in the number of obtained SNPs between using a closely related species as opposed to the species itself as reference genomes, the former yielded approximately an order of magnitude fewer SNPs compared to the latter. While the missing rate of individuals and sites of the final SNPs obtained in the two scenarios showed no significant difference. The results showed that using the reference genome of the species itself tends to be prioritized in RAD-seq studies. However, if this is unavailable, considering closely related genomes is feasible due to their wide applicability and low missing rate as alternatives. This study contributes to enrich the understanding of the impact of SNP acquisition when utilizing different reference genomes.

A roadmap of phylogenomic methods for studying polyploid plant genera

Phylogenetic inference of polyploid species is the first step towards understanding their patterns of diversification. In this paper, we review the challenges and limitations of inferring species relationships of polyploid plants using traditional phylogenetic sequencing approaches, as well as the mischaracterization of the species tree from single or multiple gene trees. We provide a roadmap to infer interspecific relationships among polyploid lineages by comparing and evaluating the application of current phylogenetic, phylogenomic, transcriptomic, and whole-genome approaches using different sequencing platforms. For polyploid species tree reconstruction, we assess the following criteria: (1) the amount of prior information or tools required to capture the genetic region(s) of interest; (2) the probability of recovering homeologs for polyploid species; and (3) the time efficiency of downstream data analysis. Moreover, we discuss bioinformatic pipelines that can reconstruct networks of polyploid species relationships. In summary, although current phylogenomic approaches have improved our understanding of reticulate species relationships in polyploid-rich genera, the difficulties of recovering reliable orthologous genes and sorting all homeologous copies for allopolyploids remain a challenge. In the future, assembled long-read sequencing data will assist the recovery and identification of multiple gene copies, which can be particularly useful for reconstructing the multiple independent origins of polyploids.

SATIN: a micro and mini satellite mining tool of total genome and coding regions with analysis of perfect repeats polymorphism in coding regions

BACKGROUND

Tandem repeats are specific sequences in genomic DNA repeated in tandem that are present in all organisms. Among the subcategories of TRs we have Satellite repeats, that is divided into macrosatellites, minisatellites, and microsatellites, being the last two of specific interest because they can identify polymorphisms between organisms due to their instability. Currently, most mining tools focus on Simple Sequence Repeats (SSR) mining, and only a few can identify SSRs in the coding regions.

RESULTS

We developed a microsatellite mining software called SATIN (Micro and Mini SATellite IdentificatioN tool) based on a new sliding window algorithm written in C and Python. It represents a new approach to SSR mining by addressing the limitations of existing tools, particularly in coding region SSR mining. SATIN is available at https://github.com/labgm/SATIN.git . It was shown to be the second fastest for perfect and compound SSR mining. It can identify SSRs from coding regions plus SSRs with motif sizes bigger than 6. Besides the SSR mining, SATIN can also analyze SSRs polymorphism on coding-regions from pre-determined groups, and identify SSRs differentially abundant among them on a per-gene basis. To validate, we analyzed SSRs from two groups of Escherichia coli (K12 and O157) and compared the results with 5 known SSRs from coding regions. SATIN identified all 5 SSRs from 237 genes with at least one SSR on it.

CONCLUSIONS

The SATIN is a novel microsatellite search software that utilizes an innovative sliding window technique based on a numerical list for repeat region search to identify perfect, and composite SSRs while generating comprehensible and analyzable outputs. It is a tool capable of using files in fasta or GenBank format as input for microsatellite mining, also being able to identify SSRs present in coding regions for GenBank files. In conclusion, we expect SATIN to help identify potential SSRs to be used as genetic markers.

Robustness in population-structure and demographic-inference results derived from the Aedes aegypti genotyping chip and whole-genome sequencing data

The mosquito Aedes aegypti is the primary vector of many human arboviruses such as dengue, yellow fever, chikungunya, and Zika, which affect millions of people worldwide. Population genetic studies on this mosquito have been important in understanding its invasion pathways and success as a vector of human disease. The Axiom aegypti1 SNP chip was developed from a sample of geographically diverse A. aegypti populations to facilitate genomic studies on this species. We evaluate the utility of the Axiom aegypti1 SNP chip for population genetics and compare it with a low-depth shotgun sequencing approach using mosquitoes from the native (Africa) and invasive ranges (outside Africa). These analyses indicate that results from the SNP chip are highly reproducible and have a higher sensitivity to capture alternative alleles than a low-coverage whole-genome sequencing approach. Although the SNP chip suffers from ascertainment bias, results from population structure, ancestry, demographic, and phylogenetic analyses using the SNP chip were congruent with those derived from low-coverage whole-genome sequencing, and consistent with previous reports on Africa and outside Africa populations using microsatellites. More importantly, we identified a subset of SNPs that can be reliably used to generate merged databases, opening the door to combined analyses. We conclude that the Axiom aegypti1 SNP chip is a convenient, more accurate, low-cost alternative to low-depth whole-genome sequencing for population genetic studies of A. aegypti that do not rely on full allelic frequency spectra. Whole-genome sequencing and SNP chip data can be easily merged, extending the usefulness of both approaches.

Genome-wide association study for yield and quality of granulated cassava processed product

The starchy storage roots of cassava are commonly processed into a variety of products, including cassava granulated processed products (gari). The commercial value of cassava roots depends on the yield and quality of processed products, directly influencing the acceptance of new varieties by farmers, processors, and consumers. This study aims to estimate genetic advance through phenotypic selection and identify genomic regions associated and candidate genes linked with gari yield and quality. Higher single nucleotide polymorphism (SNP)-based heritability estimates compared to broad-sense heritability estimates were observed for most traits highlighting the influence of genetic factors on observed variation. Using genome-wide association analysis of 188 clones, genotyped using 53,150 genome-wide SNPs, nine SNPs located on seven chromosomes were significantly associated with peel loss, gari yield, color parameters for gari and eba, bulk density, swelling index, and textural properties of eba. Future research will focus on validating and understanding the functions of identified genes and their influence on gari yield and quality traits.

Housekeeping protein-coding genes interrogated with tissue and individual variations

Housekeeping protein-coding genes are stably expressed genes in cells and tissues that are thought to be engaged in fundamental cellular biological functions. They are often utilized as normalization references in molecular biology research and are especially important in integrated bioinformatic investigations. Prior studies have examined human housekeeping protein-coding genes by analyzing various gene expression datasets. The inclusion of different tissue types significantly impacted the discovery of housekeeping genes. In this report, we investigated particularly individual human subject expression differences in protein-coding genes across different tissue types. We used GTEx V8 gene expression datasets obtained from more than 16,000 human normal tissue samples. Furthermore, the Gini index is utilized to investigate the expression variations of protein-coding genes between tissue and individual donor subjects. Housekeeping protein-coding genes found using Gini index profiles may vary depending on the tissue subtypes investigated, particularly given the diverse sample size collections across the GTEx tissue subtypes. We subsequently selected major tissues and identified subsets of housekeeping genes with stable expression levels among human donors within those tissues. In this work, we provide alternative sets of housekeeping protein-coding genes that show more consistent expression patterns in human subjects across major solid organs. Weblink: https://hpsv.ibms.sinica.edu.tw .

High-density genetic map and quantitative trait loci map of skin color in hawthorn (Crataegus pinnatifida bge. Var. major N.E.Br.)

Fruit skin color is an important trait of the hawthorn tree, which has an important influence on fruit quality. Crataegus pinnatifida Bge. var. Major N.E.Br. Is one of the most widely cultivated varieties in China and has a long history of medicinal use. In recent years, it has attracted the attention of the world due to its nutritional and medicinal values. Skin color is the focus of breeders and food processors. At present, skin color-related genes have still not been mapped. In this study, "Shandong Da Mianqiu" (♀, red skin color), "Da Huang Mianzha" (♂, yellow skin color) and 131 F1 hybrids were used to construct genetic map of hawthorn by RAD-seq, and QTL mapping was performed by combining these features with the hue angle and the observed color. In this study, 13,260 SNP was assigned to 17 linkage groups, with an integrated map covering 2,297.75 cM was constructed. A total of 5 QTLs related to hawthorn skin color were detected on LG1, LG3 and LG15. Whether hue angle or pericarp color acts as phenotype for QTL mapping, the candidate genes include bHLH086, WD repeat regions and Myb-like. bHLH, WD and Myb play an important role in the color regulation of Hawthorn skin color. These results lay a solid foundation for QTL mapping and molecular marker-assisted breeding of hawthorn.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

A novel approach of kinship determination based on the physical length of genetically shared regions of chromosomes

BACKGROUND

Determination of genetic relatedness between individuals plays a crucial role in resolving numerous civil cases involving familial relationships and in forensic investigation concerning missing persons. Short tandem repeats (STRs), known for their high degree of DNA polymorphism, have traditionally been the primary choice of DNA markers in genetic testing, but their application for kinships testing is limited to cases involving close kinship. SNPs have emerged as promising supplementary markers for kinship determination. Nevertheless, the challenging remains in discriminating between third-degree or more distant relatives, such as first cousins, using SNPs.

OBJECTIVE

To investigate a kinship analysis method for distant degree of familial relationships using high-density SNP data.

METHODS

A high-density SNP data from 337 individuals of Korean families using Affymetrix Axiom KORV1.0-96 Array was obtained for this study. SNPs were aligned by chromosomal positions, and identity-by-state (IBS) was determined, and then shared regions as consecutive SNPs with IBS of 1 or 2 were investigated. The physical lengths of these IBS segments were measured and summed them to create an Index, as a measure of kinship.

RESULTS

The kinship was determined by the physical length of shared chromosomal regions that are distinguished by each kinship. Using this method, the relationship was able be distinguished up to the fourth degree of kinship, and non-relatives were clearly distinguished from true relatives. We also found a potential for this approach to be used universally, regardless of microarray platforms for SNP genotyping and populations.

CONCLUSION

This method has a potential to determine the different degree of kinship between individuals and to distinguish non-relatives from true relatives, which can be of great help for practical applications in kinship determination.

Deciphering variations, identification of marker-trait associations and candidate genes for seed oil content under terminal heat stress in Indian mustard (Brassica juncea L. Czern & Coss) germplasm stock

This research paper investigates the variability in seed oil content (SOC) in Indian mustard (Brassica juncea L.) under terminal heat stress (THS) conditions. A genetic stock of 488 genotypes of B. juncea was evaluated over two years and grouped into five classes based on the reduction in oil content under THS compared to normal sown crop. Based on heat susceptibility index (HSI), a diverse panel of 96 genotypes was selected and evaluated under THS. Twenty-two heat-tolerant donor genotypes were identified, including introgression lines derived from B. tournefortii, B. carinata and Erucastrum cardaminoides. This study is the first to report on marker-trait associations for SOC in B. juncea under THS using a GWAS approach. Furthermore, candidate genes associated with abiotic stress tolerance and lipid metabolism were identified near the significant SNPs, emphasizing their role in SOC regulation under stress. Notable candidate genes include BjuA003240 (encoding for alcohol-forming fatty acyl-CoA reductase), BjuA003242 (involving in lipid biosynthesis), BjuA003244 (associated with mitochondrial functions and stress tolerance), and BjuA003245 (related to MYB transcription factors regulating lipid biosynthesis). This study provides valuable insights into the genetic basis of SOC variation under THS in B. juncea, highlighting potential breeding targets for improved heat stress resilience in Indian mustard cultivation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03985-w.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Effect of genotyping errors on linkage map construction based on repeated chip analysis of two recombinant inbred line populations in wheat (Triticum aestivum L.)

Linkage maps are essential for genetic mapping of phenotypic traits, gene map-based cloning, and marker-assisted selection in breeding applications. Construction of a high-quality saturated map requires high-quality genotypic data on a large number of molecular markers. Errors in genotyping cannot be completely avoided, no matter what platform is used. When genotyping error reaches a threshold level, it will seriously affect the accuracy of the constructed map and the reliability of consequent genetic studies. In this study, repeated genotyping of two recombinant inbred line (RIL) populations derived from crosses Yangxiaomai × Zhongyou 9507 and Jingshuang 16 × Bainong 64 was used to investigate the effect of genotyping errors on linkage map construction. Inconsistent data points between the two replications were regarded as genotyping errors, which were classified into three types. Genotyping errors were treated as missing values, and therefore the non-erroneous data set was generated. Firstly, linkage maps were constructed using the two replicates as well as the non-erroneous data set. Secondly, error correction methods implemented in software packages QTL IciMapping (EC) and Genotype-Corrector (GC) were applied to the two replicates. Linkage maps were therefore constructed based on the corrected genotypes and then compared with those from the non-erroneous data set. Simulation study was performed by considering different levels of genotyping errors to investigate the impact of errors and the accuracy of error correction methods. Results indicated that map length and marker order differed among the two replicates and the non-erroneous data sets in both RIL populations. For both actual and simulated populations, map length was expanded as the increase in error rate, and the correlation coefficient between linkage and physical maps became lower. Map quality can be improved by repeated genotyping and error correction algorithm. When it is impossible to genotype the whole mapping population repeatedly, 30% would be recommended in repeated genotyping. The EC method had a much lower false positive rate than did the GC method under different error rates. This study systematically expounded the impact of genotyping errors on linkage analysis, providing potential guidelines for improving the accuracy of linkage maps in the presence of genotyping errors.

Low Diversity and High Genetic Structure for Platonia insignis Mart., an Endangered Fruit Tree Species

Platonia insignis is a fruit tree native to Brazil of increasing economic importance, with its pulp trading among the highest market values. This study aimed to evaluate the structure and genomic diversity of P. insignis (bacurizeiro) accessions from six locations in the Brazilian States of Roraima, Amazonas, Pará (Amazon biome), and Maranhão (Cerrado biome). A total of 2031 SNP markers were obtained using genotyping-by-sequencing (GBS), from which 625 outlier SNPs were identified. High genetic structure was observed, with most of the genetic variability (59%) concentrated among locations, mainly between biomes (Amazon and Cerrado). A positive and significant correlation (r = 0.85; p < 0.005) was detected between genetic and geographic distances, indicating isolation by distance. The highest genetic diversity was observed for the location in the Cerrado biome (HE = 0.1746; HO = 0.2078). The locations in the Amazon biome showed low genetic diversity indexes with significant levels of inbreeding. The advance of urban areas, events of burning, and expansion of agricultural activities are most probably the main factors for the genetic diversity reduction of P. insignis. Approaches to functional analysis showed that most of the outlier loci found may be related to genes involved in cellular and metabolic processes.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Population genetics and genomics of Triatoma brasiliensis (Hemiptera, Reduviidae) in an area of high pressure of domiciliary infestation in Northeastern Brazil

Understanding the population dynamics of vectors is crucial for effective control of vector-borne diseases. In the Northeastern Brazilian semi-arid region, Triatoma brasiliensis persists as the most significant Chagas disease vector, frequently displaying recurrent domiciliary infestations. This situation raises relevant public health concerns in the municipality of Currais Novos in the state of Rio Grande do Norte. This area has experienced a high prevalence of peridomiciliary re-infestations by T. brasiliensis, coupled with elevated rates of Trypanosoma cruzi infection. Therefore, we assessed the distribution of genetic variation via mitochondrial Cytochrome b gene (MT-CYB) sequencing (n = 109) and single nucleotide polymorphisms (SNPs, n = 86) to assess the gene flow among distinct populations distributed in varied geographic spots and environments, mainly sylvatic and peridomiciliary. Insects were collected from rural communities at Currais Novos, enclosed within a 16 km radius. Sampling included 13 populations: one intradomiciliary, eight peridomiciliary, and four sylvatic. Furthermore, an external population located 220 km from Currais Novos was also included in the study. The method employed to obtain SNP information relied on ddRAD-seq genotyping-by-sequencing (GBS), enabling a genome-wide analysis to infer genetic variation. Through AMOVA analysis of MT-CYB gene variation, we identified four distinct population groups with statistical significance (FCT= 0.42; p<0.05). We identified a total of 3,013 SNPs through GBS, with 11 loci showing putative signs of being under selection. The variation based on 3,002 neutral loci evidenced low genetic structuration based on low FST values (p>0.05), indicating local panmixia. However, resampling algorithms pointed out that three samples from the external population were assigned (>98 %) in a cluster contrasting from the ones putatively under local panmixia - validating the newly applied genome-wide marker for studies on the population genetics at finer-scale resolution for T. brasiliensis. The presence of population structuring in some of the sampled points, as suggested by the mitochondrial marker, leads us to assume that infestations were probably initiated by small populations of females - demographic event poses a risk for rapid re-infestations. The local panmictic pattern revealed by the GBS marker poses a challenge for vector control measures, as re-infestation foci may be distributed over a wide geographical and ecological range. In such instances, vectors exhibit reduced susceptibility to conventional insecticide spraying operations since sylvatic populations are beyond the reach of these interventions. The pattern of infestation exhibited by T. brasiliensis necessitates integrating innovative strategies into the existing control framework, holding the potential to create a more resilient and adaptive vector control program. In our dataset, the results demonstrated that the genetic signals from both markers were complementary. Therefore, it is essential to consider the nature and inheritance pattern of each marker when inferring the pattern of re-infestations.

Resequencing Analyses Revealed Genetic Diversity and Selection Signatures during Rabbit Breeding and Improvement

Domestication has shaped the diverse characteristics of rabbits, including coat color, fur structure, body size, and various physiological traits. Utilizing whole-genome resequencing (DNBSEQ-T7), we analyzed the genetic diversity, population structure, and genomic selection across 180 rabbits from 17 distinct breeds to uncover the genetic basis of these traits. We conducted whole-genome sequencing on 17 rabbit breeds, identifying 17,430,184 high-quality SNPs and analyzing genomic diversity, patterns of genomic variation, population structure, and selection signatures related to coat color, coat structure, long hair, body size, reproductive capacity, and disease resistance. Through PCA and NJ tree analyses, distinct clusters emerged among Chinese indigenous rabbits, suggesting varied origins and domestication histories. Selective sweep testing pinpointed regions and genes linked to domestication and key morphological and economic traits, including those affecting coat color (TYR, ASIP), structure (LIPH), body size (INSIG2, GLI3), fertility (EDNRA, SRD5A2), heat stress adaptation (PLCB1), and immune response (SEC31A, CD86, LAP3). Our study identified key genomic signatures of selection related to traits such as coat color, fur structure, body size, and fertility; these findings highlight the genetic basis underlying phenotypic diversification in rabbits and have implications for breeding programs aiming to improve productive, reproductive, and adaptive traits. The detected genomic signatures of selection also provide insights into rabbit domestication and can aid conservation efforts for indigenous breeds.

Comparison of ddRADseq and EUChip60K SNP genotyping systems for population genetics and genomic selection in Eucalyptus dunnii (Maiden)

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii. Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenance-progeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positively with narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed.

Investigating the genetic basis of salt-tolerance in common bean: a genome-wide association study at the early vegetative stage

Salinity poses a significant challenge to global crop productivity, affecting approximately 20% of cultivated and 33% of irrigated farmland, and this issue is on the rise. Negative impact of salinity on plant development and metabolism leads to physiological and morphological alterations mainly due to high ion concentration in tissues and the reduced water and nutrients uptake. Common bean (Phaseolus vulgaris L.), a staple food crop accounting for a substantial portion of consumed grain legumes worldwide, is highly susceptible to salt stress resulting in noticeable reduction in dry matter gain in roots and shoots even at low salt concentrations. In this study we screened a common bean panel of diversity encompassing 192 homozygous genotypes for salt tolerance at seedling stage. Phenotypic data were leveraged to identify genomic regions involved in salt stress tolerance in the species through GWAS. We detected seven significant associations between shoot dry weight and SNP markers. The candidate genes, in linkage with the regions associated to salt tolerance or harbouring the detected SNP, showed strong homology with genes known to be involved in salt tolerance in Arabidopsis. Our findings provide valuable insights onto the genetic control of salt tolerance in common bean and represent a first contribution to address the challenge of salinity-induced yield losses in this species and poses the ground to eventually breed salt tolerant common bean varieties.

An island 'endemic' born out of hybridization between introduced lineages

Humans have profoundly impacted the distribution of plant and animal species over thousands of years. The most direct example of these effects is human-mediated movement of individuals, either through translocation of individuals within their range or through the introduction of species to new habitats. While human involvement may be suspected in species with obvious range disjunctions, it can be difficult to detect natural versus human-mediated dispersal events for populations at the edge of a species' range, and this uncertainty muddles how we understand the evolutionary history of populations and broad biogeographical patterns. Studies combining genetic data with archaeological, linguistic and historical evidence have confirmed prehistoric examples of human-mediated dispersal; however, it is unclear whether these methods can disentangle recent dispersal events, such as species translocated by European colonizers during the past 500 years. We use genomic DNA from historical museum specimens and historical records to evaluate three hypotheses regarding the timing and origin of Northern Bobwhites (Colinus virginianus) in Cuba, whose status as an endemic or introduced population has long been debated. We discovered that bobwhites from southern Mexico arrived in Cuba between the 12th and 16th centuries, followed by the subsequent introduction of bobwhites from the southeastern USA to Cuba between the 18th and 20th centuries. These dates suggest the introduction of bobwhites to Cuba was human-mediated and concomitant with Spanish colonial shipping routes between Veracruz, Mexico and Havana, Cuba during this period. Our results identify endemic Cuban bobwhites as a genetically distinct population born of hybridization between divergent, introduced lineages.