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

Premise: Detecting single-nucleotide polymorphisms (SNPs) in a cost-effective way is fundamental in any plant breeding pipeline. Here, we compare three genotyping techniques for their ability to reproduce the allele dosage of SNPs of interest in sugarcane (Saccharum spp.). Methods: To identify a reproducible technique to estimate allele dosage for the validation of SNP markers, the correlation between Flex-Seq, kompetitive allele-specific PCR (KASP), and genotyping-by-sequencing and restriction site-associated DNA sequencing (GBS+RADseq) was determined for a set of 76 SNPs. To find alternative methodologies for allele dosage estimation, the KASP and Flex-Seq techniques were compared for the same set of SNPs. For the three techniques, a population of 53 genotypes from the diverse sugarcane panel of the Centro de Investigación de la Caña de Azúcar (Cenicaña), Colombia, was selected. Results: The average Pearson correlation coefficients between GBS+RADseq and Flex-Seq, GBS+RADseq and KASP, and Flex-Seq and KASP were 0.62 ± 0.27, 0.38 ± 0.27, and 0.38 ± 0.30, respectively. Discussion: Flex-Seq reproduced the allele dosages determined using GBS+RADseq with good levels of precision because of its depth of sequencing and ability to target specific positions in the genome. Additionally, Flex-Seq outperformed KASP by allowing the conversion of a higher number of SNPs and a more accurate estimation of the allele dosage. Flex-Seq has therefore become the genotyping methodology of choice for marker validation at Cenicaña.

Premisa: Detectar polimorfismos de un único nucleótido (SNP) de forma costo­efectiva es fundamental en cualquier programa de mejoramiento genético. En este artículo nosotros comparamos tres técnicas de genotipado para medir su habilidad en reproducir las dosis alélicas de SNPs de interés en caña de azúcar (Saccharum spp.). Métodos: Para identificar una técnica reproducible para la estimación de dosis alélicas durante los pasos de validación de marcadores, la correlación entre Flex­Seq, kompetitive allele­specific PCR (KASP), y genotyping­by­sequencing and restriction site­associated DNA sequencing (GBS+RADseq) fue determinada para un set de 76 SNPs. Para identificar metodologías alternativas en la estimación de las dosis alélicas, las tecnologías KASP y Flex­Seq fueron comparadas para el mismo grupo de SNPs. Para las tres técnicas, una población de 53 genotipos fue seleccionados de la población diversa de caña de azúcar del Centro de Investigación de la Caña de Azúcar (Cenicaña), Colombia. Resultados: El promedio del coeficiente de correlación de Pearson entre GBS+RADseq y Flex­Seq, GBS+RADseq y KASP, y Flex­Seq y KASP fue de 0.62 ± 0.27, 0.38 ± 0.27, y 0.38 ± 0.30, respectivamente. Discusión: Flex­Seq reprodujo las dosis alélicas determinadas usando GBS+RADseq con buenos niveles de precisión debido a su profundidad de secuenciación y habilidad de secuenciar posiciones especificas en el genoma. Adicionalmente, Flex­Seq superó a KASP al permitir la conversión de un número mayor de SNPs y al estimar las dosis alélicas de forma más precisa. Flex­Seq por tanto se convierte en la metodología de genotipado de elección para la validación de marcadores en Cenicaña.

High-Density Genetic Map Construction and Quantitative Trait Locus Analysis of Fruit- and Oil-Related Traits in Camellia oleifera Based on Double Digest Restriction Site-Associated DNA Sequencing.

