Fetal exome sequencing for isolated increased nuchal translucency: should we be doing it?

OBJECTIVE: To evaluate the utility of prenatal exome sequencing (ES) for isolated increased nuchal translucency (NT) and to investigate factors that increase diagnostic yield. DESIGN: Retrospective analysis of data from two prospective cohort studies. SETTING: Fetal medicine centres in the UK and USA. POPULATION: Fetuses with increased NT ≥3.5 mm at 11-14 weeks of gestation recruited to the Prenatal Assessment of Genomes and Exomes (PAGE) and Columbia fetal whole exome sequencing studies (n = 213). METHODS: We grouped cases based on (1) the presence of additional structural abnormalities at presentation in the first trimester or later in pregnancy, and (2) NT measurement at presentation. We compared diagnostic rates between groups using Fisher exact test. MAIN OUTCOME MEASURES: Detection of diagnostic genetic variants considered to have caused the observed fetal structural anomaly. RESULTS: Diagnostic variants were detected in 12 (22.2%) of 54 fetuses presenting with non-isolated increased NT, 12 (32.4%) of 37 fetuses with isolated increased NT in the first trimester and additional abnormalities later in pregnancy, and 2 (1.8%) of 111 fetuses with isolated increased NT in the first trimester and no other abnormalities on subsequent scans. Diagnostic rate also increased with increasing size of NT. CONCLUSIONS: The diagnostic yield of prenatal ES is low for fetuses with isolated increased NT but significantly higher where there are additional structural anomalies. Prenatal ES may not be appropriate for truly isolated increased NT but timely, careful ultrasound scanning to identify other anomalies emerging later can direct testing to focus where there is a higher likelihood of diagnosis.

Fetal hydrops and the Incremental yield of Next-generation sequencing over standard prenatal Diagnostic testing (FIND) study: prospective cohort study and meta-analysis.

OBJECTIVE: To determine the incremental yield of exome sequencing (ES) over chromosomal microarray analysis (CMA) or karyotyping in prenatally diagnosed non-immune hydrops fetalis (NIHF). METHODS: A prospective cohort study (comprising an extended group of the Prenatal Assessment of Genomes and Exomes (PAGE) study) was performed which included 28 cases of prenatally diagnosed NIHF undergoing trio ES following negative CMA or karyotyping. These cases were combined with data from a systematic review of the literature. MEDLINE, EMBASE, CINAHL and ClinicalTrials.gov databases were searched electronically (January 2000 to October 2020) for studies reporting on the incremental yield of ES over CMA or karyotyping in fetuses with prenatally detected NIHF. Inclusion criteria for the systematic review were: (i) at least two cases of NIHF undergoing sequencing; (ii) testing initiated based on prenatal ultrasound-based phenotype; and (iii) negative CMA or karyotyping result. The incremental diagnostic yield of ES was assessed in: (i) all cases of NIHF; (ii) isolated NIHF; (iii) NIHF associated with an additional fetal structural anomaly; and (iv) NIHF according to severity (i.e. two vs three or more cavities affected). RESULTS: In the extended PAGE study cohort, the additional diagnostic yield of ES over CMA or karyotyping was 25.0% (7/28) in all NIHF cases, 21.4% (3/14) in those with isolated NIHF and 28.6% (4/14) in those with non-isolated NIHF. In the meta-analysis, the pooled incremental yield based on 21 studies (306 cases) was 29% (95% CI, 24-34%; P < 0.00001; I2 = 0%) in all NIHF, 21% (95% CI, 13-30%; P < 0.00001; I2 = 0%) in isolated NIHF and 39% (95% CI, 30-49%; P < 0.00001; I2 = 1%) in NIHF associated with an additional fetal structural anomaly. In the latter group, congenital limb contractures were the most prevalent additional structural anomaly associated with a causative pathogenic variant, occurring in 17.3% (19/110) of cases. The incremental yield did not differ significantly according to hydrops severity. The most common genetic disorders identified were RASopathies, occurring in 30.3% (27/89) of cases with a causative pathogenic variant, most frequently due to a PTPN11 variant (44.4%; 12/27). The predominant inheritance pattern in causative pathogenic variants was autosomal dominant in monoallelic disease genes (57.3%; 51/89), with most being de novo (86.3%; 44/51). CONCLUSIONS: Use of prenatal next-generation sequencing in both isolated and non-isolated NIHF should be considered in the development of clinical pathways. Given the wide range of potential syndromic diagnoses and heterogeneity in the prenatal phenotype of NIHF, exome or whole-genome sequencing may prove to be a more appropriate testing approach than a targeted gene panel testing strategy. © 2021 The Authors. Ultrasound in Obstetrics & Gynecology published by John Wiley & Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology.

