Henckelia longisepala (H. W. Li) D. J. Middleton & Mich. Möller is a rare and endangered plant species found only in Southeastern Yunnan, China, and Northern Vietnam. It is listed as a threatened species in China and recognized as a plant species with extremely small populations (PSESP), while also having high ornamental value and utilization potential. This study used ddRAD-seq technology to quantify genetic diversity and structure for 32 samples from three extant populations of H. longisepala. The H. longisepala populations were found to have low levels of genetic diversity (Ho = 0.1216, He = 0.1302, Pi = 0.1731, FIS = 0.1456), with greater genetic differentiation observed among populations (FST = 0.3225). As indicated by genetic structure and phylogenetic analyses, samples clustered into three distinct genetic groups that corresponded to geographically separate populations. MaxEnt modeling was used to identify suitable areas for H. longisepala across three time periods and two climate scenarios (SSP1-2.6, SSP5-8.5). High-suitability areas were identified in Southeastern Yunnan Province, Northern Vietnam, and Eastern Laos. Future H. longisepala distribution was predicted to remain centered in these areas, but with a decrease in the total amount of suitable habitat. The present study provides key data on H. longisepala genetic diversity, as well as a theoretical basis for the conservation, development, and utilization of its germplasm resources.

BACKGROUND: Biodiversity in freshwater ecosystems is declining due to an increased anthropogenic footprint. Freshwater crayfish are keystone species in freshwater ecosystems and play a crucial role in shaping the structure and function of their habitats. The Idle Crayfish Austropotamobius bihariensis is a native European species with a narrow distribution range, endemic to the Apuseni Mountains (Romania). Although its area is small, the populations are anthropogenically fragmented. In this context, the assessment of its conservation status is timely.
RESULTS: Using a reduced representation sequencing approach, we identified 4875 genomic SNPs from individuals belonging to 13 populations across the species distribution range. Subsequent population genomic analyses highlighted low heterozygosity levels, low number of private alleles and small effective population size. Our structuring analyses revealed that the genomic similarity of the populations is conserved within the river basins.
CONCLUSIONS: Genomic SNPs represented excellent tools to gain insights into intraspecific genomic diversity and population structure of the Idle Crayfish. Our study highlighted that the analysed populations are at risk due to their limited genetic diversity, which makes them extremely vulnerable to environmental alterations. Thus, our results emphasize the need for conservation measures and can be used as a baseline to establish species management programs.

The east coast of the Indochinese Peninsula is a well-known transition zone from subtropical to tropical systems, yet only a small number of studies have been conducted on the biogeography and phylogeography of aquatic organisms in this region. The Hau Giang medaka, Oryzias haugiangensis, was originally described from the Mekong Delta in southern Vietnam, and later reported also from southeastern Thailand, west of the Mekong Delta region. However, the species\' full geographic range and population genetic structures remain unknown. Field surveys showed a widespread distribution of this species along the east coast of the Indochinese Peninsula, as far as northern Vietnam. A mitochondrial gene phylogeny and population genetic structure analysis using genome-wide single nucleotide polymorphisms revealed that the populations of O. haugiangensis are highly structuralized along the east coast of Vietnam, with the southernmost Mekong Delta population clearly separated from three populations north of central Vietnam. Further field collections are necessary to determine the boundary between the southern and northern populations, and the presence or absence of a hybrid zone.

Despite decades of molecular research, phylogenetic relationships in Palearctic vipers (genus Vipera) still essentially rely on a few loci, such as mitochondrial barcoding genes. Here we examined the diversity and evolution of Vipera with ddRAD-seq data from 33 representative species and subspecies. Phylogenomic analyses of ∼ 1.1 Mb recovered nine major clades corresponding to known species/species complexes which are generally consistent with the mitochondrial phylogeny, albeit with a few deep discrepancies that highlight past hybridization events. The most spectacular case is the Italian-endemic V. walser, which is grouped with the alpine genetic diversity of V. berus in the nuclear tree despite carrying a divergent mitogenome related to the Caucasian V. kaznakovi complex. Clustering analyses of SNPs suggest potential admixture between diverged Iberian taxa (V. aspis zinnikeri and V. seoanei), and confirm that the Anatolian V. pontica corresponds to occasional hybrids between V. (ammodytes) meridionalis and V. kaznakovi. Finally, all analyzed lineages of the V. berus complex (including V. walser and V. barani) form vast areas of admixture and may be delimited as subspecies. Our study sets grounds for future taxonomic and phylogeographic surveys on Palearctic vipers, a group of prime interest for toxinological, ecological, biogeographic and conservation research.