Camellia oleifera, an important tree species and source of edible oil in China, has received significant attention owing to the oil's high unsaturated fatty acid content, which has benefits for human health. However, the mechanisms underlying C. oleifera yield and oil quality are largely unknown. In this study, 180 F1 progenies were obtained from two parents with obvious differences in fruit- and oil-related traits. We constructed a high-density genetic map using a double digest restriction site-associated DNA sequencing (ddRAD-Seq) strategy in C. oleifera. This map spanned 3327 cM and anchored 2780 markers in 15 linkage groups (LGs), with an average marker interval of 1.20 cM. A total of 221 quantitative trait loci (QTLs) associated with fruit- and oil-related traits were identified across three years' worth of phenotypic data. Nine QTLs were detected simultaneously in at least two different years, located on LG02, LG04, LG05, LG06, and LG11, and explained 8.5-16.6% of the phenotypic variation in the corresponding traits, respectively. Seventeen major QTLs were obtained that explained 13.0-16.6% of the phenotypic variance. Eleven and five flanking SNPs of major QTLs for fruit- and oil-related traits were detected which could be used for marker-assisted selection in C. oleifera breeding programs. Furthermore, 202 potential candidate genes in QTL regions were identified based on the collinearity of the genetic map and the C. oleifera "CON" genome. A potential regulatory network controlling fruit development and oil biosynthesis was constructed to dissect the complex mechanism of oil accumulation. The dissection of these QTLs will facilitate the gene cloning underlying lipid synthesis and increase our understanding in order to enhance C. oleifera oil yield and quality.

A mutation in POLR2A gene associated with body size traits in Dezhou donkeys revealed with GWAS.

The Dezhou donkey is a famous local donkey breed in China. The aim of the present study was to identify the genes associated with the body size traits of the Dezhou donkey and facilitate the breeding activities of the donkeys. A total of 349 donkeys from 2 generations (113 individuals in F0 and 236 in F1) were analyzed with restriction-site-associated DNA sequencing. A genome-wide association study revealed that the region between 13.7 and 15.6 Mb of chromosome 13 is significantly associated with body sizes. Candidate genes related to body size development, including POLR2A, CHRNB1, FGF11, and ZBTB4, were identified. The results of GO and KEGG analysis indicated that the genes involved in many GO terms were related to metabolic processes and developmental processes. Additionally, a T>C mutation (Chr13:14312485) was found at intron 10 of the POLR2A gene. The association analysis showed significant differences among genotypes for the size traits. The body size of the individuals with the TT genotype was significantly higher than that with the CC genotype. The results showed that the polymorphism of POLR2A has the potential to be used as a marker in the breeding programs of the Dezhou donkeys.

Body size is a crucial economic trait in donkeys, as it is closely related to meat and skin production. The aim of this study was to identify the genes and loci associated with body size traits, using the Dezhou donkey as an experimental population. The study findings make contributions to a better understanding on the molecular genetic mechanism of body size traits. The significant loci screened out in the present study may facilitate gene-assisted selection breeding and accelerate genetic selection, which is of great significance to the breeding of donkeys and the development of the donkey industry.

Population genetics and phylogeographic history of the insular lizard Podarcis lilfordi (Gunther, 1874) from the Balearic Islands based on genome-wide polymorphic data.

Islands provide a great system to explore the processes that maintain genetic diversity and promote local adaptation. We explored the genomic diversity of the Balearic lizard Podarcis lilfordi, an endemic species characterized by numerous small insular populations with large phenotypic diversity. Using the newly available genome for this species, we characterized more than 300,000 SNPs, merging genotyping-by-sequencing (GBS) data with previously published restriction site-associated DNA sequencing (RAD-Seq) data, providing a dataset of 16 island populations (191 individuals) across the range of species distribution (Menorca, Mallorca, and Cabrera). Results indicate that each islet hosts a well-differentiated population (F ST = 0.247 ± 0.09), with no recent immigration/translocation events. Contrary to expectations, most populations harbor a considerable genetic diversity (mean nucleotide diversity, P i = 0.144 ± 0.021), characterized by overall low inbreeding values (F IS < 0.1). While the genetic diversity significantly decreased with decreasing islet surface, maintenance of substantial genetic diversity even in tiny islets suggests variable selection or other mechanisms that buffer genetic drift. Maximum-likelihood tree based on concatenated SNP data confirmed the existence of the two major independent lineages of Menorca and Mallorca/Cabrera. Multiple lines of evidence, including admixture and root testing, robustly placed the origin of the species in the Mallorca Island, rather than in Menorca. Outlier analysis mainly retrieved a strong signature of genome differentiation between the two major archipelagos, especially in the sexual chromosome Z. A set of proteins were target of multiple outliers and primarily associated with binding and catalytic activity, providing interesting candidates for future selection studies. This study provides the framework to explore crucial aspects of the genetic basis of phenotypic divergence and insular adaptation.

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

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

