Targeted sequencing of the human X chromosome exome

We used a RainDance Technologies (RDT) expanded content library to enrich the human X chromosome exome (2.5 Mb) from 26 male samples followed by Illumina sequencing. Our multiplex primer library covered 98.05% of the human X chromosome exome in a single tube with 11,845 different PCR amplicons. Illumina sequencing of 24 male samples showed coverage for 97% of the targeted sequences. Sequence from 2 HapMap samples confirmed missing data rates of 2–3% at sites successfully typed by the HapMap project, with an accuracy of at least ~ 99.5% as compared to reported HapMap genotypes. Our demonstration that a RDT expanded content library can efficiently enrich and enable the routine sequencing of the human X chromosome exome suggests a wide variety of potential research and clinical applications for this platform.     

Mobilization of Stowaway-like MITEs in newly formed allohexaploid wheat species

Transposable elements (TEs) dominate the genetic capacity of most eukaryotes, especially plants, where they can account for up to 90?% of the genome, such as in wheat. The relationship between TEs and their hosts and the role of TEs in organismal biology are poorly understood. In this study, we have applied next generation sequencing, together with a transposon display technique in order to test whether a Stowaway-like MITE, termed Minos, transposes following allopolyploidization events in wheat. We have generated a 454-pyrosequencing database of Minos-specific amplicons (transposon display products) from a newly formed wheat allohexaploid and its parental lines and retrieved hundreds of novel MITE insertions in the allohexaploid. Clear mobilization of Minos was also seen by site-specific PCR analysis and sequence validation. In addition, using real-time qPCR analysis we observed an insignificant change in the relative quantity of Minos from the expected value of merging the two parental genomes, indicating that, despite its activation, no significant burst in Minos copy number can be seen in the newly formed allohexaploid. Interestingly, we found that CCGG sites surrounding Minos underwent massive hypermethylation following the allohexaploidization process. Our data suggest that MITEs have maintained their capacity for activity throughout the evolution of wheat and might be epigenetically deregulated in the first generations following allopolyploidization.     

Natural variation for fertile triploid F1 hybrid formation in allohexaploid wheat speciation

The tempo, mode, and geography of allopolyploid speciation are influenced by natural variation in the ability of parental species to express postzygotic reproductive phenotypes that affect hybrid fertility. To shed light on the impact of such natural variations, we used allohexaploid Triticum aestivum wheats’ evolution as a model and analyzed the geographic and phylogenetic distributions of Aegilops tauschii (diploid progenitor) accessions involved in the expression of abnormality and fertility in triploid F₁ hybrids with Triticum turgidum (tetraploid progenitor). Artificial-cross experiments and chloroplast-DNA-based evolutionary analyses showed that hybrid-abnormality-causing accessions had limited geographic and phylogenetic distributions, indicative that postzygotic hybridization barriers are underdeveloped between these species. In contrast, accessions that are involved with fertile triploid F₁ hybrid formation have wide geographic and phylogenetic distributions, indicative of a deep evolutionary origin. Wide-spread hybrid-fertilizing accessions support the theory that T. aestivum speciation occurred at multiple sites within the species range of Ae. tauschii, in which existing conditions enabled natural hybridization with T. turgidum. Implications of our findings on how natural variation in the ability of Ae. tauschii to express those postzygotic reproductive phenotypes diversified and contributed to the speciation of T. aestivum are discussed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes

