Diversity of human copy number variation and multicopy genes

Neural progenitor cells and their role in the development and evolutionary expansion of the neocortex

The evolutionary expansion of the mammalian brain, notably the neocortex, provides a platform for the higher cognitive abilities that characterize humans. Cortical expansion is accompanied by increased folding of the pial surface, which gives rise to a gyrencephalic (folded) rather than lissencephalic (unfolded) neocortex. This expansion reflects the prolonged and increased proliferation of neural stem and progenitor cells (NPCs). Distinct classes of NPCs can be distinguished based on either cell biological criteria (apical progenitors [APs], basal progenitors [BPs]) or lineage (primary progenitors and secondary progenitors). Cortical expansion in development and evolution is linked to an increased abundance and proliferative capacity of BPs, notably basal radial glial cells, a recently characterized type of secondary progenitor derived from apical radial glial cells, the primary progenitors. To gain insight into the molecular basis underlying the prolonged and increased proliferation of NPCs and in particular BPs, comparative genomic and transcriptomic approaches, mostly for human versus mouse, have been employed and applied to specific NPC types and subpopulations. These have revealed two principal sets of molecular changes. One concerns differences in the expression of common genes between species with different degrees of cortical expansion. The other comprises human-specific genes or genomic regulatory sequences. Various systems that allow functional testing of these genomic and gene expression differences between species have emerged, including transient and stable transgenesis, genome editing, cerebral organoids, and organotypic slice cultures. These provide future avenues for uncovering the molecular basis of cortical expansion. WIREs Dev Biol 2017, 6:e256. doi: 10.1002/wdev.256 For further resources related to this article, please visit the WIREs website.

Characterizing nuclear and mitochondrial DNA in spent embryo culture media: genetic contamination identified

OBJECTIVE

To characterize nuclear and mitochondrial DNA (mtDNA) in spent culture media from normally developing blastocysts to determine whether it could be used for noninvasive genetic assessment.

DESIGN

Prospective embryo cohort study.

SETTING

Academic center and private in vitro fertilization (IVF) clinic.

PATIENT(S)

Seventy patients undergoing intracytoplasmic sperm injection (ICSI) and 227 blastocysts.

INTERVENTION(S)

Culture media assessment, artificial blastocoele fluid collapse and DNA analysis using digital polymerase chain reaction (dPCR), long-range PCR, quantitative PCR (qPCR), and DNA fingerprinting.

MAIN OUTCOME MEASURE(S)

Presence of nuclear and mtDNA in three different commercial culture media from Vitrolife and Irvine Scientific, spent embryo media assessment at the cleavage and blastocyst stages of development, and analysis of the internal media controls for each patient that had been exposed to identical conditions as embryo media but did not come into contact with embryos.

RESULT(S)

Higher levels of nuclear and mtDNA were observed in the culture media that had been exposed to embryos compared with the internal media controls. Nuclear DNA (∼4 copies) and mtDNA (∼600 copies) could be detected in spent media, and the levels increased at the blastocyst stage. No increase in DNA was detected after artificial blastocoele fluid collapse. Mixed sex chromosome DNA was detected. This originated from contamination in the culture media and from maternal (cumulus) cells. Due to the limited amount of template, the presence of embryonic nuclear DNA could not be confirmed by DNA fingerprinting analysis.

CONCLUSION(S)

Currently DNA from culture media cannot be used for genetic assessment because embryo-associated structures release DNA into the culture medium and the DNA is of mixed origin.

Arabidopsis thaliana population analysis reveals high plasticity of the genomic region spanning MSH2, AT3G18530 and AT3G18535 genes and provides evidence for NAHR-driven recurrent CNV events occurring in this location

Background

Intraspecies copy number variations (CNVs), defined as unbalanced structural variations of specific genomic loci, ≥1 kb in size, are present in the genomes of animals and plants. A growing number of examples indicate that CNVs may have functional significance and contribute to phenotypic diversity. In the model plant Arabidopsis thaliana at least several hundred protein-coding genes might display CNV; however, locus-specific genotyping studies in this plant have not been conducted.

Results

We analyzed the natural CNVs in the region overlapping MSH2 gene that encodes the DNA mismatch repair protein, and AT3G18530 and AT3G18535 genes that encode poorly characterized proteins. By applying multiplex ligation-dependent probe amplification and droplet digital PCR we genotyped those genes in 189 A. thaliana accessions. We found that AT3G18530 and AT3G18535 were duplicated (2–14 times) in 20 and deleted in 101 accessions. MSH2 was duplicated in 12 accessions (up to 12-14 copies) but never deleted. In all but one case, the MSH2 duplications were associated with those of AT3G18530 and AT3G18535. Considering the structure of the CNVs, we distinguished 5 genotypes for this region, determined their frequency and geographical distribution. We defined the CNV breakpoints in 35 accessions with AT3G18530 and AT3G18535 deletions and tandem duplications and showed that they were reciprocal events, resulting from non-allelic homologous recombination between 99 %-identical sequences flanking these genes. The widespread geographical distribution of the deletions supported by the SNP and linkage disequilibrium analyses of the genomic sequence confirmed the recurrent nature of this CNV.

Conclusions

We characterized in detail for the first time the complex multiallelic CNV in Arabidopsis genome. The region encoding MSH2, AT3G18530 and AT3G18535 genes shows enormous variation of copy numbers among natural ecotypes, being a remarkable example of high Arabidopsis genome plasticity. We provided the molecular insight into the mechanism underlying the recurrent nature of AT3G18530-AT3G18535 duplications/deletions. We also performed the first direct comparison of the two leading experimental methods, suitable for assessing the DNA copy number status. Our comprehensive case study provides foundation information for further analyses of CNV evolution in Arabidopsis and other plants, and their possible use in plant breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3221-1) contains supplementary material, which is available to authorized users.

Genetics of Congenital Heart Disease: Past and Present

