Molecular mechanisms for genomic disorders

Long-read sequencing and structural variant characterization in 1,019 samples from the 1000 Genomes Project

Structural variants (SVs) contribute significantly to human genetic diversity and disease 1-4 . Previously, SVs have remained incompletely resolved by population genomics, with short-read sequencing facing limitations in capturing the whole spectrum of SVs at nucleotide resolution 5-7 . Here we leveraged nanopore sequencing 8 to construct an intermediate coverage resource of 1,019 long-read genomes sampled within 26 human populations from the 1000 Genomes Project. By integrating linear and graph-based approaches for SV analysis via pangenome graph-augmentation, we uncover 167,291 sequence-resolved SVs in these samples, considerably advancing SV characterization compared to population-wide short-read sequencing studies 3,4 . Our analysis details diverse SV classes-deletions, duplications, insertions, and inversions-at population-scale. LINE-1 and SVA retrotransposition activities frequently mediate transductions 9,10 of unique sequences, with both mobile element classes transducing sequences at either the 3'- or 5'-end, depending on the source element locus. Furthermore, analyses of SV breakpoint junctions suggest a continuum of homology-mediated rearrangement processes are integral to SV formation, and highlight evidence for SV recurrence involving repeat sequences. Our open-access dataset underscores the transformative impact of long-read sequencing in advancing the characterisation of polymorphic genomic architectures, and provides a resource for guiding variant prioritisation in future long-read sequencing-based disease studies.

Unravelling the Genetic Landscape of Hemiplegic Migraine: Exploring Innovative Strategies and Emerging Approaches

Migraine is a severe, debilitating neurovascular disorder. Hemiplegic migraine (HM) is a rare and debilitating neurological condition with a strong genetic basis. Sequencing technologies have improved the diagnosis and our understanding of the molecular pathophysiology of HM. Linkage analysis and sequencing studies in HM families have identified pathogenic variants in ion channels and related genes, including CACNA1A, ATP1A2, and SCN1A, that cause HM. However, approximately 75% of HM patients are negative for these mutations, indicating there are other genes involved in disease causation. In this review, we explored our current understanding of the genetics of HM. The evidence presented herein summarises the current knowledge of the genetics of HM, which can be expanded further to explain the remaining heritability of this debilitating condition. Innovative bioinformatics and computational strategies to cover the entire genetic spectrum of HM are also discussed in this review.

Multiple paralogues and recombination mechanisms drive the high incidence of 22q11.2 Deletion Syndrome

The 22q11.2 deletion syndrome (22q11.2DS) is the most common microdeletion disorder. Why the incidence of 22q11.2DS is much greater than that of other genomic disorders remains unknown. Short read sequencing cannot resolve the complex segmental duplicon structure to provide direct confirmation of the hypothesis that the rearrangements are caused by non-allelic homologous recombination between the low copy repeats on chromosome 22 (LCR22s). To enable haplotype-specific assembly and rearrangement mapping in LCR22 regions, we combined fiber-FISH optical mapping with whole genome (ultra-)long read sequencing or rearrangement-specific long-range PCR on 24 duos (22q11.2DS patient and parent-of-origin) comprising several different LCR22-mediated rearrangements. Unexpectedly, we demonstrate that not only different paralogous segmental duplicon but also palindromic AT-rich repeats (PATRR) are driving 22q11.2 rearrangements. In addition, we show the existence of two different inversion polymorphisms preceding rearrangement, and somatic mosaicism. The existence of different recombination sites and mechanisms in paralogues and PATRRs which are copy number expanding in the human population are a likely explanation for the high 22q11.2DS incidence.

A 25-year odyssey of genomic technology advances and structural variant discovery

This perspective focuses on advances in genome technology over the last 25 years and their impact on germline variant discovery within the field of human genetics. The field has witnessed tremendous technological advances from microarrays to short-read sequencing and now long-read sequencing. Each technology has provided genome-wide access to different classes of human genetic variation. We are now on the verge of comprehensive variant detection of all forms of variation for the first time with a single assay. We predict that this transition will further transform our understanding of human health and biology and, more importantly, provide novel insights into the dynamic mutational processes shaping our genomes.

Mycoplasma DnaK increases DNA copy number variants in vivo

The human microbiota affects critical cellular functions, although the responsible mechanism(s) is still poorly understood. In this regard, we previously showed that Mycoplasma fermentans DnaK, an HSP70 chaperone protein, hampers the activity of important cellular proteins responsible for DNA integrity. Here, we describe a novel DnaK knock-in mouse model generated in our laboratory to study the effect of M. fermentans DnaK expression in vivo. By using an array-based comparative genomic hybridization assay, we demonstrate that exposure to DnaK was associated with a higher number of DNA copy number variants (CNVs) indicative of unbalanced chromosomal alterations, together with reduced fertility and a high rate of fetal abnormalities. Consistent with their implication in genetic disorders, one of these CNVs caused a homozygous Grid2 deletion, resulting in an aberrant ataxic phenotype that recapitulates the extensive biallelic deletion in the Grid2 gene classified in humans as autosomal recessive spinocerebellar ataxia 18. Our data highlight a connection between components of the human urogenital tract microbiota, namely Mycoplasmas, and genetic abnormalities in the form of DNA CNVs, with obvious relevant medical, diagnostic, and therapeutic implications.

Discovery, classification, evolution and diversity of Siglecs

Immunoglobulin (Ig) superfamily proteins play diverse roles in vertebrates, including regulation of cellular responses by sensing endogenous or exogenous ligands. Siglecs are a family of glycan-recognizing proteins belonging to the Ig superfamily (i.e., I-type lectins). Siglecs are expressed on various leukocyte types and are involved in diverse aspects of immunity, including the regulation of inflammatory responses, leukocyte proliferation, host-microbe interaction, and cancer immunity. Sialoadhesin/Siglec-1, CD22/Siglec-2, and myelin-associated glycoprotein/Siglec-4 were among the first to be characterized as members of the Siglec family, and along with Siglec-15, they are relatively well-conserved among tetrapods. Conversely, CD33/Siglec-3-related Siglecs (CD33rSiglecs, so named as they show high sequence similarity with CD33/Siglec-3) are encoded in a gene cluster with many interspecies variations and even intraspecies variations within some lineages such as humans. The rapid evolution of CD33rSiglecs expressed on leukocytes involved in innate immunity likely reflects the selective pressure by pathogens that interact and possibly exploit these Siglecs. Human Siglecs have several additional unique and/or polymorphic properties as compared with closely related great apes, changes possibly related to the loss of the sialic acid Neu5Gc, another distinctly human event in sialobiology. Multiple changes in human CD33rSiglecs compared to great apes include many examples of human-specific expression in non-immune cells, coinciding with human-specific diseases involving such cell types. Some Siglec gene polymorphisms have dual consequences-beneficial in a situation but detrimental in another. The association of human Siglec gene polymorphisms with several infectious and non-infectious diseases likely reflects the ongoing competition between the host and microbial pathogens.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Mutation of Proteolipid Protein 1 Gene: From Severe Hypomyelinating Leukodystrophy to Inherited Spastic Paraplegia