BackgroundBread wheat is an allopolyploid species with a large, highly repetitive genome. To investigate the impact of selection on variants distributed among homoeologous wheat genomes and to build a foundation for understanding genotype-phenotype relationships, we performed population-scale re-sequencing of a diverse panel of wheat lines.ResultsA sample of 62 diverse lines was re-sequenced using the whole exome capture and genotyping-by-sequencing approaches. We describe the allele frequency, functional significance, and chromosomal distribution of 1.57 million single nucleotide polymorphisms and 161,719 small indels. Our results suggest that duplicated homoeologous genes are under purifying selection. We find contrasting patterns of variation and inter-variant associations among wheat genomes; this, in addition to demographic factors, could be explained by differences in the effect of directional selection on duplicated homoeologs. Only a small fraction of the homoeologous regions harboring selected variants overlapped among the wheat genomes in any given wheat line. These selected regions are enriched for loci associated with agronomic traits detected in genome-wide association studies.ConclusionsEvidence suggests that directional selection in allopolyploids rarely acted on multiple parallel advantageous mutations across homoeologous regions, likely indicating that a fitness benefit could be obtained by a mutation at any one of the homoeologs. Additional advantageous variants in other homoelogs probably either contributed little benefit, or were unavailable in populations subjected to directional selection. We hypothesize that allopolyploidy may have increased the likelihood of beneficial allele recovery by broadening the set of possible selection targets.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0606-4) contains supplementary material, which is available to authorized users.     

Multiplex single nucleotide polymorphism (SNP)-based genotyping in allohexaploid wheat using padlock probes

Single nucleotide polymorphisms are the most common polymorphism in plant and animal genomes and, as such, are the logical choice for marker-assisted selection. However, many plants are also polyploid, and marker-assisted selection can be complicated by the presence of highly similar, but non-allelic, homoeologous sequences. Despite this, there is practical and academic demand for high-throughput genotyping in several polyploid crop species, such as allohexaploid wheat. In this paper, we present such a system, which utilizes public single nucleotide polymorphisms previously identified in both agronomically important genes and in randomly selected, mapped, expressed sequence tags developed by the wheat community. To achieve relatively high levels of multiplexing, we used non-amplified genomic DNA and padlock probe pairs, together with high annealing temperatures, to differentiate between similar sequences in the wheat genome. Our results suggest that padlock probes are capable of discriminating between homoeologous sequences and hence can be used to efficiently genotype wheat varieties.     

Targeted enrichment beyond the consensus coding DNA sequence exome reveals exons with higher variant densities

Background  

Enrichment of loci by DNA hybridization-capture, followed by high-throughput sequencing, is an important tool in modern genetics. Currently, the most common targets for enrichment are the protein coding exons represented by the consensus coding DNA sequence (CCDS). The CCDS, however, excludes many actual or computationally predicted coding exons present in other databases, such as RefSeq and Vega, and non-coding functional elements such as untranslated and regulatory regions. The number of variants per base pair (variant density) and our ability to interrogate regions outside of the CCDS regions is consequently less well understood.     

Massive alterations of the methylation patterns around DNA transposons in the first four generations of a newly formed wheat allohexaploid

Rapid and reproducible genomic changes can be induced during the early stages of the life of nascent allopolyploid species. In a previous study, it was shown that following allopolyploidization, cytosine methylation changes can affect up to 11% of the wheat genome. However, the methylation patterns around transposable elements (TEs) were never studied in detail. We used transposon methylation display (TMD) to assess the methylation patterns of CCGG sites flanking three TE families (Balduin, Apollo, and Thalos) in the first four generations of a newly formed wheat allohexaploid. In addition, transposon display (TD), using a methylation-insensitive restriction enzyme, was applied to search for genomic rearrangements at the TE insertion sites. We observed that up to 54% of CCGG sites flanking the three TE families showed changes in methylation patterns in the first four generations of a newly formed wheat allohexaploid, where hypermethylation was predominant. Over 70% of the changes in TMD patterns occurred in the first two generations of the newly formed allohexaploid. Furthermore, analysis of 555 TE insertion sites by TD and 18 cases by site-specific PCR revealed a full additive pattern in the allohexaploid, an indication for lack of massive rearrangements. These data indicate that following allopolyplodization, DNA-TE insertion sites can undergo a significantly high level of methylation changes compared with methylation changes of other genomic sequences.     

Irregular patterns of transgene silencing in allohexaploid oat