Population genetics of Camellia granthamiana, an endangered plant species with extremely small populations in China.

Introduction: Camellia, the largest genus of Theaceae, is well-known for having high economic values. Camellia granthamiana demonstrates large beautiful flowers with some primitive characters, such as multiple large and persistent bracteoles and sepals, was listed as Vulnerable species on the IUCN Red List. Methods: In this study, we investigated all possible records of the species, and sampled four natural populations and five cultivated individuals. By applying shallow-genome sequencing for nine individuals and RAD-seq sequencing for all the sampled 77 individuals, we investigated population genetic diversity and population structure of the species. Results and discussion: The results showed that the population sampled from Fengkai, previously identified as C. albogigias, possessed different plastid genome from other species possibly due to plastid capture; the species possesses strong population structure possibly due to the effect of isolation by distance, habitat fragmentation, and self-crossing tendency of the species, whose effective population size declined quickly in the past 4,000 years. Nevertheless, C. granthamiana maintains a medium level of genetic diversity within population, and significant differentiation was observed among the four investigated populations, it is anticipated that more populations are expected to be found and all these extant populations should be taken into instant protection.

2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics.

Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.

Genotyping by Sequencing (GBS) for Genome-Wide SNP Identification in Plants.

Marker-assisted selection has played a pivotal role in developing several elite varieties in the past two decades. Molecular markers employed in plant breeding programs have recently shifted from microsatellites or simple sequence repeats (SSRs) to single nucleotide polymorphisms (SNPs) due to the ubiquity of SNP markers in the genome and the availability of various high-throughput SNP genotyping platforms. Rapid advances in sequencing technologies and the reduction in sequencing cost have facilitated SNP discovery in several plant species including non-model organisms with little or no genomic resources. Despite the lower cost of sequencing, genome complexity reduction approaches are still useful for SNP identification because many applications do not require every base of the genome to be sequenced. Genotyping-by-sequencing (GBS) is a quick and affordable reduced representation method that can simultaneously identify and genotype a large number of SNPs that has been successfully applied to a wide range of plant species. This chapter describes a robust two-enzyme GBS method for SNP discovery and genotyping that has been verified in non-model plant species.

Development of a High-Density Genetic Map for Muscadine Grape Using a Mapping Population from Selfing of the Perfect-Flowered Vine 'Dixie'.

Intraspecific diversity of the immune grape Muscadinia rotundifolia Michaux. can serve as a rich source of valuable resistance loci to the most widespread pathogens and pests of grapevine. While only one Run1/Rpg1 resistance locus has been introgressed from M. rotundifolia to the Vitis vinifera gene pool, a number of other genes conferring resistance to powdery mildew and downy mildew have been identified in various Muscadinia cultivars. A larger introduction of Muscadinia varieties to the European continent would greatly facilitate experiments of interspecific crosses as well as stimulate biotechnological efforts to overcome the main barrier to F1 fertility caused by the differences in chromosome number. For the successful introduction of Muscadinia into the new European environment, it is necessary to overcome the difficulties associated with the physiological characteristics of the species, such as insufficient cold tolerance and very late fruit ripening. To facilitate the further discovery of valuable loci in Muscadinia and their transfer to grapevine breeding programs, we constructed a high-density linkage map using an S1 mapping population obtained from the self-pollination of M. rotundifolia cv. Dixie maintained on the southern coast of Crimea. Using ddRADseq, 3730 SNPs were ordered across 20 linkage groups spanning 2753.6 cM of the total map length. No segregation in resistance to diseases and pests was observed among the 'Dixie' S1 population, suggesting the presence of homozygous non-segregating resistant loci in the genetic background of 'Dixie'. Markers with high segregation distortion showed a bias towards chromosomal intervals on linkage groups 10 and 20, where loci affecting the survival of 'Dixie' S1 progeny may be localized. QTLs with significant additive and dominance effects were discovered on LG14 and LG18, affecting the morphological traits associated with the vigor of growth and adaptability of young Muscadinia vines in the conditions of Crimea.