Pelizaeus–Merzbacher Disease (PMD) is an inherited leukodystrophy affecting the central nervous system (CNS)—a rare disorder that especially concerns males. Its estimated prevalence is 1.45–1.9 per 100,000 individuals in the general population. Patients affected by PMD exhibit a drastic reduction or absence of myelin sheaths in the white matter areas of the CNS. The Proteolipid Protein 1 (PLP1) gene encodes a transmembrane proteolipid protein. PLP1 is the major protein of myelin, and it plays a key role in the compaction, stabilization, and maintenance of myelin sheaths. Its function is predominant in oligodendrocyte development and axonal survival. Mutations in the PLP1 gene cause the development of a wide continuum spectrum of leukopathies from the most severe form of PMD for whom patients exhibit severe CNS hypomyelination to the relatively mild late-onset type 2 spastic paraplegia, leading to the concept of PLP1-related disorders. The genetic diversity and the biochemical complexity, along with other aspects of PMD, are discussed to reveal the obstacles that hinder the development of treatments. This review aims to provide a clinical and mechanistic overview of this spectrum of rare diseases.

Mutation spectrum of congenital heart disease in a consanguineous Turkish population

BACKGROUNDS

While many studies agree that consanguinity increases the rate of congenital heart disease (CHD), few genome analyses have been conducted with consanguineous CHD cohorts.

METHODS

We recruited 73 CHD probands from consanguineous families in Turkey and used whole-exome sequencing (WES) to identify genetic lesions in these patients.

RESULTS

On average, each patient had 6.95 rare damaging homozygous variants, 0.68 of which are loss-of-function (LoF) variants. Seven patients (9.6%) carried damaging homozygous variants in five causal CHD genes. Six of those patients exhibited laterality defects (six HTX and one D-TGA). Three additional patients (4.1%) harbored other types of CHD-associated genomic alterations, which overall explained 13.7% (10/73) of the cohort. The contribution from recessive variants in our cohort is higher than 1.8% reported from a cohort of 2871 CHD subjects where 5.6% of subjects met the criteria for consanguinity.

CONCLUSIONS

Our WES screen of a Turkish consanguineous population with structural CHD revealed its unique genetic architecture. Six of seven damaging homozygous variants in CHD causal genes occur in the setting of laterality defects implies a strong contribution from consanguinity to these defects specifically. Our study thus provided valuable information about the genetic landscape of CHD in consanguineous families in Turkey.

JAX-CNV: A Whole-genome Sequencing-based Algorithm for Copy Number Detection at Clinical Grade Level

We aimed to develop a whole-genome sequencing (WGS)-based copy number variant (CNV) calling algorithm with the potential of replacing chromosomal microarray assay (CMA) for clinical diagnosis. JAX-CNV is thus developed for CNV detection from WGS data. The performance of this CNV calling algorithm was evaluated in a blinded manner on 31 samples and compared to the 112 CNVs reported by clinically validated CMAs for these 31 samples. The result showed that JAX-CNV recalled 100% of these CNVs. Besides, JAX-CNV identified an average of 30 CNVs per individual that was an approximately seven-fold increase compared to calls of clinically validated CMAs. Experimental validation of 24 randomly selected CNVs showed one false positive, i.e., a false discovery rate (FDR) of 4.17%. A robustness test on lower-coverage data revealed a 100% sensitivity for CNVs larger than 300 kb (the current threshold for College of American Pathologists) down to 10× coverage. For CNVs larger than 50 kb, sensitivities were 100% for coverages deeper than 20×, 97% for 15×, and 95% for 10×. We developed a WGS-based CNV pipeline, including this newly developed CNV caller JAX-CNV, and found it capable of detecting CMA-reported CNVs at a sensitivity of 100% with about a FDR of 4%. We propose that JAX-CNV could be further examined in a multi-institutional study to justify the transition of first-tier genetic testing from CMAs to WGS. JAX-CNV is available at https://github.com/TheJacksonLaboratory/JAX-CNV.

Distinct sequence features underlie microdeletions and gross deletions in the human genome

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1 to 28,394,429 bp. We analyzed the DNA sequence within 1 kb of the breakpoint junctions and found that the frequencies of non‐B DNA‐forming repeats, GC‐content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ, and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) while significantly (χ², p < 2E−16) more microhomologies were identified flanking short deletions than long deletions (length >30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut‐off of 25–30 bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

Sex-specific variation in the genome-wide recombination rate

In most species that reproduce sexually, successful gametogenesis requires recombination during meiosis. The number and placement of crossovers (COs) vary among individuals, with females and males often presenting the most striking contrasts. Despite the recognition that the sexes recombine at different rates (heterochiasmy), existing data fail to answer the question of whether patterns of genetic variation in recombination rate are similar in the two sexes. To fill this gap, we measured the genome-wide recombination rate in both sexes from a panel of wild-derived inbred strains from multiple subspecies of house mice (Mus musculus) and from a few additional species of Mus. To directly compare recombination rates in females and males from the same genetic backgrounds, we applied established methods based on immunolocalization of recombination proteins to inbred strains. Our results reveal discordant patterns of genetic variation in the two sexes. Whereas male genome-wide recombination rates vary substantially among strains, female recombination rates measured in the same strains are more static. The direction of heterochiasmy varies within two subspecies, Mus musculus molossinus and Mus musculus musculus. The direction of sex differences in the length of the synaptonemal complex and CO positions is consistent across strains and does not track sex differences in genome-wide recombination rate. In males, contrasts between strains with high recombination rate and strains with low recombination rate suggest more recombination is associated with stronger CO interference and more double-strand breaks. The sex-specific patterns of genetic variation we report underscore the importance of incorporating sex differences into recombination research.

Modelling segmental duplications in the human genome

Background

Segmental duplications (SDs) are long DNA sequences that are repeated in a genome and have high sequence identity. In contrast to repetitive elements they are often unique and only sometimes have multiple copies in a genome. There are several well-studied mechanisms responsible for segmental duplications: non-allelic homologous recombination, non-homologous end joining and replication slippage. Such duplications play an important role in evolution, however, we do not have a full understanding of the dynamic properties of the duplication process.

Results

We study segmental duplications through a graph representation where nodes represent genomic regions and edges represent duplications between them. The resulting network (the SD network) is quite complex and has distinct features which allow us to make inference on the evolution of segmantal duplications. We come up with the network growth model that explains features of the SD network thus giving us insights on dynamics of segmental duplications in the human genome. Based on our analysis of genomes of other species the network growth model seems to be applicable for multiple mammalian genomes.

Conclusions

Our analysis suggests that duplication rates of genomic loci grow linearly with the number of copies of a duplicated region. Several scenarios explaining such a preferential duplication rates were suggested.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07789-7).

Double-strand breaks induce short-scale DNA replication and damage amplification in the fully grown mouse oocytes

Break-induced replication (BIR) is essential for the repair of DNA double-strand breaks (DSBs) with single ends. DSBs-induced microhomology-mediated BIR (mmBIR) and template-switching can increase the risk of complex genome rearrangement. In addition, DSBs can also induce the multi-invasion-mediated DSB amplification. The mmBIR-induced genomic rearrangement has been identified in cancer cells and patients with rare diseases. However, when and how mmBIR is initiated have not been fully and deeply studied. Furthermore, it is not well understood about the conditions for initiation of multi-invasion-mediated DSB amplification. In the G2 phase oocyte of mouse, we identified a type of short-scale BIR (ssBIR) using the DNA replication indicator 5-ethynyl-2'-deoxyuridine (EdU). These ssBIRs could only be induced in the fully grown oocytes but not the growing oocytes. If the DSB oocytes were treated with Rad51 or Chek1/2 inhibitors, both EdU signals and DSB marker γH2A.X foci would decrease. In addition, the DNA polymerase inhibitor Aphidicolin could inhibit the ssBIR and another inhibitor ddATP could reduce the number of γH2A.X foci in the DSB oocytes. In conclusion, our results showed that DNA DSBs in the fully grown oocytes can initiate ssBIR and be amplified by Rad51 or DNA replication.

