Serial sequencing of isolength RAD tags for cost-efficient genome-wide profiling of genetic and epigenetic variations

Genome-wide association analysis reveals the genetic basis of thermal tolerance in dwarf surf clam Mulinia lateralis

Recently, elevated seawater temperatures have resulted numerous adverse effects, including significant mortality among bivalves. The dwarf surf clam, Mulinia lateralis, is considered a valuable model species for bivalve research due to its rapid growth and short generation time. The successful cultivation in laboratory setting throughout its entire life cycle makes it an ideal candidate for exploring the potential mechanisms underlying bivalve responses to thermal stress. In this study, a total of 600 clams were subjected to a 17-day thermal stress experiment at a temperature of 30 °C which is the semi-lethal temperature for this species. Ninety individuals who perished initially were classified as heat-sensitive populations (HSP), while 89 individuals who survived the experiment were classified as heat-tolerant populations (HTP). Subsequently, 179 individuals were then sequenced, and 21,292 single nucleotide polymorphisms (SNPs) were genotyped for downstream analysis. The heritability estimate for survival status was found to be 0.375 ± 0.127 suggesting a genetic basis for thermal tolerance trait. Furthermore, a genome-wide association study (GWAS) identified three SNPs and 10 candidate genes associated with thermal tolerance trait in M. lateralis. These candidate genes were involved in the ETHR/EHF signaling pathway and played pivotal role in signal sensory, cell adhesion, oxidative stress, DNA damage repair, etc. Additionally, qPCR results indicated that, excluding MGAT4A, ZAN, and RFC1 genes, all others exhibited significantly higher expression in the HTP (p < 0.05), underscoring the critical involvement of the ETHR/EHF signaling pathway in M. lateralis' thermal tolerance. These results unveil the presence of standing genetic variations associated with thermal tolerance in M. lateralis, highlighting the regulatory role of the ETHR/EHF signaling pathway in the bivalve's response to thermal stress, which contribute to comprehension of the genetic basis of thermal tolerance in bivalves.

Tracking the hologenome dynamics in aquatic invertebrates by the holo-2bRAD approach

The "hologenome" concept is an increasingly popular way of thinking about microbiome-host for marine organisms. However, it is challenging to track hologenome dynamics because of the large amount of material, with tracking itself usually resulting in damage or death of the research object. Here we show the simple and efficient holo-2bRAD approach for the tracking of hologenome dynamics in marine invertebrates (i.e., scallop and shrimp) from one holo-2bRAD library. The stable performance of our approach was shown with high genotyping accuracy of 99.91% and a high correlation of r > 0.99 for the species-level profiling of microorganisms. To explore the host-microbe association underlying mass mortality events of bivalve larvae, core microbial species changed with the stages were found, and two potentially associated host SNPs were identified. Overall, our research provides a powerful tool with various advantages (e.g., cost-effective, simple, and applicable for challenging samples) in genetic, ecological, and evolutionary studies.

Analysis of Elimination Effects of Inbreeding on Genotype Frequency in Larval Stages of Chinese Shrimp

Marine animals possess genomes of considerable complexity and heterozygosity. Their unique reproductive system, characterized by high fecundity and substantial early mortality rates, increases the risk of inbreeding, potentially leading to severe inbreeding depression during various larval developmental stages. In this study, we established a set of inbred families of Fenneropenaeus chinensis, with an inbreeding coefficient of 0.25, and investigated elimination patterns and the manifestations of inbreeding depression during major larval developmental stages. Reduced-representation genome sequencing was utilized to explore the genotype frequency characteristics across two typical elimination stages. The results revealed notable mortality in hatching and metamorphosis into mysis and post-larvae stages. Inbreeding depression was also evident during these developmental stages, with depression rates of 24.36%, 29.23%, and 45.28%. Segregation analysis of SNPs indicated an important role of gametic selection before hatching, accounting for 45.95% of deviation in the zoea stage. During the zygotic selection phase of larval development, homozygote deficiency and heterozygote excess were the main selection types. Summation of the two types explained 82.31% and 89.91% of zygotic selection in the mysis and post-larvae stage, respectively. The overall distortion ratio decreased from 22.37% to 12.86% in the late developmental stage. A total of 783 loci were identified through selective sweep analysis. We also found the types of distortion at the same locus could change after the post-larvae stage. The predominant shifts included a transition of gametic selection toward normal segregation and other forms of distortion to heterozygous excess. This may be attributed to high-intensity selection on deleterious alleles and genetic hitchhiking effects. Following larval elimination, a greater proportion of heterozygous individuals were preserved. We detected an increase in genetic diversity parameters such as expected heterozygosity, observed heterozygosity, and polymorphic information content in the post-larvae stage. These findings suggest the presence of numerous recessive deleterious alleles and their linkage and suggest a major role of the partial dominance hypothesis. The results provide valuable insights into the mechanisms of inbreeding depression in marine animals and offer guidance for formulating breeding strategies in shrimp populations.

Signatures of selection in Mulinia lateralis underpinning its rapid adaptation to laboratory conditions

The dwarf surf clam, Mulinia lateralis, is considered as a model species for bivalves because of its rapid growth and short generation time. Recently, successful breeding of this species for multiple generations in our laboratory revealed its acquisition of adaptive advantages during artificial breeding. In this study, 310 individuals from five different generations were genotyped with 22,196 single nucleotide polymorphisms (SNPs) with the aim of uncovering the genetic basis of their adaptation to laboratory conditions. Results revealed that M. lateralis consistently maintained high genetic diversity across generations, characterized by high observed heterozygosity (H o: 0.2733-0.2934) and low levels of inbreeding (F is: -0.0244-0.0261). Population analysis indicated low levels of genetic differentiation among generations of M. lateralis during artificial breeding (F st <0.05). In total, 316 genomic regions exhibited divergent selection, with 168 regions under positive selection. Furthermore, 227 candidate genes were identified in the positive selection regions, which have functions including growth, stress resistance, and reproduction. Notably, certain selection signatures with significantly higher F st value were detected in genes associated with male reproduction, such as GAL3ST1, IFT88, and TSSK2, which were significantly upregulated during artificial breeding. This suggests a potential role of sperm-associated genes in the rapid evolutionary response of M. lateralis to selection in laboratory conditions. Overall, our findings highlight the phenotypic and genetic changes, as well as selection signatures, in M. lateralis during artificial breeding. This contributes to understanding their adaptation to laboratory conditions and underscores the potential for using this species to explore the adaptive evolution of bivalves.

