Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

The evolution of the human DNA replication timing program.

DNA is replicated according to a defined spatiotemporal program that is linked to both gene regulation and genome stability. The evolutionary forces that have shaped replication timing programs in eukaryotic species are largely unknown. Here, we studied the molecular causes and consequences of replication timing evolution across 94 humans, 95 chimpanzees, and 23 rhesus macaques. Replication timing differences recapitulated the species' phylogenetic tree, suggesting continuous evolution of the DNA replication timing program in primates. Hundreds of genomic regions had significant replication timing variation between humans and chimpanzees, of which 66 showed advances in replication origin firing in humans, while 57 were delayed. Genes overlapping these regions displayed correlated changes in expression levels and chromatin structure. Many human-chimpanzee variants also exhibited interindividual replication timing variation, pointing to ongoing evolution of replication timing at these loci. Association of replication timing variation with genetic variation revealed that DNA sequence evolution can explain replication timing variation between species. Taken together, DNA replication timing shows substantial and ongoing evolution in the human lineage that is driven by sequence alterations and could impact regulatory evolution at specific genomic sites.

Delayed DNA replication in haploid human embryonic stem cells.

Haploid human embryonic stem cells (ESCs) provide a powerful genetic system but diploidize at high rates. We hypothesized that diploidization results from aberrant DNA replication. To test this, we profiled DNA replication timing in isogenic haploid and diploid ESCs. The greatest difference was the earlier replication of the X Chromosome in haploids, consistent with the lack of X-Chromosome inactivation. We also identified 21 autosomal regions that had delayed replication in haploids, extending beyond the normal S phase and into G2/M. Haploid-delays comprised a unique set of quiescent genomic regions that are also underreplicated in polyploid placental cells. The same delays were observed in female ESCs with two active X Chromosomes, suggesting that increased X-Chromosome dosage may cause delayed autosomal replication. We propose that incomplete replication at the onset of mitosis could prevent cell division and result in re-entry into the cell cycle and whole genome duplication.

Accurate and Automated Genotyping of the CFTR Poly-T/TG Tract with CFTR-TIPS.

Cystic fibrosis is caused by biallelic pathogenic variants in the CFTR gene, which contains a polymorphic (TG)mTn sequence (the "poly-T/TG tract") in intron 9. While T9 and T7 alleles are benign, T5 alleles with longer TG repeats, e.g., (TG)12T5 and (TG)13T5, are clinically significant. Thus, professional medical societies currently recommend reporting the TG repeat size when T5 is detected. Sanger sequencing is a cost-effective method of genotyping the (TG)mTn tract; however, its polymorphic length substantially complicates data analysis. We developed CFTR-TIPS, a freely available web-based software tool that infers the (TG)mTn genotype from Sanger sequencing data. This tool detects the (TG)mTn tract in the chromatograms, quantifies goodness of fit with expected patterns, and visualizes the results in a graphical user interface. It is broadly compatible with any Sanger chromatogram that contains the (TG)mTn tract ± 15 bp. We evaluated CFTR-TIPS using 835 clinical samples previously analyzed in a CLIA-certified, CAP-accredited laboratory. When operated fully automatically, CFTR-TIPS achieved 99.8% concordance with our clinically validated manual workflow, while generally taking less than 10 s per sample. There were two discordant samples: one due to a co-occurring heterozygous duplication that confounded the tool and the other due to incomplete (TG)mTn tract detection in the reverse chromatogram. No clinically significant misclassifications were observed. CFTR-TIPS is a free, accurate, and rapid tool for CFTR (TG)mTn tract genotyping using cost-effective Sanger sequencing. This tool is suitable both for automated use and as an aid to manual review to enhance accuracy and reduce analysis time.

Positive and Negative Regulation of DNA Replication Initiation.

Genomic DNA is replicated every cell cycle by the programmed activation of replication origins at specific times and chromosomal locations. The factors that define the locations of replication origins and their typical activation times in eukaryotic cells are poorly understood. Previous studies highlighted the role of activating factors and epigenetic modifications in regulating replication initiation. Here, we review the role that repressive pathways - and their alleviation - play in establishing the genomic landscape of replication initiation. Several factors mediate this repression, in particular, factors associated with inactive chromatin. Repression can support organized, yet stochastic, replication initiation, and its absence could explain instances of rapid and random replication or re-replication.

SCIP: software for efficient clinical interpretation of copy number variants detected by whole-genome sequencing.