Congenital heart disease is the most common congenital anomaly, representing an important cause of infant morbidity and mortality. Congenital heart disease represents a group of heart anomalies that include septal defects, valve defects, and outflow tract anomalies. The exact genetic, epigenetic, or environmental basis of congenital heart disease remains poorly understood, although the exact mechanism is likely multifactorial. However, the development of new technologies including copy number variants, single-nucleotide polymorphism, next-generation sequencing are accelerating the detection of genetic causes of heart anomalies. Recent studies suggest a role of small non-coding RNAs, micro RNA, in congenital heart disease. The recently described epigenetic factors have also been found to contribute to cardiac morphogenesis. In this review, we present past and recent genetic discoveries in congenital heart disease.

Inter-varietal structural variation in grapevine genomes

Grapevine (Vitis vinifera L.) is one of the world's most important crop plants, which is of large economic value for fruit and wine production. There is much interest in identifying genomic variations and their functional effects on inter-varietal, phenotypic differences. Using an approach developed for the analysis of human and mammalian genomes, which combines high-throughput sequencing, array comparative genomic hybridization, fluorescent in situ hybridization and quantitative PCR, we created an inter-varietal atlas of structural variations and single nucleotide variants (SNVs) for the grapevine genome analyzing four economically and genetically relevant table grapevine varieties. We found 4.8 million SNVs and detected 8% of the grapevine genome to be affected by genomic variations. We identified more than 700 copy number variation (CNV) regions and more than 2000 genes subjected to CNV as potential candidates for phenotypic differences between varieties.

Copy Number Variants in Patients with Autism and Additional Clinical Features: Report of VIPR2 Duplication and a Novel Microduplication Syndrome

Autism is a common neurodevelopmental disorder estimated to affect 1 in 68 children. Many studies have shown the role of copy number variants (CNVs) as a major contributor in the etiology of autism with the overall detection rate of about 10-15 % and over 20 % when syndromic forms of autism exist. In this study, we used array CGH to identify CNVs in 15 Iranian patients with autism. To elevate our diagnostic yield, we selected the sporadic patients who had additional clinical features including intellectual disability (ID), craniofacial anomaly, and seizure. Six out of 15 patients showed clinically relevant CNVs including pathogenic and likely pathogenic copy number gains or losses. We report a novel gene duplication syndrome (10q21.2q21.3 microduplication) and present a new evidence for VIPR2 duplication, as a candidate gene for autism. Furthermore, we describe the first manifesting carrier female with deletion of SLC6A8 and BCAP31 genes on Xq28. Our findings suggest that there might be a higher prevalence of clinically significant CNVs in patients with autism and additional clinical manifestations. The CNV analysis in such patients could lead to the discovery of novel syndromes as well as unraveling the etiology of autism.

Interchromosomal core duplicons drive both evolutionary instability and disease susceptibility of the Chromosome 8p23.1 region

Recurrent rearrangements of Chromosome 8p23.1 are associated with congenital heart defects and developmental delay. The complexity of this region has led to inconsistencies in the current reference assembly, confounding studies of genetic variation. Using comparative sequence-based approaches, we generated a high-quality 6.3-Mbp alternate reference assembly of an inverted Chromosome 8p23.1 haplotype. Comparison with nonhuman primates reveals a 746-kbp duplicative transposition and two separate inversion events that arose in the last million years of human evolution. The breakpoints associated with these rearrangements map to an ape-specific interchromosomal core duplicon that clusters at sites of evolutionary inversion (P = 7.8 × 10⁻⁵). Refinement of microdeletion breakpoints identifies a subgroup of patients that map to the same interchromosomal core involved in the evolutionary formation of the duplication blocks. Our results define a higher-order genomic instability element that has shaped the structure of specific chromosomes during primate evolution contributing to rearrangements associated with inversion and disease.

Copy number variation contributes to cryptic genetic variation in outbreak lineages of Cryptococcus gattii from the North American Pacific Northwest

Background

Copy number variants (CNVs) are a class of structural variants (SVs) and are defined as fragments of DNA that are present at variable copy number in comparison with a reference genome. Recent advances in bioinformatics methodologies and sequencing technologies have enabled the high-resolution quantification of genome-wide CNVs. In pathogenic fungi SVs have been shown to alter gene expression, influence host specificity, and drive fungicide resistance, but little attention has focused specifically on CNVs. Using publicly available sequencing data, we identified 90 isolates across 212 Cryptococcus gattii genomes that belong to the VGII subgroups responsible for the recent deadly outbreaks in the North American Pacific Northwest. We generated CNV profiles for each sample to investigate the prevalence and function of CNV in C. gattii.

Results

We identified eight genetic clusters among publicly available Illumina whole genome sequence data from 212 C. gattii isolates through population structure analysis. Three clusters represent the VGIIa, VGIIb, and VGIIc subgroups from the North American Pacific Northwest. CNV was bioinformatically predicted and affected ~300–400 Kilobases (Kb) of the C. gattii VGII subgroup genomes. Sixty-seven loci, encompassing 58 genes, showed highly divergent patterns of copy number variation between VGII subgroups. Analysis of PFam domains within divergent CN variable genes revealed enrichment of protein domains associated with transport, cell wall organization and external encapsulating structure.

Conclusions

CNVs may contribute to pathological and phenotypic differences observed between the C. gattii VGIIa, VGIIb, and VGIIc subpopulations. Genes overlapping with population differentiated CNVs were enriched for several virulence related functional terms. These results uncover novel candidate genes to examine the genetic and functional underpinnings of C. gattii pathogenicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3044-0) contains supplementary material, which is available to authorized users.

Fine-Mapping of the 1p11.2 Breast Cancer Susceptibility Locus