Donor-derived Cell-free DNA: Advancing a Novel Assay to New Heights in Renal Transplantation

Despite advances in transplant immunosuppression, long-term renal allograft outcomes remain suboptimal because of the occurrence of rejection, recurrent disease, and interstitial fibrosis with tubular atrophy. This is largely due to limitations in our understanding of allogeneic processes coupled with inadequate surveillance strategies. The concept of donor-derived cell-free DNA as a signal of allograft stress has therefore rapidly been adopted as a noninvasive monitoring tool. Refining it for effective clinical use, however, remains an ongoing effort. Furthermore, its potential to unravel new insights in alloimmunity through novel molecular techniques is yet to be realized. This review herein summarizes current knowledge and active endeavors to optimize cell-free DNA-based diagnostic techniques for clinical use in kidney transplantation. In addition, the integration of DNA methylation and microRNA may unveil new epigenetic signatures of allograft health and is also explored in this report. Directing research initiatives toward these aspirations will not only improve diagnostic precision but may foster new paradigms in transplant immunobiology.

Genomic Variation, Evolvability, and the Paradox of Mental Illness

Twentieth-century genetics was hard put to explain the irregular behavior of neuropsychiatric disorders. Autism and schizophrenia defy a principle of natural selection; they are highly heritable but associated with low reproductive success. Nevertheless, they persist. The genetic origins of such conditions are confounded by the problem of variable expression, that is, when a given genetic aberration can lead to any one of several distinct disorders. Also, autism and schizophrenia occur on a spectrum of severity, from mild and subclinical cases to the overt and disabling. Such irregularities reflect the problem of missing heritability; although hundreds of genes may be associated with autism or schizophrenia, together they account for only a small proportion of cases. Techniques for higher resolution, genomewide analysis have begun to illuminate the irregular and unpredictable behavior of the human genome. Thus, the origins of neuropsychiatric disorders in particular and complex disease in general have been illuminated. The human genome is characterized by a high degree of structural and behavioral variability: DNA content variation, epistasis, stochasticity in gene expression, and epigenetic changes. These elements have grown more complex as evolution scaled the phylogenetic tree. They are especially pertinent to brain development and function. Genomic variability is a window on the origins of complex disease, neuropsychiatric disorders, and neurodevelopmental disorders in particular. Genomic variability, as it happens, is also the fuel of evolvability. The genomic events that presided over the evolution of the primate and hominid lineages are over-represented in patients with autism and schizophrenia, as well as intellectual disability and epilepsy. That the special qualities of the human genome that drove evolution might, in some way, contribute to neuropsychiatric disorders is a matter of no little interest.

The prenatal diagnosis and genetic counseling of chromosomal micro-duplication on 10q24.3 in a fetus: A case report and a brief review of the literature

RATIONALE

Split-hand/split-foot malformation (SHFM), also known as ectrodactyly, is a congenital limb malformation affecting the central rays of the autopod extending to syndactyly, median clefts of the hands and feet, aplasia/hypoplasia of phalanges, metacarpals and metatarsals. Duplication of this 10q24 region is associated with SHFM3. While the clinical and genetic heterogeneity of SHFM makes the prenatal diagnosis and genetic counseling more challenging and difficult.

PATIENT CONCERNS

A physically normal pregnant woman had a systemic ultrasound at the second trimester, only identified the deformity of both hands and feet on the fetus.

DIAGNOSES

The fetus was diagnosed as sporadic SHFM3.

INTERVENTIONS

After seeking advice from genetic counseling, she decided to terminate the pregnancy. The induction of infant was done after appearance of bipedal clefts, lobster-claw appearance and partial loss of phalanges and metacarpals, leaving behind 2nd finger in the left hand and the 5th in the right hand. Furthermore, collection of umbilical cord is recommended to this fetus for genome-wide detection.

OUTCOMES

An outcome of the gene detection from abortion shows that there is variation in copy number in genome of chromosome 1 and chromosome 10.

LESSONS

This case study confirms an association between SHFM3 and chromosomal micro-duplication on 10q24.3, and the extension of clinical spectrum of SHFM3. It also proposes some prenatal diagnosis and genetic counseling to help in planning and management in affected pregnancy. This will reduce the congenital and development abnormalities in birth rate, as well as relive the economic, psychological, and physical burden to the affected families.

Consequences of 22q11.2 Microdeletion on the Genome, Individual and Population Levels

Chromosomal 22q11.2 deletion syndrome (22q11.2DS) (ORPHA: 567) caused by microdeletion in chromosome 22 is the most common chromosomal microdeletion disorder in humans. Despite the same change on the genome level, like in the case of monozygotic twins, phenotypes are expressed differently in 22q11.2 deletion individuals. The rest of the genome, as well as epigenome and environmental factors, are not without influence on the variability of phenotypes. The penetrance seems to be more genotype specific than deleted locus specific. The transcript levels of deleted genes are not usually reduced by 50% as assumed due to haploinsufficiency. 22q11.2DS is often an undiagnosed condition, as each patient may have a different set out of 180 possible clinical manifestations. Diverse dysmorphic traits are present in patients from different ethnicities, which makes diagnosis even more difficult. 22q11.2 deletion syndrome serves as an example of a genetic syndrome that is not easy to manage at all stages: diagnosis, consulting and dealing with.

The repair of endo/exogenous DNA double-strand breaks and its effects on meiotic chromosome segregation in oocytes

Germ cell-derived genomic structure variants not only drive the evolution of species but also induce developmental defects in offspring. The genomic structure variants have different types, but most of them are originated from DNA double-strand breaks (DSBs). It is still not well known whether DNA DSBs exist in adult mammalian oocytes and how the growing and fully grown oocytes repair their DNA DSBs induced by endogenous or exogenous factors. In this study, we detected the endogenous DNA DSBs in the growing and fully grown mouse oocytes and found that the DNA DSBs mainly localized at the centromere-adjacent regions, which are also copy number variation hotspots. When the exogenous DNA DSBs were introduced by Etoposide, we found that Rad51-mediated homologous recombination (HR) was used to repair the broken DNA. However, the HR repair caused the chromatin intertwined and impaired the homologous chromosome segregation in oocytes. Although we had not detected the indication about HR repair of endogenous centromere-adjacent DNA DSBs, we found that Rad52 and RNA:DNA hybrids colocalized with these DNA DSBs, indicating that a Rad52-dependent DNA repair might exist in oocytes. In summary, our results not only demonstrated an association between endogenous DNA DSBs with genomic structure variants but also revealed one specific DNA DSB repair manner in oocytes.

The 22q11 low copy repeats are characterized by unprecedented size and structural variability

