Whole Exome Sequencing of Distant Relatives in Multiplex Families Implicates Rare Variants in Candidate Genes for Oral Clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.     

Exome sequencing of three cases of familial exceptional longevity

Exceptional longevity (EL) is a rare phenotype that can cluster in families, and co‐segregation of genetic variation in these families may point to candidate genes that could contribute to extended lifespan. In this study, for the first time, we have sequenced a total of seven exomes from exceptionally long‐lived siblings (probands ≥ 103 years and at least one sibling ≥ 97 years) that come from three separate families. We have focused on rare functional variants (RFVs) which have ≤ 1% minor allele frequency according to databases and that are likely to alter gene product function. Based on this, we have identified one candidate longevity gene carrying RFVs in all three families, APOB. Interestingly, APOB is a component of lipoprotein particles together with APOE, and variants in the genes encoding these two proteins have been previously associated with human longevity. Analysis of nonfamilial EL cases showed a trend, without reaching statistical significance, toward enrichment of APOB RFVs. We have also identified candidate longevity genes shared between two families (5–13) or within individual families (66–156 genes). Some of these genes have been previously linked to longevity in model organisms, such as PPARGC1A, NRG1, RAD52, RAD51, NCOR1, and ADCY5 genes. This work provides an initial catalog of genes that could contribute to exceptional familial longevity.     

Paired analysis of tumor mutation burden calculated by targeted deep sequencing panel and whole exome sequencing in non-small cell lung cancer

Owing to rapid advancements in NGS (next generation sequen-cing), genomic alteration is now considered an essential pre-dictive biomarkers that impact the treatment decision in many cases of cancer. Among the various predictive biomarkers, tumor mutation burden (TMB) was identified by NGS and was con-sidered to be useful in predicting a clinical response in cancer cases treated by immunotherapy. In this study, we directly com-pared the lab-developed-test (LDT) results by target sequencing panel, K-MASTER panel v3.0 and whole-exome sequencing (WES) to evaluate the concordance of TMB. As an initial step, the reference materials (n = 3) with known TMB status were used as an exploratory test. To validate and evaluate TMB, we used one hundred samples that were acquired from surgically resected tissues of non-small cell lung cancer (NSCLC) patients. The TMB of each sample was tested by using both LDT and WES methods, which extracted the DNA from samples at the same time. In addition, we evaluated the impact of capture re-gion, which might lead to different values of TMB; the evalu-ation of capture region was based on the size of NGS and target sequencing panels. In this pilot study, TMB was evalu-ated by LDT and WES by using duplicated reference samples; the results of TMB showed high concordance rate (R² = 0.887). This was also reflected in clinical samples (n = 100), which showed R² of 0.71. The difference between the coding sequence ratio (3.49%) and the ratio of mutations (4.8%) indicated that the LDT panel identified a relatively higher number of mutations. It was feasible to calculate TMB with LDT panel, which can be useful in clinical practice. Furthermore, a customized approach must be developed for calculating TMB, which differs according to cancer types and specific clinical settings.     

高通量测序技术在遗传性耳聋研究中的应用及研究进展

耳聋是一种常见的严重出生缺陷,阐明遗传性耳聋的致病机理不仅能够在临床上辅助诊断,为遗传咨询及耳聋预防提供依据,而且能促进人们更深入地了解耳聋的致病机制,开发新的治疗方法。随着基因组研究技术不断创新,以全基因组测序、全外显子组测序、目标区域测序为代表的高通量测序技术在遗传性耳聋研究中已得到广泛应用。本文总结了近5年全外显子组测序和目标区域测序在遗传性耳聋致病基因研究及临床分子诊断中应用及研究进展,希望能够有助于我国临床耳聋基因诊断技术的发展及诊断水平的提升。     

Sequencing of 19,219 exomes identifies a low-frequency variant in FKBP5 promoter predisposing to high myopia in a Han Chinese population

1. Download : Download high-res image (134KB)
2. Download : Download full-size image

     

Experimental studies of adaptation in Clarkia xantiana. III. Phenotypic selection across a subspecies border

Sister taxa with distinct phenotypes often occupy contrasting environments in parapatric ranges, yet we generally do not know whether trait divergence reflects spatially varying selection. We conducted a reciprocal transplant experiment to test whether selection favors “native phenotypes” in two subspecies of Clarkia xantiana (Onagraceae), an annual plant in California. For four quantitative traits that differ between subspecies, we estimated phenotypic selection in subspecies’ exclusive ranges and their contact zone in two consecutive years. We predicted that in the arid, pollinator‐scarce eastern region, selection favors phenotypes of the native subspecies parviflora: small leaves, slow leaf growth, early flowering, and diminutive flowers. In the wetter, pollinator‐rich, western range of subspecies xantiana, we expected selection for opposite phenotypes. We investigated pollinator contributions to selection by comparing naturally pollinated and pollen‐supplemented individuals. For reproductive traits and for subspecies xantiana, selection generally matched expectations. The contact zone sometimes showed distinctive selection, and in ssp. parviflora selection sometimes favored nonnative phenotypes. Pollinators influenced selection on flowering time but not on flower size. Little temporal variation in selection occurred, possibly because of plastic trait responses across years. Though there were exceptions and some causes of selection remain obscure, phenotypic differentiation between subspecies appears to reflect spatially variable selection.     

