Hiding in plain sight: Genome-wide recombination and a dynamic accessory genome drive diversity in Fusarium oxysporum f.sp. ciceris.

Understanding the origins of variation in agricultural pathogens is of fundamental interest and practical importance, especially for diseases that threaten food security. Fusarium oxysporum is among the most important of soil-borne pathogens, with a global distribution and an extensive host range. The pathogen is considered to be asexual, with horizontal transfer of chromosomes providing an analog of assortment by meiotic recombination. Here, we challenge those assumptions based on the results of population genomic analyses, describing the pathogen's diversity and inferring its origins and functional consequences in the context of a single, long-standing agricultural system. We identify simultaneously low nucleotide distance among strains, and unexpectedly high levels of genetic and genomic variability. We determine that these features arise from a combination of genome-scale recombination, best explained by widespread sexual reproduction, and presence-absence variation consistent with chromosomal rearrangement. Pangenome analyses document an accessory genome more than twice the size of the core genome, with contrasting evolutionary dynamics. The core genome is stable, with low diversity and high genetic differentiation across geographic space, while the accessory genome is paradoxically more diverse and unstable but with lower genetic differentiation and hallmarks of contemporary gene flow at local scales. We suggest a model in which episodic sexual reproduction generates haplotypes that are selected and then maintained through clone-like dynamics, followed by contemporary genomic rearrangements that reassort the accessory genome among sympatric strains. Taken together, these processes contribute unique genome content, including reassortment of virulence determinants that may explain observed variation in pathogenic potential.

Ultralow-input single-tube linked-read library method enables short-read second-generation sequencing systems to routinely generate highly accurate and economical long-range sequencing information.

Long-range sequencing information is required for haplotype phasing, de novo assembly, and structural variation detection. Current long-read sequencing technologies can provide valuable long-range information but at a high cost with low accuracy and high DNA input requirements. We have developed a single-tube Transposase Enzyme Linked Long-read Sequencing (TELL-seq) technology, which enables a low-cost, high-accuracy, and high-throughput short-read second-generation sequencer to generate over 100 kb of long-range sequencing information with as little as 0.1 ng input material. In a PCR tube, millions of clonally barcoded beads are used to uniquely barcode long DNA molecules in an open bulk reaction without dilution and compartmentation. The barcoded linked-reads are used to successfully assemble genomes ranging from microbes to human. These linked-reads also generate megabase-long phased blocks and provide a cost-effective tool for detecting structural variants in a genome, which are important to identify compound heterozygosity in recessive Mendelian diseases and discover genetic drivers and diagnostic biomarkers in cancers.

Global-level population genomics reveals differential effects of geography and phylogeny on horizontal gene transfer in soil bacteria.

Although microorganisms are known to dominate Earth's biospheres and drive biogeochemical cycling, little is known about the geographic distributions of microbial populations or the environmental factors that pattern those distributions. We used a global-level hierarchical sampling scheme to comprehensively characterize the evolutionary relationships and distributional limitations of the nitrogen-fixing bacterial symbionts of the crop chickpea, generating 1,027 draft whole-genome sequences at the level of bacterial populations, including 14 high-quality PacBio genomes from a phylogenetically representative subset. We find that diverse Mesorhizobium taxa perform symbiosis with chickpea and have largely overlapping global distributions. However, sampled locations cluster based on the phylogenetic diversity of Mesorhizobium populations, and diversity clusters correspond to edaphic and environmental factors, primarily soil type and latitude. Despite long-standing evolutionary divergence and geographic isolation, the diverse taxa observed to nodulate chickpea share a set of integrative conjugative elements (ICEs) that encode the major functions of the symbiosis. This symbiosis ICE takes 2 forms in the bacterial chromosome-tripartite and monopartite-with tripartite ICEs confined to a broadly distributed superspecies clade. The pairwise evolutionary relatedness of these elements is controlled as much by geographic distance as by the evolutionary relatedness of the background genome. In contrast, diversity in the broader gene content of Mesorhizobium genomes follows a tight linear relationship with core genome phylogenetic distance, with little detectable effect of geography. These results illustrate how geography and demography can operate differentially on the evolution of bacterial genomes and offer useful insights for the development of improved technologies for sustainable agriculture.

Genotyping and lipid profiling of 601 cultivated sunflower lines reveals novel genetic determinants of oil fatty acid content.