The Cancer Genetic Markers of Susceptibility genome-wide association study (GWAS) originally identified a single nucleotide polymorphism (SNP) rs11249433 at 1p11.2 associated with breast cancer risk. To fine-map this locus, we genotyped 92 SNPs in a 900kb region (120,505,799-121,481,132) flanking rs11249433 in 45,276 breast cancer cases and 48,998 controls of European, Asian and African ancestry from 50 studies in the Breast Cancer Association Consortium. Genotyping was done using iCOGS, a custom-built array. Due to the complicated nature of the region on chr1p11.2: 120,300,000-120,505,798, that lies near the centromere and contains seven duplicated genomic segments, we restricted analyses to 429 SNPs excluding the duplicated regions (42 genotyped and 387 imputed). Per-allelic associations with breast cancer risk were estimated using logistic regression models adjusting for study and ancestry-specific principal components. The strongest association observed was with the original identified index SNP rs11249433 (minor allele frequency (MAF) 0.402; per-allele odds ratio (OR) = 1.10, 95% confidence interval (CI) 1.08-1.13, P = 1.49 x 10-21). The association for rs11249433 was limited to ER-positive breast cancers (test for heterogeneity P≤8.41 x 10-5). Additional analyses by other tumor characteristics showed stronger associations with moderately/well differentiated tumors and tumors of lobular histology. Although no significant eQTL associations were observed, in silico analyses showed that rs11249433 was located in a region that is likely a weak enhancer/promoter. Fine-mapping analysis of the 1p11.2 breast cancer susceptibility locus confirms this region to be limited to risk to cancers that are ER-positive.

From stem cells to comparative corticogenesis: a bridge too far?

It has been hypothesized that the higher number of neurons in human cortex compared to the chimpanzee and other primate species is key to high cognitive function. Are human cortical precursors endowed with specific properties that drive greater neuronal expansion than in other non-human primates? Otani et al. 2016 addressed this issue taking advantage of comparative in vitro corticogenesis models based on human, chimpanzee and macaque pluripotent stem cells. Clonal analysis revealed a heterochrony of early developmental events possibly leading to a relatively higher expansion of human cortical precursor population. In absence of evidence going beyond putative correlation, the claim that stem cell models of cortical development indicate mechanism of cortical size regulation needs to be further examined notably with respect to in vivo observations of cortical precursor lineages.

Emergence of a Homo sapiens-specific gene family and chromosome 16p11.2 CNV susceptibility

Genetic differences that specify unique aspects of human evolution have typically been identified by comparative analyses between the genomes of humans and closely related primates, including more recently the genomes of archaic hominins. Not all regions of the genome, however, are equally amenable to such study. Recurrent copy number variation (CNV) at chromosome 16p11.2 accounts for approximately 1% of cases of autism and is mediated by a complex set of segmental duplications, many of which arose recently during human evolution. Here we reconstruct the evolutionary history of the locus and identify bolA family member 2 (BOLA2) as a gene duplicated exclusively in Homo sapiens. We estimate that a 95-kilobase-pair segment containing BOLA2 duplicated across the critical region approximately 282 thousand years ago (ka), one of the latest among a series of genomic changes that dramatically restructured the locus during hominid evolution. All humans examined carried one or more copies of the duplication, which nearly fixed early in the human lineage--a pattern unlikely to have arisen so rapidly in the absence of selection (P < 0.0097). We show that the duplication of BOLA2 led to a novel, human-specific in-frame fusion transcript and that BOLA2 copy number correlates with both RNA expression (r = 0.36) and protein level (r = 0.65), with the greatest expression difference between human and chimpanzee in experimentally derived stem cells. Analyses of 152 patients carrying a chromosome 16p11. rearrangement show that more than 96% of breakpoints occur within the H. sapiens-specific duplication. In summary, the duplicative transposition of BOLA2 at the root of the H. sapiens lineage about 282 ka simultaneously increased copy number of a gene associated with iron homeostasis and predisposed our species to recurrent rearrangements associated with disease.

A Highly Polymorphic Copy Number Variant in the NSF Gene is Associated with Cocaine Dependence

Cocaine dependence is a complex psychiatric disorder involving both genetic and environmental factors. Several neurotransmitter systems mediate cocaine’s effects, dependence and relapse, being the components of the neurotransmitter release machinery good candidates for the disorder. Previously, we identified a risk haplotype for cocaine dependence in the NSF gene, encoding the protein N-Ethylmaleimide-Sensitive Factor essential for synaptic vesicle turnover. Here we examined the possible contribution to cocaine dependence of a large copy number variant (CNV) that encompasses part of the NSF gene. We performed a case-control association study in a discovery sample (359 cases and 356 controls) and identified an association between cocaine dependence and the CNV (P = 0.013), that was confirmed in the replication sample (508 cases and 569 controls, P = 7.1e-03) and in a pooled analysis (P = 1.8e-04), with an over-representation of low number of copies in cases. Subsequently, we studied the functional impact of the CNV on gene expression and found that the levels of two NSF transcripts were significantly increased in peripheral blood mononuclear cells (PBMC) along with the number of copies of the CNV. These results, together with a previous study from our group, support the role of NSF in the susceptibility to cocaine dependence.

The genome and transcriptome of Trichormus sp. NMC-1: insights into adaptation to extreme environments on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. In this study, we generated a draft genome sequence of cyanobacteria Trichormus sp. NMC-1 in the QTP and performed whole transcriptome sequencing under low temperature to investigate the genetic mechanism by which T. sp. NMC-1 adapted to the specific environment. Its genome sequence was 5.9 Mb with a G+C content of 39.2% and encompassed a total of 5362 CDS. A phylogenomic tree indicated that this strain belongs to the Trichormus and Anabaena cluster. Genome comparison between T. sp. NMC-1 and six relatives showed that functionally unknown genes occupied a much higher proportion (28.12%) of the T. sp. NMC-1 genome. In addition, functions of specific, significant positively selected, expanded orthogroups, and differentially expressed genes involved in signal transduction, cell wall/membrane biogenesis, secondary metabolite biosynthesis, and energy production and conversion were analyzed to elucidate specific adaptation traits. Further analyses showed that the CheY-like genes, extracellular polysaccharide and mycosporine-like amino acids might play major roles in adaptation to harsh environments. Our findings indicate that sophisticated genetic mechanisms are involved in cyanobacterial adaptation to the extreme environment of the QTP.