Oral microbiota disorder in GC patients revealed by 2b-RAD-M

BACKGROUND

Microbiota alterations are linked with gastric cancer (GC). However, the relationship between the oral microbiota (especially oral fungi) and GC is not known. In this study, we aimed to apply 2b-RAD sequencing for Microbiome (2b-RAD-M) to characterize the oral microbiota in patients with GC.

METHODS

We performed 2b-RAD-M analysis on the saliva and tongue coating of GC patients and healthy controls. We carried out diversity, relative abundance, and composition analyses of saliva and tongue coating bacteria and fungi in the two groups. In addition, indicator analysis, the Gini index, and the mean decrease accuracy were used to identify oral fungal indicators of GC.

RESULTS

In this study, fungal imbalance in the saliva and tongue coating was observed in the GC group. At the species level, enriched Malassezia globosa (M. globosa) and decreased Saccharomyces cerevisiae (S. cerevisiae) were observed in saliva and tongue coating samples of the GC group. Random forest analysis indicated that M. globosa in saliva and tongue coating samples could serve as biomarkers to diagnose GC. The Gini index and mean decreases in accuracy for M. globosa in saliva and tongue coating samples were the largest. In addition, M. globosa in saliva and tongue coating samples classified GC from the control with areas under the receiver operating curve (AUCs) of 0.976 and 0.846, respectively. Further ecological analysis revealed correlations between oral bacteria and fungi.

CONCLUSION

For the first time, our data suggested that changes in oral fungi between GC patients and controls may help deepen our understanding of the complex spectrum of the different microbiotas involved in GC development. Although the cohort size was small, this study is the first to use 2b-RAD-M to reveal that oral M. globosa can be a fungal biomarker for detecting GC.

The bladder microbiome of NMIBC and MIBC patients revealed by 2bRAD-M

Background

Bladder cancer (BCa) is the most common malignancy of the urinary tract which can be divided into non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC), and their microbial differences are not fully understood. This study was conducted by performing 2bRAD sequencing for Microbiome (2bRAD-M) on NMIBC and MIBC tissue samples to investigate the microbiota differences between NMIBC and MIBC individuals.

Methods

A total of 22 patients with BCa, including 7 NMIBC and 15 MIBC, were recruited. Tumor tissues were surgically removed as samples and DNA was extracted. Type IIB restriction endonucleases were used to enzymatically cleave the microbial genome for each microbe's tag and map it to species-specific 2bRAD markers to enable qualitative and quantitative studies of microbes between MIBC and NMIBC tissues.

Results

A total of 527 species were detected. The microbial diversity of NMIBC tissues was significantly higher than that of MIBC tissues. Microbial composition of the two tumor tissues was similar, where Ralstonia_sp000620465 was the most dominant species. 4 species (Acinetobacter_guillouiae, Anoxybacillus_A_rupiensis, Brevibacillus_agri and Staphylococcus_lugdunensis) were enriched in NMIBC, while Ralstonia_mannitolilytica, Ralstonia_pickettii, and Ralstonia_sp000620465 were overrepresented in MIBC. 252 discriminatory character taxa were also revealed by linear discriminant analysis effect sizea (LEfSe). Species importance point plots identified Ralstonia_sp000620465, Cutibacterium_acnes and Ralstonia_pickettii as the three most important species between the two groups. Meanwhile, functional annotation analysis showed 3011 different COGs and 344 related signaling pathways between MIBC and NMIBC microbiome.

Conclusion

This first 2bRAD-M microbiome study on MIBC and NMIBC tissues revealed significant differences in the microbial environment between the two groups, which implies a potential association between tumor microbial dysbiosis and BCa, and provides a possible target and basis for subsequent studies on the mechanisms of BCa development and progression.

Genetic Diversity and Population Structure of Coilia nasus Revealed by 2b-RAD Sequencing

Coilia nasus is a threatened migratory species in the Yangtze River Basin. To reveal the genetic diversity of natural and farmed populations of C. nasus and the status of germplasm resources in the Yangtze River, the genetic diversity and structure of two wild populations (Yezhi Lake: YZ; Poyang Lake: PY) and two farmed populations (Zhenjiang: ZJ; Wuhan: WH) of C. nasus were analyzed using 44,718 SNPs obtained via 2b-RAD sequencing. The results indicate that both the wild and farmed populations had low genetic diversity, and germplasm resources have undergone varying degrees of degradation. Population genetic structure analyses indicated that the four populations may have come from two ancestral groups. Different amounts of gene flow were identified among WH, ZJ, and PY populations, but gene flow among YZ and other populations was low. It is speculated that the river-lake isolation of Yezhi Lake is the main cause of this phenomenon. In conclusion, this study revealed that genetic diversity reduction and germplasm resource degradation had occurred in both wild and farmed C. nasus, suggesting that conservation of its resources is of great urgency. This study provides a theoretical basis for the conservation and rational exploitation of germplasm resources for C. nasus.

Pesticide exposure affects DNA methylation patterns in natural populations of a mayfly