BACKGROUND: Sunflower is an important oilseed crop domesticated in North America approximately 4000 years ago. During the last century, oil content in sunflower was under strong selection. Further improvement of oil properties achieved by modulating its fatty acid composition is one of the main directions in modern oilseed crop breeding. RESULTS: We searched for the genetic basis of fatty acid content variation by genotyping 601 inbred sunflower lines and assessing their lipid and fatty acid composition. Our genome-wide association analysis based on the genotypes for 15,483 SNPs and the concentrations of 23 fatty acids, including minor fatty acids, revealed significant genetic associations for eleven of them. Identified genomic regions included the loci involved in rare fatty acids variation on chromosomes 3 and 14, explaining up to 34.5% of the total variation of docosanoic acid (22:0) in sunflower oil. CONCLUSIONS: This is the first large scale implementation of high-throughput lipidomic profiling to sunflower germplasm characterization. This study contributes to the genetic characterization of Russian sunflower collections, which made a substantial contribution to the development of sunflower as the oilseed crop worldwide, and provides new insights into the genetic control of oil composition that can be implemented in future studies.

Predictive genotype-phenotype relations using genetic diversity in African yam bean (Sphenostylis stenocarpa (Hochst. ex. A. Rich) Harms).

BACKGROUND: African Yam Bean (AYB) is an understudied and underutilized tuberous legume of tropical West and Central African origin. In these geographical regions, both seeds and tubers of AYB are important components of people's diets and a potential target as a nutritional security crop. The understanding of the genetic diversity among AYB accessions is thus an important component for both conservation and potential breeding programs. RESULTS: In this study, 93 AYB accessions were obtained from the International Institute of Tropical Agriculture (IITA) genebank and genotyped using 3722 SNP markers based on Restriction site-Associated DNA sequencing (RAD-Seq). Genetic data was analysed using multiple clustering methods for better understanding the distribution of genetic diversity across the population. Substantial genetic variability was observed in the present set of AYB accessions and different methodologies demonstrated that these accessions are divided into three to four main groups. The accessions were also analysed for important agronomic traits and successfully associated with their genetic clusters where great majority of accessions shared a similar phenotype. CONCLUSIONS: To our knowledge, this is the first study on predicting genotypic-phenotypic diversity relationship analysis in AYB. From a breeding perspective, we were able to identify specific diverse groups with precise phenotype such as seed or both seed and tuber yield purpose accessions. These results provide novel and important insights to support the utilization of this germplasm in AYB breeding programs.

Genomic Analysis of Vavilov's Historic Chickpea Landraces Reveals Footprints of Environmental and Human Selection.

A defining challenge of the 21st century is meeting the nutritional demands of the growing human population, under a scenario of limited land and water resources and under the specter of climate change. The Vavilov seed bank contains numerous landraces collected nearly a hundred years ago, and thus may contain 'genetic gems' with the potential to enhance modern breeding efforts. Here, we analyze 407 landraces, sampled from major historic centers of chickpea cultivation and secondary diversification. Genome-Wide Association Studies (GWAS) conducted on both phenotypic traits and bioclimatic variables at landraces sampling sites as extended phenotypes resulted in 84 GWAS hits associated to various regions. The novel haploblock-based test identified haploblocks enriched for single nucleotide polymorphisms (SNPs) associated with phenotypes and bioclimatic variables. Subsequent bi-clustering of traits sharing enriched haploblocks underscored both non-random distribution of SNPs among several haploblocks and their association with multiple traits. We hypothesize that these clusters of pleiotropic SNPs represent co-adapted genetic complexes to a range of environmental conditions that chickpea experienced during domestication and subsequent geographic radiation. Linking genetic variation to phenotypic data and a wealth of historic information preserved in historic seed banks are the keys for genome-based and environment-informed breeding intensification.

Identifying similar transcripts in a related organism from de Bruijn graphs of RNA-Seq data, with applications to the study of salt and waterlogging tolerance in Melilotus.

