Genome-wide selection signatures address trait specific candidate genes in cattle indigenous to arid regions of India

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.

Conservation genetics and potential geographic distribution modeling of Corybas taliensis, a small 'sky Island' orchid species in China

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.

Genome-wide detection of Wolbachia in natural Aedes aegypti populations using ddRAD-Seq

Background

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.

Methods

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.

Results

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.

Conclusions

Taken together, our results demonstrate the feasibility of ddRAD-Seq-based Wolbachia detection from field-collected Ae. aegypti mosquitoes.

Population genomic analyses reveal population structure and major hubs of invasive Anopheles stephensi in the Horn of Africa

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.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

Construction of a high-density genetic linkage map and QTL mapping for bioenergy-related traits in sweet sorghum [Sorghum bicolor (L.) Moench]

Sorghum is an important but arguably undervalued cereal crop, grown in large areas in Asia and Africa due to its natural resilience to drought and heat. There is growing demand for sweet sorghum as a source of bioethanol as well as food and feed. The improvement of bioenergy-related traits directly affects bioethanol production from sweet sorghum; therefore, understanding the genetic basis of these traits would enable new cultivars to be developed for bioenergy production. In order to reveal the genetic architecture behind bioenergy-related traits, we generated an F₂ population from a cross between sweet sorghum cv. 'Erdurmus' and grain sorghum cv. 'Ogretmenoglu'. This was used to construct a genetic map from SNPs discovered by double-digest restriction-site associated DNA sequencing (ddRAD-seq). F₃ lines derived from each F₂ individual were phenotyped for bioenergy-related traits in two different locations and their genotypes were analyzed with the SNPs to identify QTL regions. On chromosomes 1, 7, and 9, three major plant height (PH) QTLs (qPH1.1, qPH7.1, and qPH9.1) were identified, with phenotypic variation explained (PVE) ranging from 10.8 to 34.8%. One major QTL (qPJ6.1) on chromosome 6 was associated with the plant juice trait (PJ) and explained 35.2% of its phenotypic variation. For fresh biomass weight (FBW), four major QTLs (qFBW1.1, qFBW6.1, qFBW7.1, and qFBW9.1) were determined on chromosomes 1, 6, 7, and 9, which explained 12.3, 14.5, 10.6, and 11.9% of the phenotypic variation, respectively. Moreover, two minor QTLs (qBX3.1 and qBX7.1) of Brix (BX) were mapped on chromosomes 3 and 7, explaining 8.6 and 9.7% of the phenotypic variation, respectively. The QTLs in two clusters (qPH7.1/qBX7.1 and qPH7.1/qFBW7.1) overlapped for PH, FBW and BX. The QTL, qFBW6.1, has not been previously reported. In addition, eight SNPs were converted into cleaved amplified polymorphic sequences (CAPS) markers, which can be easily detected by agarose gel electrophoresis. These QTLs and molecular markers can be used for pyramiding and marker-assisted selection studies in sorghum, to develop advanced lines that include desirable bioenergy-related traits.

A high-quality chromosome-level genome assembly of rohu carp, Labeo rohita, and its utilization in SNP-based exploration of gene flow and sex determination

Labeo rohita (rohu) is a carp important to aquaculture in South Asia, with a production volume close to Atlantic salmon. While genetic improvements to rohu are ongoing, the genomic methods commonly used in other aquaculture improvement programs have historically been precluded in rohu, partially due to the lack of a high-quality reference genome. Here we present a high-quality de novo genome produced using a combination of next-generation sequencing technologies, resulting in a 946 Mb genome consisting of 25 chromosomes and 2,844 unplaced scaffolds. Notably, while approximately half the size of the existing genome sequence, our genome represents 97.9% of the genome size newly estimated here using flow cytometry. Sequencing from 120 individuals was used in conjunction with this genome to predict the population structure, diversity, and divergence in three major rivers (Jamuna, Padma, and Halda), in addition to infer a likely sex determination mechism in rohu. These results demonstrate the utility of the new rohu genome in modernizing some aspects of rohu genetic improvement programs.

Genomic Analyses Implicate the Amazon-Orinoco Plume as the Driver of Cryptic Speciation in a Swimming Crab

The Amazon-Orinoco plume (AOP) is the world's largest freshwater and sediment discharge into the ocean. Previous studies limited to mtDNA suggest that the swimming crab Callinectes ornatus Ordway, 1863 exists as two distinct genetic clusters separated by the AOP. However, questions concerning migration, diversification time, and species delimitation are unresolved. Densely sampling markers across the genome (SNPs) could elucidate the evolutionary processes within this species. Here, we combined mtDNA data and ddRAD-seq to explore the diversification patterns and processes within the swimming crab C. ornatus. We show great genetic differentiation between groups on the north and south sides of the plume but also signs of hybridization. Demographic modeling indicates the divergence between groups starting around 8 Mya following the AOP's formation. After a period of isolation, we detect two incidences of secondary contact with stronger migration in concordance with the North Brazil Current flow. Our results suggest speciation with gene flow explained by the interplay among the AOP, oceanographic currents, and long larval dispersal. This work represents the first investigation employing ddRAD-seq in a marine invertebrate species with distribution encompassing the north and south Atlantic and sheds light on the role of the AOP in the diversification of a marine species.

The mutual history of Schlegel's Japanese gecko (Reptilia: Squamata: Gekkonidae) and humans inscribed in genes and ancient literature

Knowing how the present distribution of organisms was formed is an essential issue in evolutionary ecology. Recently, the distribution of organisms on Earth has been significantly changed by human-mediated dispersal due to globalization. Therefore, significant attention has been paid to such processes. However, although humankind has taken considerable time to achieve modernization, the impact of ancient human activity on ecosystems has not yet been thoroughly studied. We hypothesized that ancient urban development and transitions had a non-negligible effect on species distribution. Inferring the impact of past human activity on ecosystems from ancient literature and verifying that impact by genetic analysis and human history is an effective means of tackling this problem. As geckos, a popular neighbor of human dwellings, are good material for this model, we performed this combination approach using Schlegel's Japanese gecko, Gekko japonicus. We show that G. japonicus migrated from China to the western Japanese archipelago before Christ. The gecko species dispersed itself from western to eastern the archipelago on a time scale of thousands of years. There are many synchronizations between the dispersal history of G. japonicus and the historical development of human society. It is suggested by such synchronizations that humans have influenced the distribution of G. japonicus many times throughout its dispersal history.

QTL mapping for seed morphology using the instance segmentation neural network in Lactuca spp

