The variable somatic genome

Whole-Exome Sequencing of Discordant Monozygotic Twin Families for Identification of Candidate Genes for Microtia-Atresia

Objective

We used data from twins and their families to probe the genetic factors contributing to microtia-atresia, in particular, early post-twinning variations that potentially account for the discordant phenotypes of monozygotic twin pairs.

Methods

Six families of monozygotic twins discordant for congenital microtia-atresia were recruited for study. The six patients shared a consistent clinical phenotype of unilateral microtia-atresia. Whole-exome sequencing (WES) was performed for all six twin pairs and their parents. Family segregation and multiple bioinformatics methods were applied to identify suspicious mutations in all families. Recurring mutations commonly detected in at least two families were highlighted. All variants were validated via Sanger sequencing. Gene Ontology (GO) analysis was performed to identify candidate gene sets and related pathways. Copy number variation (CNV), linkage analysis, association analysis and machine learning methods were additionally applied to isolate candidate mutations, and comparative genomics and structural modeling tools used to evaluate their potential roles in onset of microtia-atresia.

Results

Our analyses revealed 61 genes with suspected mutations associated with microtia-atresia. Five (HOXA4, MUC6, CHST15, TBX10, and AMER1) contained 7 de novo mutations that appeared in at least two families, which have been previously reported as pathogenic for other diseases. Among these, HOXA4 (c.920A>C, p.H307P) was determined as the most likely pathogenic variant for microtia-atresia. GO analysis revealed four gene sets involving 11 pathways potentially related to underlying pathogenesis of the disease. CNVs in three genes (UGT2B17, OVOS, and KATNAL2) were detected in at least two families. Linkage analysis disclosed 13 extra markers for the disease, of which two (FGFR1 and EYA1) were validated via machine learning analysis as plausible candidate genes for the disease.

Conclusion

Based on comprehensive genetic and bioinformatic analyses of WES data from six families of discordant monozygotic twins with microtia-atresia, we identified multiple candidate genes that may function in post-twinning onset of the disease. The collective findings provide novel insights into the pathogenesis of congenital microtia-atresia.

Whole-Exome Sequencing in Nine Monozygotic Discordant Twins

By definition, monozygotic (MZ) twins carry an identical set of genetic information. The observation of early post-twinning mutational events was shown to cause phenotypic discordance among MZ twin pairs. These mutational events comprise genomic alterations at different scales, ranging from single nucleotide changes to larger copy-number variations (CNVs) of varying sizes, as well as epigenetic changes. Here, we performed whole-exome sequencing (WES) in nine discordant MZ twins to identify somatic mutational events in the affected twin that might exert a dominant negative effect. Five of these MZ twin pairs were discordant for congenital heart defects (CHD), two for endocrine disorders, one for omphalocele, and one for congenital diaphragmatic hernia (CDH). Analysis of WES data from all nine MZ twin pairs using the de novo probability tool DeNovoGear detected only one apparent de novo variation in TMPRSS13 in one of the CHD-affected twins. Analysis of WES data from all nine MZ twin pairs by using standard filter criteria without the de novo probability tool DeNovoGear revealed a total of 6,657 variations in which both the twin pairs differed. After filtering for variations only present in the affected twins and absent in in-house controls, 722 variations remained. Visual inspection for read quality decreased this number to 12, present only in the affected twin. However, Sanger sequencing of the overall 13 variations failed to confirm the variation in the affected twin. These results suggest that somatic mutational events in coding regions do not seem to play a major role in the phenotypic expression of MZ discordant twin pairs.

Genomic Variants Revealed by Invariably Missing Genotypes in Nelore Cattle

High density genotyping panels have been used in a wide range of applications. From population genetics to genome-wide association studies, this technology still offers the lowest cost and the most consistent solution for generating SNP data. However, in spite of the application, part of the generated data is always discarded from final datasets based on quality control criteria used to remove unreliable markers. Some discarded data consists of markers that failed to generate genotypes, labeled as missing genotypes. A subset of missing genotypes that occur in the whole population under study may be caused by technical issues but can also be explained by the presence of genomic variations that are in the vicinity of the assayed SNP and that prevent genotyping probes from annealing. The latter case may contain relevant information because these missing genotypes might be used to identify population-specific genomic variants. In order to assess which case is more prevalent, we used Illumina HD Bovine chip genotypes from 1,709 Nelore (Bos indicus) samples. We found 3,200 missing genotypes among the whole population. NGS re-sequencing data from 8 sires were used to verify the presence of genomic variations within their flanking regions in 81.56% of these missing genotypes. Furthermore, we discovered 3,300 novel SNPs/Indels, 31% of which are located in genes that may affect traits of importance for the genetic improvement of cattle production.