BACKGROUND: A popular strategy to study alternative splicing in non-model organisms starts from sequencing the entire transcriptome, then assembling the reads by using de novo transcriptome assembly algorithms to obtain predicted transcripts. A similarity search algorithm is then applied to a related organism to infer possible function of these predicted transcripts. While some of these predictions may be inaccurate and transcripts with low coverage are often missed, we observe that it is possible to obtain a more complete set of transcripts to facilitate possible functional assignments by starting the search from the intermediate de Bruijn graph that contains all branching possibilities. RESULTS: We develop an algorithm to extract similar transcripts in a related organism by starting the search from the de Bruijn graph that represents the transcriptome instead of from predicted transcripts. We show that our algorithm is able to recover more similar transcripts than existing algorithms, with large improvements in obtaining longer transcripts and a finer resolution of isoforms. We apply our algorithm to study salt and waterlogging tolerance in two Melilotus species by constructing new RNA-Seq libraries. CONCLUSIONS: We have developed an algorithm to identify paths in the de Bruijn graph that correspond to similar transcripts in a related organism directly. Our strategy bypasses the transcript prediction step in RNA-Seq data and makes use of support from evolutionary information.

Functional Dissection of the Chickpea (Cicer arietinum L.) Stay-Green Phenotype Associated with Molecular Variation at an Ortholog of Mendel's I Gene for Cotyledon Color: Implications for Crop Production and Carotenoid Biofortification.

"Stay-green" crop phenotypes have been shown to impact drought tolerance and nutritional content of several crops. We aimed to genetically describe and functionally dissect the particular stay-green phenomenon found in chickpeas with a green cotyledon color of mature dry seed and investigate its potential use for improvement of chickpea environmental adaptations and nutritional value. We examined 40 stay-green accessions and a set of 29 BC2F4-5 stay-green introgression lines using a stay-green donor parent ICC 16340 and two Indian elite cultivars (KAK2, JGK1) as recurrent parents. Genetic studies of segregating populations indicated that the green cotyledon trait is controlled by a single recessive gene that is invariantly associated with the delayed degreening (extended chlorophyll retention). We found that the chickpea ortholog of Mendel's I locus of garden pea, encoding a SGR protein as very likely to underlie the persistently green cotyledon color phenotype of chickpea. Further sequence characterization of this chickpea ortholog CaStGR1 (CaStGR1, for carietinum stay-green gene 1) revealed the presence of five different molecular variants (alleles), each of which is likely a loss-of-function of the chickpea protein (CaStGR1) involved in chlorophyll catabolism. We tested the wild type and green cotyledon lines for components of adaptations to dry environments and traits linked to agronomic performance in different experimental systems and different levels of water availability. We found that the plant processes linked to disrupted CaStGR1 gene did not functionality affect transpiration efficiency or water usage. Photosynthetic pigments in grains, including provitaminogenic carotenoids important for human nutrition, were 2-3-fold higher in the stay-green type. Agronomic performance did not appear to be correlated with the presence/absence of the stay-green allele. We conclude that allelic variation in chickpea CaStGR1 does not compromise traits linked to environmental adaptation and agronomic performance, and is a promising genetic technology for biofortification of provitaminogenic carotenoids in chickpea.

Genomic exaptation enables Lasius niger adaptation to urban environments.

BACKGROUND: The world is rapidly urbanizing, and only a subset of species are able to succeed in stressful city environments. Efficient genome-enabled stress response appears to be a likely prerequisite for urban adaptation. Despite the important role ants play in the ecosytem, only the genomes of ~13 have been sequenced so far. Here, we present the draft genome assembly of the black garden ant Lasius niger - the most successful urban inhabitant of all ants - and we compare it with the genomes of other ant species, including the closely related Camponotus floridanus. RESULTS: Sequences from 272 M Illumina reads were assembled into 41,406 contigs with total length of 245 MB, and N50 of 16,382 bp, similar to other ant genome assemblies enabling comparative genomic analysis. Remarkably, the predicted proteome of L. niger is significantly enriched relative to other ant genomes in terms of abundance of domains involved in nucleic acid binding, DNA repair, and nucleotidyl transferase activity, reflecting transposable element proliferation and a likely genomic response. With respect to environmental stress, we note a proliferation of various detoxification genes, including glutatione-S-transferases and those in the cytochrome P450 families. Notably, the CYP9 family is highly expanded with 19 complete and 21 nearly complete members - over twice as many compared to other ants. This family exhibits the signatures of strong directional selection, with eleven positively selected positions in ligand-binding pockets of enzymes. Gene family contraction was detected for several components of the olfactory system, accompanied by instances of both directional selection and relaxation. CONCLUSIONS: Our results suggest that the success of L. niger in urbanized areas may be the result of fortuitous coincidence of several factors, including the expansion of the CYP9 cytochrome family due to coevolution with parasitic fungi, the diversification of DNA repair systems as an answer to proliferation of retroelements, and the reduction of olfactory system and behavioral preadaptations from non-territorial subdominant life strategies found in natural environments. Diversification of cytochromes and DNA repair systems along with reduced odorant communication are the basis of L. niger pollutant resistance and polyphagy, while non-territorial and mobilization strategies allows more efficient exploitation of large but patchy food sources.