Wild species of lettuce (Lactuca sp.) are thought to have first been domesticated for oilseed contents to provide seed oil for human consumption. Although seed morphology is an important trait contributing to oilseed in lettuce, the underlying genetic mechanisms remain elusive. Since lettuce seeds are small, a manual phenotypic determination required for a genetic dissection of such traits is challenging. In this study, we built and applied an instance segmentation-based seed morphology quantification pipeline to measure traits in seeds generated from a cross between the domesticated oilseed type cultivar 'Oilseed' and the wild species 'UenoyamaMaruba' in an automated manner. Quantitative trait locus (QTL) mapping following ddRAD-seq revealed 11 QTLs linked to 7 seed traits (area, width, length, length-to-width ratio, eccentricity, perimeter length, and circularity). Remarkably, the three QTLs with the highest LOD scores, qLWR-3.1, qECC-3.1, and qCIR-3.1, for length-to-width ratio, eccentricity, and circularity, respectively, mapped to linkage group 3 (LG3) around 161.5 to 214.6 Mb, a region previously reported to be associated with domestication traits from wild species. These results suggest that the oilseed cultivar harbors genes acquired during domestication to control seed shape in this genomic region. This study also provides genetic evidence that domestication arose, at least in part, by selection for the oilseed type from wild species and demonstrates the effectiveness of image-based phenotyping to accelerate discoveries of the genetic basis for small morphological features such as seed size and shape.

Genetic mapping of QTLs for drought tolerance in chickpea (Cicer arietinum L.)

Chickpea yield is severely affected by drought stress, which is a complex quantitative trait regulated by multiple small-effect genes. Identifying genomic regions associated with drought tolerance component traits may increase our understanding of drought tolerance mechanisms and assist in the development of drought-tolerant varieties. Here, a total of 187 F8 recombinant inbred lines (RILs) developed from an interspecific cross between drought-tolerant genotype GPF 2 (Cicer arietinum) and drought-sensitive accession ILWC 292 (C. reticulatum) were evaluated to identify quantitative trait loci (QTLs) associated with drought tolerance component traits. A total of 21 traits, including 12 morpho-physiological traits and nine root-related traits, were studied under rainfed and irrigated conditions. Composite interval mapping identified 31 QTLs at Ludhiana and 23 QTLs at Faridkot locations for morphological and physiological traits, and seven QTLs were identified for root-related traits. QTL analysis identified eight consensus QTLs for six traits and five QTL clusters containing QTLs for multiple traits on linkage groups CaLG04 and CaLG06. The identified major QTLs and genomic regions associated with drought tolerance component traits can be introgressed into elite cultivars using genomics-assisted breeding to enhance drought tolerance in chickpea.

Resolving species-level diversity of Beringiana and Sinanodonta mussels (Bivalvia: Unionidae) in the Japanese archipelago using genome-wide data

Accurate species identification is of primary importance in ecology and evolutionary biology. For a long time, the unionid mussels Beringiana and Sinanodonta have puzzled researchers trying to unravel their diversity because of their poorly discernible morphologies. A recent study conducted species delineation of unionid mussels based on mitochondrial DNA variation, opening up a new avenue to grasp species diversity of the mussels. However, mtDNA-based classification may not align with species boundaries because mtDNA is prone to introgression and incomplete lineage sorting that cause discordance between species affiliation and gene phylogeny. In this study, we evaluated the validity of the mtDNA-based classification of unionid mussels Beringiana and Sinanodonta in Japan using mitochondrial sequence data, double digest restriction site-associated DNA library (ddRAD) sequencing, and morphological data. We found significant inconsistencies in the mitochondrial and nuclear DNA phylogenies, casting doubt on the reliability of the mtDNA-based classification in this group. In addition, nuclear DNA phylogeny revealed that there are at least two unionid lineages hidden in the mtDNA phylogeny. Although molecular dating technique indicates that Beringiana and Sinanodonta diverged >35 million years ago, their shell morphologies are often indistinguishable. Specifically, morphological analyses exhibited the parallel appearance of nearly identical ball-like shell forms in the two genera in Lake Biwa, which further complicates species identification and the morphological evolution of unionid mussels. Our study adds to a growing body of literature that accurate species identification of unionid mussels is difficult when using morphological characters alone. Although mtDNA-based classification is a simple and convenient way to classify unionid mussels, considerable caution is warranted for its application in ecological and evolutionary studies.

Extensive polyploid clonality was a successful strategy for seagrass to expand into a newly submerged environment

Polyploidy has the potential to allow organisms to outcompete their diploid progenitor(s) and occupy new environments. Shark Bay, Western Australia, is a World Heritage Area dominated by temperate seagrass meadows including Poseidon's ribbon weed, Posidonia australis. This seagrass is at the northern extent of its natural geographic range and experiences extremes in temperature and salinity. Our genomic and cytogenetic assessments of 10 meadows identified geographically restricted, diploid clones (2n = 20) in a single location, and a single widespread, high-heterozygosity, polyploid clone (2n = 40) in all other locations. The polyploid clone spanned at least 180 km, making it the largest known example of a clone in any environment on earth. Whole-genome duplication through polyploidy, combined with clonality, may have provided the mechanism for P. australis to expand into new habitats and adapt to new environments that became increasingly stressful for its diploid progenitor(s). The new polyploid clone probably formed in shallow waters after the inundation of Shark Bay less than 8500 years ago and subsequently expanded via vegetative growth into newly submerged habitats.

Phylogeography of a canopy-forming kelp, Eisenia bicyclis (Laminariales, Phaeophyceae), based on a genome-wide sequencing analysis

Analyses of phylogeographic patterns and genetic diversity provide fundamental information for the management and conservation of species. However, little is published about these patterns in Japanese kelp species. In this study, we conducted phylogeographic analyses of a canopy-forming kelp, Eisenia bicyclis, based on genome-wide SNPs identified by ddRAD-seq. We obtained 1,299 SNPs for 76 samples from nine localities across the distribution. STRUCTURE, NeighborNet, and discriminant analysis of principal components consistently showed high genetic differentiation among the Eastern Pacific, Central Pacific, and Sea of Japan coastal regions. Relatively strong gene flow was detected only within populations in the Eastern Pacific and in the Sea of Japan. Genetic diversity and genetic uniqueness were high in the Central Pacific and low in the Sea of Japan. These results suggest that there were at least three independent refugia corresponding to the three regions during the Last Glacial Maximum (LGM). Furthermore, relatively larger populations in the Central Pacific and smaller populations in the Sea of Japan have been maintained in the demographic history from before the LGM to the present. These phylogeographic histories were supported by an Approximate Bayesian Computation analysis. From a conservation genetics perspective, the loss of southern populations in the Central Pacific would greatly reduce the total genetic diversity of the species. Southern populations in the Sea of Japan, which have relatively low genetic diversity, may be highly vulnerable to environmental change, such as heat waves and increased feeding. Therefore, careful monitoring and conservation are needed in the two regions.

Subtle limits to connectivity revealed by outlier loci within two divergent metapopulations of the deep-sea hydrothermal gastropod Ifremeria nautilei