Anchored pseudo-de novo assembly of human genomes identifies extensive sequence variation from unmapped sequence reads

The human genome reference (HGR) completion marked the genomics era beginning, yet despite its utility universal application is limited by the small number of individuals used in its development. This is highlighted by the presence of high-quality sequence reads failing to map within the HGR. Sequences failing to map generally represent 2-5 % of total reads, which may harbor regions that would enhance our understanding of population variation, evolution, and disease. Alternatively, complete de novo assemblies can be created, but these effectively ignore the groundwork of the HGR. In an effort to find a middle ground, we developed a bioinformatic pipeline that maps paired-end reads to the HGR as separate single reads, exports unmappable reads, de novo assembles these reads per individual and then combines assemblies into a secondary reference assembly used for comparative analysis. Using 45 diverse 1000 Genomes Project individuals, we identified 351,361 contigs covering 195.5 Mb of sequence unincorporated in GRCh38. 30,879 contigs are represented in multiple individuals with ~40 % showing high sequence complexity. Genomic coordinates were generated for 99.9 %, with 52.5 % exhibiting high-quality mapping scores. Comparative genomic analyses with archaic humans and primates revealed significant sequence alignments and comparisons with model organism RefSeq gene datasets identified novel human genes. If incorporated, these sequences will expand the HGR, but more importantly our data highlight that with this method low coverage (~10-20×) next-generation sequencing can still be used to identify novel unmapped sequences to explore biological functions contributing to human phenotypic variation, disease and functionality for personal genomic medicine.

SRGAP2 and Its Human-Specific Paralog Co-Regulate the Development of Excitatory and Inhibitory Synapses

The proper function of neural circuits requires spatially and temporally balanced development of excitatory and inhibitory synapses. However, the molecular mechanisms coordinating excitatory and inhibitory synaptogenesis remain unknown. Here we demonstrate that SRGAP2A and its human-specific paralog SRGAP2C co-regulate the development of excitatory and inhibitory synapses in cortical pyramidal neurons in vivo. SRGAP2A promotes synaptic maturation, and ultimately the synaptic accumulation of AMPA and GABAA receptors, by interacting with key components of both excitatory and inhibitory postsynaptic scaffolds, Homer and Gephyrin. Furthermore, SRGAP2A limits the density of both types of synapses via its Rac1-GAP activity. SRGAP2C inhibits all identified functions of SRGAP2A, protracting the maturation and increasing the density of excitatory and inhibitory synapses. Our results uncover a molecular mechanism coordinating critical features of synaptic development and suggest that human-specific duplication of SRGAP2 might have contributed to the emergence of unique traits of human neurons while preserving the excitation/inhibition balance.

Copy Number Variation in TAS2R Bitter Taste Receptor Genes: Structure, Origin, and Population Genetics

Bitter taste receptor genes (TAS2Rs) harbor extensive diversity, which is broadly distributed across human populations and strongly associated with taste response phenotypes. The majority of TAS2R variation is composed of single-nucleotide polymorphisms. However, 2 closely positioned loci at 12p13, TAS2R43 and -45, harbor high-frequency deletion (Δ) alleles in which genomic segments are absent, resulting in copy number variation (CNV). To resolve their chromosomal structure and organization, we generated maps using long-range contig alignments and local sequencing across the TAS2R43-45 region. These revealed that the deletion alleles (43Δ and 45Δ) are 37.8 and 32.2kb in length, respectively and span the complete coding region of each gene (~1kb) along with extensive up- and downstream flanking sequence, producing separate CNVs at the 2 loci. Comparisons with a chimpanzee genome, which contained intact homologs of TAS2R43, -45, and nearby TAS2Rs, indicated that the deletions evolved recently, through unequal recombination in a cluster of closely related loci. Population genetic analyses in 946 subjects from 52 worldwide populations revealed that copy number ranged from 0 to 2 at both TAS2R43 and TAS2R45, with 43Δ and 45Δ occurring at high global frequencies (0.33 and 0.18). Estimated recombination rates between the loci were low (ρ = 2.7×10(-4); r = 6.6×10(-9)) and linkage disequilibrium was high (D' = 1.0), consistent with their adjacent genomic positioning and recent origin. Geographic variation pointed to an African origin for the deletions. However, no signatures of natural selection were found in population structure or integrated haplotype scores spanning the region, suggesting that patterns of diversity at TAS2R43 and -45 are primarily due to genetic drift.

LILRA6 copy number variation correlates with susceptibility to atopic dermatitis

Leukocyte immunoglobulin-like receptors (LILR) are expressed mostly on myelomonocytic cells where they are mediators of immunological tolerance. Two LILR genes, LILRA3 and LILRA6, exhibit marked copy number variation. We assessed the contribution of these genes to atopic dermatitis (AD) by analysing transmission in 378 AD families. The data indicated that copies of LILRA6 were over-transmitted to affected patients. They are consistent with a contribution of LILR genes to AD. They could affect the equilibrium between activating and inhibitory signals in the immune response.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits

Background

Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance.

Results

Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively.

Conclusions

Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2752-9) contains supplementary material, which is available to authorized users.

A High-Resolution SNP Array-Based Linkage Map Anchors a New Domestic Cat Draft Genome Assembly and Provides Detailed Patterns of Recombination