Genome-wide SNPs reveal novel patterns of spatial genetic structure in Aedes albopictus (Diptera Culicidae) population in China.

Introduction: Since the second half of the 20th century, Aedes albopictus, a vector for more than 20 arboviruses, has spread worldwide. Aedes albopictus is the main vector of infectious diseases transmitted by Aedes mosquitoes in China, and it has caused concerns regarding public health. A comprehensive understanding of the spatial genetic structure of this vector species at a genomic level is essential for effective vector control and the prevention of vector-borne diseases. Methods: During 2016-2018, adult female Ae. albopictus mosquitoes were collected from eight different geographical locations across China. Restriction site-associated DNA sequencing (RAD-seq) was used for high-throughput identification of single nucleotide polymorphisms (SNPs) and genotyping of the Ae. albopictus population. The spatial genetic structure was analyzed and compared to those exhibited by mitochondrial cytochrome c oxidase subunit 1 (cox1) and microsatellites in the Ae. albopictus population. Results: A total of 9,103 genome-wide SNP loci in 101 specimens and 32 haplotypes of cox1 in 231 specimens were identified in the samples from eight locations in China. Principal component analysis revealed that samples from Lingshui and Zhanjiang were more genetically different than those from the other locations. The SNPs provided a better resolution and stronger signals for novel spatial population genetic structures than those from the cox1 data and a set of previously genotyped microsatellites. The fixation indexes from the SNP dataset showed shallow but significant genetic differentiation in the population. The Mantel test indicated a positive correlation between genetic distance and geographical distance. However, the asymmetric gene flow was detected among the populations, and it was higher from south to north and west to east than in the opposite directions. Conclusions: The genome-wide SNPs revealed seven gene pools and fine spatial genetic structure of the Ae. albopictus population in China. The RAD-seq approach has great potential to increase our understanding of the spatial dynamics of population spread and establishment, which will help us to design new strategies for controlling vectors and mosquito-borne diseases.

Genome-Wide Population Structure and Selection Signatures of Yunling Goat Based on RAD-seq.

Animal diseases impose a huge burden on the countries where diseases are endemic. Conventional control strategies of vaccines and veterinary drugs are to control diseases from a pharmaceutical perspective. Another alternative approach is using pre-existing genetic disease resistance or tolerance. We know that the Yunling goat is an excellent local breed from Yunnan, southwestern China, which has characteristics of strong disease resistance and remarkable adaptability. However, genetic information about the selection signatures of Yunling goats is limited. We reasoned that the genes underlying the observed difference in disease resistance might be identified by investigating selection signatures between two different goat breeds. Herein, we selected the Nubian goat as the reference group to perform the population structure and selection signature analysis by using RAD-seq technology. The results showed that two goat breeds were divided into two clusters, but there also existed gene flow. We used Fst (F-statistics) and π (pi/θπ) methods to carry out selection signature analysis. Eight selected regions and 91 candidate genes were identified, in which some genes such as DOK2, TIMM17A, MAVS, and DOCK8 related to disease and immunity and some genes such as SPEFI, CDC25B, and MIR103 were associated with reproduction. Four GO (Gene Ontology) terms (GO:0010591, GO:001601, GO:0038023, and GO:0017166) were associated with cell migration, signal transduction, and immune responses. The KEGG (Kyoto Encyclopedia of Genes and Genomes) signaling pathways were mainly associated with immune responses, inflammatory responses, and stress reactions. This study preliminarily revealed the genetic basis of strong disease resistance and adaptability of Yunling goats. It provides a theoretical basis for the subsequent genetic breeding of disease resistance of goats.

A snapshot of progenitor-derivative speciation in Iberodes (Boraginaceae).