Hydrothermal vents form archipelagos of ephemeral deep-sea habitats that raise interesting questions about the evolution and dynamics of the associated endemic fauna, constantly subject to extinction-recolonization processes. These metal-rich environments are coveted for the mineral resources they harbor, thus raising recent conservation concerns. The evolutionary fate and demographic resilience of hydrothermal species strongly depend on the degree of connectivity among and within their fragmented metapopulations. In the deep sea, however, assessing connectivity is difficult and usually requires indirect genetic approaches. Improved detection of fine-scale genetic connectivity is now possible based on genome-wide screening for genetic differentiation. Here, we explored population connectivity in the hydrothermal vent snail Ifremeria nautilei across its species range encompassing five distinct back-arc basins in the Southwest Pacific. The global analysis, based on 10 570 single nucleotide polymorphism (SNP) markers derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), depicted two semi-isolated and homogeneous genetic clusters. Demo-genetic modeling suggests that these two groups began to diverge about 70 000 generations ago, but continue to exhibit weak and slightly asymmetrical gene flow. Furthermore, a careful analysis of outlier loci showed subtle limitations to connectivity between neighboring basins within both groups. This finding indicates that migration is not strong enough to totally counterbalance drift or local selection, hence questioning the potential for demographic resilience at this latter geographical scale. These results illustrate the potential of large genomic datasets to understand fine-scale connectivity patterns in hydrothermal vents and the deep sea.

Genomic Data Reveals Population Genetic and Demographic History of Magnolia fistulosa (Magnoliaceae), a Plant Species With Extremely Small Populations in Yunnan Province, China

Elucidating the genetic background of threatened species is fundamental to their management and conservation, and investigating the demographic history of these species is helpful in the determination of the threats facing them. The woody species of the genus Magnolia (Magnoliaceae) have high economic, scientific and ecological values. Although nearly half of all Magnolia species have been evaluated as threatened, to date there has been no population genetic study employing Next Generation Sequencing (NGS) technology in this genus. In the present study, we investigate the conservation genomics of Magnolia fistulosa, a threatened species endemic to the limestone area along the Sino-Vietnamese border, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. To increase the reliability of our statistical inferences, we employed two approaches, Stacks and ipyrad, for SNP calling. A total of 15,272 and 18,960, respectively, putatively neutral SNPs were generated by Stacks and ipyrad. Relatively high genetic diversity and large population divergence were detected in M. fistulosa. Although higher absolute values were calculated using the ipyrad data set, the two data sets showed the same trends in genetic diversity (π, H e), population differentiation (F ST) and inbreeding coefficients (F IS). A change in the effective population size of M. fistulosa within the last 1 Ma was detected, including a population decline about 0.5-0.8 Ma ago, a bottleneck event about 0.2-0.3 Ma ago, population fluctuations during the last glacial stage, and the recovery of effective population size after the last glacial maximum. Our findings not only lay the foundation for the future conservation of this species, but also provide new insights into the evolutionary history of the genus Magnolia in southeastern Yunnan, China.

Geographic patterns of genomic diversity and structure in the C4 grass Panicum hallii across its natural distribution

Geographic patterns of within-species genomic diversity are shaped by evolutionary processes, life history and historical and contemporary factors. New genomic approaches can be used to infer the influence of such factors on the current distribution of infraspecific lineages. In this study, we evaluated the genomic and morphological diversity as well as the genetic structure of the C₄ grass Panicum hallii across its complex natural distribution in North America. We sampled extensively across the natural range of P. hallii in Mexico and the USA to generate double-digestion restriction-associated DNA (ddRAD) sequence data for 423 individuals from 118 localities. We used these individuals to study the divergence between the two varieties of P. hallii, P. hallii var. filipes and P. hallii var. hallii as well as the genetic diversity and structure within these groups. We also examined the possibility of admixture in the geographically sympatric zone shared by both varieties, and assessed distribution shifts related with past climatic fluctuations. There is strong genetic and morphological divergence between the varieties and consistent genetic structure defining seven genetic clusters that follow major ecoregions across the range. South Texas constitutes a hotspot of genetic diversity with the co-occurrence of all genetic clusters and admixture between the two varieties. It is likely a recolonization and convergence point of populations that previously diverged in isolation during fragmentation events following glaciation periods.

The natural hybridization between species Ligularia nelumbifolia and Cremanthodium stenoglossum (Senecioneae, Asteraceae) suggests underdeveloped reproductive isolation and ambiguous intergeneric boundary

Natural hybridization is frequent in plants and is considered an important factor facilitating speciation. The natural intergeneric hybridization between Ligularia and Cremanthodium was previously confirmed using a couple of DNA markers. However, the mechanism of this intergeneric hybridization and the role of reproductive isolation in the process of hybridization remain unclear. Here we used double-digest restriction site-associated DNA sequencing (ddRAD-seq) to further quantify the occurrence of hybridization, the genetic structure of the hybrid population and the role of reproductive isolation between Ligularia nelumbifolia and Cremanthodium stenoglossum. The results based on the ddRAD-seq SNP data sets indicated that hybridization between L. nelumbifolia and C. stenoglossum was restricted to F₁s, and no gene introgression was identified between these two species. STRUCTURE analysis and maximum likelihood (ML) tree results showed a slightly larger genetic contribution of L. nelumbifolia to putative hybrid F₁s. We deduced that the reproductive isolation between these two parent species is not well-developed but still strong enough to maintain the genetic integrity of the species, and that their F₁s are sterile or with low fertility. Given the poorly resolved phylogenetic relationship between Ligularia and Cremanthodium, the occurrence of natural hybridization between L. nelumbifolia and C. stenoglossum may provide new insights into the re-circumscription and re-delimitation of these two genera.

A high-density SNP-based genetic map and several economic traits-related loci in Pelteobagrus vachelli

BACKGROUND

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection in aquaculture species. Pelteobagrus vachelli is a very popular commercial species in Asia. However, some specific characters hindered achievement of the traditional selective breeding based on phenotypes, such as lack of large-scale genomic resource and short of markers tightly associated with growth, sex determination and hypoxia tolerance related traits.

RESULTS

By making use of 5059 ddRAD markers in P. vachelli, a high-resolution genetic linkage map was successfully constructed. The map' length was 4047.01 cM by using an interval of 0.11 cm, which is an average marker standard. Comparative genome mapping revealed that a high proportion (83.2%) of markers with a one-to-one correspondence were observed between P. vachelli and P. fulvidraco. Based on the genetic map, 8 significant genome-wide QTLs for 4 weight, 1 body proportion, 2 sex determination, and 1 hypoxia tolerance related traits were detected on 4 LGs. Some SNPs from these significant genome-wide QTLs were observably associated with these phenotypic traits in other individuals by Kompetitive Allele Specific PCR. In addition, two candidate genes for weight, Sipa1 and HSD11B2, were differentially expressed between fast-, medium- and slow-growing P. vachelli. Sema7a, associated with hypoxia tolerance, was induced after hypoxia exposure and reoxygenation.

CONCLUSIONS

We mapped a set of suggestive and significant QTLs as well as candidate genes for 12 growth, 1 sex determination and 1 hypoxia tolerance related traits based on a high-density genetic linkage map by making use of SNP markers for P. fulvidraco. Our results have offered a valuable method about the much more efficient production of all-male, fast growth and hypoxia tolerance P. vachelli for the aquaculture industry.

Conservation Genomics of a Threatened Rhododendron: Contrasting Patterns of Population Structure Revealed From Neutral and Selected SNPs