Low copy repeats (LCRs) are recognized as a significant source of genomic instability, driving genome variability and evolution. The Chromosome 22 LCRs (LCR22s) mediate nonallelic homologous recombination (NAHR) leading to the 22q11 deletion syndrome (22q11DS). However, LCR22s are among the most complex regions in the genome, and their structure remains unresolved. The difficulty in generating accurate maps of LCR22s has also hindered localization of the deletion end points in 22q11DS patients. Using fiber FISH and Bionano optical mapping, we assembled LCR22 alleles in 187 cell lines. Our analysis uncovered an unprecedented level of variation in LCR22s, including LCR22A alleles ranging in size from 250 to 2000 kb. Further, the incidence of various LCR22 alleles varied within different populations. Additionally, the analysis of LCR22s in 22q11DS patients and their parents enabled further refinement of the rearrangement site within LCR22A and -D, which flank the 22q11 deletion. The NAHR site was localized to a 160-kb paralog shared between the LCR22A and -D in seven 22q11DS patients. Thus, we present the most comprehensive map of LCR22 variation to date. This will greatly facilitate the investigation of the role of LCR variation as a driver of 22q11 rearrangements and the phenotypic variability among 22q11DS patients.

A Deep Learning Approach for Detecting Copy Number Variation in Next-Generation Sequencing Data

Copy number variants (CNV) are associated with phenotypic variation in several species. However, properly detecting changes in copy numbers of sequences remains a difficult problem, especially in lower quality or lower coverage next-generation sequencing data. Here, inspired by recent applications of machine learning in genomics, we describe a method to detect duplications and deletions in short-read sequencing data. In low coverage data, machine learning appears to be more powerful in the detection of CNVs than the gold-standard methods of coverage estimation alone, and of equal power in high coverage data. We also demonstrate how replicating training sets allows a more precise detection of CNVs, even identifying novel CNVs in two genomes previously surveyed thoroughly for CNVs using long read data.

Genome-wide recombination map construction from single individuals using linked-read sequencing

Meiotic recombination rates vary across the genome, often involving localized crossover hotspots and coldspots. Studying the molecular basis and mechanisms underlying this variation has been challenging due to the high cost and effort required to construct individualized genome-wide maps of recombination crossovers. Here we introduce a new method, called ReMIX, to detect crossovers from gamete DNA of a single individual using Illumina sequencing of 10X Genomics linked-read libraries. ReMIX reconstructs haplotypes and identifies the valuable rare molecules spanning crossover breakpoints, allowing quantification of the genomic location and intensity of meiotic recombination. Using a single mouse and stickleback fish, we demonstrate how ReMIX faithfully recovers recombination hotspots and landscapes that have previously been built using hundreds of offspring. ReMIX provides a high-resolution, high-throughput, and low-cost approach to quantify recombination variation across the genome, providing an exciting opportunity to study recombination among multiple individuals in diverse organisms.

Variation of recombination rates within genomes has important implications in genetics and evolution. Here, the authors develop a method for building genome-wide recombination maps from single individuals using linked-read sequencing data, and report its application in mouse and stickleback fish.

A first genetic portrait of synaptonemal complex variation

The synaptonemal complex (SC) is a proteinaceous scaffold required for synapsis and recombination between homologous chromosomes during meiosis. Although the SC has been linked to differences in genome-wide crossover rates, the genetic basis of standing variation in SC structure remains unknown. To investigate the possibility that recombination evolves through changes to the SC, we characterized the genetic architecture of SC divergence on two evolutionary timescales. Applying a novel digital image analysis technique to spermatocyte spreads, we measured total SC length in 9,532 spermatocytes from recombinant offspring of wild-derived mouse strains with differences in this fundamental meiotic trait. Using this large dataset, we identified the first known genomic regions involved in the evolution of SC length. Distinct loci affect total SC length divergence between and within subspecies, with the X chromosome contributing to both. Joint genetic analysis of MLH1 foci—immunofluorescent markers of crossovers—from the same spermatocytes revealed that two of the identified loci also confer differences in the genome-wide recombination rate. Causal mediation analysis suggested that one pleiotropic locus acts early in meiosis to designate crossovers prior to SC assembly, whereas a second locus primarily shapes crossover number through its effect on SC length. One genomic interval shapes the relationship between SC length and recombination rate, likely modulating the strength of crossover interference. Our findings pinpoint SC formation as a key step in the evolution of recombination and demonstrate the power of genetic mapping on standing variation in the context of the recombination pathway.

During the first stages of meiosis, the chromosome axes are organized along a protein scaffold in preparation for recombination and their subsequent segregation. This scaffold, known as the synaptonemal complex (SC), is critical for the regular progression of recombination. A complex relationship exists between the organization of the SC, the frequency of recombination, and the likelihood of improper chromosome segregation. In this study, we investigate the genetics of synaptonemal complex variation in the house mouse and connect it with variation in the rate of recombination. We found five loci and several compelling candidate genes responsible for the evolution of synaptonemal complex length within and between mouse subspecies. Several of these loci also affect recombination rate, and our joint analyses of the phenotypes suggest an order by which their effects manifest within the recombination pathway. Our results show that evolution of SC length is crucial to recombination rate divergence. Our work here also demonstrates that genetic analysis of additional meiotic phenotypes can help explain the evolution of recombination, a fundamental evolutionary force.

Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode

Soybean cyst nematode (SCN) is the most devastating plant-parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1-b locus, while 'Peking' requires rhg1-a and Rhg4 for SCN resistance. In the present study, whole genome re-sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG-Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and 'Cloud'. At least 5.6 copies of the PI88788-type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a 'Peking'-type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad-based resistance to SCN and providing new insight into epistasis, haplotype-compatibility, CNV, promoter variation and its impact on broad-based disease resistance in plants.

Conservation of the genome-wide recombination rate in white-footed mice

Despite being linked to the fundamental processes of chromosome segregation and offspring diversification, meiotic recombination rates vary within and between species. Recent years have seen progress in quantifying recombination rate evolution across multiple temporal and genomic scales. Nevertheless, the level of variation in recombination rate within wild populations-a key determinant of evolution in this trait-remains poorly documented on the genomic scale. To address this notable gap, we used immunofluorescent cytology to quantify genome-wide recombination rates in males from a wild population of the white-footed mouse, Peromyscus leucopus. For comparison, we measured recombination rates in a second population of male P. leucopus raised in the laboratory and in male deer mice from the subspecies Peromyscus maniculatus bairdii. Although we found differences between individuals in the genome-wide recombination rate, levels of variation were low-within populations, between populations, and between species. Quantification of synaptonemal complex length and crossover positions along chromosome 1 using a novel automated approach also revealed conservation in broad-scale crossover patterning, including strong crossover interference. We propose stabilizing selection targeting recombination or correlated processes as the explanation for these patterns.

Evolutionary Story of the Low/Medium-Affinity IgG Fc Receptor Gene Cluster

Low/intermediate affinity Fc-gamma receptors (FcγR) are crucial for the recognition of immune complexes and IgG-sensitized microorganisms by phagocytic and cytotoxic effector cells. In all mammalian species studied so far, their genes are clustered in a single locus. However, this locus differs between humans and mice, both in the number of genes and the structure/function of the encoded receptors. We show that murine fcgr3 evolved through several steps into FCGR2A, its ortholog, which is specific to primates. One of these steps was the insertion of a retroviral element bringing a new intracellular exon comprising a non-canonical ITAM motif. We also show that the fcgr3-hspa6-fcgr4-fcgr2b module in mammals that has evolved in a FCGR2A-HSPA6-FCGR4-FCGR2B module in primates, was subsequently duplicated in apes through a Non-Allelic Homologous Recombination (NAHR), giving birth to FCGR2C, a hybrid gene between FCGR2B and FCGR2A. The FCGR4 duplication, which occurred simultaneously, eventually resulted in the emergence of FCGR3B, while FCGR3A remained the true FCGR4 ortholog. FCGR2C and FCGR3B, markers of this NAHR, are present in gorillas and chimpanzees, whereas they are absent in orangutans and more distant primates, such as gibbons and macaques. These data need to be taken into account when testing IgG-based therapies in animal species.