High-resolution genetic and physical maps are invaluable tools for building accurate genome assemblies, and interpreting results of genome-wide association studies (GWAS). Previous genetic and physical maps anchored good quality draft assemblies of the domestic cat genome, enabling the discovery of numerous genes underlying hereditary disease and phenotypes of interest to the biomedical science and breeding communities. However, these maps lacked sufficient marker density to order thousands of shorter scaffolds in earlier assemblies, which instead relied heavily on comparative mapping with related species. A high-resolution map would aid in validating and ordering chromosome scaffolds from existing and new genome assemblies. Here, we describe a high-resolution genetic linkage map of the domestic cat genome based on genotyping 453 domestic cats from several multi-generational pedigrees on the Illumina 63K SNP array. The final maps include 58,055 SNP markers placed relative to 6637 markers with unique positions, distributed across all autosomes and the X chromosome. Our final sex-averaged maps span a total autosomal length of 4464 cM, the longest described linkage map for any mammal, confirming length estimates from a previous microsatellite-based map. The linkage map was used to order and orient the scaffolds from a substantially more contiguous domestic cat genome assembly (Felis catus v8.0), which incorporated ∼20 × coverage of Illumina fragment reads. The new genome assembly shows substantial improvements in contiguity, with a nearly fourfold increase in N50 scaffold size to 18 Mb. We use this map to report probable structural errors in previous maps and assemblies, and to describe features of the recombination landscape, including a massive (∼50 Mb) recombination desert (of virtually zero recombination) on the X chromosome that parallels a similar desert on the porcine X chromosome in both size and physical location.

Low AMY1 Gene Copy Number Is Associated with Increased Body Mass Index in Prepubertal Boys

Background

Genome-wide association studies have identified more than 60 single nucleotide polymorphisms associated with Body Mass Index (BMI). Additional genetic variants, such as copy number variations (CNV), have also been investigated in relation to BMI. Recently, the highly polymorphic CNV in the salivary amylase (AMY1) gene, encoding an enzyme implicated in the first step of starch digestion, has been associated with obesity in adults and children. We assessed the potential association between AMY1 copy number and a wide range of BMI in a population of Italian school-children.

Methods

744 children (354 boys, 390 girls, mean age (±SD): 8.4±1.4years) underwent anthropometric assessments (height, weight) and collection of saliva samples for DNA extraction. AMY1 copies were evaluated by quantitative PCR.

Results

A significant increase of BMI z-score by decreasing AMY1 copy number was observed in boys (β: -0.117, p = 0.033), but not in girls. Similarly, waist circumference (β: -0.155, p = 0.003, adjusted for age) was negatively influenced by AMY1 copy number in boys. Boys with 8 or more AMY1 copy numbers presented a significant lower BMI z-score (p = 0.04) and waist circumference (p = 0.01) when compared to boys with less than 8 copy numbers.

Conclusions

In this pediatric-only, population-based study, a lower AMY1 copy number emerged to be associated with increased BMI in boys. These data confirm previous findings from adult studies and support a potential role of a higher copy number of the salivary AMY1 gene in protecting from excess weight gain.

Diversity and population-genetic properties of copy number variations and multicopy genes in cattle

The diversity and population genetics of copy number variation (CNV) in domesticated animals are not well understood. In this study, we analysed 75 genomes of major taurine and indicine cattle breeds (including Angus, Brahman, Gir, Holstein, Jersey, Limousin, Nelore, and Romagnola), sequenced to 11-fold coverage to identify 1,853 non-redundant CNV regions. Supported by high validation rates in array comparative genomic hybridization (CGH) and qPCR experiments, these CNV regions accounted for 3.1% (87.5 Mb) of the cattle reference genome, representing a significant increase over previous estimates of the area of the genome that is copy number variable (∼2%). Further population genetics and evolutionary genomics analyses based on these CNVs revealed the population structures of the cattle taurine and indicine breeds and uncovered potential diversely selected CNVs near important functional genes, including AOX1, ASZ1, GAT, GLYAT, and KRTAP9-1. Additionally, 121 CNV gene regions were found to be either breed specific or differentially variable across breeds, such as RICTOR in dairy breeds and PNPLA3 in beef breeds. In contrast, clusters of the PRP and PAG genes were found to be duplicated in all sequenced animals, suggesting that subfunctionalization, neofunctionalization, or overdominance play roles in diversifying those fertility-related genes. These CNV results provide a new glimpse into the diverse selection histories of cattle breeds and a basis for correlating structural variation with complex traits in the future.

Long-read sequence assembly of the gorilla genome

Accurate sequence and assembly of genomes is a critical first step for studies of genetic variation. We generated a high-quality assembly of the gorilla genome using single-molecule, real-time sequence technology and a string graph de novo assembly algorithm. The new assembly improves contiguity by two to three orders of magnitude with respect to previously released assemblies, recovering 87% of missing reference exons and incomplete gene models. Although regions of large, high-identity segmental duplications remain largely unresolved, this comprehensive assembly provides new biological insight into genetic diversity, structural variation, gene loss, and representation of repeat structures within the gorilla genome. The approach provides a path forward for the routine assembly of mammalian genomes at a level approaching that of the current quality of the human genome.

Population-genetic properties of differentiated copy number variations in cattle

While single nucleotide polymorphism (SNP) is typically the variant of choice for population genetics, copy number variation (CNV) which comprises insertion, deletion and duplication of genomic sequence, is an informative type of genetic variation. CNVs have been shown to be both common in mammals and important for understanding the relationship between genotype and phenotype. However, CNV differentiation, selection and its population genetic properties are not well understood across diverse populations. We performed a population genetics survey based on CNVs derived from the BovineHD SNP array data of eight distinct cattle breeds. We generated high resolution results that show geographical patterns of variations and genome-wide admixture proportions within and among breeds. Similar to the previous SNP-based studies, our CNV-based results displayed a strong correlation of population structure and geographical location. By conducting three pairwise comparisons among European taurine, African taurine, and indicine groups, we further identified 78 unique CNV regions that were highly differentiated, some of which might be due to selection. These CNV regions overlapped with genes involved in traits related to parasite resistance, immunity response, body size, fertility, and milk production. Our results characterize CNV diversity among cattle populations and provide a list of lineage-differentiated CNVs.

Copy number of the Adenomatous Polyposis Coli gene is not always neutral in sporadic colorectal cancers with loss of heterozygosity for the gene

Background

Changes in the number of alleles of a chromosome may have an impact upon gene expression. Loss of heterozygosity (LOH) indicates that one allele of a gene has been lost, and knowing the exact copy number of the gene would indicate whether duplication of the remaining allele has occurred. We were interested to determine the copy number of the Adenomatous Polyposis Coli (APC) gene in sporadic colorectal cancers with LOH.