Though it is well-acknowledged that next generation sequencing (NGS) technologies can provide further insights into plant conservation management than traditional molecular markers, studies employing NGS to address conservation genomics and subsequent conservation strategies for threatened plants are still rare. Rhododendron is the largest genus of woody plants in China, and many species are threatened, however, to date there has been no conservation genetic research using NGS in this genus. In the present study, we investigated the conservation genetics of R. cyanocarpum, a threatened species endemic to the Cangshan Mountains in Yunnan, China, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. Due to the availability of sufficient SNPs, we were able to distinguish between neutral and putatively selected SNPs and were able to further investigate the genetic diversity, population structure, and differentiation in R. cyanocarpum, as well as make an estimation of its demographic history. A total of 6,584 SNPs were obtained, of which 5,729 were neutral (detected using Tajima’s D). In terms of the 5,729 neutral SNPs, R. cyanocarpum had a higher genetic diversity (π = 0.0702 ± 0.0017, H_e = 0.0675 ± 0.0016) than other plant species assessed using Rad-seq methods, while population differentiation (F_st from 0.0314 to 0.0452) was weak. Interestingly, contrasting patterns of population structure were revealed from all neutral and selected SNPs, with distinct genetic clusters forming for all SNPs and neutral SNPs, but no distinct subgroups for selected ones. Moreover, we were able to detect changes in effective population size (N_e) of R. cyanocarpum from 150,000 years ago, including a bottleneck event ca. 60,000 years ago, followed by recovery of N_e over a short period, and a subsequent gradual decline in N_e to date. Implications for conserving R. cyanocarpum based on these main results are then discussed.

Development of High-Density Genetic Linkage Maps and Identification of Loci for Chestnut Gall Wasp Resistance in Castanea spp

Castanea sativa is an important multipurpose species in Europe for nut and timber production as well as for its role in the landscape and in the forest ecosystem. This species has low tolerance to chestnut gall wasp (Dryocosmus kuriphilus Yasumatsu), which is a pest that was accidentally introduced into Europe in early 2000 and devastated forest and orchard trees. Resistance to the gall wasp was found in the hybrid cultivar ‘Bouche de Bétizac’ (C. sativa × C. crenata) and studied by developing genetic linkage maps using a population derived from a cross between ‘Bouche de Bétizac’ and the susceptible cultivar ‘Madonna’ (C. sativa). The high-density genetic maps were constructed using double-digest restriction site-associated DNA-seq and simple sequence repeat markers. The map of ‘Bouche de Bétizac’ consisted of 1459 loci and spanned 809.6 cM; the map of ‘Madonna’ consisted of 1089 loci and spanned 753.3 cM. In both maps, 12 linkage groups were identified. A single major QTL was recognized on the ‘Bouche de Bétizac’ map, explaining up to 67–69% of the phenotypic variance of the resistance trait (Rdk1). The Rdk1 quantitative trait loci (QTL) region included 11 scaffolds and two candidate genes putatively involved in the resistance response were identified. This study will contribute to C. sativa breeding programs and to the study of Rdk1 genes.

ddRAD-seq derived genome-wide SNPs, high density linkage map and QTLs for fruit quality traits in strawberry (Fragaria x ananassa)

Understanding the genetic determinants are essential for improving the fruit quality traits of strawberry. In this study, we focused on mapping the loci for fruit-length (FL), -diameter (FD), -weight (FW) and -soluble solid content (SSC) using the genome-wide single nucleotide polymorphisms (SNPs) identified via ddRAD-sequencing of the F₁ population raised from Maehyang (♀) X Festival (♂). A total of 12,698 high quality SNPs were identified of which 1554 SNPs that showed significant Mendelian segregation (p < 0.05) were mapped to 53 linkage groups (LG) spanning a total of 2937.93 cM with an average marker density of 2.14 cM/locus. Six QTLs for FL and four QTLs for each of FD, FW and SSC were identified that explained 24–35%, 21–42%, 24–54% and 23–50% of overall phenotypic variations, respectively. The genes that lie within these QTL regions were extracted and discussed thoroughly. In addition, a high resolution melting marker (MF154) were designed based on the SNP A1⁷²³G of the UDP-glucose 4-epimerase GEPI48-like gene FAN_iscf00021287. The marker detected the high vs low sugar containing F₁ plants and commercial cultivars with 81.39% and 86.95% detection accuracy, respectively. These SNPs, linkage map, QTLs and candidate genes will be helpful in understanding and improving the fruit quality traits of strawberry.

Resolving complicated relationships of the Panax bipinnatifidus complex in southwestern China by RAD-seq data

The P. binpinnatifidus complex included most of the Panax species distributed in Sino-Himalaya regions except for P. pseudoginseng, P. stipuleanatus and P. notoginseng. However, the delimitation and identification of these taxa within the species complex are very difficult due to the existence of morphological intermediates, and their evolutionary relationships remain unresolved despite several studies have been carried out based on traditional DNA markers. The taxonomic uncertainty hinders the identification, conservation and exploration of these wild populations of Panax. To study this species complex, we employed ddRAD-seq data of these taxa from 18 different localities of southwestern China, using two RAD analysis pipelines, STACKS and pyRAD. Based on the results of phylogenetic analysis, the species complex was divided into four clades with high supports, which largely agreed with morphologically described species. Two clades, corresponding to P. vietnamensis and P. zingiberensis, respectively, were sister groups, indicating that these two species had a closer genetic relationship; the third clade was consisted of samples with bamboo-like rhizomes named as P. wangianus clade, and the fourth one with moniliform rhizomes was named as P. bipinnatifidus clade. The population genetic structure analysis and D-statistics test showed the localized admixture among these species, which indicated that introgression had occurred among the related lineages continuously distributed in southeastern Yunnan and adjacent regions.

Ryanodine receptor mutations (G4946E and I4790K) differentially responsible for diamide insecticide resistance in diamondback moth, Plutella xylostella L

This study examined diamondback moth (Plutella xylostella) strains showing high-level resistance to cyantraniliprole (KA17 strain) and to flubendiamide and chlorantraniliprole (KU13 strain). The LC₅₀ value of the KA17 strain against cyantraniliprole was ca. 100-fold higher than that of the KU13 strain. The KA17 strain also exhibited high-level resistance to chlorantraniliprole and flubendiamide equivalent to those in the KU13 strain. The KU13 strain showed a higher LC₅₀ value against cyantraniliprole than the susceptible strains. However, the LC₅₀ value of the KU13 strain against cyantraniliprole was below the agriculturally recommended concentration. Subsequent QTL analysis using ddRAD-seq identified the resistance responsible regions of the KA17 and KU13 strains with different diamide resistance profiles. Ryanodine receptor (RyR) gene was included in the identified regions. Single nucleotide polymorphism calling in the RyR gene using RNA-seq found previously reported G4946E (amino acid mutation from glycine to glutamic acid at amino acid position 4946) and novel I4790K (amino acid mutation from isoleucine to lysine at amino acid position 4790) mutations, respectively, in the RyR of the KU13 and KA17 strains. Functional significance of I4790K in the resistance was confirmed in calcium imaging of the human embryonic kidney 293T cell line expressing Bombyx mori RyR with the mutation. This reporting is the first describing I4790K as a fundamental mechanism responsible for the resistance to the diamides including cyantraniliprole. From this study, we also report up-regulated expression of some degradation enzymes and that of the RyR gene in the KA17 and KU13 strains based on results of RNA-seq data analysis.