Copy number variants (CNVs) represent major etiologic factors in rare genetic diseases. Current clinical CNV interpretation workflows require extensive back-and-forth with multiple tools and databases. This increases complexity and time burden, potentially resulting in missed genetic diagnoses. We present the Suite for CNV Interpretation and Prioritization (SCIP), a software package for the clinical interpretation of CNVs detected by whole-genome sequencing (WGS). The SCIP Visualization Module near-instantaneously displays all information necessary for CNV interpretation (variant quality, population frequency, inheritance pattern, and clinical relevance) on a single page-supported by modules providing variant filtration and prioritization. SCIP was comprehensively evaluated using WGS data from 1027 families with congenital cardiac disease and/or autism spectrum disorder, containing 187 pathogenic or likely pathogenic (P/LP) CNVs identified in previous curations. SCIP was efficient in filtration and prioritization: a median of just two CNVs per case were selected for review, yet it captured all P/LP findings (92.5% of which ranked 1st). SCIP was also able to identify one pathogenic CNV previously missed. SCIP was benchmarked against AnnotSV and a spreadsheet-based manual workflow and performed superiorly than both. In conclusion, SCIP is a novel software package for efficient clinical CNV interpretation, substantially faster and more accurate than previous tools (available at https://github.com/qd29/SCIP , a video tutorial series is available at https://bit.ly/SCIPVideos ).

Mutation Rate Variability across Human Y-Chromosome Haplogroups.

A common assumption in dating patrilineal events using Y-chromosome sequencing data is that the Y-chromosome mutation rate is invariant across haplogroups. Previous studies revealed interhaplogroup heterogeneity in phylogenetic branch length. Whether this heterogeneity is caused by interhaplogroup mutation rate variation or nongenetic confounders remains unknown. Here, we analyzed whole-genome sequences from cultured cells derived from >1,700 males. We confirmed the presence of branch length heterogeneity. We demonstrate that sex-chromosome mutations that appear within cell lines, which likely occurred somatically or in vitro (and are thus not influenced by nongenetic confounders) are informative for germline mutational processes. Using within-cell-line mutations, we computed a relative Y-chromosome somatic mutation rate, and uncovered substantial variation (up to 83.3%) in this proxy for germline mutation rate among haplogroups. This rate positively correlates with phylogenetic branch length, indicating that interhaplogroup mutation rate variation is a likely cause of branch length heterogeneity.

Ancestral Origins and Genetic History of Tibetan Highlanders.

The origin of Tibetans remains one of the most contentious puzzles in history, anthropology, and genetics. Analyses of deeply sequenced (30×-60×) genomes of 38 Tibetan highlanders and 39 Han Chinese lowlanders, together with available data on archaic and modern humans, allow us to comprehensively characterize the ancestral makeup of Tibetans and uncover their origins. Non-modern human sequences compose ∼6% of the Tibetan gene pool and form unique haplotypes in some genomic regions, where Denisovan-like, Neanderthal-like, ancient-Siberian-like, and unknown ancestries are entangled and elevated. The shared ancestry of Tibetan-enriched sequences dates back to ∼62,000-38,000 years ago, predating the Last Glacial Maximum (LGM) and representing early colonization of the plateau. Nonetheless, most of the Tibetan gene pool is of modern human origin and diverged from that of Han Chinese ∼15,000 to ∼9,000 years ago, which can be largely attributed to post-LGM arrivals. Analysis of ∼200 contemporary populations showed that Tibetans share ancestry with populations from East Asia (∼82%), Central Asia and Siberia (∼11%), South Asia (∼6%), and western Eurasia and Oceania (∼1%). Our results support that Tibetans arose from a mixture of multiple ancestral gene pools but that their origins are much more complicated and ancient than previously suspected. We provide compelling evidence of the co-existence of Paleolithic and Neolithic ancestries in the Tibetan gene pool, indicating a genetic continuity between pre-historical highland-foragers and present-day Tibetans. In particular, highly differentiated sequences harbored in highlanders' genomes were most likely inherited from pre-LGM settlers of multiple ancestral origins (SUNDer) and maintained in high frequency by natural selection.

Reintroduction of a Homocysteine Level-Associated Allele into East Asians by Neanderthal Introgression.

In this study, we present an analysis of Neanderthal introgression at the dipeptidase 1 gene, DPEP1. A Neanderthal origin for the putative introgressive haplotypes was demonstrated using an established three-step approach. This introgression was under positive natural selection, reached a frequency of >50%, and introduced a homocysteine level- and pigmentation-associated allele (rs460879-T) into East Asians. However, the same allele was also found in non-East Asians, but not from Neanderthal introgression. It is likely that rs460879-T was lost in East Asians and was reintroduced subsequently through Neanderthal introgression. Our findings suggest that Neanderthal introgression could reintroduce an important previously existing allele into populations where the allele had been lost. This study sheds new light on understanding the contribution of Neanderthal introgression to the adaptation of non-Africans.

Neanderthal introgression at chromosome 3p21.31 was under positive natural selection in East Asians.

Studies of the Neanderthal and Denisovan genomes demonstrate archaic hominin introgression in Eurasians. Here, we present evidence of Neanderthal introgression within the chromosome 3p21.31 region, occurring with a high frequency in East Asians (ranging from 49.4% to 66.5%) and at a low frequency in Europeans. We also detected a signal of strong positive selection in this region only in East Asians. Our data indicate that likely candidate targets of selection include rs12488302-T and its associated alleles--among which four are nonsynonymous, including rs35455589-G in HYAL2, a gene related to the cellular response to ultraviolet-B irradiation. Furthermore, suggestive evidence supports latitude-dependent selection, implicating a role of ultraviolet-B. Interestingly, the distribution of rs35455589-G suggests that this allele was lost during the exodus of ancestors of modern Eurasians from Africa and reintroduced to Eurasians from Neanderthals.

Neanderthal origin of the haplotypes carrying the functional variant Val92Met in the MC1R in modern humans.

Skin color is one of the most visible and important phenotypes of modern humans. Melanocyte-stimulating hormone and its receptor played an important role in regulating skin color. In this article, we present evidence of Neanderthal introgression encompassing the melanocyte-stimulating hormone receptor gene MC1R. The haplotypes from Neanderthal introgression diverged with the Altai Neanderthal 103.3 ka, which postdates the anatomically modern human-Neanderthal divergence. We further discovered that all of the putative Neanderthal introgressive haplotypes carry the Val92Met variant, a loss-of-function variant in MC1R that is associated with multiple dermatological traits including skin color and photoaging. Frequency of this Neanderthal introgression is low in Europeans (∼5%), moderate in continental East Asians (∼30%), and high in Taiwanese aborigines (60-70%). As the putative Neanderthal introgressive haplotypes carry a loss-of-function variant that could alter the function of MC1R and is associated with multiple traits related to skin color, we speculate that the Neanderthal introgression may have played an important role in the local adaptation of Eurasians to sunlight intensity.

Using whole genome sequence findings to assess gene-disease causality in cardiomyopathy and arrhythmia patients.

Aim: The genetic etiologies of cardiomyopathies and arrhythmias have not been fully elucidated. Materials & methods: Research findings from genome analyses in a cardiomyopathy and arrhythmia cohort were gathered. Gene-disease relationships from two databases were compared with patient phenotypes. A literature review was conducted for genes with limited evidence. Results: Of 43 genes with candidate findings from 18 cases, 23.3% of genes had never been curated, 15.0% were curated for cardiomyopathies, 16.7% for arrhythmias and 31.3% for other conditions. 25.5% of candidate findings were curated for the patient's specific phenotype with 11.8% having definitive evidence. MYH6 and TPCN1 were flagged for recuration. Conclusion: Findings from genome sequencing in disease cohorts may be useful to guide gene-curation efforts.

Germline whole genome sequencing in adults with multiple primary tumors.

Multiple primary tumors (MPTs) are a harbinger of hereditary cancer syndromes. Affected individuals often fit genetic testing criteria for a number of hereditary cancer genes and undergo multigene panel testing. Other genomic testing options, such as whole exome (WES) and whole genome sequencing (WGS) are available, but the utility of these genomic approaches as a second-tier test for those with uninformative multigene panel testing has not been explored. Here, we report our germline sequencing results from WGS in 9 patients with MPTs who had non-informative multigene panel testing. Following germline WGS, sequence (agnostic or 735 selected genes) and copy number variant (CNV) analysis was performed according to the American College of Medical Genetics (ACMG) standards and guidelines for interpreting sequence variants and reporting CNVs. In this cohort, WGS, as a second-tier test, did not identify additional pathogenic or likely pathogenic variants in cancer predisposition genes. Although we identified a CHEK2 likely pathogenic variant and a MUTYH pathogenic variant, both were previously identified in the multigene panels and were not explanatory for the presented type of tumors. CNV analysis also failed to identify any pathogenic or likely pathogenic variants in cancer predisposition genes. In summary, after multigene panel testing, WGS did not reveal any additional pathogenic variants in patients with MPTs. Our study, based on a small cohort of patients with MPT, suggests that germline gene panel testing may be sufficient to investigate these cases. Future studies with larger sample sizes may further elucidate the additional utility of WGS in MPTs.

The genetic architecture of DNA replication timing in human pluripotent stem cells.

DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome's replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) - sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

Growth performance and stress responses of larval Mississippi Paddlefish Polyodon spathula to hypoxia under different diet treatments.

A growth trial was conducted to detect the effects of different diets on the growth performance and hypoxia adaptation capacity of Mississippi Paddlefish (Polyodon spathula) larvae. The larvae were fed with live food, formulated diets, and 1/2 live food with 1/2 formulated diets. After a 15-d growth trial, final body weight and total body length were measured, and five larvae from each dietary group were subjected to 1 h of hypoxia treatment. Serum total antioxidant capacity (T-AOC), serum superoxide dismutase (SOD), and liver malondialdehyde (MDA) were measured. Final body weight and weight gain of the fish fed live food were significantly higher than the values for the other two groups. Total body length of the fish fed live food and 1/2 live food with 1/2 formulated diets exhibited no significant difference. After hypoxia treatment, serum T-AOC and SOD activities of the fish fed formulated diets were significantly lower than those of the other two groups. Liver MDA content of the fish fed with live food was significantly higher than that of the other two groups. In conclusion, larval paddlefish fed with an appropriate proportion of live food and formulated diets exhibit improved adaptive capacity to hypoxia.

Comment on "Phonemic diversity supports a serial founder effect model of language expansion from Africa".

Atkinson (Reports, 15 April 2011, p. 346) reported a declined trend of phonemic diversity from Africa that indicated the African exodus of modern languages. However, his claim was only supported when the phonemic diversities were binned into three or five levels. Analyses using raw data without simplification suggest a decline from central Asia rather than from Africa.