COngenital heart disease and the Diagnostic yield with Exome sequencing (CODE) study: prospective cohort study and systematic review.

OBJECTIVE: To determine the incremental yield of antenatal exome sequencing (ES) over chromosomal microarray analysis (CMA) or conventional karyotyping in prenatally diagnosed congenital heart disease (CHD). METHODS: A prospective cohort study of 197 trios undergoing ES following CMA or karyotyping owing to CHD identified prenatally and a systematic review of the literature were performed. MEDLINE, EMBASE, CINAHL and ClinicalTrials.gov (January 2000 to October 2019) databases were searched electronically for studies reporting on the diagnostic yield of ES in prenatally diagnosed CHD. Selected studies included those with more than three cases, with initiation of testing based upon prenatal phenotype only and that included cases in which CMA or karyotyping was negative. The incremental diagnostic yield of ES was assessed in: (1) all cases of CHD; (2) isolated CHD; (3) CHD associated with extracardiac anomaly (ECA); and (4) CHD according to phenotypic subgroup. RESULTS: In our cohort, ES had an additional diagnostic yield in all CHD, isolated CHD and CHD associated with ECA of 12.7% (25/197), 11.5% (14/122) and 14.7% (11/75), respectively (P = 0.81). The corresponding pooled incremental yields from 18 studies (encompassing 636 CHD cases) included in the systematic review were 21% (95% CI, 15-27%), 11% (95% CI, 7-15%) and 37% (95% CI, 18-56%), respectively. The results did not differ significantly when subanalysis was limited to studies including more than 20 cases, except for CHD associated with ECA, in which the incremental yield was greater (49% (95% CI, 17-80%)). In cases of CHD associated with ECA in the primary analysis, the most common extracardiac anomalies associated with a pathogenic variant were those affecting the genitourinary system (23/52 (44.2%)). The greatest incremental yield was in cardiac shunt lesions (41% (95% CI, 19-63%)), followed by right-sided lesions (26% (95% CI, 9-43%)). In the majority (68/96 (70.8%)) of instances, pathogenic variants occurred de novo and in autosomal dominant (monoallelic) disease genes. The most common (19/96 (19.8%)) monogenic syndrome identified was Kabuki syndrome. CONCLUSIONS: There is an apparent incremental yield of prenatal ES in CHD. While the greatest yield is in CHD associated with ECA, consideration could also be given to performing ES in the presence of an isolated cardiac abnormality. A policy of routine application of ES would require the adoption of robust bioinformatic, clinical and ethical pathways. Copyright © 2020 ISUOG. Published by John Wiley & Sons Ltd.

Clinically-relevant postzygotic mosaicism in parents and children with developmental disorders in trio exome sequencing data.

Mosaic genetic variants can have major clinical impact. We systematically analyse trio exome sequence data from 4,293 probands from the DDD Study with severe developmental disorders for pathogenic postzygotic mosaicism (PZM) in the child or a clinically-unaffected parent, and use ultrahigh-depth sequencing to validate candidate mosaic variants. We observe that levels of mosaicism for small genetic variants are usually equivalent in both saliva and blood and ~3% of causative de novo mutations exhibit PZM; this is an important observation, as the sibling recurrence risk is extremely low. We identify parental PZM in 21 trios (0.5% of trios), resulting in a substantially increased sibling recurrence risk in future pregnancies. Together, these forms of mosaicism account for 40 (1%) diagnoses in our cohort. Likely child-PZM mutations occur equally on both parental haplotypes, and the penetrance of detectable mosaic pathogenic variants overall is likely to be less than half that of constitutive variants.

Haploid chromosomes in molecular ecology: lessons from the human Y.

We review the potential use of haploid chromosomes in molecular ecology, using recent work on the human Y chromosome as a paradigm. Chromosomal sex-determination systems, and hence constitutively haploid chromosomes, which escape from recombination over much of their length, have evolved multiple times in the animal kingdom. In mammals, where males are the heterogametic sex, the patrilineal Y chromosome represents a paternal counterpart to mitochondrial DNA. Work on the human Y chromosome has shown it to contain the same range of polymorphic markers as the rest of the nuclear genome and these have rendered it the most informative haplotypic system in the human genome. Examples from research on the human Y chromosome are used to illustrate the common interests of anthropologists and ecologists in investigating the genetic impact of sex-specific behaviours and dispersals, as well as patterns of global diversity. We present some methodologies for extracting information from these uniquely informative yet under-utilized loci.