Chemical pollutants derived from agricultural activities represent a major threat to freshwater biota. Despite growing evidence involving epigenetic processes, such as DNA methylation, in response to pesticide contamination in agroecosystems, research on wild populations of non-model species remains scarce, particularly for endemic freshwater arthropods. Using the MethylRAD method, this study investigates whether exposure to pesticide contamination in natural populations of the endemic mayfly A. torrens produces genome wide changes in levels of DNA methylation. From a total of 1,377,147 MethylRAD markers produced from 285 specimens collected at 30 different study sites along the Limarí watershed of north-central Chile, six showed significant differential methylation between populations exposed and unexposed to pesticides. In all cases the effect of pesticides was positive, independent and stronger than the effects detected for other spatial and environmental factors. Only one candidate marker appeared correlated significantly with additional variables, nitrate and calcium levels, which also reflects the impact of agrichemicals and could additionally suggest, to a lower extent, antagonistic effects of mineral salts concentration for this specific marker. These results suggest that the effect of pesticide exposure on methylation levels is apparent at these six MethylRAD markers in A. torrens populations. Such data is challenging to obtain in natural populations and is, for the most part, lacking in ecotoxicological studies. Our study shows that DNA methylation processes are involved in the response to pesticide contamination in populations of the mayfly A. torrens in their natural habitat, and provides new evidence regarding the impact of pesticide contamination and agricultural activities on the endemic fauna of lotic ecosystems.

Screening and characterization of sex-specific sequences through 2b-RAD sequencing in American shad (Alosa sapidissima)

American shad (Alosa sapidissima), introduced from the United States, has become one of the most expensive farmed fish in the aquatic product market of China. The shad reveals significant sexual dimorphism in growth and behaviors. For the study, five male-specific tags were identified in two-generation breeding populations of Alosa sapidissima and were verified by PCR amplification. Averages of 10,245,091 and 8,685,704 raw and enzyme reads were obtained by high-throughput sequencing of the 2b-RAD library, respectively. 301,022 unique tags were obtained from the sequences of twenty samples with sequencing depths of 0 to 500. Finally, 274,324 special tags and 29,327 SNPs were selected with a sequencing depth of 3 to 500. Eleven preliminary screening male-specific tags and three male heterogametic SNP loci were isolated. After verification by PCR amplification, five male-specific sequences of 27 bp located on chromosome 3 were screened out. Chromosome 3 could be assumed to be the sex chromosome of Alosa sapidissima. Sex-specific markers will provide invaluable and systematic animal germplasm resources to allow for the precise identification of neo-males for the all-female breeding of Alosa sapidissima in commercial aquaculture.

Comparative phylogeography of two commensal rat species (Rattus tanezumi and Rattus norvegicus) in China: Insights from mitochondrial DNA, microsatellite, and 2b-RAD data

Rattus norvegicus and Rattus tanezumi are dominant species of Chinese house rats, but the colonization and demographic history of two species in China have not been thoroughly explored. Phylogenetic analyses with mitochondrial DNA including 486 individuals from 31 localities revealed that R. norvegicus is widely distributed in China, R. tanezumi is mainly distributed in southern China with currently invading northward; northeast China was the natal region of R. norvegicus, while the spread of R. tanezumi in China most likely started from the southeast coast. A total of 123 individuals from 18 localities were subjected to 2b-RAD analyses. In neighbor-joining tree, individuals of R. tanezumi grouped into geographic-specific branches, and populations from southeast coast were ancestral groups, which confirmed the colonization route from southeast coast to central and western China. However, individuals of R. norvegicus were generally grouped into two clusters instead of geographic-specific branches. One cluster comprised inland populations, and another cluster included both southeast coast and inland populations, which indicated that spread history of R. norvegicus in China was complex; in addition to on-land colonization, shipping transportation also have played great roles. ADMIXTURE and principal component analyses provided further supports for the colonization history. Demographic analyses revealed that climate changes at ~40,000 to 18,000 years ago and ~4000 years ago had led to population declines of both species; the R. norvegicus declined rapidly while the population of R. tanezumi continuously expanded since ~1500 years ago, indicating the importance of interspecies' competition in their population size changes. Our study provided a valuable framework for further investigation on phylogeography of two species in China.

The renal pelvis urobiome in the unilateral kidney stone patients revealed by 2bRAD-M

Background

The pathogenesis of kidney stone disease (KSD) is not fully understood, and potential contributing factors remain to be explored. Several studies have revealed that the urinary microbiome (urobiome) of stone formers was distinct from that of healthy individuals using 16S rRNA gene sequencing, most of which only provided microbial identification at the genus level. 2bRAD sequencing for Microbiome (2bRAD-M) is a novel sequencing technique that enables accurate characterization of the low-biomass microbiome at the species resolution. We aimed to apply 2bRAD-M to profile the renal pelvis urobiome of unilateral kidney stone patients and compared the urobiome with and without stone(s).

Method

A total of 30 patients with unilateral stones were recruited, and their renal pelvis urine from both sides was collected. A ureteroscope was inserted into the renal pelvis with stone(s) and a ureteral catheter was placed into the ureteroscope to collect renal pelvis urine. This procedure was repeated again with new devices to collect the urine of the other side. 2bRAD-M was performed to characterize the renal pelvis urobiome of unilateral stone formers to explore whether microbial differences existed between the stone side and the non-stone side.

Results

The microbial community composition of the stone side was similar to that of the non-stone side. Paired comparison showed that Corynebacterium was increased and Prevotella and Lactobacillus were decreased in the stone side. Four species (Prevotella bivia, Lactobacillus iners, Corynebacterium aurimucosum, and Pseudomonas sp_286) were overrepresented in the non-stone side. 24 differential taxa were also identified between two groups by linear discriminant analysis effect size (LEfSe). Extensive and close connections among genera and species were observed in the correlation analysis. Moreover, a random forest classifier was constructed using specific enriched species, which can distinguish the stone side from the non-stone side with an accuracy of 71.2%.

Conclusion

This first 2bRAD-M microbiome survey gave an important hint towards the potential role of urinary dysbiosis in KSD and provided a better understanding of mechanism of stone formation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03639-6.

Integrated Methylome and Transcriptome Analysis Provides Insights into the DNA Methylation Underlying the Mechanism of Cytoplasmic Male Sterility in Kenaf (Hibiscus cannabinus L.)