Traditional classification of speciation modes has focused on physical barriers to gene flow. Allopatric speciation with complete reproductive isolation is viewed as the most common mechanism of speciation. Parapatry and sympatry, by contrast, entail speciation in the face of ongoing gene flow, making them more difficult to detect. The genus Iberodes (Boraginaceae, NW Europe) comprises five species with contrasting morphological traits, habitats and species distributions. Based on the predominance of narrow and geographically distant endemic species, we hypothesized that geographical barriers were responsible for most speciation events in Iberodes. We undertook an integrative study including: (i) phylogenomics through restriction-site-associated DNA sequencing (RAD-seq), (ii) genetic structure analyses, (iii) demographic modelling, (iv) morphometrics, and (v) climatic niche modelling and niche overlap analysis. The results revealed a history of recurrent progenitor-derivative speciation manifested by a paraphyletic pattern of nested species differentiation. Budding speciation mediated by ecological differentiation is suggested for the coastal lineage, deriving from the inland widespread Iberodes linifolia during the Late Pliocene. Meanwhile, geographical isolation followed by niche shifts are suggested for the more recent differentiation of the coastland taxa. Our work provides a model for distinguishing speciation via ecological differentiation of peripheral, narrowly endemic I. kuzinskyanae and I. littoralis from a widespread extant ancestor, I. linifolia. Ultimately, our results illustrate a case of Pliocene speciation in the probable absence of geographical barriers and get away from the traditional cladistic perspective of speciation as producing two species from an extinct ancestor, thus reminding us that phylogenetic trees tell only part of the story.

Rescued back from extinction in the wild: past, present and future of the genetics of the Arabian oryx in Oman.

The Arabian oryx was the first species to be rescued from extinction in the wild by the concerted efforts of captive programmes in zoos and private collections around the world. Reintroduction efforts have used two main sources: the 'World Herd', established at the Phoenix Zoo, and private collections in Saudi Arabia. The breeding programme at the Al-Wusta Wildlife Reserve (WWR) in Oman has played a central role in the rescue of the oryx. Individuals from the 'World Herd' and the United Arab Emirates have been the main source for the WWR programme. However, no breeding strategies accounting for genetic diversity have been implemented. To address this, we investigated the diversity of the WWR population and historical samples using mitochondrial DNA (mtDNA) and single nucleotide polymorphisms (SNPs). We found individuals at WWR contain 58% of the total mtDNA diversity observed globally. Inference of ancestry and spatial patterns of SNP variation shows the presence of three ancestral sources and three different groups of individuals. Similar levels of diversity and low inbreeding were observed between groups. We identified individuals and groups that could most effectively contribute to maximizing genetic diversity. Our results will be valuable to guide breeding and reintroduction programmes at WWR.

Little hope for the polyploid endemic Pyrenean Larkspur (Delphinium montanum): Evidences from population genomics and Ecological Niche Modeling.

Species endemic to restricted geographical ranges represent a particular conservation issue, be it for their heritage interest. In a context of global change, this is particularly the case for plants which belong to high-mountain ecosystems and, because of their ecological requirements, are doomed to survive or disappear on their "sky islands". The Pyrenean Larkspur (Delphinium montanum, Ranunculaceae) is endemic to the Eastern part of the Pyrenees (France and Spain). It is now only observable at a dozen of localities and some populations show signs of decline, such as a recurrent lack of flowering. Implementing population genomics approach (e.g., RAD-seq like) is particularly useful to understand genomic patterns of diversity and differentiation in order to provide recommendations in term of conservation. However, it remains challenging for species such as D. montanum that are autotetraploid with a large genome size (1C-value >10 pg) as most methods currently available were developed for diploid species. A Bayesian framework able to call genotypes with uncertainty allowed us to assess genetic diversity and population structure in this system. Our results show evidence for inbreeding (mean G IS = 0.361) within all the populations and substantial population structure (mean G ST = 0.403) at the metapopulation level. In addition to a lack of connectivity between populations, spatial projections of Ecological Niche Modeling (ENM) analyses under different climatic scenarios predict a dramatic decrease of suitable habitat for D. montanum in the future. Based on these results, we discuss the relevance and feasibility of different conservation measures.

Population genomics of free-ranging Great Plains white-tailed and mule deer reflects a long history of interspecific hybridization.