Quantifying effects of biodiversity on ecosystem functioning across times and places

Biodiversity loss decreases ecosystem functioning at the local scales at which species interact, but it remains unclear how biodiversity loss affects ecosystem functioning at the larger scales of space and time that are most relevant to biodiversity conservation and policy. Theory predicts that additional insurance effects of biodiversity on ecosystem functioning could emerge across time and space if species respond asynchronously to environmental variation and if species become increasingly dominant when and where they are most productive. Even if only a few dominant species maintain ecosystem functioning within a particular time and place, ecosystem functioning may be enhanced by many different species across many times and places (β‐diversity). Here, we develop and apply a new approach to estimate these previously unquantified insurance effects of biodiversity on ecosystem functioning that arise due to species turnover across times and places. In a long‐term (18‐year) grassland plant diversity experiment, we find that total insurance effects are positive in sign and substantial in magnitude, amounting to 19% of the net biodiversity effect, mostly due to temporal insurance effects. Species loss can therefore reduce ecosystem functioning both locally and by eliminating species that would otherwise enhance ecosystem functioning across temporally fluctuating and spatially heterogeneous environments.     

Detection of Rare Genomic Variants from Pooled Sequencing Using SPLINTER

As DNA sequencing technology has markedly advanced in recent years², it has become increasingly evident that the amount of genetic variation between any two individuals is greater than previously thought³. In contrast, array-based genotyping has failed to identify a significant contribution of common sequence variants to the phenotypic variability of common disease^4,5. Taken together, these observations have led to the evolution of the Common Disease / Rare Variant hypothesis suggesting that the majority of the "missing heritability" in common and complex phenotypes is instead due to an individual''s personal profile of rare or private DNA variants^6-8. However, characterizing how rare variation impacts complex phenotypes requires the analysis of many affected individuals at many genomic loci, and is ideally compared to a similar survey in an unaffected cohort. Despite the sequencing power offered by today''s platforms, a population-based survey of many genomic loci and the subsequent computational analysis required remains prohibitive for many investigators.To address this need, we have developed a pooled sequencing approach^1,9 and a novel software package¹ for highly accurate rare variant detection from the resulting data. The ability to pool genomes from entire populations of affected individuals and survey the degree of genetic variation at multiple targeted regions in a single sequencing library provides excellent cost and time savings to traditional single-sample sequencing methodology. With a mean sequencing coverage per allele of 25-fold, our custom algorithm, SPLINTER, uses an internal variant calling control strategy to call insertions, deletions and substitutions up to four base pairs in length with high sensitivity and specificity from pools of up to 1 mutant allele in 500 individuals. Here we describe the method for preparing the pooled sequencing library followed by step-by-step instructions on how to use the SPLINTER package for pooled sequencing analysis (http://www.ibridgenetwork.org/wustl/splinter). We show a comparison between pooled sequencing of 947 individuals, all of whom also underwent genome-wide array, at over 20kb of sequencing per person. Concordance between genotyping of tagged and novel variants called in the pooled sample were excellent. This method can be easily scaled up to any number of genomic loci and any number of individuals. By incorporating the internal positive and negative amplicon controls at ratios that mimic the population under study, the algorithm can be calibrated for optimal performance. This strategy can also be modified for use with hybridization capture or individual-specific barcodes and can be applied to the sequencing of naturally heterogeneous samples, such as tumor DNA.     

Variation in complex courtship traits across a hybrid zone between grasshopper species of the Chorthippus albomarginatus group

Closely related grasshopper species of the Chorthippus albomarginatus group are notable for their extremely complex courtship songs, accompanied by a visual display. Two species of this group, Ch. albomarginatus and Ch. oschei, were previously shown to hybridize in a wide mosaic hybrid zone in Ukraine and Moldova. In this paper, variation in five courtship song characters, one character of visual display and the number of stridulatory pegs were analysed across the hybrid zone to estimate selection against hybrids and strength of assortative mating. Comparison of cline width and position across the hybrid zone showed concordant and coincident clines in four traits, such as three song characters and one morphological character, and discordant and non‐coincident clines in two other song characters and the character of visual display. Concordance of clines in different characters suggests an equal strength of selection acting on underlying loci. Increase of variance and covariance between phenotypic traits at the cline centre could more likely result from assortative mating than from selection against hybrids. Most pairwise cases showed the highest covariance for the oschei‐like, than for the albomarginatus‐like hybrid populations. This indicates that introgression of the oschei genes into the albomarginatus genome is stronger than vice versa, and may be evidence of the movement of the hybrid zone in favour of Ch. albomarginatus. Analysis of associations between phenotype and local vegetation showed that mosaic structure of the hybrid zone is explained to a great extent by habitat–phenotype associations. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 275–291.