CNVmap: A Method and Software To Detect and Map Copy Number Variants from Segregation Data

Single nucleotide polymorphisms (SNPs) are used widely for detecting quantitative trait loci, or for searching for causal variants of diseases. Nevertheless, structural variations such as copy-number variants (CNVs) represent a large part of natural genetic diversity, and contribute significantly to trait variation. Numerous methods and softwares based on different technologies (amplicons, CGH, tiling, or SNP arrays, or sequencing) have already been developed to detect CNVs, but they bypass a wealth of information such as genotyping data from segregating populations, produced, e.g., for QTL mapping. Here, we propose an original method to both detect and genetically map CNVs using mapping panels. Specifically, we exploit the apparent heterozygous state of duplicated loci: peaks in appropriately defined genome-wide allelic profiles provide highly specific signatures that identify the nature and position of the CNVs. Our original method and software can detect and map automatically up to 33 different predefined types of CNVs based on segregation data only. We validate this approach on simulated and experimental biparental mapping panels in two maize populations and one wheat population. Most of the events found correspond to having just one extra copy in one of the parental lines, but the corresponding allelic value can be that of either parent. We also find cases with two or more additional copies, especially in wheat, where these copies locate to homeologues. More generally, our computational tool can be used to give additional value, at no cost, to many datasets produced over the past decade from genetic mapping panels.

CNVmap: a method and software to detect and map copy number variants from segregation data

Single nucleotide polymorphisms (SNPs) are widely used for detecting quantitative trait loci or for searching for causal variants of diseases. Nevertheless, structural variations such as copy-number variants (CNVs) represent a large part of natural genetic diversity and contribute significantly to trait variation. Over the past decade, numerous methods and softwares have been developed to detect CNVs. Such approaches are based on exploiting sequencing data or SNP arrays, but they bypass a wealth of information such as genotyping data from segregating populations, produced e.g. for QTL mapping. Here we propose an original method to both detect and genetically map CNVs using mapping panels. Specifically, we exploit the apparent heterozygous state of duplicated loci: peaks in appropriately defined genome-wide allelic profiles provide highly specific signatures that identify the nature and position of the CNVs. Our original method and software can detect and map automatically up to 33 different predefined types of CNVs based on segregation data only. We validate this approach on simulated and experimental bi-parental mapping panels in two maize and one wheat populations. Most of the events found correspond to having just one extra copy in one of the parental lines but the corresponding allelic value can be that of either parent. We also find cases with two or more additional copies, especially in wheat where these copies locate to homeologues. More generally, our computational tool can be used to give additional value, at no cost, to many datasets produced over the past decade from genetic mapping panels.

Landscape of gene transposition-duplication within the Brassicaceae family

ABSTRACTWe developed the CLfinder-OrthNet pipeline that detects co-linearity in gene arrangement among multiple closely related genomes; find ortholog groups; and encodes the evolutionary history of each ortholog group into a representative network (OrthNet). Using a search based on network topology, out of a total of 17,432 OrthNets in six Brassicaceae genomes, we identified 1,394 that included gene transposition-duplication (tr-d) events in one or more genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. They also had a higher chance of being expressed and derived from older tr-d events. We also found an enrichment, compared to random chance, of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3. The CLfinder-OrthNet pipeline provides a flexible and a modular toolkit to compare gene order, encode and visualize evolutionary paths among orthologs as networks, and identify all gene loci that share the same evolutionary history using network topology searches.Funding source: This work was supported by National Science Foundation (MCB 1616827) and the Next Generation BioGreen21 Program (PJ011379) of the Rural Development Administration, Republic of Korea.Online-only Supplementary materials includes supplementary text (S1-S10), methods (M1-M4), figures (S1-S7), and tables (S1-S3), in two PDF files, one for text and methods and the other for figures and tables. Additionally, Supplementary Dataset S1 is available at the Figshare repository (https://doi.org/10.6084/m9.figshare.5825937) and Dataset S2 and S3 as separate Excel files.

Ribosomal DNA locus variation and REMAP analysis of the diploid and triploid complexes of Lilium lancifolium

Lilium lancifolium Thunb. (2n = 2x = 24) is a cytologically conspicuous species with both diploids and triploids in nature. Cytological and molecular genetic analyses were carried out in both diploids and triploids that were collected from 55 geographical locations in Korea, Japan, and China. While the 5S rRNA gene loci were located at duplicated loci on the long arm of chromosome 2, the 45S rRNA gene loci were present in chromosomes 1, 2, 4, 6, 7, and 11. While the loci on chromosomes 1 and 7 were constant, the loci on chromosomes 2, 4, 6, 7, and 11 were variable in some plants so that the L. lancifolium accessions were grouped into 7 cytotypes in diploids and 12 cytotypes in triploids. REMAP marker analysis revealed that the diploids were classified into seven clusters, and the triploids were classified into a large cluster. Geographic, cytological, and genetic differentiations were not related in both the diploid and triploid accessions of L. lancifolium. Thus, current genetic variations occurred prior to the geographic differentiation in both diploids and triploids, and the 45S rDNA cytotype variations occurred after geographic differentiation in the current habitats of L. lancifolium.