Geographical, linguistic, and cultural influences on genetic diversity: Y-chromosomal distribution in Northern European populations.

We analyzed 10 Y-chromosomal binary markers in 363 males from 8 populations in Northern Europe and 5 Y microsatellites in 346 of these individuals. These populations can be grouped according to cultural, linguistic, or geographical criteria, and the groupings are different in each case. We can therefore ask which criterion best corresponds to the distribution of genetic variation. In an AMOVA analysis using the binary markers, 13% of the Y variation was found between populations, indicating a high level of differentiation within this small area. No significant difference was seen between the traditionally nomadic Saami and the neighboring, historically farming, populations. When the populations were divided into Uralic speakers and Indo-European speakers, 8% of the variation was found between groups, but when they were divided according to geographical location, 14% of the variation was between groups. Geographical factors have thus been the most important in limiting gene flow between these populations, but linguistic differences have also been important in the east.

Patterns of inter- and intra-group genetic diversity in the Vlax Roma as revealed by Y chromosome and mitochondrial DNA lineages.

Previous genetic studies, supported by linguistic and historical data, suggest that the European Roma, comprising a large number of socially divergent endogamous groups, may be a complex conglomerate of founder populations. The boundaries and characteristics of such founder populations and their relationship to the currently existing social stratification of the Roma have not been investigated. This study is an attempt to address the issues of common vs independent origins and the history of population fissioning in three Romani groups that are well defined and strictly endogamous relative to each other. According to linguistic classifications, these groups belong to the Vlax Roma, who account for a large proportion of the European Romani population. The analysis of mtDNA sequence variation has shown that a large proportion of maternal lineages are common to the three groups. The study of a set of Y chromosome markers of different mutability has revealed that over 70% of males belong to a single lineage that appears unique to the Roma and presents with closely related microsatellite haplotypes and MSY1 codes. The study unambiguously points to the common origins of the three Vlax groups and the recent nature of the population fissions, and provides preliminary evidence of limited genetic diversity in this young founder population.

Gene conversion homogenizes the CMT1A paralogous repeats.

BACKGROUND: Non-allelic homologous recombination between paralogous repeats is increasingly being recognized as a major mechanism causing both pathogenic microdeletions and duplications, and structural polymorphism in the human genome. It has recently been shown empirically that gene conversion can homogenize such repeats, resulting in longer stretches of absolute identity that may increase the rate of non-allelic homologous recombination. RESULTS: Here, a statistical test to detect gene conversion between pairs of non-coding sequences is presented. It is shown that the 24 kb Charcot-Marie-Tooth type 1A paralogous repeats (CMT1A-REPs) exhibit the imprint of gene conversion processes whilst control orthologous sequences do not. In addition, Monte Carlo simulations of the evolutionary divergence of the CMT1A-REPs, incorporating two alternative models for gene conversion, generate repeats that are statistically indistinguishable from the observed repeats. Bounds are placed on the rate of these conversion processes, with central values of 1.3 x 10(-4) and 5.1 x 10(-5) per generation for the alternative models. CONCLUSIONS: This evidence presented here suggests that gene conversion may have played an important role in the evolution of the CMT1A-REP paralogous repeats. The rates of these processes are such that it is probable that homogenized CMT1A-REPs are polymorphic within modern populations. Gene conversion processes are similarly likely to play an important role in the evolution of other segmental duplications and may influence the rate of non-allelic homologous recombination between them.

Y-chromosomal diversity in Europe is clinal and influenced primarily by geography, rather than by language.

Clinal patterns of autosomal genetic diversity within Europe have been interpreted in previous studies in terms of a Neolithic demic diffusion model for the spread of agriculture; in contrast, studies using mtDNA have traced many founding lineages to the Paleolithic and have not shown strongly clinal variation. We have used 11 human Y-chromosomal biallelic polymorphisms, defining 10 haplogroups, to analyze a sample of 3,616 Y chromosomes belonging to 47 European and circum-European populations. Patterns of geographic differentiation are highly nonrandom, and, when they are assessed using spatial autocorrelation analysis, they show significant clines for five of six haplogroups analyzed. Clines for two haplogroups, representing 45% of the chromosomes, are continentwide and consistent with the demic diffusion hypothesis. Clines for three other haplogroups each have different foci and are more regionally restricted and are likely to reflect distinct population movements, including one from north of the Black Sea. Principal-components analysis suggests that populations are related primarily on the basis of geography, rather than on the basis of linguistic affinity. This is confirmed in Mantel tests, which show a strong and highly significant partial correlation between genetics and geography but a low, nonsignificant partial correlation between genetics and language. Genetic-barrier analysis also indicates the primacy of geography in the shaping of patterns of variation. These patterns retain a strong signal of expansion from the Near East but also suggest that the demographic history of Europe has been complex and influenced by other major population movements, as well as by linguistic and geographic heterogeneities and the effects of drift.