Genotyping for Species Identification and Diversity Assessment Using Double-Digest Restriction Site-Associated DNA Sequencing (ddRAD-Seq)

Genotyping-by-sequencing (GBS) is a powerful approach for studying the genetic diversity of legume species. By using restriction enzymes or other methods to generate a reduced representation of the genome for sequencing, GBS can provide genome-wide single nucleotide polymorphisms (SNP) for diversity analysis at high throughput and low cost. Here we describe a novel double-digest restriction site-associated DNA sequencing (ddRAD-seq) approach. We also describe the downstream bioinformatic analysis of the sequencing data, including alignment to a reference genome, de novo assembly, SNP calling, phylogenetic analysis, and structure analysis.

Phylogenomic analyses reveal intractable evolutionary history of a temperate bamboo genus (Poaceae: Bambusoideae)

Shibataea is a genus of temperate bamboos (Poaceae: Bambusoideae) endemic to China, but little is known about its phylogenetic position and interspecific relationships. To elucidate the phylogenetic relationship of the bamboo genus Shibataea, we performed genome-scale phylogenetic analysis of all seven species and one variety of the genus using double digest restriction-site associated DNA sequencing (ddRAD-seq) and whole plastid genomes generated using genome skimming. Our phylogenomic analyses based on ddRAD-seq and plastome data congruently recovered Shibataea as monophyletic. The nuclear data resolved S. hispida as the earliest diverged species, followed by S. chinensis, while the rest of Shibataea can be further divided into two clades. However, the plastid and nuclear topologies conflict significantly. By comparing the results of network analysis and topologies reconstructed from different datasets, we identify S. kumasasa as the most admixed species, which may be caused by incomplete lineage sorting (ILS) or interspecific gene flow with four sympatric species. This study highlights the power of ddRAD and plastome data in resolving complex relationships in the intractable bamboo genus.

Genome-Wide Association Study Reveals Candidate Genes for Flowering Time Variation in Common Bean (Phaseolus vulgaris L.)

The common bean is one of the most important staples in many areas of the world. Extensive phenotypic and genetic characterization of unexplored bean germplasm are still needed to unlock the breeding potential of this crop. Dissecting genetic control of flowering time is of pivotal importance to foster common bean breeding and to develop new varieties able to adapt to changing climatic conditions. Indeed, flowering time strongly affects yield and plant adaptation ability. The aim of this study was to investigate the genetic control of days to flowering using a whole genome association approach on a panel of 192 highly homozygous common bean genotypes purposely developed from landraces using Single Seed Descent. The phenotypic characterization was carried out at two experimental sites throughout two growing seasons, using a randomized partially replicated experimental design. The same plant material was genotyped using double digest Restriction-site Associated DNA sequencing producing, after a strict quality control, a dataset of about 50 k Single Nucleotide Polymorphisms (SNPs). The Genome-Wide Association Study revealed significant and meaningful associations between days to flowering and several SNP markers; seven genes are proposed as the best candidates to explain the detected associations.

QTL Mapping for Red Blotches in Malaysia Red Tilapia (Oreochromis spp.)

Body color is an interesting economic trait in fish. Red tilapia with red blotches may decrease its commercial values. Conventional selection of pure red color lines is a time-consuming and labor-intensive process. To accelerate selection of pure lines through marker-assisted selection, in this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) technology was applied to genotype a full-sib mapping family of Malaysia red tilapia (Oreochromis spp.) (N = 192). Genome-wide significant quantitative trait locus (QTL)-controlling red blotches were mapped onto two chromosomes (chrLG5 and chrLG15) explaining 9.7% and 8.2% of phenotypic variances by a genome-wide association study (GWAS) and linkage-based QTL mapping. Six SNPs from the chromosome chrLG5 (four), chrLG15 (one), and unplaced supercontig GL831288-1 (one) were significantly associated to the red blotch trait in GWAS analysis. We developed nine microsatellite markers and validated significant correlations between genotypes and blotch data (p < 0.05). Our study laid a foundation for exploring a genetic mechanism of body colors and carrying out genetic improvement for color quality in tilapia.

Phylogenomics clarifies biogeographic and evolutionary history, and conservation status of West Indian tremblers and thrashers (Aves: Mimidae)

The West Indian avifauna has provided fundamental insights into island biogeography, taxon cycles, and the evolution of avian behavior. Our interpretations, however, should rely on robust hypotheses of evolutionary relationships and consistent conclusions about taxonomic status in groups with many endemic island populations. Here we present a phylogenetic study of the West Indian thrashers, tremblers, and allies, an assemblage of at least 5 species found on 29 islands, including what is considered the Lesser Antilles' only avian radiation. We improve on previous phylogenetic studies of this group by using double-digest restriction site-associated DNA sequencing (ddRAD-seq) to broadly sample loci scattered across the nuclear genome. A variety of analyses, based on either nucleotide variation in 2223 loci recovered in all samples or at 13,282 loci confidently scored as present or absent in all samples, converged on a single well-supported phylogenetic hypothesis. Results indicate that the resident West Indian taxa form a monophyletic group, exclusive of the Neotropical-Nearctic migratory Gray Catbird Dumetella carolinensis, which breeds in North America; this outcome differs from earlier studies suggesting that Gray Catbird was nested within a clade of island resident species. Thus, our findings imply a single colonization of the West Indies without the need to invoke a subsequent 'reverse colonization' of the mainland by West Indian taxa. Additionally, our study is the first to sample both endemic subspecies of the endangered White-breasted Thrasher Ramphocinclus brachyurus. We find that these subspecies have a long history of evolutionary independence with no evidence of gene flow, and are as genetically divergent from each other as other genera in the group. These findings support recognition of R. brachyurus (restricted to Martinique) and the Saint Lucia Thrasher R. sanctaeluciae as two distinct, single-island endemic species, and indicate the need to re-evaluate conservation plans for these taxa. Our results demonstrate the utility of phylogenomic datasets for generating robust systematic hypotheses.

Mapping of a novel clubroot resistance QTL using ddRAD-seq in Chinese cabbage (Brassica rapa L.)

BACKGROUND

Plasmodiophora brassicae is a soil-borne plant pathogen that causes clubroot disease, which results in crop yield loss in cultivated Brassica species. Here, we investigated whether a quantitative trait locus (QTL) in B. rapa might confer resistance to a Korean P. brassicae pathotype isolate, Seosan. We crossed resistant and susceptible parental lines and analyzed the segregation pattern in a F₂ population of 348 lines. We identified and mapped a novel clubroot resistance QTL using the same mapping population that included susceptible Chinese cabbage and resistant turnip lines. Forty-five resistant and 45 susceptible F₂ lines along with their parental lines were used for double digest restriction site-associated DNA sequencing (ddRAD-seq). High resolution melting (HRM)-based validation of SNP positions was conducted to confirm the novel locus.