Molluscs have undergone many transitions between separate sexes and hermaphroditism, which is of interest in studying the evolution of sex determination and differentiation. Here, we combined multi-locus genotypes obtained from restriction site-associated DNA (RAD) sequencing with anatomical observations of the gonads of three deep-sea hydrothermal vent gastropods of the genus Alviniconcha living in the southwest Pacific. We found that all three species (Alviniconcha boucheti, Alviniconcha strummeri, and Alviniconcha kojimai) share the same male-heterogametic XY sex-determination system but that the gonads of XX A. kojimai individuals are invaded by a variable proportion of male reproductive tissue. The identification of Y-specific RAD loci (found only in A. boucheti) and the phylogenetic analysis of three sex-linked loci shared by all species suggested that X-Y recombination has evolved differently within each species. This situation of three species showing variation in gonadal development around a common sex-determination system provides new insights into the reproductive mode of poorly known deep-sea species and opens up an opportunity to study the evolution of recombination suppression on sex chromosomes and its association with mixed or transitory sexual systems.

Self-compatibility is a highly desirable trait for pear breeding programs. Our breeding program previously developed a novel self-compatible pollen-part Japanese pear mutant (Pyrus pyrifolia Nakai), \'415-1\', by using γ-irradiated pollen. \'415-1\' carries the S-genotype S4dS5S5, with \"d\" indicating a duplication of S 5 responsible for breakdown of self-incompatibility. Until now, the size and inheritance of the duplicated segment was undetermined, and a reliable detection method was lacking. Here, we examined genome duplications and their inheritance in 140 F1 seedlings resulting from a cross between \'515-20\' (S1S3) and \'415-1\'. Amplicon sequencing of S-RNase and SFBB18 clearly detected S-haplotype duplications in the seedlings. Intriguingly, 30 partially triploid seedlings including genotypes S1S4dS5, S3S4dS5, S1S5dS5, S3S5dS5, and S3S4dS4 were detected among the 140 seedlings. Depth-of-coverage analysis using ddRAD-seq showed that the duplications in those individuals were limited to chromosome 17. Further analysis through resequencing confirmed an 11-Mb chromosome duplication spanning the middle to the end of chromosome 17. The duplicated segment remained consistent in size across generations. The presence of an S3S4dS4 seedling provided evidence for recombination between the duplicated S5 segment and the original S4haplotype, suggesting that the duplicated segment can pair with other parts of chromosome 17. This research provides valuable insights for improving pear breeding programs using partially triploid individuals.

BACKGROUND: Corybas taliensis is an endemic species of sky islands in China. Its habitat is fragile and unstable, and it is likely that the species is threatened. However, it is difficult to determine the conservation priority or unit without knowing the genetic background and the overall distribution of this species. In this study, we used double digest restriction-site associated DNA-sequencing (ddRAD-seq) to investigate the conservation genomics of C. taliensis. At the same time, we modeled the extent of suitable habitat for C. taliensis in present and future (2030 and 2090) habitat using the maximum-entropy (MaxEnt) model.
RESULTS: The results suggested that the related C. fanjingshanensis belongs to C. taliensis and should not be considered a separate species. All the sampling locations were divided into three genetic groups: the Sichuan & Guizhou population (SG population), the Hengduan Mountains population (HD population) and Himalayan population (HM population), and we found that there was complex gene flow between the sampling locations of HD population. MT was distinct genetically from the other sampling locations due to the unique environment in Motuo. The genetic diversity (π, He) of C. taliensis was relatively high, but its contemporary effective population size (Ne) was small. C. taliensis might be currently affected by inbreeding depression, although its large population density may be able to reduce the effect of this. The predicted areas of suitable habitat currently found in higher mountains will not change significantly in the future, and these suitable habitats are predicted to spread to other higher mountains under future climate change. However, suitable habitat in relatively low altitude areas may disappear in the future. This suggests that C. taliensis will be caught in a \'summit trap\' in low altitude areas, however, in contrast, the high altitude of the Himalaya and the Hengduan Mountains are predicted to act as \'biological refuges\' for C. taliensis in the future.
CONCLUSIONS: These results not only provide a new understanding of the genetic background and potential resource distribution of C. taliensis, but also lay the foundation for its conservation and management.