Development of F1 hybrid population and the high-density linkage map for European aspen (Populus tremula L.) using RADseq technology.

BACKGROUND: Restriction-site associated DNA sequencing (RADseq) technology was recently employed to identify a large number of single nucleotide polymorphisms (SNP) for linkage mapping of a North American and Eastern Asian Populus species. However, there is also the need for high-density genetic linkage maps for the European aspen (P. tremula) as a tool for further mapping of quantitative trait loci (QTLs) and marker-assisted selection of the Populus species native to Europe. RESULTS: We established a hybrid F1 population from the cross of two aspen parental genotypes diverged in their phenological and morphological traits. We performed RADseq of 122 F1 progenies and two parents yielding 15,732 high-quality SNPs that were successfully identified using the reference genome of P. trichocarpa. 2055 SNPs were employed for the construction of maternal and paternal linkage maps. The maternal linkage map was assembled with 1000 SNPs, containing 19 linkage groups and spanning 3054.9 cM of the genome, with an average distance of 3.05 cM between adjacent markers. The paternal map consisted of 1055 SNPs and the same number of linkage groups with a total length of 3090.56 cM and average interval distance of 2.93 cM. The linkage maps were employed for QTL mapping of one-year-old seedlings height variation. The most significant QTL (LOD = 5.73) was localized to LG5 (96.94 cM) of the male linkage map, explaining 18% of the phenotypic variation. CONCLUSIONS: The set of 15,732 SNPs polymorphic in aspen and high-density genetic linkage maps constructed for the P. tremula intra-specific cross will provide a valuable source for QTL mapping and identification of candidate genes facilitating marker-assisted selection in European aspen.

Whole exome sequencing of wild-derived inbred strains of mice improves power to link phenotype and genotype.

The house mouse is a powerful model to dissect the genetic basis of phenotypic variation, and serves as a model to study human diseases. Despite a wealth of discoveries, most classical laboratory strains have captured only a small fraction of genetic variation known to segregate in their wild progenitors, and existing strains are often related to each other in complex ways. Inbred strains of mice independently derived from natural populations have the potential to increase power in genetic studies with the addition of novel genetic variation. Here, we perform exome-enrichment and high-throughput sequencing (~8× coverage) of 26 wild-derived strains known in the mouse research community as the "Montpellier strains." We identified 1.46 million SNPs in our dataset, approximately 19% of which have not been detected from other inbred strains. This novel genetic variation is expected to contribute to phenotypic variation, as they include 18,496 nonsynonymous variants and 262 early stop codons. Simulations demonstrate that the higher density of genetic variation in the Montpellier strains provides increased power for quantitative genetic studies. Inasmuch as the power to connect genotype to phenotype depends on genetic variation, it is important to incorporate these additional genetic strains into future research programs.

Races of the Celery Pathogen Fusarium oxysporum f. sp. apii Are Polyphyletic.

Fusarium oxysporum species complex (FOSC) isolates were obtained from celery with symptoms of Fusarium yellows between 1993 and 2013 primarily in California. Virulence tests and a two-gene dataset from 174 isolates indicated that virulent isolates collected before 2013 were a highly clonal population of F. oxysporum f. sp. apii race 2. In 2013, new highly virulent clonal isolates, designated race 4, were discovered in production fields in Camarillo, California. Long-read Illumina data were used to analyze 16 isolates: six race 2, one of each from races 1, 3, and 4, and seven genetically diverse FOSC that were isolated from symptomatic celery but are nonpathogenic on this host. Analyses of a 10-gene dataset comprising 38 kb indicated that F. oxysporum f. sp. apii is polyphyletic; race 2 is nested within clade 3, whereas the evolutionary origins of races 1, 3, and 4 are within clade 2. Based on 6,898 single nucleotide polymorphisms from the core FOSC genome, race 3 and the new highly virulent race 4 are highly similar with Nei's Da = 0.0019, suggesting that F. oxysporum f. sp. apii race 4 evolved from race 3. Next generation sequences were used to develop PCR primers that allow rapid diagnosis of races 2 and 4 in planta.

Sex-specific signaling in the blood-brain barrier is required for male courtship in Drosophila.