A PHLDB1 variant associated with the nonfunctional pituitary adenoma

PURPOSE

Previous studies have revealed that PHLDB1 single-nucleotide polymorphisms (SNPs) are associated with glioma risk. Nonetheless, the association between PHLDB1 SNPs and the risk of pituitary adenoma has not been studied. The present study evaluated the association of PHLDB1 SNPs with the risk of pituitary adenomas.

METHODS

We genotyped 27 PHLDB1 tagging and exon SNPs in a case-control study that included 148 patients who got a diagnosis of nonfunctional pituitary adenoma (NFPA) and 375 normal controls within the Korean population. Statistical analyses of the association between PHLDB1 SNPs and the NFPA risk were conducted using logistic regression.

RESULTS

We detected an association between a PHLDB1 SNP and the risk of NFPA in the Korean population. Rs67307131 in intron 2 was significantly associated with NFPA (odds ratio [OR] = 2.15, 95% confidence interval [CI] 1.44-3.20; P = 0.0002 in the dominant model). In the referent analysis, a higher OR and stronger association (lower P value) were observed among patients with the "C/T" genotype (OR = 2.39, 95% CI 1.60-3.58; P = 0.00002). In a functional analysis with a SNP annotation tool, this SNP was predicted to be a CpG site and copy number variant; these properties are associated with susceptibility to diseases.

CONCLUSIONS

Our findings suggest that genetic variation of PHLDB1 may be associated with the risk of NFPA. This is the first report of an association between PHLDB1 variants and NFPA. Further research is needed to confirm the impact of this SNP on NFPA susceptibility.

Pelizaeus-Merzbacher Disease: Molecular and Cellular Pathologies and Associated Phenotypes

Pelizaeus-Merzbacher disease (PMD) represents a group of disorders known as hypomyelinating leukodystrophies, which are characterized by abnormal development and maintenance of myelin in the central nervous system. PMD is caused by different types of mutations in the proteolipid protein 1 (PLP1) gene, which encodes a major myelin membrane lipoprotein. These mutations in the PLP1 gene result in distinct cellular and molecular pathologies and a spectrum of clinical phenotypes. In this chapter, I discuss the historical aspects and current understanding of the mechanisms underlying how different PLP1 mutations disrupt the normal process of myelination and result in PMD and other disorders.

A Genomic Region Containing REC8 and RNF212B Is Associated with Individual Recombination Rate Variation in a Wild Population of Red Deer (Cervus elaphus)

Recombination is a fundamental feature of sexual reproduction, ensuring proper disjunction, preventing mutation accumulation and generating new allelic combinations upon which selection can act. However it is also mutagenic, and breaks up favorable allelic combinations previously built up by selection. Identifying the genetic drivers of recombination rate variation is a key step in understanding the causes and consequences of this variation, how loci associated with recombination are evolving and how they affect the potential of a population to respond to selection. However, to date, few studies have examined the genetic architecture of recombination rate variation in natural populations. Here, we use pedigree data from ∼ 2,600 individuals genotyped at ∼ 38,000 SNPs to investigate the genetic architecture of individual autosomal recombination rate in a wild population of red deer (Cervus elaphus). Female red deer exhibited a higher mean and phenotypic variance in autosomal crossover counts (ACC). Animal models fitting genomic relatedness matrices showed that ACC was heritable in females ([Formula: see text] = 0.12) but not in males. A regional heritability mapping approach showed that almost all heritable variation in female ACC was explained by a genomic region on deer linkage group 12 containing the candidate loci REC8 and RNF212B, with an additional region on linkage group 32 containing TOP2B approaching genome-wide significance. The REC8/RNF212B region and its paralogue RNF212 have been associated with recombination in cattle, mice, humans and sheep. Our findings suggest that mammalian recombination rates have a relatively conserved genetic architecture in both domesticated and wild systems, and provide a foundation for understanding the association between recombination loci and individual fitness within this population.

Molecular characterization of HDAC8 deletions in individuals with atypical Cornelia de Lange syndrome

Cornelia de Lange syndrome (CdLS) is a rare neurodevelopmental syndrome for which mutations in five causative genes that encode (SMC1A, SMC3, RAD21) or regulate (NIPBL, HDAC8) the cohesin complex, account for ~70% of cases. Herein we report on four female Subjects who were found to carry novel intragenic deletions in HDAC8. In one case, the deletion was found in mosaic state and it was determined to be present in ~38% of blood lymphocytes and in nearly all cells of a buccal sample. All deletions, for which parental blood samples were available, were shown to have arisen de novo. X-chromosome inactivation studies demonstrated marked skewing, suggesting strong selection against the mutated HDAC8 allele. Based on an investigation of the deletion breakpoints, we hypothesize that microhomology-mediated replicative mechanisms may be implicated in the formation of some of these rearrangements. This study broadens the mutational spectrum of HDAC8, provides the first description of a causative HDAC8 somatic mutation and increases the knowledge on possible mutational mechanisms underlying copy number variations in HDAC8. Moreover our findings highlight the clinical utility of considering copy number analysis in HDAC8 as well as the analysis on DNA from more than one tissue as an indispensable part of the routine molecular diagnosis of individuals with CdLS or CdLS-overlapping features.

Pyruvate Kinase and Fcγ Receptor Gene Copy Numbers Associated With Malaria Phenotypes

Genetic factors are associated with susceptibility to many infectious diseases and may be determinants of clinical progression. Gene copy number variation (CNV) has been shown to be associated with phenotypes of numerous diseases, including malaria. We quantified gene copy numbers of the pyruvate kinase, liver, and red blood cell (PKLR) gene as well as of the Fcγ receptor 2A and Fcγ receptor 2C (FCGR2A, FCGR2C) and Fcγ receptor 3 (FCGR3) genes using real-time quantitative polymerase chain reaction (RT-qPCR) assays in Gabonese children with severe (n = 184) or and mild (n = 189) malaria and in healthy Gabonese and white individuals (n = 76 each). The means of PKLR, FCGR2A, FCGR2C, and FCGR3 copy numbers were significantly higher among children with severe malaria compared to those with mild malaria (P < .002), indicating that a surplus of copies of those genes is significantly associated with malaria severity. Copy numbers of the FCGR2A and FCGR2C genes were significantly lower (P = .005) in Gabonese individuals compared with white individuals. In conclusion, CNV of the PKLR, FCGR2A, FCGR2C, and FCGR3 genes is associated with malaria severity, and our results provide evidence for a role of CNV in host responses to malaria.

Genetics of Genome-Wide Recombination Rate Evolution in Mice from an Isolated Island

Recombination rate is a heritable quantitative trait that evolves despite the fundamentally conserved role that recombination plays in meiosis. Differences in recombination rate can alter the landscape of the genome and the genetic diversity of populations. Yet our understanding of the genetic basis of recombination rate evolution in nature remains limited. We used wild house mice (Mus musculus domesticus) from Gough Island (GI), which diverged recently from their mainland counterparts, to characterize the genetics of recombination rate evolution. We quantified genome-wide autosomal recombination rates by immunofluorescence cytology in spermatocytes from 240 F2 males generated from intercrosses between GI-derived mice and the wild-derived inbred strain WSB/EiJ. We identified four quantitative trait loci (QTL) responsible for inter-F2 variation in this trait, the strongest of which had effects that opposed the direction of the parental trait differences. Candidate genes and mutations for these QTL were identified by overlapping the detected intervals with whole-genome sequencing data and publicly available transcriptomic profiles from spermatocytes. Combined with existing studies, our findings suggest that genome-wide recombination rate divergence is not directional and its evolution within and between subspecies proceeds from distinct genetic loci.