An irregular pattern of transgene silencing was revealed in expression and inheritance studies conducted over multiple generations following transgene introduction by microprojectile bombardment of allohexaploid cultivated oat (Avena sativa L.). Expression of two transgenes, bar and uidA, delivered on the same plasmid was investigated in 23 transgenic oat lines. Twenty-one transgenic lines, each derived from an independently selected transformed tissue culture, showed expression of both bar and uidA while two lines expressed only bar. The relationship of the transgenic phenotypes to the presence of the transgenes in the study was determined using (1) phenotypic scoring combined with Southern blot analyses of progeny, (2) coexpression of the two transgenic phenotypes since the two transgenes always cosegregated, and (3) reactivation of a transgenic phenotype in self-pollinated progenies of transgenic plants that did not exhibit a transgenic phenotype. Transgene silencing was observed in 19 of the 23 transgenic lines and resulted in distorted segregation of transgenic phenotypes in 10 lines. Silencing and inheritance distortions were irregular and unpredictable. They were often reversible in a subsequent generation of self-pollinated progeny and abnormally segregating progenies were as likely to trace back to parents that exhibited normal segregation in a previous generation as to parents showing segregation distortions. Possible causes of the irregular patterns of transgene silencing are discussed.     

Targeted/exome sequencing identified mutations in ten Chinese patients diagnosed with Noonan syndrome and related disorders

Background

Noonan syndrome (NS) and Noonan syndrome with multiple lentigines (NSML) are autosomal dominant developmental disorders. NS and NSML are caused by abnormalities in genes that encode proteins related to the RAS-MAPK pathway, including PTPN11, RAF1, BRAF, and MAP2K. In this study, we diagnosed ten NS or NSML patients via targeted sequencing or whole exome sequencing (TS/WES).

Methods

TS/WES was performed to identify mutations in ten Chinese patients who exhibited the following manifestations: potential facial dysmorphisms, short stature, congenital heart defects, and developmental delay. Sanger sequencing was used to confirm the suspected pathological variants in the patients and their family members.

Results

TS/WES revealed three mutations in the PTPN11 gene, three mutations in RAF1 gene, and four mutations in BRAF gene in the NS and NSML patients who were previously diagnosed based on the abovementioned clinical features. All the identified mutations were determined to be de novo mutations. However, two patients who carried the same mutation in the RAF1 gene presented different clinical features. One patient with multiple lentigines was diagnosed with NSML, while the other patient without lentigines was diagnosed with NS. In addition, a patient who carried a hotspot mutation in the BRAF gene was diagnosed with NS instead of cardiofaciocutaneous syndrome (CFCS).

Conclusions

TS/WES has emerged as a useful tool for definitive diagnosis and accurate genetic counseling of atypical cases. In this study, we analyzed ten Chinese patients diagnosed with NS and related disorders and identified their correspondingPTPN11, RAF1, and BRAF mutations. Among the target genes, BRAF showed the same degree of correlation with NS incidence as that of PTPN11 or RAF1.

     

Targeted and efficient transfer of value-added genes into a wheat variety

With an objective to optimize an approach to transfer value-added genes to a wheat variety while maintaining and improving agronomic performance, two alleles (Als1 and Als2) with mutations in the acetolactate synthase (ALS) gene located on the long arm of wheat chromosomes 6B and 6D providing tolerance to imazamox herbicide were transferred to Eltan, a popular soft white common winter wheat cultivar in the Pacific Northwest (PNW), USA. Four-step marker-assisted background selection and marker assisted forward breeding approaches were used to develop a wheat variety carrying two genes (Als1 and Als2) for imazamox tolerance along with improvements in many other agronomic traits. Screening of 1600 BC₁ plants for imazamox tolerance identified 378 plants that were further screened with markers to identify seven plants that were used to make a population of 1400 BC₂ plants, and the selection cycle was repeated. Progeny of 17 selected BC₂F₁ plants was evaluated for various agronomic and quality parameters to select 12 plants that were increased for field testing. Field evaluation of these lines conducted over 58 location-years along with evaluation in the greenhouse/growth chamber led to the selection of a line “WA8143” carrying the two genes for imazamox tolerance that yielded >3% higher than Eltan did. With 96.8% similarity to the recurrent parent, WA8143 also showed a better disease resistance package and grain quality needed in a successful Pacific Northwest wheat variety and was subsequently released for cultivation under the name of “Curiosity CL+.”     