RESULTS

A 3:1 ratio was observed for resistant: susceptible genotypes, which is in accordance with Mendelian segregation. ddRAD-seq identified a new locus, CRs, on chromosome A08 that was different from the clubroot resistance (CR) locus, Crr1. HRM analysis validated SNP positions and constricted CRs region. Four out of seventeen single nucleotide polymorphisms (SNPs) positions were within a 0.8-Mb region that included three NBS-LRR candidate genes but not Crr1.

CONCLUSION

The newly identified CRs locus is a novel clubroot resistance locus, as the cultivar Akimeki bears the previously known Crr1 locus but remains susceptible to the Seosan isolate. These results could be exploited to develop molecular markers to detect Seosan-resistant genotypes and develop resistant Chinese cabbage cultivars.

ddRAD-seq based phylogeographic study of Sargassum thunbergii (Phaeophyceae, Heterokonta) around Japanese coast

Using genome-wide SNP data obtained from high-throughput techniques based on double digest restriction site-associated DNA sequencing (ddRAD-seq), we elucidated the migration history and genetic diversity of the Japanese population of the ecologically important brown seaweed Sargassum thunbergii (Mertens ex Roth) Kuntze. STRUCTURE and NeighborNet analyses showed a clear genetic differentiation among populations of four geographic regions: Kyushu (POP1); Sea of Japan (POP2); Hokkaido and Tohoku (POP3); and Pacific coast from Kyushu to Kanto (POP4). Approximate Bayesian Computation (ABC) analysis indicated that POP4 diverged first, followed by the separation between POP2 (the largest effective population size) and POP3; POP1 was the last to form, shaped by the mixture of POP2 (73%) and POP4 (27%). High genetic diversity was detected in POP1 and POP2, whereas low genetic diversity was detected in POP3 and POP4. These results indicated that S. thunbergii populations of Kyushu and the Sea of Japan might have been maintained as large and stable populations gathered different lineages from China, Korea and/or Japan.

Assessment of genetic diversity in Coho salmon (Oncorhynchus kisutch) populations with no family records using ddRAD-seq

OBJECTIVE

Selective breeding for desirable traits is becoming popular in aquaculture. In Miyagi prefecture, Japan, a selectively bred population of Coho salmon (Oncorhynchus kisutch) has been established with the original, randomly breeding population maintained separately. Since they have been bred without family records, the genetic diversity within these populations remains unknown. In this study, we estimated the genetic diversity and key quantitative genetic parameters such as heritability and genomic breeding value for body size traits by means of genomic best linear unbiased prediction to assess the genetic health of these populations.

RESULTS

Ninety-nine and 83 females from the selective and random groups, respectively, were genotyped at 2350 putative SNPs by means of double digest restriction associated DNA sequencing. The genetic diversity in the selectively bred group was low, as were the estimated heritability and prediction accuracy for length and weight (h² = 0.26-0.28; accuracy = 0.34), compared to the randomly bred group (h² = 0.50-0.60; accuracy = 0.51-0.54). Although the tested sample size was small, these results suggest that further selection is difficult for the selectively bred population, while there is some potential for the randomly bred group, especially with the aid of genomic information.

Development of Diagnostic SNP Markers To Monitor Hybridization Between Sika Deer (Cervus nippon) and Wapiti (Cervus elaphus)

Sika deer (Cervus Nippon) and wapiti (Cervus elaphus) are closely related species and their hybridization can result in significant allele-shift of their gene pool. Additive genetic effects and putative heterotic effects of their hybridization on growth performance could confer considerable economic advantage in deer farming. Here, we used double-digest restriction site-associated DNA sequencing technology (ddRAD-seq) and detected ∼320,000 genome-wide SNPs from 30 captive individuals: 7 sika deer, 6 wapiti and 17 F1 hybrids (reciprocal cross). By screening observed heterozygosity of each SNP across four taxonomic groups, we report for the first time a resource of 2,015 putative diagnostic SNP markers (species-specific SNPs for sika deer and wapiti), which can be used to design tools for assessing or monitoring the degree of hybridization between sika deer and wapiti. These ddRAD-seq data and SNP datasets are also valuable resources for genome-wide studies, including trait discovery for breeders of domestic deer.

Phylogenomic history of enigmatic pygmy perches: implications for biogeography, taxonomy and conservation

Pygmy perches (Percichthyidae) are a group of poorly dispersing freshwater fishes that have a puzzling biogeographic disjunction across southern Australia. Current understanding of pygmy perch phylogenetic relationships suggests past east-west migrations across a vast expanse of now arid habitat in central southern Australia, a region lacking contemporary rivers. Pygmy perches also represent a threatened group with confusing taxonomy and potentially cryptic species diversity. Here, we present the first study of the evolutionary history of pygmy perches based on genome-wide information. Data from 13 991 ddRAD loci and a concatenated sequence of 1 075 734 bp were generated for all currently described and potentially cryptic species. Phylogenetic relationships, biogeographic history and cryptic diversification were inferred using a framework that combines phylogenomics, species delimitation and estimation of divergence times. The genome-wide phylogeny clarified the biogeographic history of pygmy perches, demonstrating multiple east-west events of divergence within the group across the Australian continent. These results also resolved discordance between nuclear and mitochondrial data from a previous study. In addition, we propose three cryptic species within a southwestern species complex. The finding of potentially new species demonstrates that pygmy perches may be even more susceptible to ecological and demographic threats than previously thought. Our results have substantial implications for improving conservation legislation of pygmy perch lineages, especially in southwestern Western Australia.

A comparison of different methods for preserving plant molecular materials and the effect of degraded DNA on ddRAD sequencing

Obtaining high-quality plant materials for experiments is challenging for many research projects. Therefore, it is of special importance to determine the best method for preserving biological macromolecules like DNA, which degrade over time. Although some research has demonstrated that DNA degradation has little effect on traditional molecular markers, the effects of DNA degradation on ddRAD-seq, a popular reduced-representation sequencing technology, have not been adequately investigated. In this study, we first chose six woody bamboo species (Bambusoideae, Poaceae) to explore appropriate methods for preserving molecular materials with two DNA extraction approaches. Then we sequenced twenty-one bamboos and examined the effects of DNA quality on data generation using the ddRAD-seq technique (MiddRAD-seq). Finally, we reconstructed phylogenies of twenty woody bamboo species. We found that the integrity of dry-powdered DNA was preserved longer than that of TE-dissolved DNA, regardless of whether the DNA was extracted by a modified CTAB protocol or DNAsecure plant kit. The ddRAD-seq data were robust, except when DNA was severely degraded. In addition, we resolved the phylogenetic positions of the sampled Phyllostachys spp. Our results suggest that dry-powdered DNA is the most appropriate preservation method for plant molecular materials. Furthermore, a moderate level of DNA degradation has little effect on reduced representation sequencing techniques represented by ddRAD-seq.

Genetic architecture of traits associated with reproductive barriers in Silene: Coupling, sex chromosomes and variation