Cytoplasmic male sterility (CMS) is widely exploited in hybrid seed production. Kenaf is an important fiber crop with high heterosis. The molecular mechanism of kenaf CMS remains unclear, particularly in terms of DNA methylation. Here, using the anthers of a kenaf CMS line (P3A) and its maintainer line (P3B), comparative physiological, DNA methylation, and transcriptome analyses were performed. The results showed that P3A had considerably lower levels of IAA, ABA, photosynthetic products and ATP contents than P3B. DNA methylome analysis revealed 650 differentially methylated genes (DMGs) with 313 up- and 337 down methylated, and transcriptome analysis revealed 1788 differentially expressed genes (DEGs) with 558 up- and 1230 downregulated genes in P3A compared with P3B. Moreover, 45 genes were characterized as both DEGs and DMGs, including AUX,CYP, BGL3B, SUS6, AGL30 and MYB21. Many DEGs may be regulated by related DMGs based on methylome and transcriptome studies. These DEGs were involved in carbon metabolism, plant hormone signal transduction, the TCA cycle and the MAPK signaling pathway and were shown to be important for CMS in kenaf. These results provide new insights into the epigenetic mechanism of CMS in kenaf and other crops.

Genetic Diversity and Genetic Differentiation of Populations of Golden-Backed Carp (Cyprinus carpio var. Jinbei) in Traditional Rice Fields in Guizhou, China

Simple Summary

Rice-fish farming refers to the cultivation of fish in rice fields. China’s Guizhou Province has a long history (over 1000 years) of Rice-fish farming. Carp (golden-backed carp) that have been cultivated for generations in rice fields have evolved to be extremely adaptable to that environment. Unfortunately, farmers and enterprises have found that the feeding ability and growth rate of golden-backed carp are becoming increasingly weaker and slower, respectively. Therefore, the present study carried out to examine the population genetics of golden-backed, and the results indicated that the status of germplasm resources of this species is not very good. The aim of this study was to further elucidate the status of golden-backed carp germplasm resources in this province. In conclusion, this work can provide reliable fundamental data for the future conservation and breeding of golden-backed carp.

Abstract

The aim of this study was to assess the current status of the germplasm resources of golden-backed carp (Cyprinus carpio var. Jinbei) cultured in paddy fields in Guizhou Province, China. Five populations of golden-backed carp in Liping County, Jinping County, Huangping County, Congjiang County and Duyun City in Guizhou Province were subjected to high-throughput sequencing by 2b-RAD technology, and their genetic diversity and genetic differentiation were analysed. Based on sequencing, 44,896 SNP loci were obtained, and all five population genetic diversity indicators showed low diversity. In the NJ tree, the Congjiang and Liping populations were mixed together, and the other three groups formed a cluster. A cross-validation error box plot and pong cluster plot were constructed to show the K value results. When K = 1, the cross-validation error rate was the lowest. Principal component analysis showed that the Duyun population formed a group separate from the group comprising the other four populations. The genetic differentiation index and genetic distances between the Duyun population and the remaining four populations were greater than 0.05, indicating population differentiation. The genetic diversity of the five populations of golden-backed carp in Guizhou Province was low, the genetic differentiation of the Duyun population was the most significant, and the Duyun population was separate from the other four groups.

Genome-wide analysis of methylation in rat fetal heart under hyperglycemia by methylation-dependent restriction site–associated DNA sequencing

Diabetes mellitus causes an increased incidence of congenital heart malformations. However, the pathogenesis and potential epigenetic mechanism involved in this process are unclear. In this study, we used MethylRAD sequencing to compare changes in methylation levels in the genomic landscapes in the fetal heart in a rat model of hyperglycemia. Our results showed that methylation of CCGG/CCNGG sites were mostly enriched in intergenic regions, followed by intron, exon, upstream and the 5′ and 3′ untranslated regions. qRT-PCR results confirmed the MethylRAD sequencing findings, suggesting that abnormal CCGG/CCNGG methylation in the upstream region regulated gene expression. The differential methylation genes (DMGs) based on the CCGG and CCNGG sites in the upstream region were examined by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. Gene Ontology indicated that the CCGG-based DMGs involved in biological process and function were mainly related to transcription and co-SMAD binding. The CCNGG-based DMGs were mainly related to transcription and cytokine-mediated signaling pathways. Kyoto Encyclopedia of Genes and Genomes analysis indicated that CCGG-based DMGs were mainly involved in the Wnt signaling and TGF-β signaling pathways. CCNGG-based DMGs were involved in the TNF signaling and apoptosis pathways. These genes may play dominant roles in cardiomyocyte apoptosis and heart disease and require further study. These genes may also serve as potential molecular targets or diagnostic biomarkers for heart malformations under hyperglycemia.

IIb-RAD-sequencing coupled with random forest classification indicates regional population structuring and sex-specific differentiation in salmon lice (Lepeophtheirus salmonis)

The aquaculture industry has been dealing with salmon lice problems forming serious threats to salmonid farming. Several treatment approaches have been used to control the parasite. Treatment effectiveness must be optimized, and the systematic genetic differences between subpopulations must be studied to monitor louse species and enhance targeted control measures. We have used IIb-RAD sequencing in tandem with a random forest classification algorithm to detect the regional genetic structure of the Norwegian salmon lice and identify important markers for sex differentiation of this species. We identified 19,428 single nucleotide polymorphisms (SNPs) from 95 individuals of salmon lice. These SNPs, however, were not able to distinguish the differential structure of lice populations. Using the random forest algorithm, we selected 91 SNPs important for geographical classification and 14 SNPs important for sex classification. The geographically important SNP data substantially improved the genetic understanding of the population structure and classified regional demographic clusters along the Norwegian coast. We also uncovered SNP markers that could help determine the sex of the salmon louse. A large portion of the SNPs identified to be under directional selection was also ranked highly important by random forest. According to our findings, there is a regional population structure of salmon lice associated with the geographical location along the Norwegian coastline.

A High-Density Genetic Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Red Swamp Crayfish (Procambarus clarkii)