Methods

We selected 38 carcinomas with LOH for the APC gene region of chromosome 5, as determined by amplification of the CA repeat region within the D5S346 loci. The copy number status of APC was ascertained using the SALSA® MLPA® P043-B1 APC Kit. LOH for the DCC gene, KRAS gene mutation, and microsatellite instability were also evaluated for each tumor, utilizing standard polymerase chain reaction methods.

Results

No tumor demonstrated microsatellite instability. LOH of the DCC gene was also present in 33 of 36 (91.7 %) informative tumors. A KRAS gene mutation was present in 16 of the 38 (42.1 %) tumors. Twenty-four (63.2 %) of the tumors were copy number neutral, 10 (26.3 %) tumors demonstrated major loss, while two (5.3 %) showed partial loss. Two tumors (5.3 %) had copy number gain.

Conclusions

Results of APC and DCC LOH, KRAS and microsatellite instability indicate our colorectal cancer cases were typical of sporadic cancers following the ‘chromosomal instability’ pathway. The majority of our colorectal carcinomas with LOH for APC gene are copy number neutral. However, one-third of our cases showed copy number loss, suggesting that duplication of the remaining allele is not required for the development of a colorectal carcinoma.

Feature co-localization landscape of the human genome

Although feature co-localizations could serve as useful guide-posts to genome architecture, a comprehensive and quantitative feature co-localization map of the human genome has been lacking. Herein we show that, in contrast to the conventional bipartite division of genomic sequences into genic and inter-genic regions, pairwise co-localizations of forty-two genomic features in the twenty-two autosomes based on 50-kb to 2,000-kb sequence windows indicate a tripartite zonal architecture comprising Genic zones enriched with gene-related features and Alu-elements; Proximal zones enriched with MIR- and L2-elements, transcription-factor-binding-sites (TFBSs), and conserved-indels (CIDs); and Distal zones enriched with L1-elements. Co-localizations between single-nucleotide-polymorphisms (SNPs) and copy-number-variations (CNVs) reveal a fraction of sequence windows displaying steeply enhanced levels of SNPs, CNVs and recombination rates that point to active adaptive evolution in such pathways as immune response, sensory perceptions, and cognition. The strongest positive co-localization observed between TFBSs and CIDs suggests a regulatory role of CIDs in cooperation with TFBSs. The positive co-localizations of cancer somatic CNVs (CNVT) with all Proximal zone and most Genic zone features, in contrast to the distinctly more restricted co-localizations exhibited by germline CNVs (CNVG), reveal disparate distributions of CNVTs and CNVGs indicative of dissimilarity in their underlying mechanisms.

Design and Implementation of the International Genetics and Translational Research in Transplantation Network

Genetic association studies of transplantation outcomes have been hampered by small samples and highly complex multifactorial phenotypes, hindering investigations of the genetic architecture of a range of comorbidities which significantly impact graft and recipient life expectancy. We describe here the rationale and design of the International Genetics & Translational Research in Transplantation Network. The network comprises 22 studies to date, including 16494 transplant recipients and 11669 donors, of whom more than 5000 are of non-European ancestry, all of whom have existing genomewide genotype data sets.

iGeneTRAiN is a consortium that has genome-wide genotype datasets. These genomic data allows robust statistically analysis of genetic associations that impact graft and patients variables such as, such as: graft survival, acute rejection, new onset of diabetes after transplantation, and delayed graft kidney function. Supplemental digital content is available in the text.

Genome Sequencing of Autism-Affected Families Reveals Disruption of Putative Noncoding Regulatory DNA

We performed whole-genome sequencing (WGS) of 208 genomes from 53 families affected by simplex autism. For the majority of these families, no copy-number variant (CNV) or candidate de novo gene-disruptive single-nucleotide variant (SNV) had been detected by microarray or whole-exome sequencing (WES). We integrated multiple CNV and SNV analyses and extensive experimental validation to identify additional candidate mutations in eight families. We report that compared to control individuals, probands showed a significant (p = 0.03) enrichment of de novo and private disruptive mutations within fetal CNS DNase I hypersensitive sites (i.e., putative regulatory regions). This effect was only observed within 50 kb of genes that have been previously associated with autism risk, including genes where dosage sensitivity has already been established by recurrent disruptive de novo protein-coding mutations (ARID1B, SCN2A, NR3C2, PRKCA, and DSCAM). In addition, we provide evidence of gene-disruptive CNVs (in DISC1, WNT7A, RBFOX1, and MBD5), as well as smaller de novo CNVs and exon-specific SNVs missed by exome sequencing in neurodevelopmental genes (e.g., CANX, SAE1, and PIK3CA). Our results suggest that the detection of smaller, often multiple CNVs affecting putative regulatory elements might help explain additional risk of simplex autism.

Database resources of the National Center for Biotechnology Information

The National Center for Biotechnology Information (NCBI) provides a large suite of online resources for biological information and data, including the GenBank(®) nucleic acid sequence database and the PubMed database of citations and abstracts for published life science journals. Additional NCBI resources focus on literature (PubMed Central (PMC), Bookshelf and PubReader), health (ClinVar, dbGaP, dbMHC, the Genetic Testing Registry, HIV-1/Human Protein Interaction Database and MedGen), genomes (BioProject, Assembly, Genome, BioSample, dbSNP, dbVar, Epigenomics, the Map Viewer, Nucleotide, Probe, RefSeq, Sequence Read Archive, the Taxonomy Browser and the Trace Archive), genes (Gene, Gene Expression Omnibus (GEO), HomoloGene, PopSet and UniGene), proteins (Protein, the Conserved Domain Database (CDD), COBALT, Conserved Domain Architecture Retrieval Tool (CDART), the Molecular Modeling Database (MMDB) and Protein Clusters) and chemicals (Biosystems and the PubChem suite of small molecule databases). The Entrez system provides search and retrieval operations for most of these databases. Augmenting many of the web applications are custom implementations of the BLAST program optimized to search specialized datasets. All of these resources can be accessed through the NCBI home page at www.ncbi.nlm.nih.gov.