Hybridization is a natural process at species-range boundaries that may variably promote the speciation process or break down species barriers but minimally will influence management outcomes of distinct populations. White-tailed deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) have broad and overlapping distributions in North America and a recognized capacity for interspecific hybridization. In response to contemporary environmental change to any of one or multiple still-unknown factors, mule deer range is contracting westward accompanied by a westward expansion of white-tailed deer, leading to increasing interactions, opportunities for gene flow, and associated conservation implications. To quantify genetic diversity, phylogenomic structure, and dynamics of hybridization in sympatric populations of white-tailed and mule deer, we used mitochondrial cytochrome b data coupled with SNP loci discovered with double-digest restriction site-associated DNA sequencing. We recovered 25,018 SNPs across 92 deer samples from both species, collected from two regions of western Kansas. Eight individuals with unambiguous external morphology representing both species were of hybrid origin (8.7%), and represented the product of multi-generational backcrossing. Mitochondrial data showed both ancient and recent directional discordance with morphological species assignments, reflecting a legacy of mule deer males mating with white-tailed deer females. Mule deer had lower genetic diversity than white-tailed deer, and both mitochondrial and nuclear data suggest contemporary mule deer effective population decline. Landscape genetic analyses show relative isolation between the two study regions for white-tailed deer, but greater connectivity among mule deer, with predominant movement from north to south. Collectively, our results suggest a long history of gene flow between these species in the Great Plains and hint at evolutionary processes that purge incompatible functional genomic elements as a result of hybridization. Surviving hybrids evidently may be reproductive, but with unknown consequences for the future integrity of these species, population trajectories, or relative susceptibility to emerging pathogens.

Application of Restriction Site-Associated DNA Sequencing (RAD-Seq) for Copy Number Variation and Triploidy Detection in Human.

At present, low-pass whole-genome sequencing (WGS) is frequently used in clinical research and in the screening of copy number variations (CNVs). However, there are still some challenges in the detection of triploids. Restriction site-associated DNA sequencing (RAD-Seq) technology is a reduced-representation genome sequencing technology developed based on next-generation sequencing. Here, we verified whether RAD-Seq could be employed to detect CNVs and triploids. In this study, genomic DNA of 11 samples was extracted employing a routine method and used to build libraries. Five cell lines of known karyotypes and 6 triploid abortion tissue samples were included for RAD-Seq testing. The triploid samples were confirmed by STR analysis and also tested by low-pass WGS. The accuracy and efficiency of detecting CNVs and triploids by RAD-Seq were then assessed, compared with low-pass WGS. In our results, RAD-Seq detected 11 out of 11 (100%) chromosomal abnormalities, including 4 deletions and 1 aneuploidy in the purchased cell lines and all triploid samples. By contrast, these triploids were missed by low-pass WGS. Furthermore, RAD-Seq showed a higher resolution and more accurate allele frequency in the detection of triploids than low-pass WGS. Our study shows that, compared with low-pass WGS, RAD-Seq has relatively higher accuracy in CNV detection at a similar cost and is capable of identifying triploids. Therefore, the application of this technique in medical genetics has a significant potential value.

Genome-Wide SNPs Provide Insights on the Cryptic Genetic Structure and Signatures of Climate Adaption in Amorphophallus albus Germplasms.

Domesticated species represent unique systems in which the evolutionary genomic consequences of intensive selective breeding and adaptation can be thoroughly investigated. Amorphophallus albus occurs naturally and is in cultivation throughout the downstream region of the Jinshajiang River in Southwest China. This species is characterised by high konjac glucomannan content, and has been cultivated in China for nearly 2,000 years. To study genetic differentiation and local adaption of A. albus, we sampled 13 distinct local cultivated populations of this species. Restriction site-associated DNA sequencing was conducted with 87 samples, resulting in 24,225 SNPs. The population structure analyses suggest two main genetic groups: one in the relatively upstream region, and one downstream. We found evidence of additional sub-structure within the upstream group, demonstrating the statistical power of genomic SNPs in discovering subtle genetic structure. The environmental and geographic factors were all identified as significant in shaping the genetic differentiation of this species. Notably, the proportion of environmental factors was larger than geographic factors in influencing the population genetic patterns of A. albus. We also discovered loci that were associated with local adaptation. These findings will help us understand the genetic differentiation of this newly domesticated species, thereby informing future breeding programs of A. albus.