Red swamp crayfish (Procambarus clarkii) is a commercially important species in global aquaculture and most successfully invasive freshwater shrimp in China. In order to determine the genetic basis of growth- and sex-related traits, a high-density genetic linkage map was constructed using 2b-RAD sequencing technology in a full-sib family. The consensus map contains 4,878 SNP markers assigned to 94 linkage groups (LGs) and spanned 6,157.737 cM with an average marker interval of 1.26 cM and 96.93% genome coverage. The quantitative trait locus (QTL) mapping for growth and sex traits was performed for the first time. QTL mapping uncovers 28 QTLs for growth-related traits in nine LGs, explaining 7.9-14.4% of the phenotypic variation, and identifies some potential candidate growth-related genes such as mih, lamr, golgb1, nurf301, and tbcd1 within the QTL intervals. A single major locus for sex determination was revealed in LG20 that explains 59.3-63.7% of the phenotypic variations. Some candidate sex-related genes, such as vps4bl, ssrf, and acot1, were identified in the QTL intervals and found to be differentially expressed in the muscle tissues between the females and the males. Furthermore, the identified SNPs were revealed to be female heterozygotes, suggesting that red swamp crayfish might have the female heterogametic ZZ/ZW sex determination system. The present study provides a valuable resource for marker-assisted selection and genetic improvement and for further genetic and genomic research in red swamp crayfish.

Genetic diversity of Chinese laboratory macaques based on 2b-RAD simplified genome sequencing

BACKGROUND

Currently, Chinese laboratory macaques are widely used in biomedical research. Correspondingly, clarity regarding the genetic diversity of Chinese laboratory macaques is important for both vendors and users.

METHODS

Genome sequences of 55 laboratory macaques (40 cynomolgus macaques and 15 rhesus macaques) bred in South China were analyzed using 2b-RAD simplified genome sequencing.

RESULTS

A total of 115,681 single-nucleotide polymorphisms (SNPs) were found that were distributed in 21 chromosomes and an unplaced scaffold. Genetic diversity indices varied across populations and exhibited low values. The results of principal coordinate analysis (PCA) were consistent with those of the arithmetic mean (UPGMA) clustered tree and supported the structure analysis, demonstrating that the genetic differentiation in rhesus macaques was higher than that in cynomolgus macaques. Introgressive hybridization with the Chinese rhesus macaque was supported in more than 80% (32/40) of cynomolgus macaques.

CONCLUSIONS

Chinese laboratory macaques had relatively low genetic diversity at the genomic level, and genetic differentiation in Chinese rhesus macaques was higher than in cynomolgus macaques. The genome of cynomolgus macaques has been shaped by introgression after hybridization with the Chinese rhesus macaques.

Species-resolved sequencing of low-biomass or degraded microbiomes using 2bRAD-M

Microbiome samples with low microbial biomass or severe DNA degradation remain challenging for amplicon-based or whole-metagenome sequencing approaches. Here, we introduce 2bRAD-M, a highly reduced and cost-effective strategy which only sequences ~ 1% of metagenome and can simultaneously produce species-level bacterial, archaeal, and fungal profiles. 2bRAD-M can accurately generate species-level taxonomic profiles for otherwise hard-to-sequence samples with merely 1 pg of total DNA, high host DNA contamination, or severely fragmented DNA from degraded samples. Tests of 2bRAD-M on various stool, skin, environmental, and clinical FFPE samples suggest a successful reconstruction of comprehensive, high-resolution microbial profiles.

Linkage mapping and QTL analysis of growth traits in Rhopilema esculentum

R. esculentum is a popular seafood in Asian countries and an economic marine fishery resource in China. However, the genetic linkage map and growth-related molecular markers are still lacking, hindering marker assisted selection (MAS) for genetic improvement of R. esculentum. Therefore, we firstly used 2b-restriction site-associated DNA (2b-RAD) method to sequence 152 R. esculentum specimens and obtained 9100 single nucleotide polymorphism (SNP) markers. A 1456.34 cM linkage map was constructed using 2508 SNP markers with an average interval of 0.58 cM. Then, six quantitative trait loci (QTLs) for umbrella diameter and body weight were detected by QTL analysis based on the new linkage map. The six QTLs are located on four linkage groups (LGs), LG4, LG13, LG14 and LG15, explaining 9.4% to 13.4% of the phenotypic variation. Finally, 27 candidate genes in QTLs regions of LG 14 and 15 were found associated with growth and one gene named RE13670 (sushi, von Willebrand factor type A, EGF and pentraxin domain-containing protein 1-like) may play an important role in controlling the growth of R. esculentum. This study provides valuable information for investigating the growth mechanism and MAS breeding in R. esculentum.

Characterization of novel genotyping-by-sequencing (GBS)-based simple sequence repeats (SSRs) and their application for population genomics of Capoeta aculeata (Valenciennes, 1844)

BACKGROUND

The species Capoeta aculeata (Valenciennes, 1844) is one of the most important freshwater species endemic to Iran. However, the investigation of the population genetic structure of this species is limited by the low number of molecular markers currently described.

METHODS AND RESULTS

In this study, we implemented next generation sequencing technology to identify polymorphic microsatellite markers and investigate the population genetic structure of C. aculeata sampled from three geographical sites in Iran. We characterized and developed 36 novel polymorphic microsatellite markers and these loci were examined in 120 individuals from three populations occurring in the Zagros basin. The average number of alleles per locus varied from 1.7 to 16 (average = 7.89). The results showed that, the polymorphism information content (PIC) of these simple sequence repeat (SSR) loci varied from 0.254 to 0.888. The observed heterozygosity (H_O) per locus ranged from 0.170 to 0.881, while the expected heterozygosity (H_E) per locus was from 0.170 to 0.881. Among these SSR loci, 20 loci deviated significantly from the Hardy-Weinberg equilibrium after Bonferroni correction (p < 0.05).

CONCLUSIONS

These microsatellite markers could provide a valuable tool for future population and conservation genetics studies of C. aculeate and other closely related species.

Bacterial diversity in ballast water and sediments revealed by 2b-RAD sequencing