Genomic divergence and brain evolution: How regulatory DNA influences development of the cerebral cortex

The cerebral cortex controls our most distinguishing higher cognitive functions. Human-specific gene expression differences are abundant in the cerebral cortex, yet we have only begun to understand how these variations impact brain function. This review discusses the current evidence linking non-coding regulatory DNA changes, including enhancers, with neocortical evolution. Functional interrogation using animal models reveals converging roles for our genome in key aspects of cortical development including progenitor cell cycle and neuronal signaling. New technologies, including iPS cells and organoids, offer potential alternatives to modeling evolutionary modifications in a relevant species context. Several diseases rooted in the cerebral cortex uniquely manifest in humans compared to other primates, thus highlighting the importance of understanding human brain differences. Future studies of regulatory loci, including those implicated in disease, will collectively help elucidate key cellular and genetic mechanisms underlying our distinguishing cognitive traits.

Image processing for optical mapping

Optical Mapping is an established single-molecule, whole-genome analysis system, which has been used to gain a comprehensive understanding of genomic structure and to study structural variation of complex genomes. A critical component of Optical Mapping system is the image processing module, which extracts single molecule restriction maps from image datasets of immobilized, restriction digested and fluorescently stained large DNA molecules. In this review, we describe robust and efficient image processing techniques to process these massive datasets and extract accurate restriction maps in the presence of noise, ambiguity and confounding artifacts. We also highlight a few applications of the Optical Mapping system.

Reassessing the Evolutionary History of the 17q21 Inversion Polymorphism

A polymorphic inversion that lies on chromosome 17q21 comprises two major haplotype families (H1 and H2) that not only differ in orientation but also in copy-number. Although the processes driving the spread of the inversion-associated lineage (H2) in humans remain unclear, a selective advantage has been proposed for one of its subtypes. Here, we genotyped a large panel of individuals from previously overlooked populations using a custom array with a unique panel of H2-specific single nucleotide polymorphisms and found a patchy distribution of H2 haplotypes in Africa, with North Africans displaying a higher frequency of inverted subtypes, when compared with Sub-Saharan groups. Interestingly, North African H2s were found to be closer to "non-African" chromosomes further supporting that these populations may have diverged more recently from groups outside Africa. Our results uncovered higher diversity within the H2 family than previously described, weakening the hypothesis of a strong selective sweep on all inverted chromosomes and suggesting a rather complex evolutionary history at this locus.

Determining multiallelic complex copy number and sequence variation from high coverage exome sequencing data

BACKGROUND

Copy number variation (CNV) is a major component of genomic variation, yet methods to accurately type genomic CNV lag behind methods that type single nucleotide variation. High-throughput sequencing can contribute to these methods by using sequence read depth, which takes the number of reads that map to a given part of the reference genome as a proxy for copy number of that region, and compares across samples. Furthermore, high-throughput sequencing also provides information on the sequence differences between copies within and between individuals.

METHODS

In this study we use high-coverage phase 3 exome sequences of the 1000 Genomes project to infer diploid copy number of the beta-defensin genomic region, a well-studied CNV that carries several beta-defensin genes involved in the antimicrobial response, signalling, and fertility. We also use these data to call sequence variants, a particular challenge given the multicopy nature of the region.

RESULTS

We confidently call copy number and sequence variation of the beta-defensin genes on 1285 samples from 26 global populations, validate copy number using Nanostring nCounter and triplex paralogue ratio test data. We use the copy number calls to verify the genomic extent of the CNV and validate sequence calls using analysis of cloned PCR products. We identify novel variation, mostly individually rare, predicted to alter amino-acid sequence in the beta-defensin genes. Such novel variants may alter antimicrobial properties or have off-target receptor interactions, and may contribute to individuality in immunological response and fertility.

CONCLUSIONS

Given that 81% of identified sequence variants were not previously in dbSNP, we show that sequence variation in multiallelic CNVs represent an unappreciated source of genomic diversity.

Genetic variation and the de novo assembly of human genomes

The discovery of genetic variation and the assembly of genome sequences are both inextricably linked to advances in DNA-sequencing technology. Short-read massively parallel sequencing has revolutionized our ability to discover genetic variation but is insufficient to generate high-quality genome assemblies or resolve most structural variation. Full resolution of variation is only guaranteed by complete de novo assembly of a genome. Here, we review approaches to genome assembly, the nature of gaps or missing sequences, and biases in the assembly process. We describe the challenges of generating a complete de novo genome assembly using current technologies and the impact that being able to perfectly sequence the genome would have on understanding human disease and evolution. Finally, we summarize recent technological advances that improve both contiguity and accuracy and emphasize the importance of complete de novo assembly as opposed to read mapping as the primary means to understanding the full range of human genetic variation.

An integrated map of structural variation in 2,504 human genomes

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

Chromosomal Rearrangements as Barriers to Genetic Homogenization between Archaic and Modern Humans

Chromosomal rearrangements, which shuffle DNA throughout the genome, are an important source of divergence across taxa. Using a paired-end read approach with Illumina sequence data for archaic humans, I identify changes in genome structure that occurred recently in human evolution. Hundreds of rearrangements indicate genomic trafficking between the sex chromosomes and autosomes, raising the possibility of sex-specific changes. Additionally, genes adjacent to genome structure changes in Neanderthals are associated with testis-specific expression, consistent with evolutionary theory that new genes commonly form with expression in the testes. I identify one case of new-gene creation through transposition from the Y chromosome to chromosome 10 that combines the 5'-end of the testis-specific gene Fank1 with previously untranscribed sequence. This new transcript experienced copy number expansion in archaic genomes, indicating rapid genomic change. Among rearrangements identified in Neanderthals, 13% are transposition of selfish genetic elements, whereas 32% appear to be ectopic exchange between repeats. In Denisovan, the pattern is similar but numbers are significantly higher with 18% of rearrangements reflecting transposition and 40% ectopic exchange between distantly related repeats. There is an excess of divergent rearrangements relative to polymorphism in Denisovan, which might result from nonuniform rates of mutation, possibly reflecting a burst of transposable element activity in the lineage that led to Denisovan. Finally, loci containing genome structure changes show diminished rates of introgression from Neanderthals into modern humans, consistent with the hypothesis that rearrangements serve as barriers to gene flow during hybridization. Together, these results suggest that this previously unidentified source of genomic variation has important biological consequences in human evolution.