The evolution of reproductive barriers and their underlying genetic architecture is of central importance for the formation of new species. Reproductive barriers can be controlled either by few large-effect loci suggesting strong selection on key traits, or by many small-effect loci, consistent with gradual divergence or with selection on polygenic or multiple traits. Genetic coupling between reproductive barrier loci further promotes divergence, particularly divergence with ongoing gene flow. In this study, we investigated the genetic architectures of ten morphological, phenological and life history traits associated with reproductive barriers between the hybridizing sister species Silene dioica and S. latifolia; both are dioecious with XY-sex determination. We used quantitative trait locus (QTL) mapping in two reciprocal F₂ crosses. One to six QTLs per trait, including nine major QTLs (PVE > 20%), were detected on 11 of the 12 linkage groups. We found strong evidence for coupling of QTLs for uncorrelated traits and for an important role of sex chromosomes in the genetic architectures of reproductive barrier traits. Unexpectedly, QTLs detected in the two F₂ crosses differed largely, despite limited phenotypic differences between them and sufficient statistical power. The widely dispersed genetic architectures of traits associated with reproductive barriers suggest gradual divergence or multifarious selection. Coupling of the underlying QTLs likely promoted divergence with gene flow in this system. The low congruence of QTLs between the two crosses further points to variable and possibly redundant genetic architectures of traits associated with reproductive barriers, with important implications for the evolutionary dynamics of divergence and speciation.

Genomic data reveal cryptic lineage diversification and introgression in Californian golden cup oaks (section Protobalanus)

Here we study hybridization, introgression and lineage diversification in the widely distributed canyon live oak (Quercus chrysolepis) and the relict island oak (Q. tomentella), two Californian golden cup oaks with an intriguing biogeographical history. We employed restriction-site-associated DNA sequencing and integrated phylogenomic and population genomic analyses to study hybridization and reconstruct the evolutionary past of these taxa. Our analyses revealed the presence of two cryptic lineages within Q. chrysolepis. One of these lineages shares its most recent common ancestor with Q. tomentella, supporting the paraphyly of Q. chrysolepis. The split of these lineages was estimated to take place during the late Pliocene or the early Pleistocene, a time corresponding well with the common presence of Q. tomentella in the fossil records of continental California. Analyses also revealed historical hybridization among lineages, high introgression from Q. tomentella into Q. chrysolepis in their current area of sympatry, and widespread admixture between the two lineages of Q. chrysolepis in contact zones. Our results support that the two lineages of Q. chrysolepis behave as a single functional species phenotypically and ecologically well differentiated from Q. tomentella, a situation that can be only accommodated considering hybridization and speciation as a continuum with diffuse limits.

Molecular evidence for the first records of facultative parthenogenesis in elapid snakes

Parthenogenesis is a form of asexual reproduction by which embryos develop from unfertilized eggs. Parthenogenesis occurs in reptiles; however, it is not yet known to occur in the widespread elapid snakes (Elapidae), which include well-known taxa such as cobras, mambas, taipans and sea snakes. Here, we describe the production of viable parthenogens in two species of Australo-Papuan elapids with divergent reproductive modes: the oviparous coastal/Papuan taipan (Oxyuranus scutellatus) and the viviparous southern death adder (Acanthophis antarcticus). Analyses of nuclear SNP data excluded paternity for putative fathers and convincingly demonstrated asexual reproduction, thus representing the first evidence of facultative parthenogenesis in Elapidae. Our finding has broad implications for understanding the evolution of reproductive diversity in snakes, as well as managing the conservation of genetic diversity in wild and captive populations.

Ecological disturbance influences adaptive divergence despite high gene flow in golden perch (Macquaria ambigua): Implications for management and resilience to climate change

Populations that are adaptively divergent but maintain high gene flow may have greater resilience to environmental change as gene flow allows the spread of alleles that have already been tested elsewhere. In addition, populations naturally subjected to ecological disturbance may already hold resilience to future environmental change. Confirming this necessitates ecological genomic studies of high dispersal, generalist species. Here we perform one such study on golden perch (Macquaria ambigua) in the Murray-Darling Basin (MDB), Australia, using a genome-wide SNP data set. The MDB spans across arid to wet and temperate to subtropical environments, with low to high ecological disturbance in the form of low to high hydrological variability. We found high gene flow across the basin and three populations with low neutral differentiation. Genotype-environment association analyses detected adaptive divergence predominantly linked to an arid region with highly variable riverine flow, and candidate loci included functions related to fat storage, stress and molecular or tissue repair. The high connectivity of golden perch in the MDB will likely allow locally adaptive traits in its most arid and hydrologically variable environment to spread and be selected in localities that are predicted to become arid and hydrologically variable in future climates. High connectivity in golden perch is likely due to their generalist life history and efforts of fisheries management. Our study adds to growing evidence of adaptation in the face of gene flow and highlights the importance of considering ecological disturbance and adaptive divergence in biodiversity management.

Genome-Wide SNP Discovery and Analysis of Genetic Diversity in Farmed Sika Deer (Cervus nippon) in Northeast China Using Double-Digest Restriction Site-Associated DNA Sequencing

Sika deer are an economically valuable species owing to their use in traditional Chinese medicine, particularly their velvet antlers. Sika deer in northeast China are mostly farmed in enclosure. Therefore, genetic management of farmed sika deer would benefit from detailed knowledge of their genetic diversity. In this study, we generated over 1.45 billion high-quality paired-end reads (288 Gbp) across 42 unrelated individuals using double-digest restriction site-associated DNA sequencing (ddRAD-seq). A total of 96,188 (29.63%) putative biallelic SNP loci were identified with an average sequencing depth of 23×. Based on the analysis, we found that the majority of the loci had a deficit of heterozygotes (F_IS >0) and low values of H_obs, which could be due to inbreeding and Wahlund effects. We also developed a collection of high-quality SNP probes that will likely be useful in a variety of applications in genotyping for cervid species in the future.

Rapid genetic and morphologic divergence between captive and wild populations of the endangered Leon Springs pupfish, Cyprinodon bovinus

The Leon Springs pupfish (Cyprinodon bovinus) is an endangered species currently restricted to a single desert spring and a separate captive habitat in southwestern North America. Following establishment of the captive population from wild stock in 1976, the wild population has undergone natural population size fluctuations, intentional culling to purge genetic contamination from an invasive congener (Cyprinodon variegatus) and augmentation/replacement of wild fish from the captive stock. A severe population decline following the most recent introduction of captive fish prompted us to examine whether the captive and wild populations have differentiated during the short time they have been isolated from one another. If so, the development of divergent genetic and/or morphologic traits between populations could contribute to a diminished ability of fish from one location to thrive in the other. Examination of genomewide single nucleotide polymorphisms and morphologic variation revealed no evidence of residual C. variegatus characteristics in contemporary C. bovinus samples. However, significant genetic and morphologic differentiation was detected between the wild and captive populations, some of which might reflect local adaptation. Our results indicate that genetic and physical characteristics can diverge rapidly between isolated subdivisions of managed populations, potentially compromising the value of captive stock for future supplementation efforts. In the case of C. bovinus, our findings underscore the need to periodically inoculate the captive population with wild genetic material to help mitigate genetic, and potentially morphologic, divergence between them and also highlight the utility of parallel morphologic and genomic evaluation to inform conservation management planning.