Next generation sequencing: implications in personalized medicine and pharmacogenomics

A breakthrough in next generation sequencing (NGS) in the last decade provided an unprecedented opportunity to investigate genetic variations in humans and their roles in health and disease. NGS offers regional genomic sequencing such as whole exome sequencing of coding regions of all genes, as well as whole genome sequencing. RNA-seq offers sequencing of the entire transcriptome and ChIP-seq allows for sequencing the epigenetic architecture of the genome. Identifying genetic variations in individuals can be used to predict disease risk, with the potential to halt or retard disease progression. NGS can also be used to predict the response to or adverse effects of drugs or to calculate appropriate drug dosage. Such a personalized medicine also provides the possibility to treat diseases based on the genetic makeup of the patient. Here, we review the basics of NGS technologies and their application in human diseases to foster human healthcare and personalized medicine.

Inherited human IRAK-1 deficiency selectively impairs TLR signaling in fibroblasts

Most members of the Toll-like receptor (TLR) and interleukin-1 receptor (IL-1R) families transduce signals via a canonical pathway involving the MyD88 adapter and the interleukin-1 receptor-associated kinase (IRAK) complex. This complex contains four molecules, including at least two (IRAK-1 and IRAK-4) active kinases. In mice and humans, deficiencies of IRAK-4 or MyD88 abolish most TLR (except for TLR3 and some TLR4) and IL-1R signaling in both leukocytes and fibroblasts. TLR and IL-1R responses are weak but not abolished in mice lacking IRAK-1, whereas the role of IRAK-1 in humans remains unclear. We describe here a boy with X-linked MECP2 deficiency-related syndrome due to a large de novo Xq28 chromosomal deletion encompassing both MECP2 and IRAK1 Like many boys with MECP2 null mutations, this child died very early, at the age of 7 mo. Unlike most IRAK-4- or MyD88-deficient patients, he did not suffer from invasive bacterial diseases during his short life. The IRAK-1 protein was completely absent from the patient's fibroblasts, which responded very poorly to all TLR2/6 (PAM₂CSK₄, LTA, FSL-1), TLR1/2 (PAM₃CSK₄), and TLR4 (LPS, MPLA) agonists tested but had almost unimpaired responses to IL-1β. By contrast, the patient's peripheral blood mononuclear cells responded normally to all TLR1/2, TLR2/6, TLR4, TLR7, and TLR8 (R848) agonists tested, and to IL-1β. The death of this child precluded long-term evaluations of the clinical consequences of inherited IRAK-1 deficiency. However, these findings suggest that human IRAK-1 is essential downstream from TLRs but not IL-1Rs in fibroblasts, whereas it plays a redundant role downstream from both TLRs and IL-1Rs in leukocytes.

Doublet-Mediated DNA Rearrangement-A Novel and Potentially Underestimated Mechanism for the Formation of Recurrent Pathogenic Deletions

Deletions and duplications of genomic DNA contribute to evolution, phenotypic diversity, and human disease. The underlying mechanisms are incompletely understood. We identified deletions of exon 10 of the SPAST gene in two unrelated families with hereditary spastic paraplegia. We excluded a founder event, but observed that the breakpoints map to identical repeat regions. These regions likely represent an intragenic "doublet," that is, an enigmatic class of local duplications. The fusion sequences for both deletions are compatible with recombination-based as well as with replication-based mechanisms. Searching the literature, we identified a partial SLC24A4 deletion that involved two copies of another doublet, and was likely formed in an analogous way. Comparing the SPAST and the SLC24A4 doublets with doublets identified previously suggested that many additional doublets have a high potential for triggering rearrangements. Considering that doublets are still being formed in the human genome, and that they likely create high local instability, we suggest that a two-step mechanism consisting of doublet generation and subsequent doublet-mediated deletion/duplication may underlie certain copy-number changes for which other mechanisms are currently assumed. Further studies are necessary to delineate the significance of the thus-far understudied doublets for the formation of copy-number variation.

Population Parameters Underlying an Ongoing Soft Sweep in Southeast Asian Malaria Parasites

Multiple kelch13 alleles conferring artemisinin resistance (ART-R) are currently spreading through Southeast Asian malaria parasite populations, providing a unique opportunity to observe an ongoing soft selective sweep, investigate why resistance alleles have evolved multiple times and determine fundamental population genetic parameters for Plasmodium. We sequenced kelch13 (n = 1,876), genotyped 75 flanking SNPs, and measured clearance rate (n = 3,552) in parasite infections from Western Thailand (2001–2014). We describe 32 independent coding mutations including common mutations outside the kelch13 propeller associated with significant reductions in clearance rate. Mutations were first observed in 2003 and rose to 90% by 2014, consistent with a selection coefficient of ∼0.079. ART-R allele diversity rose until 2012 and then dropped as one allele (C580Y) spread to high frequency. The frequency with which adaptive alleles arise is determined by the rate of mutation and the population size. Two factors drive this soft sweep: (1) multiple kelch13 amino-acid mutations confer resistance providing a large mutational target—we estimate the target is 87–163 bp. (2) The population mutation parameter (Θ = 2N_eμ) can be estimated from the frequency distribution of ART-R alleles and is ∼5.69, suggesting that short term effective population size is 88 thousand to 1.2 million. This is 52–705 times greater than N_e estimated from fluctuation in allele frequencies, suggesting that we have previously underestimated the capacity for adaptive evolution in Plasmodium. Our central conclusions are that retrospective studies may underestimate the complexity of selective events and the N_e relevant for adaptation for malaria is considerably higher than previously estimated.

Effect of manipulating recombination rates on response to selection in livestock breeding programs

Background

In this work, we performed simulations to explore the potential of manipulating recombination rates to increase response to selection in livestock breeding programs.

Methods

We carried out ten replicates of several scenarios that followed a common overall structure but differed in the average rate of recombination along the genome (expressed as the length of a chromosome in Morgan), the genetic architecture of the trait under selection, and the selection intensity under truncation selection (expressed as the proportion of males selected). Recombination rates were defined by simulating nine different chromosome lengths: 0.10, 0.25, 0.50, 1, 2, 5, 10, 15 and 20 Morgan, respectively. One Morgan was considered to be the typical chromosome length for current livestock species. The genetic architecture was defined by the number of quantitative trait variants (QTV) that affected the trait under selection. Either a large (10,000) or a small (1000 or 500) number of QTV was simulated. Finally, the proportions of males selected under truncation selection as sires for the next generation were equal to 1.2, 2.4, 5, or 10 %.

Results

Increasing recombination rate increased the overall response to selection and decreased the loss of genetic variance. The difference in cumulative response between low and high recombination rates increased over generations. At low recombination rates, cumulative response to selection tended to asymptote sooner and the genetic variance was completely eroded. If the trait under selection was affected by few QTV, differences between low and high recombination rates still existed, but the selection limit was reached at all rates of recombination.

Conclusions