Wolbachia, an endosymbiotic bacterium, is globally used to control arboviruses because of its ability to block arboviral replication and manipulate the reproduction of Wolbachia host, Aedes aegypti. Polymerase chain reaction (PCR)-based Wolbachia detection has been recently reported from natural Ae. aegypti populations. However, due to the technical limitations of PCR, such as primer incompatibility, PCR-based assays are not sufficiently reliable or accurate. In this study, we examined double digestion restriction site-associated DNA sequencing (ddRAD-Seq) efficiency and limitations in Wolbachia detection and quantification in field-collected Ae. aegypti natural populations in Metro Manila, the Philippines, compared with PCR-based assays.
A total of 217 individuals Ae. aegypti were collected from Metropolitan Manila, Philippines. We separated it into 14 populations consisting of 7 female and male populations. We constructed a library for pool ddRAD-Seq per population and also screened for Wolbachia by PCR assays using wsp and 16S rRNA. Wolbachia density per population were measured using RPS17 as the housekeeping gene.
From 146,239,637 sequence reads obtained, 26,299 and 43,778 reads were mapped across the entire Wolbachia genome (with the wAlbA and wAlbB strains, respectively), suggesting that ddRAD-Seq complements PCR assays and supports more reliable Wolbachia detection from a genome-wide perspective. The number of reads mapped to the Wolbachia genome per population positively correlated with the number of Wolbachia-infected individuals per population based on PCR assays and the relative density of Wolbachia in the Ae. aegypti populations based on qPCR, suggesting ddRAD-Seq-based semi-quantification of Wolbachia by ddRAD-Seq. Male Ae. aegypti exhibited more reads mapped to the Wolbachia genome than females, suggesting higher Wolbachia prevalence rates in their case. We detected 150 single nucleotide polymorphism loci across the Wolbachia genome, allowing for more accurate the detection of four strains: wPip, wRi, TRS of Brugia malayi, and wMel.
Taken together, our results demonstrate the feasibility of ddRAD-Seq-based Wolbachia detection from field-collected Ae. aegypti mosquitoes.

The peculiarity of Indian cattle lies in milk quality, resistance to diseases and stressors as well as adaptability. The investigation addressed selection signatures in Gir and Tharparkar cattle, belonging to arid ecotypes of India. Double digest restriction-site associated DNA sequencing (ddRAD-seq) yielded nearly 26 million high-quality reads from unrelated seven Gir and seven Tharparkar cows. In all, 19,127 high-quality SNPs were processed for selection signature analysis. An approach involving within-population composite likelihood ratio (CLR) statistics and between-population FST statistics was used to capture selection signatures within and between the breeds, respectively. A total of 191 selection signatures were addressed using CLR and FST approaches. Selection signatures overlapping 86 and 73 genes were detected as Gir- and Tharparkar-specific, respectively. Notably, genes related to production (CACNA1D, GHRHR), reproduction (ESR1, RBMS3), immunity (NOSTRIN, IL12B) and adaptation (ADAM22, ASL) were annotated to selection signatures. Gene pathway analysis revealed genes in insulin/IGF pathway for milk production, gonadotropin releasing hormone pathway for reproduction, Wnt signalling pathway and chemokine and cytokine signalling pathway for adaptation. This is the first study where selection signatures are identified using ddRAD-seq in indicine cattle breeds. The study shall help in conservation and leveraging genetic improvements in Gir and Tharparkar cattle.

Anopheles stephensi invasion in the Horn of Africa (HoA) poses a substantial risk of increased malaria disease burden in the region. An understanding of the history of introduction(s), establishment(s) and potential A. stephensi sources in the HoA is needed to predict future expansions and establish where they may be effectively controlled. To this end, we take a landscape genomic approach to assess A. stephensi origins and spread throughout the HoA, information essential for vector control. Specifically, we assayed 2070 genome-wide single nucleotide polymorphisms across 214 samples spanning 13 populations of A. stephensi from Ethiopia and Somaliland collected in 2018 and 2020, respectively. Principal component and genetic ancestry analyses revealed clustering that followed an isolation-by-distance pattern, with genetic divergence among the Ethiopian samples significantly correlating with geographical distance. Additionally, genetic relatedness was observed between the northeastern and east central Ethiopian A. stephensi populations and the Somaliland A. stephensi populations. These results reveal population differentiation and genetic connectivity within HoA A. stephensi populations. Furthermore, based on genetic network analysis, we uncovered that Dire Dawa, the site of a spring 2022 malaria outbreak, was one of the major hubs from which sequential founder events occurred in the rest of the eastern Ethiopian region. These findings can be useful for the selection of sites for heightened control to prevent future malaria outbreaks. Finally, we did not detect significant genotype-environmental associations, potentially due to the recency of their colonization and/or other anthropogenic factors leading to the initial spread and establishment of A. stephensi. Our study highlights how coupling genomic data at landscape levels can shed light into even ongoing invasions.