Genetic diversity of Atlantic Bluefin tuna in the Mediterranean Sea: insights from genome-wide SNPs and microsatellites

Background

Elucidating the patterns of the Atlantic Bluefin tuna [ABFT, Thunnus thynnus (Linnaeus, 1758)] population structure constitutes a challenging task of great importance. Most of the unique challenges stem from its biology, as well as the attributes of the marine realm in which it disperses. Accurate information is urgently needed for stock assessment, and the identification of critical features to the persistence and adaptation of populations in order to formulate and adopt effective strategies for ABFT conservation and management. Conclusions of a great number of ABFT genetic studies on the Mediterranean Sea stock structure are rather controversial and not yet conclusive. In this study, ABFT genomic diversity was investigated in the Mediterranean Sea, which is the most important area for the species’ reproduction.

Results

Analyzing genome-wide SNPs and microsatellites from ABFT samples collected throughout the Mediterranean Sea did not provide strong evidence of genetic structure, pointing towards the existence of a single panmictic unit. An alternative view would recognize a failure to reject the null hypothesis of a panmictic unit as an effect of the study’s sampling design, the type of markers used, and the effectiveness/suitability of analysis methods in respect to the species biological characteristics or any combination of the above.

Conclusions

Unravelling the drivers of ABFT population diversity would require the consideration of important aspects of the species spawning behavior for the determination of the appropriate sampling design. Novel approaches and methods of analysis that will bring together experts in genetics/-omics, ecology and oceanography are deemed necessary. Analyzing ABFT genetic data under the discipline of seascape genetics could provide the analysis framework under which major abiotic and biotic forces controlling ABFT recruitment could be identified, elucidating the complicated population dynamics of the species, while multiple and continuous fisheries monitoring should in all cases be considered as a prerequisite in order to achieve efficient and long-term ABFT conservation.

Application of a dense genetic map for assessment of genomic responses to selection and inbreeding in Heliothis virescens

Adaptation of pest species to laboratory conditions and selection for resistance to toxins in the laboratory are expected to cause inbreeding and genetic bottlenecks that reduce genetic variation. Heliothis virescens, a major cotton pest, has been colonized in the laboratory many times, and a few laboratory colonies have been selected for Bacillus thuringiensis (Bt) resistance. We developed 350-bp double-digest restriction-site associated DNA-sequencing (ddRAD-seq) molecular markers to examine and compare changes in genetic variation associated with laboratory adaptation, artificial selection and inbreeding in this nonmodel insect species. We found that allelic and nucleotide diversity declined dramatically in laboratory-reared H. virescens as compared with field-collected populations. The declines were primarily a result of the loss of low frequency alleles present in field-collected H. virescens. A further, albeit modest decline in genetic diversity was observed in a Bt-selected population. The greatest decline was seen in H. virescens that were sib-mated for 10 generations, in which more than 80% of loci were fixed for a single allele. To determine which regions of the genome were resistant to fixation in our sib-mated line, we generated a dense intraspecific linkage map containing three PCR-based and 659 ddRAD-seq markers. Markers that retained polymorphism were observed in small clusters spread over multiple linkage groups, but this clustering was not statistically significant. Overall, we have confirmed and extended the general expectations for reduced genetic diversity in laboratory colonies, provided tools for further genomic analyses and produced highly homozygous genomic DNA for future whole genome sequencing of H. virescens.

ddRAD-seq phylogenetics based on nucleotide, indel, and presence-absence polymorphisms: Analyses of two avian genera with contrasting histories

Genotype-by-sequencing (GBS) methods have revolutionized the field of molecular ecology, but their application in molecular phylogenetics remains somewhat limited. In addition, most phylogenetic studies based on large GBS data sets have relied on analyses of concatenated data rather than species tree methods that explicitly account for genealogical stochasticity among loci. We explored the utility of "double-digest" restriction site-associated DNA sequencing (ddRAD-seq) for phylogenetic analyses of the Lagonosticta firefinches (family Estrildidae) and the Vidua brood parasitic finches (family Viduidae). As expected, the number of homologous loci shared among samples was negatively correlated with genetic distance due to the accumulation of restriction site polymorphisms. Nonetheless, for each genus, we obtained data sets of ∼3000 loci shared in common among all samples, including a more distantly related outgroup taxon. For all samples combined, we obtained >1000 homologous loci despite ∼20my divergence between estrildid and parasitic finches. In addition to nucleotide polymorphisms, the ddRAD-seq data yielded large sets of indel and locus presence-absence polymorphisms, all of which had higher consistency indices than mtDNA sequence data in the context of concatenated parsimony analyses. Species tree methods, using individual gene trees or single nucleotide polymorphisms as input, generated results broadly consistent with analyses of concatenated data, particularly for Lagonosticta, which appears to have a well resolved, bifurcating history. Results for Vidua were also generally consistent across methods and data sets, although nodal support and results from different species tree methods were more variable. Lower gene tree congruence in Vidua is likely the result of its unique evolutionary history, which includes rapid speciation by host shift and occasional hybridization and introgression due to incomplete reproductive isolation. We conclude that ddRAD-seq is a cost-effective method for generating robust phylogenetic data sets, particularly for analyses of closely related species and genera.

Genetic mapping of horizontal stripes in Lake Victoria cichlid fishes: benefits and pitfalls of using RAD markers for dense linkage mapping

The genetic dissection of naturally occurring phenotypes sheds light on many fundamental and longstanding questions in speciation and adaptation and is a central research topic in evolutionary biology. Until recently, forward-genetic approaches were virtually impossible to apply to nonmodel organisms, but the development of next-generation sequencing techniques eases this difficulty. Here, we use the ddRAD-seq method to map a colour trait with a known adaptive function in cichlid fishes, well-known textbook examples for rapid rates of speciation and astonishing phenotypic diversification. A suite of phenotypic key innovations is related to speciation and adaptation in cichlids, among which body coloration features prominently. The focal trait of this study, horizontal stripes, evolved in parallel in several cichlid radiations and is associated with piscivorous foraging behaviour. We conducted interspecific crosses between Haplochromis sauvagei and H. nyererei and constructed a linkage map with 867 SNP markers distributed on 22 linkage groups and total size of 1130.63 cM. Lateral stripes are inherited as a Mendelian trait and map to a single genomic interval that harbours a paralog of a gene with known function in stripe patterning. Dorsolateral and mid-lateral stripes were always coinherited and are thus under the same genetic control. Additionally, we directly quantify the genotyping error rates in RAD markers and offer guidelines for identifying and dealing with errors. Uncritical marker selection was found to severely impact linkage map construction. Fortunately, by applying appropriate quality control steps, a genotyping accuracy of >99.9% can be reached, thus allowing for efficient linkage mapping of evolutionarily relevant traits.