Targeted modification of wheat grain protein to reduce the content of celiac causing epitopes

The prolamin peptides in wheat gluten and in the homologous storage proteins of barley and rye cause painful chronic erasure of microvilli of the small intestine epithelium in celiac patients. If untreated, it can lead to chronic diarrhea, abdominal distension, osteoporosis, weight-loss due to malabsorption of nutrients, and anemia. In addition to congenital cases, life-long exposure to gluten proteins in bread and pasta can also induce development of celiac sprue in adults. To date, the only effective treatment is life-long strict abstinence from the staple food grains. Complete exclusion of dietary gluten is, however, difficult due to use of wheat in many foods, incomplete labeling and social constraints. Thus, finding alternative therapies for this most common foodborne disease remained an active area of research, which has led to many suggestions in last few years. The pros and cons associated with these therapies were reviewed in the present communication. As different celiac patients are immunogenic to different members of the undigestible proline/glutamine rich peptides of ～149 gliadins and low molecular weight glutenin subunits as well as the six high molecular weight glutenin subunits, an exhaustive digestion of the immunogenic peptides in the stomach, duodenum, jejunum, and ileum of celiacs is required. In view of the above, we evaluated the capacity of cereal grains to synthesize and store the enzymes prolyl endopeptidase from Flavobacterium meningosepticum and the barley cysteine endoprotease B2, which in combination are capable of detoxifying immunogenic gluten peptides in a novel treatment of celiac disease.     

Simultaneous gene editing of three homoeoalleles in self-incompatible allohexaploid grasses

Clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) system has been widely used for precise gene editing in plants. However, simultaneous gene editing of multiple homoeoalleles remains challenging, especially in self-incompatible polyploid plants. Here, we simultaneously introduced targeted mutations in all three homoeoalleles of two genes in the self-incompatible allohexaploid tall fescue, using both CRISPR/Cas9 and LbCas12a (LbCpf1) systems. Loss-of-function mutants of FaPDS exhibited albino leaves, while knockout of FaHSP17.9 resulted in impaired heat resistance in T0 generation of tall fescue. Moreover, these mutations were inheritable. Our findings demonstrate the feasibility of generating loss-of-function mutants in T0 generation polyploid perennial grasses using CRISPR/Cas systems.     

Targeted high throughput sequencing of a cancer-related exome subset by specific sequence capture with a fully automated microarray platform

Sequence capture methods for targeted next generation sequencing promise to massively reduce cost of genomics projects compared to untargeted sequencing. However, evaluated capture methods specifically dedicated to biologically relevant genomic regions are rare. Whole exome capture has been shown to be a powerful tool to discover the genetic origin of disease and provides a reduction in target size and thus calculative sequencing capacity of > 90-fold compared to untargeted whole genome sequencing. For further cost reduction, a valuable complementing approach is the analysis of smaller, relevant gene subsets but involving large cohorts of samples. However, effective adjustment of target sizes and sample numbers is hampered by the limited scalability of enrichment systems. We report a highly scalable and automated method to capture a 480 Kb exome subset of 115 cancer-related genes using microfluidic DNA arrays. The arrays are adaptable from 125 Kb to 1 Mb target size and/or one to eight samples without barcoding strategies, representing a further 26 – 270-fold reduction of calculative sequencing capacity compared to whole exome sequencing. Illumina GAII analysis of a HapMap genome enriched for this exome subset revealed a completeness of > 96%. Uniformity was such that > 68% of exons had at least half the median depth of coverage. An analysis of reference SNPs revealed a sensitivity of up to 93% and a specificity of 98.2% or higher.     

Reducing the exome search space for Mendelian diseases using genetic linkage analysis of exome genotypes

Many exome sequencing studies of Mendelian disorders fail to optimally exploit family information. Classical genetic linkage analysis is an effective method for eliminating a large fraction of the candidate causal variants discovered, even in small families that lack a unique linkage peak. We demonstrate that accurate genetic linkage mapping can be performed using SNP genotypes extracted from exome data, removing the need for separate array-based genotyping. We provide software to facilitate such analyses.