Soluble circulating proteins play an important role in the regulation of mating behavior in Drosophila melanogaster. However, how these factors signal through the blood-brain barrier (bbb) to interact with the sex-specific brain circuits that control courtship is unknown. Here we show that male identity of the blood-brain barrier is necessary and that male-specific factors in the bbb are physiologically required for normal male courtship behavior. Feminization of the bbb of adult males significantly reduces male courtship. We show that the bbb-specific G-protein coupled receptor moody and bbb-specific Go signaling in adult males are necessary for normal courtship. These data identify sex-specific factors and signaling processes in the bbb as important regulators of male mating behavior.

CYCP2;1 integrates genetic and nutritional information to promote meristem cell division in Arabidopsis.

In higher plants, cell cycle activation in the meristems at germination is essential for the initiation of post-embryonic development. We previously identified the signaling pathways of homeobox transcription factor STIMPY and metabolic sugars as two interacting branches of the regulatory network that is responsible for activating meristematic tissue proliferation in Arabidopsis. In this study, we found that CYCP2;1 is both a direct target of STIMPY transcriptional activation and an early responder to sugar signals. Genetic and molecular studies show that CYCP2;1 physically interacts with three of the five mitotic CDKs in Arabidopsis, and is required for the G2 to M transition during meristem activation. Taken together, our results suggest that CYCP2;1 acts as a permissive control of cell cycle progression during seedling establishment by directly linking genetic control and nutritional cues with the activity of the core cell cycle machinery.

Genetic variation of transgenerational plasticity of offspring germination in response to salinity stress and the seed transcriptome of Medicago truncatula.

BACKGROUND: Transgenerational plasticity provides phenotypic variation that contributes to adaptation. For plants, the timing of seed germination is critical for offspring survival in stressful environments, as germination timing can alter the environmental conditions a seedling experiences. Stored seed transcripts are important determinants of seed germination, but have not previously been linked with transgenerational plasticity of germination behavior. In this study we used RNAseq and growth chamber experiments of the model legume M. trucantula to test whether parental exposure to salinity stress influences the expression of stored seed transcripts and early offspring traits and test for genetic variation. RESULTS: We detected genotype-dependent parental environmental effects (transgenerational plasticity) on the expression levels of stored seed transcripts, seed size, and germination behavior of four M. truncatula genotypes. More than 50% of the transcripts detected in the mature, ungerminated seed transcriptome were annotated as regulating seed germination, some of which are involved in abiotic stress response and post-embryonic development. Some genotypes showed increased seed size in response to parental exposure to salinity stress, but no parental environmental influence on germination timing. In contrast, other genotypes showed no seed size differences across contrasting parental conditions but displayed transgenerational plasticity for germimation timing, with significantly delayed germination in saline conditions when parental plants were exposed to salinity. In genotypes that show significant transgenerational plastic germination response, we found significant coexpression networks derived from salt responsive transcripts involved in post-transcriptional regulation of the germination pathway. Consistent with the delayed germination response to saline conditions in these genotypes, we found genes associated with dormancy and up-regulation of abscisic acid (ABA). CONCLUSIONS: Our results demonstrate genetic variation in transgenerational plasticity within M. truncatula and show that parental exposure to salinity stress influences the expression of stored seed transcripts, seed weight, and germination behavior. Furthermore, we show that the parental environment influences gene expression to modulate biological pathways that are likely responsible for offspring germination responses to salinity stress.

Evolutionary genomics of Culex pipiens: global and local adaptations associated with climate, life-history traits and anthropogenic factors.

We present the first genome-wide study of recent evolution in Culex pipiens species complex focusing on the genomic extent, functional targets and likely causes of global and local adaptations. We resequenced pooled samples of six populations of C. pipiens and two populations of the outgroup Culex torrentium. We used principal component analysis to systematically study differential natural selection across populations and developed a phylogenetic scanning method to analyse admixture without haplotype data. We found evidence for the prominent role of geographical distribution in shaping population structure and specifying patterns of genomic selection. Multiple adaptive events, involving genes implicated with autogeny, diapause and insecticide resistance were limited to specific populations. We estimate that about 5-20% of the genes (including several histone genes) and almost half of the annotated pathways were undergoing selective sweeps in each population. The high occurrence of sweeps in non-genic regions and in chromatin remodelling genes indicated the adaptive importance of gene expression changes. We hypothesize that global adaptive processes in the C. pipiens complex are potentially associated with South to North range expansion, requiring adjustments in chromatin conformation. Strong local signature of adaptation and emergence of hybrid bridge vectors necessitate genomic assessment of populations before specifying control agents.