Higher recombination rates can enhance the efficiency of breeding programs to turn genetic variation into response to selection. However, to increase response to selection significantly, the recombination rate would need to be increased 10- or 20-fold. The biological feasibility and consequences of such large increases in recombination rates are unknown.

Electronic supplementary material

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

Conserved Genetic Architecture Underlying Individual Recombination Rate Variation in a Wild Population of Soay Sheep (Ovis aries)

Meiotic recombination breaks down linkage disequilibrium (LD) and forms new haplotypes, meaning that it is an important driver of diversity in eukaryotic genomes. Understanding the causes of variation in recombination rate is important in interpreting and predicting evolutionary phenomena and in understanding the potential of a population to respond to selection. However, despite attention in model systems, there remains little data on how recombination rate varies at the individual level in natural populations. Here we used extensive pedigree and high-density SNP information in a wild population of Soay sheep (Ovis aries) to investigate the genetic architecture of individual autosomal recombination rates. Individual rates were high relative to other mammal systems and were higher in males than in females (autosomal map lengths of 3748 and 2860 cM, respectively). The heritability of autosomal recombination rate was low but significant in both sexes (h² = 0.16 and 0.12 in females and males, respectively). In females, 46.7% of the heritable variation was explained by a subtelomeric region on chromosome 6; a genome-wide association study showed the strongest associations at locus RNF212, with further associations observed at a nearby ∼374-kb region of complete LD containing three additional candidate loci, CPLX1, GAK, and PCGF3. A second region on chromosome 7 containing REC8 and RNF212B explained 26.2% of the heritable variation in recombination rate in both sexes. Comparative analyses with 40 other sheep breeds showed that haplotypes associated with recombination rates are both old and globally distributed. Both regions have been implicated in rate variation in mice, cattle, and humans, suggesting a common genetic architecture of recombination rate variation in mammals.

Type 2 diabetes mellitus disease risk genes identified by genome wide copy number variation scan in normal populations

AIMS

To identify the role of copy number variations (CNVs) on disease risk genes and its effect on disease phenotypes in type 2 diabetes mellitus (T2DM) in 12 random populations using high throughput arrays.

METHODS

CNV analysis was carried out on a total of 1715 individuals from 12 populations, from ArrayExpress Archive of the European Bioinformatics Institute along with our subjects using Affymetrix Genome Wide SNP 6.0 array. CNV effect on T2DM genes were analyzed using several bioinformatics tools and a molecular protein interaction network was constructed to identify the disease mechanism altered by the CNVs.

RESULTS

Analysis showed 34.4% of the total population to be under CNV burden for T2DM, with 83 disease causal and associated genes being under CNV influence. Hotspots were identified on chromosomes 22, 12, 6, 19 and 11.Overlap studies with case cohorts revealed significant disease risk genes such as EGFR, E2F1, PPP1R3A, HLA and TSPAN8.

CONCLUSIONS

CNVs play a significant role in predisposing T2DM in normal cohorts and contribute to the phenotypic effects. Thus, CNVs should be considered as one of the major contributors in predisposition of the disease.

A monoallelic-to-biallelic T-cell transcriptional switch regulates GATA3 abundance

Ku et al. show that loss of one Gata3 allele leads to diminished expansion of immature T cells as well as aberrant induction of myeloid transcription factor PU.1. Gata3 is monoallelically expressed in hematopoietic stem cells and early T-cell progenitors. Half of the developing cells switch to biallelic Gata3 transcription abruptly at midthymopoiesis.

Protein abundance must be precisely regulated throughout life, and nowhere is the stringency of this requirement more evident than during T-cell development: A twofold increase in the abundance of transcription factor GATA3 results in thymic lymphoma, while reduced GATA3 leads to diminished T-cell production. GATA3 haploinsufficiency also causes human HDR (hypoparathyroidism, deafness, and renal dysplasia) syndrome, often accompanied by immunodeficiency. Here we show that loss of one Gata3 allele leads to diminished expansion (and compromised development) of immature T cells as well as aberrant induction of myeloid transcription factor PU.1. This effect is at least in part mediated transcriptionally: We discovered that Gata3 is monoallelically expressed in a parent of origin-independent manner in hematopoietic stem cells and early T-cell progenitors. Curiously, half of the developing cells switch to biallelic Gata3 transcription abruptly at midthymopoiesis. We show that the monoallelic-to-biallelic transcriptional switch is stably maintained and therefore is not a stochastic phenomenon. This unique mechanism, if adopted by other regulatory genes, may provide new biological insights into the rather prevalent phenomenon of monoallelic expression of autosomal genes as well as into the variably penetrant pathophysiological spectrum of phenotypes observed in many human syndromes that are due to haploinsufficiency of the affected gene.

Meiosis evolves: adaptation to external and internal environments

306 I. 306 II. 307 III. 312 IV. 317 V. 318 319 References 319 SUMMARY: Meiosis is essential for the fertility of most eukaryotes and its structures and progression are conserved across kingdoms. Yet many of its core proteins show evidence of rapid or adaptive evolution. What drives the evolution of meiosis proteins? How can constrained meiotic processes be modified in response to challenges without compromising their essential functions? In surveying the literature, we found evidence of two especially potent challenges to meiotic chromosome segregation that probably necessitate adaptive evolutionary responses: whole-genome duplication and abiotic environment, especially temperature. Evolutionary solutions to both kinds of challenge are likely to involve modification of homologous recombination and synapsis, probably via adjustments of core structural components important in meiosis I. Synthesizing these findings with broader patterns of meiosis gene evolution suggests that the structural components of meiosis coevolve as adaptive modules that may change in primary sequence and function while maintaining three-dimensional structures and protein interactions. The often sharp divergence of these genes among species probably reflects periodic modification of entire multiprotein complexes driven by genomic or environmental changes. We suggest that the pressures that cause meiosis to evolve to maintain fertility may cause pleiotropic alterations of global crossover rates. We highlight several important areas for future research.

Early evolutionary history and genomic features of gene duplicates in the human genome

BACKGROUND

Human gene duplicates have been the focus of intense research since the development of array-based and targeted next-generation sequencing approaches in the last decade. These studies have primarily concentrated on determining the extant copy-number variation from a population-genomic perspective but lack a robust evolutionary framework to elucidate the early structural and genomic characteristics of gene duplicates at emergence and their subsequent evolution with increasing age.

RESULTS

We analyzed 184 gene duplicate pairs comprising small gene families in the draft human genome with 10% or less synonymous sequence divergence. Human gene duplicates primarily originate from DNA-mediated events, taking up genomic residence as intrachromosomal copies in direct or inverse orientation. The distribution of paralogs on autosomes follows random expectations in contrast to their significant enrichment on the sex chromosomes. Furthermore, human gene duplicates exhibit a skewed gradient of distribution along the chromosomal length with significant clustering in pericentromeric regions. Surprisingly, despite the large average length of human genes, the majority of extant duplicates (83%) are complete duplicates, wherein the entire ORF of the ancestral copy was duplicated. The preponderance of complete duplicates is in accord with an extremely large median duplication span of 36 kb, which enhances the probability of capturing ancestral ORFs in their entirety. With increasing evolutionary age, human paralogs exhibit declines in (i) the frequency of intrachromosomal paralogs, and (ii) the proportion of complete duplicates. These changes may reflect lower survival rates of certain classes of duplicates and/or the role of purifying selection. Duplications arising from RNA-mediated events comprise a small fraction (11.4%) of all human paralogs and are more numerous in older evolutionary cohorts of duplicates.