Despite that the ballast water management (BWM) convention has come into force to prevent the spread of harmful aquatic organisms, to date, very few bacteria can be identified through microbial culture method. In this study, we explored a reduced-representation sequencing of 2b-RAD approach to investigate the bacterial diversity in ballast water and sediments (BWS). Our results indicated a large amount of bacteria species (1496) detected in BWS up to now, including 13 pathogens that are seriously concerning in marine environment and aquaculture like the most harmful Vibrio harveyi and Aurantimonas coralicida. We showed that the ballast water had relative lower species, which was dominated by Proteobacteria. In contrast, the sediments had richer species, which was dominated by Bacteroidetes. Although BWS differed significantly in species composition, sediments shared most of the concerned pathogens with ballast water, highlighting the importance of sediment management. In conclusion, 2b-RAD sequencing shows promise in future BWM.

Genome-Wide Association Study Reveals PC4 as the Candidate Gene for Thermal Tolerance in Bay Scallop (Argopecten irradians irradians)

The increasing sea temperature caused by global warming has resulted in severe mortalities in maricultural scallops. Therefore, improving thermal tolerance has become an active research area in the scallop farming industry. Bay scallop (Argopecten irradians irradians) was introduced into China in 1982 and has developed into a vast aquaculture industry in northern China. To date, genetic studies on thermal tolerance in bay scallops are limited, and no systematic screening of thermal tolerance-related loci or genes has been conducted in this species. In the present study, we conducted a genome-wide association study (GWAS) for thermal tolerance using the Arrhenius break temperature (ABT) indicators of 435 bay scallops and 38,011 single nucleotide polymorphism (SNP) markers. The GWAS identified 1,906 significant thermal tolerance-associated SNPs located in 16 chromosomes of bay scallop. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that 638 genes were enriched in 42 GO terms, while 549 annotated genes were enriched in aggregation pathways. Additionally, the SNP (15-5091-20379557-1) with the lowest P value was located in the transcriptional coactivator p15 (PC4) gene, which is involved in regulating DNA damage repair and stabilizing genome functions. Further analysis in another population identified two new thermal tolerance-associated SNPs in the first coding sequence of PC4 in bay scallops (AiPC4). Moreover, AiPC4 expression levels were significantly correlated (r = 0.675–0.962; P < 0.05) with the ABT values of the examined bay scallops. Our data suggest that AiPC4 might be a positive regulator of thermal tolerance and a potential candidate gene for molecular breeding in bay scallop aiming at thermal tolerance improvement.

Conservation Genomics of Wild Red Sage (Salvia miltiorrhiza) and Its Endangered Relatives in China: Population Structure and Interspecific Relationships Revealed From 2b-RAD Data

Red sage (Salvia miltiorrhiza) is a widely used medicinal plant for treatment of cardiovascular and cerebrovascular diseases. Because of excessive excavation by huge market demand and habitat loss by human activities, the wild population resources of S. miltiorrhiza have reduced drastically in recent years. Meanwhile, population status of two closely related species S. bowleyana and S. paramiltiorrhiza were in a trend of decreasing due to their potential replacement of S. miltiorrhiza. Particularly, S. paramiltiorrhiza was threatened and endemic to a small region in eastern China. However, to date there has been no conservation genetic research reported for wild S. miltiorrhiza population and its endangered relatives. Assess the wild germplasm diversity for S. miltiorrhiza and its related species would provide fundamental genetic background for cultivation and molecular breeding of this medicinally important species. In the present study, we investigated the genetic diversity, population structure, and intra/inter-specific differentiation of S. miltiorrhiza and above two relatives using 2b-RAD genome-wide genotyping method. By investigating 81 individuals of S. miltiorrhiza, 55 individuals of S. bowleyana and 15 individuals of S. paramiltiorrhiza from 23 locations in China, we obtained 23,928 SNPs in total. A comparatively high genetic diversity was observed in S. miltiorrhiza (π = 0.0788, H e = 0.0783 ± 0.0007). The observed and expected heterozygosity in populations of these three species ranged from 0.0297 to 0.1481 and 0.0251 to 0.831, respectively. Two major lineage groups were detected in the examined S. miltiorrhiza populations. The results indicated that Dabie Mountain as a genetic diversity center of S. miltiorrhiza and possible complex inter-specific genetic exchange/hybridization occurred between S. miltiorrhiza and the two relatives. We suggest that strategic conservation and germplasm preservation should be considered not only for wild populations of S. miltiorrhiza, but also for its related S. bowleyana and S. paramiltiorrhiza.

The Genetic Component of Seagrass Restoration: What We Know and the Way Forwards

Seagrasses are marine flowering plants providing key ecological services and functions in coasts and estuaries across the globe. Increased environmental changes fueled by human activities are affecting their existence, compromising natural habitats and ecosystems’ biodiversity and functioning. In this context, restoration of disturbed seagrass environments has become a worldwide priority to reverse ecosystem degradation and to recover ecosystem functionality and associated services. Despite the proven importance of genetic research to perform successful restoration projects, this aspect has often been overlooked in seagrass restoration. Here, we aimed to provide a comprehensive perspective of genetic aspects related to seagrass restoration. To this end, we first reviewed the importance of studying the genetic diversity and population structure of target seagrass populations; then, we discussed the pros and cons of different approaches used to restore and/or reinforce degraded populations. In general, the collection of genetic information and the development of connectivity maps are critical steps for any seagrass restoration activity. Traditionally, the selection of donor population preferred the use of local gene pools, thought to be the best adapted to current conditions. However, in the face of rapid ocean changes, alternative approaches such as the use of climate-adjusted or admixture genotypes might provide more sustainable options to secure the survival of restored meadows. Also, we discussed different transplantation strategies applied in seagrasses and emphasized the importance of long-term seagrass monitoring in restoration. The newly developed information on epigenetics as well as the application of assisted evolution strategies were also explored. Finally, a view of legal and ethical issues related to national and international restoration management is included, highlighting improvements and potential new directions to integrate with the genetic assessment. We concluded that a good restoration effort should incorporate: (1) a good understanding of the genetic structure of both donors and populations being restored; (2) the analysis of local environmental conditions and disturbances that affect the site to be restored; (3) the analysis of local adaptation constraints influencing the performances of donor populations and native plants; (4) the integration of distribution/connectivity maps with genetic information and environmental factors relative to the target seagrass populations; (5) the planning of long-term monitoring programs to assess the performance of the restored populations. The inclusion of epigenetic knowledge and the development of assisted evolution programs are strongly hoped for the future.