Identification of QTL related to anther color and hull color by RAD sequencing in a RIL population of Setaria italica.

BACKGROUND: Foxtail millet (Setaria italica) is one of the oldest domesticated crops and has been considered as an ideal model plant for C4 grasses. It has abundant type of anther and hull colors which is not only a most intuitive morphological marker for color selection in seed production, but also has very important biological significance for the study of molecular mechanism of regulating the synthesis and metabolism of flavonoids and lignin. However, only a few genetic studies have been reported for anther color and hull color in foxtail millet. RESULTS: Quantitative trait loci (QTL) analysis for anther color and hull color was conducted using 400 F6 and F7 recombinant inbreed lines (RILs) derived from a cross between parents Yugu18 and Jigu19. Using restriction-site associated DNA sequencing, 43,001 single-nucleotide polymorphisms (SNPs) and 3,022 indels were identified between both the parents and the RILs. A total of 1,304 bin markers developed from the SNPs and indels were used to construct a genetic map that spanned 2196 cM of the foxtail millet genome with an average of 1.68 cM/bin. Combined with this genetic map and the phenotypic data observed in two locations for two years, two QTL located on chromosome 6 (Chr6) in a 1.215-Mb interval (33,627,819-34,877,940 bp) for anther color (yellow - white) and three QTL located on Chr1 in a 6.23-Mb interval (1-6,229,734 bp) for hull color (gold-reddish brown) were detected. To narrow the QTL regions identified from the genetic map and QTL analysis, we developed a new method named "inconsistent rate analysis" and efficiently narrowed the QTL regions of anther color into a 60-kb interval (34.13-34.19 Mb) in Chr6, and narrowed the QTL regions of hull color into 70-kb (5.43-5.50 Mb) and 30-kb (5.69-5.72 Mb) intervals in Chr1. Two genes (Seita.6G228600.v2.2 and Seita.6G228700.v2.2) and a cinnamyl alcohol dehydrogenase (CAD) gene (Seita.1G057300.v2.2) with amino acid changes between the parents detected by whole-genome resequencing were identified as candidate genes for anther and hull color, respectively. CONCLUSIONS: This work presents the related QTL and candidate genes of anther and hull color in foxtail millet and developed a new method named inconsistent rate analysis to detect the chromosome fragments linked with the quality trait in RILs. This is the first study of the QTL related to hull color in foxtail millet and clarifying that the CAD gene (Seita.1G057300.v2.2) is the key gene responsible for this trait. It lays the foundation for further cloning of the functional genes and provides a powerful tool to detect the chromosome fragments linked with quality traits in RILs.

New perspective on the genetic structure and habitat adaptation of Pampus minor off the coast of China based on RAD-seq.

Understanding the genetic diversity and habitat adaptation patterns of different geographic populations of marine organisms is essential for the rational maintenance, development, and utilization of their resources. Pampus minor Liu and Li 1998 is an economically valuable marine fish species. To determine the population genetic structure and elucidate the genetic mechanism underlying the habitat adaptation of this species, restriction site-associated DNA sequencing (RAD-seq) was used to scan the whole genomes of three P. minor putative populations along the coast of China for single-nucleotide polymorphism (SNPs) and outlier SNPs. Our population genetic structure analysis based on 2388 SNPs and 731 outlier SNPs throughout the genome revealed no significant genetic differentiation among the three populations. Results suggested that the life-cycle characteristics of P. minor, its relatively large population sizes, and ocean current transport might have shaped its current genetic pattern. The annotation information of both assembled sequences (which included outlier SNPs) and candidate loci associated with adaptations indicated that genes involved in many processes, including ion exchange, osmotic pressure regulation, metabolism, and the immune response, have been very important in the adaptations of P. minor to its heterogeneous habitats. Against the background of increased human activities, increased fishing intensity, and destruction of marine habitats, the results of this study provide basic genetic information for the accurate division of protection units of P. minor.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population.

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.