Global patterns of large copy number variations in the human genome reveal complexity in chromosome organization

Global patterns of copy number variations (CNVs) in chromosomes are required to understand the dynamics of genome organization and complexity. For this study, analysis was performed using the Affymetrix Genome-Wide Human SNP Array 6.0 chip and CytoScan High-Density arrays. We identified a total of 44 109 CNVs from 1715 genomes with a mean of 25 CNVs in an individual, which established the first drafts of population-specific CNV maps providing a rationale for prioritizing chromosomal regions. About 19 905 ancient CNVs were identified across all chromosomes and populations at varying frequencies. CNV count, and sometimes CNV size, contributed to the bulk CNV size of the chromosome. Population specific lengthening and shortening of chromosomal length was observed. Sex bias for CNV presence was largely dependent on ethnicity. Lower CNV inheritance rate was observed for India, compared to YRI and CEU. A total of 33 candidate CNV hotspots from 5382 copy number (CN) variable region (CNVR) clusters were identified. Population specific CNV distribution patterns in p and q arms disturbed the assumption that CNV counts in the p arm are less common compared to long arms, and the CNV occurrence and distribution in chromosomes is length independent. This study unraveled the force of independent evolutionary dynamics on genome organization and complexity across chromosomes and populations.

Genomic approaches to studying human-specific developmental traits

Changes in developmental regulatory programs drive both disease and phenotypic differences among species. Linking human-specific traits to alterations in development is challenging, because we have lacked the tools to assay and manipulate regulatory networks in human and primate embryonic cells. This field was transformed by the sequencing of hundreds of genomes – human and non-human – that can be compared to discover the regulatory machinery of genes involved in human development. This approach has identified thousands of human-specific genome alterations in developmental genes and their regulatory regions. With recent advances in stem cell techniques, genome engineering, and genomics, we can now test these sequences for effects on developmental gene regulation and downstream phenotypes in human cells and tissues.

Global diversity, population stratification, and selection of human copy-number variation

In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide-variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Read clouds uncover variation in complex regions of the human genome

Although an increasing amount of human genetic variation is being identified and recorded, determining variants within repeated sequences of the human genome remains a challenge. Most population and genome-wide association studies have therefore been unable to consider variation in these regions. Core to the problem is the lack of a sequencing technology that produces reads with sufficient length and accuracy to enable unique mapping. Here, we present a novel methodology of using read clouds, obtained by accurate short-read sequencing of DNA derived from long fragment libraries, to confidently align short reads within repeat regions and enable accurate variant discovery. Our novel algorithm, Random Field Aligner (RFA), captures the relationships among the short reads governed by the long read process via a Markov Random Field. We utilized a modified version of the Illumina TruSeq synthetic long-read protocol, which yielded shallow-sequenced read clouds. We test RFA through extensive simulations and apply it to discover variants on the NA12878 human sample, for which shallow TruSeq read cloud sequencing data are available, and on an invasive breast carcinoma genome that we sequenced using the same method. We demonstrate that RFA facilitates accurate recovery of variation in 155 Mb of the human genome, including 94% of 67 Mb of segmental duplication sequence and 96% of 11 Mb of transcribed sequence, that are currently hidden from short-read technologies.

A genome-wide scan for copy number variations using high-density single nucleotide polymorphism array in Simmental cattle

Copy number variations (CNVs) have recently been identified as promising sources of genetic variation, complementary to single nucleotide polymorphisms (SNPs). As a result, detection of CNVs has attracted a great deal of attention. In this study, we performed genome-wide CNV detection using Illumina Bovine HD BeadChip (770k) data on 792 Simmental cattle. A total of 263 CNV regions (CNVRs) were identified, which included 137 losses, 102 gains and 24 regions classified as both loss and gain, covering 35.48 Mb (1.41%) of the bovine genome. The length of these CNVRs ranged from 10.18 kb to 1.76 Mb, with an average length of 134.78 kb and a median length of 61.95 kb. In 136 of these regions, a total of 313 genes were identified related to biological functions such as transmembrane activity and olfactory transduction activity. To validate the results, we performed quantitative PCR to detect nine randomly selected CNVRs and successfully confirmed seven (77.6%) of them. Our results present a map of cattle CNVs derived from high-density SNP data, which expands the current CNV map of the cattle genome and provides useful information for investigation of genomic structural variation in cattle.

All-or-(N)One - an epistemological characterization of the human tumorigenic neuronal paralogous FAM72 gene loci

FAM72 is a novel neuronal progenitor cell (NPC) self-renewal supporting protein expressed under physiological conditions at low levels in other tissues. Accumulating data indicate the potential pivotal tumourigenic effects of FAM72. Our in silico human genome-wide analysis (GWA) revealed that the FAM72 gene family consists of four human-specific paralogous members, all of which are located on chromosome (chr) 1. Unique asymmetric FAM72 segmental gene duplications are most likely to have occurred in conjunction with the paired genomic neighbour SRGAP2 (SLIT-ROBO Rho GTPase activating protein), as both genes have four paralogues in humans but only one vertebra-emerging orthologue in all other species. No species with two or three FAM72/SRGAP2 gene pairs could be identified, and the four exclusively human-defining ohnologues, with different mutation patterns in Homo neanderthalensis and Denisova hominin, may remain under epigenetic control through long non-coding (lnc) RNAs.