Exploring the genetic basis of fatty liver development in geese

Although geese possess an adaptive physiological capacity for lipid storage, few candidate genes contributing to this ability are characterised. By comparing the genomes of individuals with extremely high and low fatty liver weights (FLW), candidate genes were identified, including ARAP2, GABRE, and IL6. Single-nucleotide polymorphisms in or near these genes were significantly (p < 0.05) associated with carcass traits (FLW) and biochemical indexes (very-low-density lipoprotein and N-terminal procollagen III), suggesting contribution to trait variation. A common variant at the 5'-end of LCORL explained ~ 18% and ~ 26% of the phenotypic variance in body weight with/without overfeeding and had significant effects on FLW (p < 0.01). ZFF36L1, ARHGEF1 and IQCJ, involved in bile acid metabolism, blood pressure, and lipid concentration modulation, were also identified. The presence of highly divergent haplotypes within these genes suggested involvement in protection against negative effects from excessive lipids in the liver or circulatory system. Based on this and transcriptomic data, we concluded that geese hepatosteatosis results from severe imbalance between lipid accumulation and secretion, comparable to human non-alcohol fatty liver disease but involving other genes. Our results provided valuable insights into the genesis of geese fatty liver and detected potential target genes for treatment of lipid-related diseases.

First High-Density Linkage Map and QTL Fine Mapping for Growth-Related Traits of Spotted Sea bass (Lateolabrax maculatus)

Possessing powerful adaptive capacity and a pleasant taste, spotted sea bass (Lateolabrax maculatus) has a broad natural distribution and is one of the most popular mariculture fish in China. However, the genetic improvement program for this fish is still in its infancy. Growth is the most economically important trait and is controlled by quantitative trait loci (QTL); thus, the identification of QTLs and genetic markers for growth-related traits is an essential step for the establishment of marker-assisted selection (MAS) breeding programs. In this study, we report the first high-density linkage map of spotted sea bass constructed by sequencing 333 F1 generation individuals in a full-sib family using 2b-RAD technology. A total of 6883 SNP markers were anchored onto 24 linkage groups, spanning 2189.96 cM with an average marker interval of 0.33 cM. Twenty-four growth-related QTLs, including 13 QTLs for body weight and 11 QTLs for body length, were successfully detected, with phenotypic variance explained (PVE) ranging from 5.1 to 8.6%. Thirty potential candidate growth-related genes surrounding the associated SNPs were involved in cell adhesion, cell proliferation, cytoskeleton reorganization, calcium channels, and neuromodulation. Notably, the fgfr4 gene was detected in the most significant QTL; this gene plays a pivotal role in myogenesis and bone growth. The results of this study may facilitate marker-assisted selection for breeding populations and establish the foundation for further genomic and genetic studies investigating spotted sea bass.

High accuracy of pooled DNA genotyping by 2b-RAD sequencing in the Pacific white shrimp, Litopenaeus vannamei

Using pooled DNA genotyping to estimate the proportional contributions from multiple families in a pooled sample is of particular interest for selective breeding in aquaculture. We compared different pooled libraries with separate 2b-RAD sequencing of Litopenaeus vannamei individuals to assess the effect of different population structures (different numbers of individuals and families) on pooled DNA sequencing, the accuracy of parent sequencing of the DNA pools and the effect of SNP numbers on pooled DNA sequencing. We demonstrated that small pooled DNA genotyping of up to 53 individuals by 2b-RAD sequencing could provide a highly accurate assessment of population allele frequencies. The accuracy increased as the number of individuals and families increased. The allele frequencies of the parents from each pool were highly correlated with those of the pools or the corresponding individuals in the pool. We chose 500-28,000 SNPs to test the effect of SNP number on the accuracy of pooled sequencing, and no linear relationship was found between them. When the SNP number was fixed, increasing the number of individuals in the mixed pool resulted in higher accuracy of each pooled genotyping. Our data confirmed that pooled DNA genotyping by 2b-RAD sequencing could achieve higher accuracy than that of individual-based genotyping. The results will provide important information for shrimp breeding programs.

Three Novel Players: PTK2B, SYK, and TNFRSF21 Were Identified to Be Involved in the Regulation of Bovine Mastitis Susceptibility via GWAS and Post-transcriptional Analysis

Bovine mastitis is a common inflammatory disease caused by multiple factors in early lactation or dry period. Genome wide association studies (GWAS) can provide a convenient and effective strategy for understanding the biological basis of mastitis and better prevention. 2b-RADseq is a high-throughput sequencing technique that offers a powerful method for genome-wide genetic marker development and genotyping. In this study, single nucleotide polymorphisms (SNPs) of the immune-regulated gene correlative with mastitis were screened and identified by two stage association analysis via GWAS-2b-RADseq in Chinese Holstein cows. We have screened 10,058 high quality SNPs from 7,957,920 tags and calculated their allele frequencies. Twenty-seven significant SNPs were co-labeled in two GWAS analysis models [Bayesian (P < 0.001) and Logistic regression (P < 0.01)], and only three SNPs (rs75762330, C > T, PIC = 0.2999; rs88640083, A > G, PIC = 0.1676; rs20438858, G > A, PIC = 0.3366) were annotated to immune-regulated genes (PTK2B, SYK, and TNFRSF21). Identified three SNPs are located in non-coding regions with low or moderate genetic polymorphisms. However, independent sample population validation (Case-control study) data showed that three important SNPs (rs75762330, P < 0.025, OR > 1; rs88640083, P < 0.005, OR > 1; rs20438858, P < 0.001, OR < 1) were significantly associated with clinical mastitis trait. Importantly, PTK2B and SYK expression was down-regulated in both peripheral blood leukocytes (PBLs) of clinical mastitis cows and in vitro LPS (E. coli)-stimulated bovine mammary epithelial cells, while TNFRSF21 was up-regulated. Under the same conditions, expression of Toll-like receptor 4 (TLR4), AKT1, and pro-inflammatory factors (IL-1β and IL-8) were also up-regulated. Interestingly, network analysis indicated that PTK2B and SYK are co-expressed in innate immune signaling pathway of Chinese Holstein. Taken together, these results provided strong evidence for the study of SNPs in bovine mastitis, and revealed the role of SYK, PTK2B, and TNFRSF21 in bovine mastitis susceptibility/tolerance.