The ecological genomic basis of salinity adaptation in Tunisian Medicago truncatula.

BACKGROUND: As our world becomes warmer, agriculture is increasingly impacted by rising soil salinity and understanding plant adaptation to salt stress can help enable effective crop breeding. Salt tolerance is a complex plant phenotype and we know little about the pathways utilized by naturally tolerant plants. Legumes are important species in agricultural and natural ecosystems, since they engage in symbiotic nitrogen-fixation, but are especially vulnerable to salinity stress. RESULTS: Our studies of the model legume Medicago truncatula in field and greenhouse settings demonstrate that Tunisian populations are locally adapted to saline soils at the metapopulation level and that saline origin genotypes are less impacted by salt than non-saline origin genotypes; these populations thus likely contain adaptively diverged alleles. Whole genome resequencing of 39 wild accessions reveals ongoing migration and candidate genomic regions that assort non-randomly with soil salinity. Consistent with natural selection acting at these sites, saline alleles are typically rare in the range-wide species' gene pool and are also typically derived relative to the sister species M. littoralis. Candidate regions for adaptation contain genes that regulate physiological acclimation to salt stress, such as abscisic acid and jasmonic acid signaling, including a novel salt-tolerance candidate orthologous to the uncharacterized gene AtCIPK21. Unexpectedly, these regions also contain biotic stress genes and flowering time pathway genes. We show that flowering time is differentiated between saline and non-saline populations and may allow salt stress escape. CONCLUSIONS: This work nominates multiple potential pathways of adaptation to naturally stressful environments in a model legume. These candidates point to the importance of both tolerance and avoidance in natural legume populations. We have uncovered several promising targets that could be used to breed for enhanced salt tolerance in crop legumes to enhance food security in an era of increasing soil salinization.

Genomic changes under rapid evolution: selection for parasitoid resistance.

In this study, we characterize changes in the genome during a swift evolutionary adaptation, by combining experimental selection with high-throughput sequencing. We imposed strong experimental selection on an ecologically relevant trait, parasitoid resistance in Drosophila melanogaster against Asobara tabida. Replicated selection lines rapidly evolved towards enhanced immunity. Larval survival after parasitization increased twofold after just five generations of selection. Whole-genome sequencing revealed that the fast and strong selection response in innate immunity produced multiple, highly localized genomic changes. We identified narrow genomic regions carrying a significant signature of selection, which were present across all chromosomes and covered in total less than 5% of the whole D. melanogaster genome. We identified segregating sites with highly significant changes in frequency between control and selection lines that fell within these narrow 'selected regions'. These segregating sites were associated with 42 genes that constitute possible targets of selection. A region on chromosome 2R was highly enriched in significant segregating sites and may be of major effect on parasitoid defence. The high genetic variability and small linkage blocks in our base population are likely responsible for allowing this complex trait to evolve without causing widespread erosive effects in the genome, even under such a fast and strong selective regime.

Targeted phasing of 2-200 kilobase DNA fragments with a short-read sequencer and a single-tube linked-read library method.

In the human genome, heterozygous sites refer to genomic positions with a different allele or nucleotide variant on the maternal and paternal chromosomes. Resolving these allelic differences by chromosomal copy, also known as phasing, is achievable on a short-read sequencer when using a library preparation method that captures long-range genomic information. TELL-Seq is a library preparation that captures long-range genomic information with the aid of molecular identifiers (barcodes). The same barcode is used to tag the reads derived from the same long DNA fragment within a range of up to 200 kilobases (kb), generating linked-reads. This strategy can be used to phase an entire genome. Here, we introduce a TELL-Seq protocol developed for targeted applications, enabling the phasing of enriched loci of varying sizes, purity levels, and heterozygosity. To validate this protocol, we phased 2-200 kb loci enriched with different methods: CRISPR/Cas9-mediated excision coupled with pulse-field electrophoresis for the longest fragments, CRISPR/Cas9-mediated protection from exonuclease digestion for mid-size fragments, and long PCR for the shortest fragments. All selected loci have known clinical relevance: BRCA1, BRCA2, MLH1, MSH2, MSH6, APC, PMS2, SCN5A-SCN10A, and PKI3CA. Collectively, the analyses show that TELL-Seq can accurately phase 2-200 kb targets using a short-read sequencer.