Divergent outcomes of intrachromosomal recombination on the human Y chromosome: male infertility and recurrent polymorphism.

The Y chromosome provides a unique opportunity to study mutational processes within the human genome, decoupled from the confounding effects of interchromosomal recombination. It has been suggested that the increased density of certain dispersed repeats on the Y could account for the high frequency of causative microdeletions relative to single nucleotide mutations in infertile males. Previously we localised breakpoints of an AZFa microdeletion close to two highly homologous complete human endogenous retroviral sequences (HERV), separated by 700 kb. Here we show, by sequencing across the breakpoint, that the microdeletion occurs in register within a highly homologous segment between the HERVs. Furthermore, we show that recurrent double crossovers have occurred between the HERVs, resulting in the loss of a 1.5 kb insertion from one HERV, an event underlying the first ever Y chromosomal polymorphism described, the 12f2 deletion. This event produces a substantially longer segment of absolute homology and as such may result in increased predisposition to further intrachromosomal recombination. Intrachromosomal crosstalk between these two HERV sequences can thus result in either homogenizing sequence conversion or a microdeletion causing male infertility. This represents a major subclass of AZFa deletions.

Recent male-mediated gene flow over a linguistic barrier in Iberia, suggested by analysis of a Y-chromosomal DNA polymorphism.

We have examined the worldwide distribution of a Y-chromosomal base-substitution polymorphism, the T/C transition at SRY-2627, where the T allele defines haplogroup 22; sequencing of primate homologues shows that the ancestral state cannot be determined unambiguously but is probably the C allele. Of 1,191 human Y chromosomes analyzed, 33 belong to haplogroup 22. Twenty-nine come from Iberia, and the highest frequencies are in Basques (11%; n=117) and Catalans (22%; n=32). Microsatellite and minisatellite (MSY1) diversity analysis shows that non-Iberian haplogroup-22 chromosomes are not significantly different from Iberian ones. The simplest interpretation of these data is that haplogroup 22 arose in Iberia and that non-Iberian cases reflect Iberian emigrants. Several different methods were used to date the origin of the polymorphism: microsatellite data gave ages of 1,650, 2,700, 3,100, or 3,450 years, and MSY1 gave ages of 1,000, 2,300, or 2,650 years, although 95% confidence intervals on all of these figures are wide. The age of the split between Basque and Catalan haplogroup-22 chromosomes was calculated as only 20% of the age of the lineage as a whole. This study thus provides evidence for direct or indirect gene flow over the substantial linguistic barrier between the Indo-European and non-Indo-European-speaking populations of the Catalans and the Basques, during the past few thousand years.

European Y-chromosomal lineages in Polynesians: a contrast to the population structure revealed by mtDNA.

We have used Y-chromosomal polymorphisms to trace paternal lineages in Polynesians by use of samples previously typed for mtDNA variants. A genealogical approach utilizing hierarchical analysis of eight rare-event biallelic polymorphisms, seven microsatellite loci, and internal structural analysis of the hypervariable minisatellite, MSY1, has been used to define three major paternal-lineage clusters in Polynesians. Two of these clusters, both defined by novel MSY1 modular structures and representing 55% of the Polynesians studied, are also found in coastal Papua New Guinea. Reduced Polynesian diversity, relative to that in Melanesians, is illustrated by the presence of several examples of identical MSY1 codes and microsatellite haplotypes within these lineage clusters in Polynesians. The complete lack of Y chromosomes having the M4 base substitution in Polynesians, despite their prevalence (64%) in Melanesians, may also be a result of the multiple bottleneck events during the colonization of this region of the world. The origin of the M4 mutation has been dated by use of two independent methods based on microsatellite-haplotype and minisatellite-code diversity. Because of the wide confidence limits on the mutation rates of these loci, the M4 mutation cannot be conclusively dated relative to the colonization of Polynesia, 3,000 years ago. The other major lineage cluster found in Polynesians, defined by a base substitution at the 92R7 locus, represents 27% of the Polynesians studied and, most probably, originates in Europe. This is the first Y-chromosomal evidence of major European admixture with indigenous Polynesian populations and contrasts sharply with the picture given by mtDNA evidence.