Contrasting fine-scale genetic structure of two sympatric clonal plants in an alpine swampy meadow featured by tussocks

Tussocks are unique vegetation structures in wetlands. Many tussock species mainly reproduce by clonal growth, resulting in genetically identical offspring distributed in various spatial patterns. These fine-scale patterns could influence mating patterns and thus the long-term evolution of wetland plants. Here, we contribute the first genetic and clonal structures of two key species in alpine wetlands on the Qinghai–Tibet Plateau, Kobresia tibetica and Blysmus sinocompressus, using > 5000 SNPs identified by 2b-RAD sequencing. The tussock-building species, K. tibetica, has a phalanx (clumping) growth form, but different genets could co-occur within the tussocks, indicating that it is not proper to treat a tussock as one genetic individual. Phalanx growth does not necessarily lead to increased inbreeding in K. tibetica. B. sinocompressus has a guerilla (spreading) growth form, with the largest detected clone size being 18.32 m, but genets at the local scale tend to be inbred offspring. Our results highlight that the combination of clone expansion and seedling recruitment facilitates the contemporary advantage of B. sinocompressus, but its evolutionary potential is limited by the input genetic load of the original genets. The tussocks of K. tibetica are more diverse and a valuable genetic legacy of former well-developed wet meadows, and they are worthy of conservation attention.

Predicting Growth Traits with Genomic Selection Methods in Zhikong Scallop (Chlamys farreri)

Selective breeding is a common and effective approach for genetic improvement of aquaculture stocks with parental selection as the key factor. Genomic selection (GS) has been proposed as a promising tool to facilitate selective breeding. Here, we evaluated the predictability of four GS methods in Zhikong scallop (Chlamys farreri) through real dataset analyses of four economical traits (e.g., shell length, shell height, shell width, and whole weight). Our analysis revealed that different GS models exhibited variable performance in prediction accuracy depending on genetic and statistical factors, but non-parametric method, including reproducing kernel Hilbert spaces regression (RKHS) and sparse neural networks (SNN), generally outperformed parametric linear method, such as genomic best linear unbiased prediction (GBLUP) and BayesB. Furthermore, we demonstrated that the predictability relied mainly on the heritability regardless of GS methods. The size of training population and marker density also had considerable effects on the predictive performance. In practice, increasing the training population size could better improve the genomic prediction than raising the marker density. This study is the first to apply non-linear model and neural networks for GS in scallop and should be valuable to help develop strategies for aquaculture breeding programs.

HD-Marker: a highly multiplexed and flexible approach for targeted genotyping of more than 10,000 genes in a single-tube assay

Targeted genotyping of transcriptome-scale genetic markers is highly attractive for genetic, ecological, and evolutionary studies, but achieving this goal in a cost-effective manner remains a major challenge, especially for laboratories working on nonmodel organisms. Here, we develop a high-throughput, sequencing-based GoldenGate approach (called HD-Marker), which addresses the array-related issues of original GoldenGate methodology and allows for highly multiplexed and flexible targeted genotyping of more than 12,000 loci in a single-tube assay (in contrast to fewer than 3100 in the original GoldenGate assay). We perform extensive analyses to demonstrate the power and performance of HD-Marker on various multiplex levels (296, 795, 1293, and 12,472 genic SNPs) across two sequencing platforms in two nonmodel species (the scallops Chlamys farreri and Patinopecten yessoensis), with extremely high capture rate (98%-99%) and genotyping accuracy (97%-99%). We also demonstrate the potential of HD-Marker for high-throughput targeted genotyping of alternative marker types (e.g., microsatellites and indels). With its remarkable cost-effectiveness (as low as $0.002 per genotype) and high flexibility in choice of multiplex levels and marker types, HD-Marker provides a highly attractive tool over array-based platforms for fulfilling genome/transcriptome-wide targeted genotyping applications, especially in nonmodel organisms.

Whole-Genome Restriction Mapping by "Subhaploid"-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding

Assembly of complex genomes using short reads remains a major challenge, which usually yields highly fragmented assemblies. Generation of ultradense linkage maps is promising for anchoring such assemblies, but traditional linkage mapping methods are hindered by the infrequency and unevenness of meiotic recombination that limit attainable map resolution. Here we develop a sequencing-based "in vitro" linkage mapping approach (called RadMap), where chromosome breakage and segregation are realized by generating hundreds of "subhaploid" fosmid/bacterial-artificial-chromosome clone pools, and by restriction site-associated DNA sequencing of these clone pools to produce an ultradense whole-genome restriction map to facilitate genome scaffolding. A bootstrap-based minimum spanning tree algorithm is developed for grouping and ordering of genome-wide markers and is implemented in a user-friendly, integrated software package (AMMO). We perform extensive analyses to validate the power and accuracy of our approach in the model plant Arabidopsis thaliana and human. We also demonstrate the utility of RadMap for enhancing the contiguity of a variety of whole-genome shotgun assemblies generated using either short Illumina reads (300 bp) or long PacBio reads (6-14 kb), with up to 15-fold improvement of N50 (∼816 kb-3.7 Mb) and high scaffolding accuracy (98.1-98.5%). RadMap outperforms BioNano and Hi-C when input assembly is highly fragmented (contig N50 = 54 kb). RadMap can capture wide-range contiguity information and provide an efficient and flexible tool for high-resolution physical mapping and scaffolding of highly fragmented assemblies.