CONCLUSIONS

The degree of structural resemblance, genomic location and duplication span appear to influence the long-term maintenance of paralogs in the human genome. The median duplication span in the human genome far exceeds that in C. elegans and yeast and likely contributes to the high prevalence of complete duplicates relative to structurally heterogeneous duplicates (partial and chimeric). The relative roles of regulatory sequence versus exon-intron structure changes in the acquisition of novel function by human paralogs remains to be determined.

MicroRNA-650 in a copy number-variable region regulates the production of interleukin 6 in human osteosarcoma cells

Copy number variation is a well-known genetic variation. microRNAs (miRNAs/miRs) are non-coding RNAs that mediate gene expression by regulating target mRNAs. In the present study, copy number deletions encompassing miRNA coding regions were investigated to determine the association between the deletion of miRNA and its phenotypic effects. A total of 38 human miRNAs in copy number variants were identified and miR-650, which is functional in the human osteosarcoma MG-63 cell line, was selected. Overexpression of miR-650 decreased the expression of inhibitor of growth family member 4 (ING4) in the MG-63 cells and increased interleukin (IL)6 transcription, as well as IL6 secretion in IL1B-stimulated cells. Furthermore, miR-650 downregulated the amount of nuclear factor of κ light polypeptide gene enhancer in B cells inhibitor α and increased the transcriptional activity of nuclear factor (NF)κB. Downregulation of ING4 also increased the production of IL6, similar to miR-650 overexpression. Taken together, these data indicate that miR-650 plays a significant role in the production of IL6 by regulating ING4 expression and NFκB signaling in IL1B-stimulated MG-63 osteosarcoma cells.

A decade of structural variants: description, history and methods to detect structural variation

In the past decade, the view on genomic structural variation (SV) has been changed completely. SVs, previously considered rare events, are now recognized as the largest source of interindividual genetic variation affecting more bases than single nucleotide polymorphisms, variable number of tandem repeats and other small genetic variants. They have also been shown to play a role in phenotypic variation and in disease. In this review, the authors will provide an introduction to SV; a short historical perspective on the research of this source of genomic variation; a description of the types of structural variants, and on how they may have arisen; and an overview on methods of detecting structural variants, focusing on the analysis of high-throughput sequencing data.

Interpreting human genetic variation with in vivo zebrafish assays

Rapid advances and cost erosion in exome and genome analysis of patients with both rare and common genetic disorders have accelerated gene discovery and illuminated fundamental biological mechanisms. The thrill of discovery has been accompanied, however, with the sobering appreciation that human genomes are burdened with a large number of rare and ultra rare variants, thereby posing a significant challenge in dissecting both the effect of such alleles on protein function and also the biological relevance of these events to patient pathology. In an effort to develop model systems that are able to generate surrogates of human pathologies, a powerful suite of tools have been developed in zebrafish, taking advantage of the relatively small (compared to invertebrate models) evolutionary distance of that genome to humans, the orthology of several organs and signaling processes, and the suitability of this organism for medium and high throughput phenotypic screening. Here we will review the use of this model organism in dissecting human genetic disorders; we will highlight how diverse strategies have informed disease causality and genetic architecture; and we will discuss relative strengths and limitations of these approaches in the context of medical genome sequencing. This article is part of a Special Issue entitled: From Genome to Function.

Resolving genomic disorder-associated breakpoints within segmental DNA duplications using massively parallel sequencing

The most common recurrent copy-number variants associated with autism, developmental delay and epilepsy are flanked by segmental duplications. Complete genetic characterization of these events is challenging because their breakpoints often occur within high-identity, copy-number polymorphic paralogous sequences that cannot be specifically assayed using hybridization-based methods. Here we provide a protocol for breakpoint resolution with sequence-level precision. Massively parallel sequencing is performed on libraries generated from haplotype-resolved chromosomes, genomic DNA or molecular inversion probe (MIP)-captured breakpoint-informative regions harboring paralog-distinguishing variants. Quantification of sequencing depth over informative sites enables breakpoint localization, typically within several kilobases to tens of kilobases. Depending on the approach used, the sequencing platform, and the accuracy and completeness of the reference genome sequence, this protocol takes from a few days to several months to complete. Once established for a specific genomic disorder, it is possible to process thousands of DNA samples within as little as 3-4 weeks.

In vitro microbial culture models and their application in drug development

Drug development faces its nemesis in the form of drug resistance. The rate of bacterial resistance to antibiotics, or tumor resistance to chemotherapy decisively depends on the surrounding heterogeneous tissue. However, in vitro drug testing is almost exclusively done in well stirred, homogeneous environments. Recent advancements in microfluidics and microfabrication introduce opportunities to develop in vitro culture models that mimic the complex in vivo tissue environment. In this review, we will first discuss the design principles underlying such models. Then we will demonstrate two types of microfluidic devices that combine stressor gradients, cell motility, large population of competing/cooperative cells and time varying dosage of drugs. By incorporating ideas from how natural selection and evolution move drug resistance forward, we show that drug resistance can occur at much greater rates than in well-stirred environments. Finally, we will discuss the future direction of in vitro microbial culture models and how to extend the lessons learned from microbial systems to eukaryotic cells.

Common copy number variation detection from multiple sequenced samples

Common copy number variations (CNVs) are small regions of genomic variations at the same loci across multiple samples, which can be detected with high resolution from next-generation sequencing (NGS) technique. Multiple sequencing data samples are often available from genomic studies; examples include sequences from multiple platforms and sequences from multiple individuals. By integrating complementary information from multiple data samples, detection power can be potentially improved. However, most of current CNV detection methods often process an individual sequence sample, or two samples in an abnormal versus matched normal study; researches on detecting common CNVs across multiple samples have been very limited but are much needed. In this paper, we propose a novel method to detect common CNVs from multiple sequencing samples by exploiting the concurrency of genomic variations in read depth signals derived from multiple NGS data. We use a penalized sparse regression model to fit multiple read depth profiles, based on which common CNV identification is formulated as a change-point detection problem. Finally, we validate the proposed method on both simulation and real data, showing that it can give both higher detection power and better break point estimation over several published CNV detection methods.

Copy-number changes in evolution: rates, fitness effects and adaptive significance

Gene copy-number differences due to gene duplications and deletions are rampant in natural populations and play a crucial role in the evolution of genome complexity. Per-locus analyses of gene duplication rates in the pre-genomic era revealed that gene duplication rates are much higher than the per nucleotide substitution rate. Analyses of gene duplication and deletion rates in mutation accumulation lines of model organisms have revealed that these high rates of copy-number mutations occur at a genome-wide scale. Furthermore, comparisons of the spontaneous duplication and deletion rates to copy-number polymorphism data and bioinformatic-based estimates of duplication rates from sequenced genomes suggest that the vast majority of gene duplications are detrimental and removed by natural selection. The rate at which new gene copies appear in populations greatly influences their evolutionary dynamics and standing gene copy-number variation in populations. The opportunity for mutations that result in the maintenance of duplicate copies, either through neofunctionalization or subfunctionalization, also depends on the equilibrium frequency of additional gene copies in the population, and hence on the spontaneous gene duplication (and loss) rate. The duplication rate may therefore have profound effects on the role of adaptation in the evolution of duplicated genes as well as important consequences for the evolutionary potential of organisms. We further discuss the broad ramifications of this standing gene copy-number variation on fitness and adaptive potential from a population-genetic and genome-wide perspective.