Known Mutations at the Cause of Alpha-1 Antitrypsin Deficiency an Updated Overview of SERPINA1 Variation Spectrum

Abstract

Alpha-1-Antitrypsin deficiency (AATD), caused by SERPINA1 mutations, is one of the most prevalent Mendelian disorders among individuals of European descend. However, this condition, which is characterized by reduced serum levels of alpha-1-antitrypsin (AAT) and associated with increased risks of pulmonary emphysema and liver disease in both children and adults, remains frequently underdiagnosed. AATD clinical manifestations are often correlated with two pathogenic variants, the Z allele (p.Glu342Lys) and the S allele (p.Glu264Val), which can be combined in severe ZZ or moderate SZ risk genotypes. Yet, screenings of AATD cases and large sequencing efforts carried out in both control and disease populations are disclosing outstanding numbers of rare SERPINA1 variants (>500), including many pathogenic and other likely deleterious mutations. Generally speaking, pathogenic variants can be subdivided into either loss- or gain-of-function according to their pathophysiological effects. In AATD, the loss-of-function is correlated with an uncontrolled activity of elastase by its natural inhibitor, the AAT. This phenomenon can result from the absence of circulating AAT (null alleles), poor AAT secretion from hepatocytes (deficiency alleles) or even from a modified inhibitory activity (dysfunctional alleles). On the other hand, the gain-of-function is connected with the formation of AAT polymers and their switching on of cellular stress and inflammatory responses (deficiency alleles). Less frequently, the gain-of-function is related to a modified protease affinity (dysfunctional alleles). Here, we revisit SERPINA1 mutation spectrum, its origins and population history with a greater emphasis on variants fitting the aforementioned processes of AATD pathogenesis. Those were selected based on their clinical significance and wider geographic distribution. Moreover, we also provide some directions for future studies of AATD clinically heterogeneity and comprehensive diagnosis.

Keywords:

• SERPINA1 variants
• Z allele
• S allele
• M_Malton allele and Q_Ourém allele

Introduction

Alpha-1-antitrypsin deficiency (AATD) is an autosomal codominant disorder caused by multiple mutations affecting SERPINA1 gene. Although in most cases these are translated into a functional protein present in the bloodstream at lower concentrations, some mutations can also lead to dysfunctional molecules or null variants.^1–3 The major clinical implications of AATD are an increased risk of chronic obstructive pulmonary disease (COPD) and emphysema in adults, and hepatic illness during childhood or later in adult life. Nonetheless, AATD has been associated with disorders affecting not only the respiratory system but also other organs, such as bronchiectasis, anti-neutrophil cytoplasmic antibody (ANCA) – associated vasculitis, and panniculitis.^1,4–7 AATD has been typically correlated with European populations, where its prevalence can reach up to 1:2000 subjects.^8,9 However, this genetic condition can be virtually detected anywhere due to a worldwide dispersal of two more prevalent deficiency alleles, the S and Z alleles, and a large spectrum of rare pathogenic mutations, which are typically more geographically confined.^2,10–13

The AAT Molecule and Pathophysiology of AATD

Alpha-1-Antitrypsin (AAT) is the major protease inhibitor in the serum and the archetypical member of the SERerine Protease INhibitors (SERPIN) superfamily. Therefore, it is often referred to as SERPINA1 – SERPIN clade A member.^14,15 AAT is a monomeric glycoprotein that once in the bloodstream weighs approximately 52 KDa and is composed of 394 amino acid residues. However, prior to its secretion, the AAT molecule is post-translationally modified by the removal of its signal peptide, composed of 24 amino acids, and by the linkage of three carbohydrate side chains to asparagine residues located at positions 46, 83 and 247.^5,7,16

AAT is primarily synthetized at high levels by hepatocytes, but also to a lower extent by many other tissues, such as whole blood innate immunity cells (neutrophils and monocytes), epithelial lung cells, the kidney and the skin (Figure 1A).^4,17–19 At the structural level, AAT is a globular protein packed into three β-sheets (A-C) and nine α-helices (A-I), showing a solvent-exposed reactive center loop (RCL; Figure 1B) that determines its inhibitory affinity towards distinct proteases.^3,20 Precisely, it is the Met358-Ser359 bond (P1-P1ʹ; Figure 1B) located within the RCL that by mimicking a protease substrate, works like a specific bite in a kind of molecular mouse-trap.^3,20 This remarkable feature of AAT can only be accomplished by the ability of SERPINs to undergo major conformational rearrangements upon RCL cleavage, in which proteases are moved into the opposite pole of the SERPIN molecule. Thereafter, the structure of the protease is disturbed and permanently inactivated in a stable complex with its inhibitor (Figure 1C).^3,20 Although AAT is mainly recognized as a potent inhibitor of neutrophil elastase, which is released in the lung during infectious and inflammatory states, it is also capable of regulating other proteases, including proteinase 3, kallikrein 7, matriptase, and caspase-3.^4,21 Importantly, AAT is also an acute phase molecule with multiple immunomodulatory and anti-inflammatory functions, such as in the control and release of cytokines, the recruitment of innate immunity cells and in responses to diverse pathogens.^4,21,22

Figure 1 (A) SERPINA1 expression across different tissues. The graph was built in GTEx portal (https://gtexportal.org/home/) using the RNAseq data from multiple individuals generated by the Genotype-Tissue Expression (GTEx) project. (B) Three-dimensional structure of an active AAT molecule (PDB code: 1QLP). The breach, shutter and gate domains are highlighted by dark blue circles. The reactive center loop (RCL) and Met-Ser bond (P1-P1ʹ) therein are indicated. (C) Three-dimensional structure of an AAT molecule bond to trypsin in a protease-inhibitor complex (PDB code: 1EXZ). The inserted β-strand of the RCL and the Met358 and Ser359 residues from the cleaved P1-P1ʹ are indicated. α-helixes are shown in light blue, β-strands in red and loops in magenta. The three-dimensional images were generated using PyMOL.¹¹²

Broadly speaking, the clinical manifestations of the AATD condition are partially attributed to the low levels of AAT and an unbalanced activity of proteases, which can lead to tissue degradation through the cleavage of extracellular matrix proteins as described for pulmonary emphysema.^9,19 In general, a serum concentration of 11 uM (57 mg/dl) is considered as the minimum protective threshold against the lung proteolytic attack.^1,5 Nevertheless, some authors already questioned its universal application given that some AATD patients with values above the proposed limit may still exhibit respiratory symptoms and benefit from specific clinical interventions.^23,24 On the other hand, AATD complications can be also correlated with the AAT conformational flexibility and its propensity to polymerize when misfolded. The toxic effects of AAT polymers inside the endoplasmic reticulum (ER) have long been recognized as a primary cause of the hepatic disease associated with AATD. More recently, the role of AAT polymers in the pathophysiology of pulmonary emphysema was also acknowledged after their detection in the lungs.^9,19 The underlying mechanisms of panniculitis and ANCA vasculitis linked to AATD are not yet fully elucidated, but those are probably correlated with both processes of unbalanced proteolysis and polymer accumulation.^1,6

The SERPINA1 Gene and Its Variants

AAT is encoded by SERPINA1 (SERPIN family A member 1) gene, which is located at chromosomal region 14q32.13, and spans for about 13.9 kb (GRCh38.p13/hg38 genomic coordinates chr14:94,376,747–94,390,654). SERPINA1 is organized into seven exons of which the first three are untranslated regions (UTR), often labeled as IA, IB and IC, and the remaining four comprise the coding region and the terminal UTR (exons II–V). The 5ʹ UTR contains distinct transcription initiation sites being IA and IB exons reported as alternatively used by monocytes, and IC specifically employed in the gene expression by hepatocytes.^6,7,25

So far, more than 120 SERPINA1 mutations have been reported in the literature, of which near 40% are described to cause some type of AATD.² Despite the vast majority of pathogenic mutations is localized within the coding region (II–V exon) and engages either single nucleotide substitutions or small insertion and deletions (1 to 3 nucleotides), there are also a few reports of AATD causing variants affecting splice sites, as well as large deletions eliminating partially or entirely the sequence of SERPINA1 gene (Table 1).²

Table 1 List of Known SERPINA1 Mutations Defining Different AAT Protein Alleles. Mutations are Listed According to Their Order of Appearance in Gene (Exon IA to Exon V) or Protein Sequence (Amino Acid −24 to 394)

SERPINA1 Mutation (5ʹ→3ʹ)	Marker ID	Allele Name^a (Molecular Background)	DNA Sequence^b	Population Group and Frequency (gnomAD Exomes)^c	Mutation Effect^d
SERPINA1 deletion		Q0Riedenburg			AATD: protein absence (LoF)
g.5307_5308ins8bp		Q0Savannah			AATD: protein absence (LoF)
c.-5+1G > A	rs775786225	Q0Porto	GT→AT		AATD: protein absence (LoF)
c.-5+2dupT		Q0Madrid (M3)	GT^T		AATD: protein absence (LoF)
		Q0Faro (M1Val)
Del 17Kb including exons II–V		Q0Isola di Procida			AATD: protein absence (LoF)
c.646+1G>T	rs751235320	Q0West	AG→AT		AATD: protein absence (LoF)
c.647–1delG	rs1555368758	Q0Bonny blue	AgGA		AATD: protein absence (LoF)
p.Ser-19Leu	rs140814100	ZWrexham (Z)	TCG→TTG	AFR - 0.02%; EAS - 0.05%; NFE <0.01%; SAS <0.01%	None
p.Asp2Ala	rs199422212	VMunich	GAT→GCT		None
p.Ser14Phe	rs745463238	SDonosti (S)	TCC→TTC	AMR <0.01%; NFE <0.01%	None
p.His15Leu		M1Bruxelles	CAC→CTC		None
p.Ala34Thr	rs149319176	M5Karlsruhe	GCC→ACC	AFR - 0.01%; AMR <0.01%; NFE <0.01%; SAS <0.01%	None
p.Tyr38Ter	rs762321137	Q0Knowloon	TAC→TAA	SAS <0.01%	AATD: protein absence (LoF)
p.Arg39Cys	rs28931570	I	CGC→TGC	AFR - 0.05%; AMR - 0.04%; EAS - 0.02%; FIN <0.01%; NFE - 0.21%; OTH - 0.07%; SAS <0.01%	AATD: protein deficiency (LoF and GoF)
p.Arg39His	rs764726147	MRouen	CGC→CAC	EAS - 0.01%; NFE<0.01%	Unknown
p.Leu41Pro	rs28931569	MProcida	CTG→CCG	NFE <0.01%; OTH - 0.03%	AATD: protein deficiency (LoF)
p.Leu41Trpfs*15		MVarallo	30bp del replaced by a new 22bp sequence		Uncertain AATD: protein deficiency or absence (LoF)
p.Ser45Phe	rs199687431	M6Bonn	TCC→TTC	FIN <0.01%; NFE - 0.01%	None
p.Ser47Arg	rs11575873	SRoubaix	AGC→C GC	AMR - 0.06%; NFE - 0.01%; OTH - 0.13%; SAS <0.01%	None
p.Phe52del	rs775982338	MMalton (M2)	Atc tTC TTC	AFR <0.01%;	AATD: protein deficiency (LoF and GoF)
		MPalermo (M1Val)	TCC, ATC TTC	AMR - 0.03%;
		MNichinan (V)	Ttc tCC or ATC	EAS - 0.04%;
		Q0LaPalma (S)	TTC ttc TCC	NFE - 0.01%
p.Ser53Phe	rs55819880	Siiyama	TCC→TTC		AATD: protein deficiency (LoF and GoF)
p.Ala60Thr	rs111850950	M6Passau	GCC→ ACC	AFR <0.01%; AMR - 0.02%; FIN - 0.01%; NFE - 0.05%; OTH - 0.03%	None
p.Gly67Glu	rs28931568	Mmineral Springs	GGG→GAG		AATD: protein deficiency (LoF)
p.His73Metfs*7	rs1057516212	Q0Casablanca (M2)	GAC ACt cac GAT		AATD: protein absence (LoF)
		Q0Lille (Z)
p.Thr68Ile	rs1490133295	Q0Lisbon	ACC→ATC	AMR - 0.02%	AATD: Protein absence (LoF)
p.Glu75Val		Trento	GGA→GTA		AATD: protein deficiency (LoF)
p.Thr85Met	rs199422213	ZBristol (M1Val)	ACG→ATG	EAS <0.01%; NFE <0.01%; SAS <0.01%	AATD: protein deficiency (LoF and GoF)
p.Pro88Thr	rs886044322	M5Berlin	CCG → ACG		None
p.Ile92Asn	rs28931572	Q0Ludwisghafen	ATC→AAC		AATD: protein absence (LoF)
p.Gly95Val	rs749295615	M1Saint-rambert	GGC→GTC	SAS <0.01%	None
p.Arg101His	rs709932	M2 (M3)	CGT→CAT	AFR - 3.32%; AMR - 10.79%; ASJ - 16.21%; EAS - 20.66%; FIN - 12.06%; NFE - 16.4%; OTH - 16.61%; SAS - 26.9%	None
		M4 (M1Val)
p.Thr102Profs*10		Q0Soest	aCC CTC		AATD: protein absence (LoF)
p.Gly115Ser	rs11558261	Q0Devon (Unknown)	GGC→AGC	NFE <0.01%; SAS - 0.07%	AATD: protein absence (LoF)
		Q0Newport (Z)
p.Glu122Lys	rs537285845	WSaint-Avre	GAG→AAG	AFR - 0.02%; NFE <0.01%; OTH - 0.02%; SAS <0.01%	Unknown
p.Leu126Arg		WVernaison	CTA→CGA		AATD: protein deficiency (LoF)
p. Tyr138Cys		Q0Vigo (S)	TAC→TGC		AATD: protein absence (LoF and likely GoF)
p.Gly148Arg	rs112030253	V (M1Val)	GGG→AGG	AFR - 0.17%; AMR - 0.02%; EAS <0.01%; NFE - 0.04%; OTH - 0.05%	None
		MNichinan (p.Phe52del)
p.Gly148Trp	rs112030253	M2Obernburg	GGG→TGG	AFR - 0.01%; AMR <0.01%; FIN <0.01%; NFE - 0.13%; OTH - 0.13%; SAS - 0.24%	None
p. Glu151Lys	rs149770048	M1Cadiz	GAA→AAA	AMR <0.01%; NFE <0.01%; SAS <0.01%	None
p.Lys154Asn		Queen´s	AAG→AAC		AATD: protein deficiency (LoF and GoF)
p.Gln156Ter	rs864622051	Q0Chillichote	CAG→TAG		AATD: protein absence (LoF)
p.Gln156Glu	rs864622051	LFrankfurt	CAG→GAG		None
p.Asp159Asn	rs759578830	OThonon-les-bains	GAT→AAT	NFE <0.01%; SAS - 0.03%	None
p.Tyr160Ter	rs267606950	Q0Granite Falls	TAc GTG		AATD: protein absence (LoF)
p.Tyr160Ter	rs199422210	Q0Bredevoor (Unknown)	TAC→TAG	NFE <0.01%	AATD: protein absence (LoF)
		Q0Amersfoort (M1Ala)
p.Glu162Gly	rs1035710135	PGaia	GAG→GGG		AATD: protein deficiency (LoF)
p.Lys163Ter		Q0Saint-Etienne	AAG→TAG		AATD: protein absence (LoF)
p.Thr180Serfs*11	rs921982028	Q0Cork	Aca GTT	AMR <0.01%	AATD: protein absence (LoF)
p.Trp194Ter	rs1445192595	Q0Trastevere	TGG→TGA		AATD: protein absence (LoF)
p.Glu204Lys	rs199422208	X	GAG→AAG	AFR <0.01%; NFE <0.01%	None
p.Val213Ala	rs6647	M1	GTG→GCG	AFR - 55.85%; AMR - 9.41%; ASJ - 25.89%; EAS - 2.13%; FIN - 29.86%; NFE - 21.76%; OTH - 21.16%; SAS – 16%	None
p.Lys217Ter	rs199422211	Q0Bellingham	AAG→TAG	AFR <0.01%; NFE <0.01%	AATD: protein absence (LoF)
p.Met221Thr	rs766260108	PLoyettes	ATG→ACG	NFE <0.01%	AATD: protein deficiency- uncertain and reduced inhibitory activity (LoF)
p.Met221Ile	rs864622053	PSolaize	ATG→ATA		AATD: protein deficiency-and reduced inhibitory activity (LoF)
p.Arg223Cys	rs28929470	F	CGT→TGT	AFR - 0.12%; AMR - 0.01%; ASJ - 0.04%; FIN - 0.05%; NFE - 0.43%; OTH - 0.15%	None - reduced inhibitory activity (LoF)
Gly225Arg	rs764220898	Pbrescia	GGC→CGC	NFE <0.01%	AATD: protein deficiency (LoF and GoF)
p.Val239Cysfs*3	rs1200349975	Q0Perugia	TGG gTG		AATD: protein absence (LoF)
p.Asn247Asp	rs755851961	XCuris	AAT→GAT	NFE <0.01%	None
p.Asp256Val	rs121912714	PDuarte (M4)	GAT→GTT	AMR - 0.05%; ASJ - 0.33%; NFE - 0.05%; OTH - 0.05%; SAS <0.01%	AATD: protein deficiency (LoF)
		PLowell (M1Val)
		YBarcelona (p.Pro391His)
p.Glu257Ter		Q0Brescia	GAG→TAG		AATD: protein absence (Lof)
p.Lys259Ile	rs864622044	MPisa	AAA→ATA	NFE <0.01%	AATD: protein deficiency (LoF and GoF)
p.Lys259Ter	rs1566753359	Q0Cairo	AAA→TAA	NFE <0.01%	AATD: protein absence (LoF)
p.Lys259-Glu264delTer		Q0Milano	17 bp del		AATD: protein absence (LoF)
p.Leu263Pro	rs1450795004	Q0Gaia	CTG→CCG	SAS <0.01%	AATD: protein absence (LoF)
p.Glu264Val	rs17580	S (M1Val)	GGA→GTA	AFR - 0.79%; AMR - 3.01%; ASJ - 1.37%; EAS <0.01%; FIN - 0.85%; NFE - 3.65%; OTH - 3.42%	AATD: protein deficiency (LoF and GoF)
		T (M4 and M3)
p.His269Gln	rs141095970	M1Lille	CAC→CAA	AFR - 0.2%; AMR <0.01%; NFE <0.01%; OTH - 0.02%	None
p.Arg281Lysfs*17		Q0Oliveira do Douro	GAC Aga AGG		AATD: protein absence (LoF)
p.Arg281del	rs748777702	IEuskadi	GAC aga AGG	AMR <0.01%	None
c.917 + 1G > A	rs750779440	Q0Achicourt	GT→AT	SAS <0.01%	AATD: protein absence (LoF)
c.918–1G > A	rs764950047	Q0Saint-Avold	AG→AA		AATD: protein absence (LoF)
p.Ala284Ser	rs141620200	M1Lyon	GCC→TCC	AFR - 0.06%; AMR <0.01%; ASJ - 0.09%; FIN - 0.21%; NFE - 0.37%; OTH - 0.23%	None
p.Try297Ter		Q0Torino	TAT→TAA		AATD: protein absence (LoF)
p.Tyr297Cys	rs774775536	M1Brest	TAT→TGT		None
p.Gln305Ter		Q0Cosenza	CAA→TAA		AATD: protein absence (LoF)
p.Ser319ArgfsTer16	rs1057519610	Q0Hong Kong	CTC tcC or Ctc TCC		AATD: protein absence (LoF)
p.Gly320Glu		Q0New Hope	GGG→GAG		AATD: protein absence (LoF)
p.Leu327Argfs*12	rs72555374	Q0Pordenone	CtG AAG		AATD: protein absence (LoF)
p.Ser330Phe	rs201788603	SMunich	TCC→TTC	AMR <0.01%; NFE - 0.03%; SAS <0.01%	None
c.1065 + 1G > A	rs781591420	Q0Amiens	GT→AT	NFE <0.01%	AATD: protein absence (LoF)
p.Val333Met	rs373630097	Not attributed	GTG→ATG	AFR <0.01%; AMR <0.01%; NFE <0.01%; OTH - 0.07%; SAS - 0.34%	AATD: protein deficiency (LoF and GoF)
p.His334Asp		King’s	CAT→GAT		AATD: protein deficiency (LoF and GoF)
p.His334Gln		M1Cremeaux	CAT→CAA		AATD: protein deficiency (LoF)
p.Ala336Thr	rs1802959	WBethesda	GCT→ACT	NFE <0.01%; OTH - 0.02%; SAS <0.01%	AATD: protein deficiency (LoF)
p.Asp341Asn	rs143370956	PDonauwoerth (M1Val)	GAC→AAC	AFR - 0.15%; AMR - 0.01%; SAS <0.01%	None
		PSaint albans (Z)
p.Glu342Lys	rs28929474	Z (M1Ala)	GAG→AAG	AFR - 0.29%; AMR - 0.38%; ASJ - 0.84%; FIN - 1.78%; NFE - 1.84%; OTH - 0.88%; SAS <0.01%	AATD: protein deficiency (LoF and GoF)
		ZAugsburg (M2)
p.Leu353Phefs*24	rs763023697	Q0Mattawa (M1Val)	TTT^T	AMR <0.01%; NFE <0.01%	AATD: protein absence (LoF)
		Q0Ourem (M3)
p.Met358Arg	rs121912713	Pittsburgh	ATG→AGG		None - Altered inhibitory activity (GoF)
p.Pro362Thr	rs12233	LOffenbach	CCC→ACC	AFR <0.01%; AMR <0.01%; EAS - 0.01%; NFE <0.01%; SAS - 0.01%	None
p.Pro362His	rs569384943	São Tomé	CCC→CAC		None
p.Glu363ArgfsTer11	rs764325655	Q0Bolton	CCc GAG GTC	NFE <0.01%	AATD: protein absence (LoF)
p.Glu363ArgfsTer14	rs764325655	Q0Clayton (M1Val)	CCC^C	NFE <0.01%	AATD: protein absence (LoF)
		Q0Saarbruecken (M1Ala)
p.Glu363Lys	rs121912712	XChristchurch	GAG→AAG	AFR <0.01%; EAS - 0.14%; FIN - 0.26%; NFE <0.01%; OTH - 0.02%; SAS - 0.02%	None
p.Lys368Glu		ETaurisano	AAA→GAA		AATD: protein deficiency (LoF and GoF)
p.Pro369Leu	rs199422209	MHeerlen (M1Ala)	CCC→CTC	NFE - 0.01%	AATD: protein absence (LoF)
		MHeerlen (M1Val)
p.Pro369Ser	rs61761869	MWurzburg	CCC→TCC	AFR <0.01%; AMR - 0.02%; FIN <0.01%; NFE - 0.05%; OTH - 0.05%; SAS <0.01	AATD: protein deficiency (LoF and GoF)
p.Phe370Leufs*4		Q0Dublin	TTt GTC		AATD: protein absence (LoF)
p.Met374Leufs*19		Q0Vila Real	TTA atg aTT		AATD: protein absence (LoF)
p.Glu376Asp	rs1303	M3	GAA→GAC	AFR - 10.96%; AMR - 32.25%; ASJ - 26.61%; EAS - 28.4%; FIN - 23.23%; NFE - 25.65%; OTH - 27.49%; SAS - 42.52%	None
p.Val389del	rs760849035	Q0Montluel	GTG gtg AAT	NFE <0.01%	AATD: protein absence (LoF)
p.Pro391His		YBarcelona (PLowell)	CCC→CAC		AATD: protein deficiency (LoF and GoF)
		YOrzinuovi (M1Val)
p.Pro391Thr		Q0Aachen (S)	CCC→ACC		AATD: protein absence (LoF and likely GoF)
p.Lys394Ter		GSaint-sorlin	AAA→TAA		Unknown

Notes: ^aThe list of known alleles was retrieved from AATD specialized literature.^{3,13,70–72,82–85,89,90,92,94,99,104–111} Molecular background is provided to discriminate alleles carrying the same mutation. ^bNucleotide substitutions are underlined, deleted bases are presented in lower case and inserted bases are preceded by an ^ symbol. ^cGnomAD groups where the mutation was detected and corresponding frequencies. No value is shown in the absence of the mutation in the population group. For variants present at very low-frequencies the values are summarized as <0.01%. ^dMutations effects are reported according to their impact in bloodstream AAT concentrations and classified into loss of function (LoF; for deficiency, null and dysfunctional alleles with reduced inhibitory activity) and/or gain of function (GoF; for increased polymerization susceptibility and altered inhibitory activity alleles).

Abbreviations: AFR, African/African American; AMR, Latino/Admixed American; ASJ, Ashkenazi Jewish (ASJ); EAS, East Asian; FIN, European–Finnish; NFE, European–non-Finnish; OTH, Other; SAS, South Asian.

To date, sequencing efforts are targeting the exome/genome screening of thousands of individuals from ethnically diverse healthy (controls; eg, gnomAD – The Genome Aggregation Database) and diseased populations (eg, TOPMed – The Trans-Omics for Precision Medicine). These are uncovering an unprecedented number of very low-frequency variants (<0.01%), many of them with predicted negative effects in AAT function, if considering distinct bioinformatics pathogenicity scores (Supplementary Table 1).^2,12 This indicates that we are only aware of the tip of the iceberg regarding AATD causing mutations. In this respect and as discussed in the following sections, whereas polymorphic variants (>1%) and several low-frequency variants tend to be dispersed worldwide, or at least found across wider geographical regions, other rarer mutations are more likely to be restricted to specific populations. While low-frequency variants will probably be detected in control populations of genomic databases, population-specific pathogenic variants are likely to be absent from those datasets and identified mainly through AATD clinical cases, or among related disease groups such as COPD subjects included in TOPMED. From a clinical standpoint, this highlights the importance of not neglecting SERPINA1 rare variants while performing AATD genetic testing, especially when subjects show symptoms and reduced AAT serum levels.

AATD and AAT Variation as a Classical Polymorphism

AATD was firstly discovered by Laurell and Eriksson more than 50 years ago, in 1963, when they noticed the absence of a protein band corresponding to AAT in the electrophoresis strips of plasma samples from several pulmonary emphysema cases.^26,27 Other separating systems with improved band resolution were later developed, and used in multiple studies to assess the genetic polymorphism of blood proteins, including AAT. In general, the alleles of those proteins were labeled as F (fast) and S (slow) according to their electrophoretic mobility in starch, agarose or polyacrylamide gels.^28,29 For AAT in particular, this classic nomenclature of protein variants is still maintained nowadays. Precisely, it relies on the band patterns achieved upon an electrophoretic separation in an isoelectric focusing (IEF) gel with a narrow pH gradient (pH 4–5). There, alleles are classified from “A – Z”, being “A” the faster ones (closer to the anode – the positive end of the IEF gel), “M” those found at the middle region, and “Z” the slowest migrating alleles (Figure 2).^28,30

Figure 2 AAT phenotypes as obtained by IEF in polyacrylamide gels after Coomassie blue staining (left side) and schematic representation of major diagnostic bands of AAT alleles.

Another typical feature of this terminology is correlated with the “Q0” usage to denote null alleles: variants in which no band can be detected in IEF gels. Furthermore, this nomenclature can also admit a linked city name to discriminate alleles with similar electrophoretic mobility, but characterized by distinct mutations (eg, M_Malton, Z_Augsburg or Q0_Ourém). This is more often employed in a rare variant classification where the city name corresponds to the birth-place of the oldest carrier with the identified mutation.²⁸ Currently, the most updated genomic databases follow the Human Genome Variation Society (HGVS) guidelines,³¹ in which the first amino acid is always the initiation methionine. Still, in the AATD field, it is a common rule to employ the residue numbering according to the mature protein (ie, without the first 24 residues of the signal peptide). Thus, for consistency with the former literature, we will use this same labeling of AAT residues, and whenever available we will provide their genomic reference codes.

Common and Rare AAT Variants

Like other genetic variations, AAT alleles can be divided not only according to their association with the AATD condition but also depending on their frequencies across distinct human populations. If using 1% as the threshold for common variants (former assumption of polymorphism), the spectrum of AAT mutations is quite limited, comprising only five variable positions in which just two are linked to AATD. Precisely, these comprise p.Arg101His (rs709932), p.Ala213Val (rs6647), p.Glu264Val (rs17580), p.Glu342Lys (rs28929474) and p.Glu376Asp (rs1303) amino acid replacements that altogether define seven prevalent protein alleles: M1 (Ala or Val), M2, M3, M4, S and Z (Table 2).^28,32 Whereas M1 to M4 are all considered normal alleles placed in the central area of IEF gels; S and Z, are both deficiency alleles associated with fainter bands, located closer to the cathode (Figure 2). Here, the Z allele defined by the p.Glu342Lys (rs28929474) mutation, which replaces a negatively charged residue with a positive one, displays the largest separation from the M zone in IEF gels and also the fainter protein bands (Figure 2). On the other hand, the S allele characterized by a p.Glu264Val (rs17580) mutation changes a negatively charged residue by a neutral one, which results in a protein band with intermediate intensity focusing between M and Z zones (Figure 2).

Table 2 Common SERPINA1 Variants Defining the Most Prevalent AAT Protein Alleles

AAT Allele	Variable Residues					Serum Concentration
	101 rs709932	213 rs6647	264 rs17580	342 rs28929474	376 rs1303
M1Ala	Arg	Ala	Glu	Glu	Glu	100%
M1Val	Arg	Val	Glu	Glu	Glu	100%
M3	Arg	Val	Glu	Glu	Asp	100%
M2	His	Val	Glu	Glu	Asp	100%
M4	His	Val	Glu	Glu	Glu	100%
S	Arg	Val	Val	Glu	Glu	50–60%
Z	Arg	Ala	Glu	Lys	Glu	15%

Note: Alterations in amino acid sequence are shown in italics.

Along with these major protein alleles, there are a large plethora of other variants segregating in human populations. In general, those are found at much lower frequencies (<1%) and thus often referred to as rare AAT alleles, which comprise variants associated with regular AAT activity and serum concentrations, as well as many pathogenic alleles (Table 1).

In overview, common and rare pathogenic variants (or mutations) can be subdivided into two non-mutually exclusively categories according to their negative effects in AAT function and/or activity:

1. Deficiency alleles: This category includes all alleles causing a significant reduction in AAT serum levels. Some authors suggested a 20 uM threshold to consider an allele as causing AATD.⁵ Null alleles (Q0), showing only residual AAT concentrations and no detectable protein band in IEF gels, can also be included in this category. In other words, this group comprises all alleles or mutations that can be connected with a loss-of-function classification.^8,9,19

2. Dysfunctional alleles: In a strict sense, this category includes only those alleles described to affect the inhibitory activity of AAT, in a qualitative or quantitative manner. Nonetheless, in a broader interpretation of dysfunction and fitting the currently used gain-of-function definition, variants increasing AAT susceptibility to polymerization may also fall within this group.^8,9,19

In the following sections, we will revisit some of the collected data about the most relevant pathogenic mutations, by considering not only their clinical importance and functional effects in AAT but also the prevalence among studied AATD cohorts and control populations from the gnomAD database.³³

The Z Allele and p.Glu342Lys Mutation

The Z allele is characterized by a p.Glu342Lys (rs28929474) mutation in an M1Ala background allele (Table 2).^28,32,34 This replacement of a glutamic acid by a lysine residue at a key region of the SERPINs structure – the breach (upper extremity of the shutter domain; Figure 1B and C; Figure 3) – has been experimentally demonstrated to change the conformational properties of AAT. Briefly, the p.Glu342Lys mutation generates the synthesis of unstable monomers, in which the breach and the shutter are opened. This structural modification predisposes the Z allele to dimerize in a process that can be explained by three alternative mechanisms: 1) the loop-sheet rearrangement, where the RCL alone occupies the shutter of another AAT molecule; 2) the β-hairpin configuration, where the RCL together with a contiguous β-strand forms a β-hairpin structure that is inserted into the core of an acceptor AAT molecule; and 3) the C-terminal swap, where the donor AAT undergoes an RCL self-insertion coupled with an abnormal folding of the C-terminal region, which in turn is incorporated in the following AAT molecule. Indeed, it is the steady state of this dimerization process that when carried in a sequential manner leads to the production of the elongated AAT polymers.^{9,18,19,35,36} According to a recent study by Faull et al (2020), it is the latter model of the C-terminal swap that best suits the structural information and high-resolution images of ex vivo polymers obtained from several ZZ liver explants. In addition, when compared with the other mechanisms, the C-terminal swap also predicates for an increased stability and irreversibility of Z allele conformational changes.³⁵

Figure 3 Three-dimensional structure of an active AAT molecule (PDB code: 1QLP) displaying the mutated residues in selected pathogenic AAT alleles. Amino acids affected in variants with increased polymerization susceptibility are shown in red, in dysfunctional alleles are highlighted in blue, in null alleles are displayed in green and combining polymerogenic and null alleles are presented in orange. The three-dimensional images were generated using PyMOL.¹¹²

Notably, although Z polymers are described to aggregate primarily in the ER of hepatocytes, they are suggested to also accumulate in other cell types, such as pneumocytes.^8,9,19,37 Furthermore, Z polymers were found extracellular as well, in the lung interstitium where these are supposed to have a pro-inflammatory activity with repercussions in the progression of emphysema.^9,18,19

This increased susceptibility of the Z allele to polymerize underlies the dual feature of p.Glu342Lys as a loss- and gain-of-function mutation. If, on one hand, it reduces the secretion of AAT into the bloodstream by 85% thus, enhancing the risk of pulmonary disease due to the unbalanced activity of proteases (loss-of-function); on the other hand, the formation of AAT granules that escape the intracellular pathways of protein degradation induces a stress response that may culminate in cell death and fibrosis (gain-of-function).^8,9,19,37 Notably, this represents a critical step not only in liver disease among ZZ individuals^9,38 but also in alveolar cell loss in emphysema, where extracellular polymers found in the interstitium were correlated with another gain-of-function effect that amplifies the damage to lung parenchyma.¹⁹ It is through a chemotactic action that the extracellular polymers stimulate the recruitment and degranulation of neutrophils in the alveolar membranes, maximizing proteolytic activity and spreading the focus of inflammation in the lower lung.^8,9,19,37 On top of this major mechanism of AATD pathophysiology, the p.Glu342Lys mutation per se is also known to reduce the inhibitory capacity of the AAT molecule further enhancing its loss-of-function effects.^3,39

From an epidemiological perspective, the Z allele is more prevalent in Scandinavia and Baltic countries, where its frequency is estimated to vary between 2% and 4%. Accordingly, the proportion of severe AATD cases caused by the ZZ genotype is also relatively high considering their population sizes – about 2000 to 4500 individuals.^10,11,40 The elevated prevalence of the Z allele in Northern Europe, together with its wide distribution across other regions, namely in Southern Europe and Mediterranean areas, is on the genesis of a Viking migration (800–1100 A.D) theory for its dispersal.^28,41 Additionally, hypotheses for a selective advantage of the Z allele were also advanced. For instance, taking into account that Z polymers magnify the inflammatory response of an individual, it could thus have been beneficial in fighting pathogens during a pre-antibiotic era.⁴¹ Another scenario was a possible increase in the reproductive success of ZZ subjects. Although an AAT function in the regulation of acrosin proteolytic activity and sperm-egg fertilization might appear as a valid argument for an improved fitness of the Z allele, this selective hypothesis was later rejected by family-based transmission studies.

More recently, the Z allele was proposed to have been positively selected by height and/or weight adaptive traits to colder climates or conversely, to represent an example of a balanced polymorphism correlated with better lung function values (FEV1: Forced Expiratory Volume in 1 s and FVC: Forced Vital Capacity) in heterozygotes.⁴² Yet, the evidence of extended Z haplotypes collected by North et al 2016, as well as the former arguments for a selective advantage of the p.Glu342Lys, may be questioned by similar findings obtained for other genetics variants. For example, the 32 bp deletion in C-C chemokine receptor 5 gene (CCR5-del32) that confers an effective resistance to HIV infection and is present in North-Europeans at a 16% frequency was proposed to have been selected during the Middle ages, due to a protective effect against the bubonic plague and/or smallpox.^43,44 Nevertheless, CCR5-del32 long-range haplotypes did not differ from other variants found in the same frequency bins, suggesting a likely neutral evolution for this variant.^45–47

Interestingly, the age estimates of CCR5-del32 mutation fall within a 3000–5000 yrs. old range,^45–47 which closely overlaps with the origins of Z allele calculated to be around 2000 to 6000 yrs. ago.^2,34,48,49 For the Z allele, this age variability has been correlated with a large haplotype diversity described in Southern European regions compared to Northern ones.⁴⁹ Indeed, this south-north cline of genetic diversity is not restricted to SERPINA1 or CCR5 mutations, but rather a global European trend as recently shown by genome-wide data. To date, it is well accepted that the genetic variability within Europe mirrors the local geography, and the observed gradients are likely to result from several population movements towards the north, like the ones associated with the expansion of agriculturalists during Neolithic times.⁵⁰ Furthermore, the larger sustained population sizes in Southern European regions, because of their improved long-term survival conditions, are likely to have contributed as well to their larger genetic diversity.⁵⁰ In summary, the updated knowledge of human genetic variation discards, not only a Scandinavian origin of the Z allele but also any linked benefit to its carries. Taking into account the limited mortality of ZZ subjects in their early life, and that major clinical complications occur after post-reproductive ages (>40 yrs. old), often connected with cigarette smoking, AATD is likely to represent a modern challenge to human health.

The S Allele and the p.Glu264Val Mutation

The S allele is defined by a p.Glu264Val (rs17580) mutation in an M1Val background allele (Table 2).^16,28,32 The pathogenic consequences of this allele result from the loss of a hydrogen bond linking the glutamic acid to a tyrosine residue located in the opposed 38 position (Figure 3). This alteration was found to disturb the shutter domain in its periphery affecting the stability of the AAT molecule. An increased turnover of the S protein was initially advanced as the source for its reduced secretion at 50–60% of the normal level.⁵¹ However, the p.Glu264Val mutation has been shown to turn AAT susceptible to polymerization, and to cause its retention in the ER, although at a lesser extent than the Z allele.^51–53 This finding provides an explanation for the S allele being less frequently associated with detectable granules in hepatocytes, a condition that may be shifted in SZ subjects, given that the two variants are capable of forming heteropolymers.^51,53,54 According to recent co-expression assays of Z and S alleles, the first molecules slow down and reduce the intracellular degradation of the misfolded S ones without compromising their secretion.^53,54 Indeed, this phenomenon of heteropolymerization affords a better understanding of the incidence of hepatic disease among SZ subjects, particularly when there is a second hit caused by other risk factors, such as infection, fat liver or alcoholism.^53,55

Oddly, the risk of lung disease among SZ individuals is not yet fully elucidated. This results from some contradictory studies concerning the association of the SZ genotype with pulmonary emphysema and COPD, but also from the small numbers of individuals analyzed so far. Notwithstanding, in smokers, this same genotype seems to increase the risk of COPD, although the disease itself does not differ from what is observed in non-AATD cases.^55–59 On the other hand, the proportion of SZ patients showing serum levels below the protective threshold (11 uM) is limited to about 16%^23,60 and to our knowledge, there are no data concerning the heteropolymerization in the lung, either intra- or extracellularly. At this point, it may be difficult to accurately address the loss- and gain-of-function effects of p.Glu264Val mutation in COPD and emphysema development, and if a second “hit” is required in SZ subjects similar to what happens in hepatic disease. Importantly, this genotype is several times more prevalent worldwide than the severe risk ZZ genotype, especially within Europe where numbers are estimated to vary within the following ranges: 600,000–700,000 for SZ and 80,000–125,000 for ZZ individuals.^10,40,61 Considering that AATD persists as an underdiagnosed condition, SZ subjects could still represent a vulnerable group to respiratory illnesses.

The S allele is more common in Iberia, where frequencies as high as 15% were described in both North and West Iberian regions.^49,62 Previous works also report an increased prevalence of the S allele in Madeira island (Portugal) with an 18% frequency, and in a few African countries, like Angola and Namibia with a prevalence of 18.8% and 14.6%, respectively.^10,11,63 Still, those reports are likely to be correlated with the past demography of the studied groups and their recognized Iberian heritage.^64,65 In this respect, the most updated and unbiased surveys of human genetic variation carried out in African and African-American populations show only minor frequencies of the S allele (<1%; Table 1).

Similar to the Z allele, the p.Glu264Val mutation was proposed to have arisen in the region where the S allele is more common to date – in the Iberian Peninsula. In this case, the 8000–14,000 yrs. age estimates place the origin of the S allele in the Late Paleolithic after the last glacial maximum, in a period when modern humans are thought to have left their refugees in Southern Mediterranean regions.^49,50 Later, Neolithic migrations might have facilitated a first dispersal of the S allele throughout Eurasia and North Africa. More recently, the Iberian diaspora initiated in the XV century with the maritime expansion could have contributed to the spread of the S allele in sub-Saharan Africa and in Central and South America.^50,64,66

The M_Malton, M_Palermo, M_Nichinan, Q0_LaPalma and S_Iiyama Alleles

The five variants M_Malton, M_Palermo, M_Nichinan, Q0_LaPalma and S_Iiyama are connected by their alterations of the contiguous 51–53 residues, where the Phe-Phe-Ser sequence is fundamental to maintain a stable conformation of the shutter domain (Figure 3; Table 1).^67–72 The deletion of a phenylalanine p.Phe52del (rs775982338) found in M_Malton, M_Palermo or M_Nichinan and the replacement of serine by a phenylalanine p.Ser53Phe (rs55819880) as observed in S_Iiyama, cause in a similar way to the Z allele the opening of the breach/shutter, enhancing the susceptibility of these variants to polymerize. However, whereas the S_Iiyama allele is an extremely rare variant identified in a few families from Japan,^67,68,73 the other alleles characterized by the p.Phe52del are widely distributed across different human groups as indicated by the genomics studies of large control populations, where this mutation was found at low or very frequencies (<0.01–0.04%) (Table 1). The discrimination of M_Malton, M_Palermo or M_Nichinan alleles is only achieved by the analysis of their molecular background. The M_Malton has a p.Phe52del mutation linked to an M2 allele and the M_Palermo to an M1Val instead.^69–71 The M_Nichinan allele, on the other hand, has a p.Phe52del in an M1Val background together with another rare mutation that does not lead to AATD, the p.Gly148Arg (V allele – rs112030253).⁶⁸ Finally, in Q0_LaPalma allele, the p.Phe52del is in phase with the p.Glu264Val mutation, which typically defines the S variant. Interestingly, this combination of pathogenic mutations appears to not only impair AAT secretion, as well as its polymerization. No evidence of hepatic dysfunction was detected in Q0_LaPalma carriers.⁷²

Whereas M_Malton and M_Palermo alleles have been identified in several AATD cases across different European and North African countries,^13,74,75 the M_Nichinan was only reported in AATD subjects from Japan.⁷³ So far, only a few studies performed a clinical evaluation of M_Malton carriers, homozygous and heterozygous, confirming the association of p.Phe52del with both pulmonary and hepatic diseases.^76–78 Interestingly, one of those studies reports a case of end-stage liver disease in the presence of a non-deficiency allele, in an M3M_Malton genotype, suggesting a high pathogenicity and heteropolymerization potential of p.Phe52del variant.^53,79 According to a recent work, M_Malton granules found in the liver of AATD patients contain calcium and are mineralized. These modifications were hypothesized to upregulate ER chaperones fostering cell necrosis in a distinct molecular process from one of the Z and S_Iiyama alleles.⁸⁰

The haplotype analysis of M_Malton and M_Palermo alleles disclosed their independent origin, indicating that the p.Phe52del mutation occurred at least twice.² In addition, such data were used to estimate their time to the most common recent ancestor, which in both instances is close to 5000 years old.² This date is consistent with a European and Mediterranean distribution, as well as, a spread associated with Neolithic migrations.⁵⁰ Although the molecular base linked to p.Phe52del is rarely acknowledged, thus not allowing the discrimination of M_Malton and M_Palermo alleles, their association with an M2 or M1Val is not expected to further modify protein structure or disease risk.

Other Rare Variants Leading to AAT Polymerization

There are at least 10 other variants with polymerogenic potential, if taking into account their association with hepatic disease and/or evidence from experimental studies.^12,18 Among those only King’s (p.His334Asp),⁸¹ E_Taurisano (p.Lys368Glu)⁸² and Y_Orzinuovi (p.Pro391His)⁸³ are rare variants absent from public databases (Table 1). Conversely, the I (p.Arg39Cys; rs28931570)⁷¹ and M_Wurzburg (p.Pro369Ser; rs61761869)⁸⁴ alleles are the most prevalent of those low-frequency variants within Europe (0.2% - 0.05%, respectively); the p.Val333Met (rs373630097)¹² is another mutation observed mainly in South Asians at 0.3% (Table 1). Importantly, to the best of our knowledge, only the I and M_Wurzburg alleles have their polymerogenic properties carefully and experimentally validated.^53,54 Whereas the I variant is a functional equivalent of the S allele (the residue 39 faces the 264 one, Figure 3) interfering as well in the periphery of the shutter domain, the M_Wurzburg is closely located to the C-terminal region of the AAT molecule, and described to affect the stability of another key SERPIN domain – the gate.^2,54,83 Moreover, the amino acid replacements of I and M_Wurzburg may also contribute to their ER accumulation by increasing AAT hydrophilicity, although at a lesser extent than Z or M_Malton, given that their association with hepatic disease is dependent on a stronger polymerogenic allele.^53,54,83 Although the low prevalence of these alleles may prevent an accurate evaluation of their impact in AATD manifestations, their detection among symptomatic cases supports the hypothesis of some associated risk of liver and lung diseases.

Dysfunctional AAT Alleles

The best example of a dysfunctional allele is illustrated by the Pittsburgh, which is characterized by a p.Met358Arg substitution (rs121912713), altering the P1 residue of the reactive bond located within the RCL (Figure 1B and Figure 3).⁸⁵ This change causes the AAT molecule to shift its protease affinity towards thrombin instead of its regular anti-elastase inhibitory activity. Interestingly, the clinical outcomes associated with the Pittsburgh variant are not connected with the loss of AAT function, but rather with an abnormal control of coagulation cascades. Indeed, this mutation was first described in a young boy who died of a fatal bleeding disorder, in which the AAT properties as an acute-phase reactant after trauma were proposed to increase its new anti-coagulation potential, consequently causing the patient severe hemorrhages.⁸⁵ This variant is extremely rare and it has been reported in few isolated cases and in a single family, where it always arose as a de novo mutation.⁸⁶ Unsurprisingly, this variant is also absent from control populations currently available through large genomic databases (Table 1).

The F allele is another well-known example of a dysfunctional variant. In this case, if compared with a normal M allele, the F variant displays a considerably reduced inhibitory activity towards neutrophil elastase but also proteinase 3.^39,87 This functional impaired variant results from a p.Arg223Cys (rs28929470) substitution that in spite of being located outside of the RCL is believed to alter the architecture of P1-P1ʹ region.^39,88 This allele has been proposed to increase the risk of emphysema; however, this may only apply to FZ heterozygous, given that no evidence of elastase-induced damage was detected in an FF homozygous patient.³⁹ Importantly, the F allele is a relatively uncommon variant that is found mainly in Europeans, reaching its highest frequency in Sweden (0.9% in gnomAD). Considering that there is also some heterogeneity in the clinical manifestations of AATD cases, the study of additional patients would be fundamental to clarify whether FF subjects are susceptible to lung disease, or not. Still, in no circumstance, the respiratory manifestations would be a consequence of the lower AAT concentrations given that the secretion of F allele is within normal ranges.^39,88 There are recent reports of additional dysfunctional variants but their rarity has prevented so far, a detailed evaluation of their association with AATD disease.⁸⁹

Null Alleles

In most instances, null alleles (Q0) are caused by the premature termination of protein sequence, due to either a nonsense mutation, or the introduction of small insertions and deletions altering the open reading frame. Nonetheless, in a few instances, AAT null alleles can result from standard amino acid substitutions where a misfolded protein is synthetized and rapidly degraded by the ER-associated machinery.^2,90 The M_Heerlen (p.Pro369Leu; rs199422209) and Q0_Lisbon (p.Thr68Ile; rs1490133295) alleles are two such examples, in which a full protein is produced but not detected intracellularly as granules and only poorly secreted from cells, thus being nearly absent from the serum.^71,84,91 In addition, these are likely to be the most prevalent of the null alleles, in populations of European descent and native American admixed (approximately 0.01% and 0.02%, respectively; Table 1), as suggested also by studies of AATD patients.^2,13,74,92 In South Asians, a similar situation occurs with the p.Gly115Ser (rs11558261) replacement, which defines the Q0_Devon and Q0_Newport alleles (0.07%, Table 1).

Splice site changes are another source of null alleles. Interestingly, some of the few cases reported in the literature engage exon IC donor site (Table 1) and an impairment of SERPINA1 expression restricted to hepatocytes.^2,93,94 Given that other cell types use upstream transcription initiation sites located in exons IA and IB, some AAT will be synthetized, but those molecules are far from reestablishing the normal serum levels secreted by the liver. The Q0_Porto (rs775786225), Q0_Faro and Q0_Madrid alleles belong to this category of splice site mutations and were proposed to produce abnormal transcripts targeted by mechanisms of nonsense-mediated mRNA decay (NMD).^2,93,94 The same NMD phenomenon is likely to also be involved in the processing of most truncated variants preventing the accumulation of unusually long or short transcripts within the cells.^2,90,95

Independently of the underlying molecular mechanisms of null alleles, those are by default clinically connected with the loss-of-function effects, and subsequently, with AATD pulmonary manifestations. Here, it is important to strengthen that they will be strictly correlated with the unopposed activity of neutrophil elastase and other proteases. Once there is no formation of AAT polymers, the progression of pulmonary emphysema cannot be exacerbated by a polymer amplified pro-inflammatory response.

The Example of Q0_Ourém, Q0_Mattawa and the p.Leu353Phefs*24 Mutation

The Q0_Ourém, which is characterized by an insertion of a thymine at 353 codon (p.Leu353Phefs*24; rs763023697) in an M3 background allele, is perhaps one of the most prevalent truncated variants of SERPINA1 (Table 1).^93,96 Although it has been identified in only three individuals from control populations (2 Europeans and 1 Latino American), in the Iberian Peninsula it has been associated with several unrelated cases of severe AATD, including several homozygous patients.^74,96–98 The haplotype analysis of Q0_Ourém alleles identified in Portuguese patients confirmed not only a single origin of p.Leu353Phefs*24 mutation but also allowed to estimate its age in 650 yrs.⁹⁶ Remarkably, this timeframe concurs with the major settlements of Ourém village, in central Portugal, which took place during the XIV-XV centuries. Considering that most cases cluster in Ourém and other neighboring areas, a founder effect offers a plausible explanation for its spread in Portugal.⁹⁶

Notably, the p.Leu353Phefs*24 mutation was originally identified in an M1Val background and named as Q0_Mattawa according to the birthplace of the index case, in Ontario, Canada.⁹⁹ This suggested that the p.Leu353Phefs*24 mutation could have occurred at least twice. However, there is recent evidence for the segregation of Q0_Ourém allele in two unrelated subjects from Quebec city population (Quebec, Canada).⁹⁷ Although an independent origin of the two alleles may now be questionable considering the proximity of Mattawa and Quebec cities (700 km) in Canada, the hypothesis may only be ruled out by an updated characterization of Q0_Mattawa haplotypes.

A Comprehensive Analysis of AATD Causing Mutations

The mutational spectrum of SERPINA1 was proposed to be composed of a limited repertoire of genetic variants, in which specific sites were more prone to accumulate events of homoplasy.^2,100 Although this conjecture may still hold true for some pathogenic variants (eg, p.Phe52del, p.Glu264Val and p.Glu342Lys; Table 1), the improvements achieved over the last years in AATD testing, as well as the efforts made in the exome/genome sequencing of thousands of individuals, uncovered unprecedented numbers of SERPINA1 variants present in different populations at low and very low frequencies (<1% and <0.01%).^2,12 In the absence of clinical or experimental data addressing the functional impact of those variants, there are currently multiple predictive scores of variant pathogenicity that offer an alternative bioinformatics approach to evaluate disease-causing mutations.¹⁰¹ In the context of AATD, a recent study has already analyzed the accuracy of those predictive scores to infer the pathogenicity of known and novel AAT mutations.¹² The authors confirmed not only the scores feasibility in the classification of recognized AATD causing variants, as well as, in the detection of new ones later validated by functional assays.¹² Currently, the compiling of variants identified mainly through the screening of AATD patients (Table 1) with those reported in the gnomAD database³³ linked to their pathogenicity scores for a series of predictive algorithms (Supplementary Table 1), shows that SERPINA1 has a vast number of mutations spreading throughout its coding region and splice sites, many of them with expected negative consequences. In overview, SERPINA1 variants exceed the 500, in which more than a third likely have implications for AATD according to clinical data and/or computational predictive tools.^2,12,89

Still and as mentioned above, there are some variants fitting the concept of homoplasy,^2,100 for which the existence of haplotype data is fundamental to address their underlying molecular mechanism. The M_Malton and M_Palermo (p.Phe52del in M2 and M1Val, respectively) alleles occurring in a short repeat sequence (Table 1), are one of the few examples with sustained evidence of mutational recurrence.² The T allele is another known case of homoplasy (p.Glu264Val in an M4⁷¹ or M3 background² instead of the common M1Val) possibly resulting from either crossover or mutation recurrence events.² The Z_Augsburg (p.Glu342Lys in M2 instead of M1Ala), the M_Heerlen (p.Pro369Leu in M1Ala or M1Val), the M_Nichinan (p.Phe52del in a V allele) among other alleles (Table 1) represent additional examples, in which the molecular events behind their origin were not yet elucidated.^2,68,91,102

Although the AATD condition is excepted to be mainly caused by pathogenic mutations, we cannot rule out the influence of their extended genomic background in symptoms variability and clinical outcomes. In this respect, it was demonstrated that the p.Ala213Val variant modifies AAT interactions with neutrophil elastase and lipoproteins,¹⁰³ and on the other hand, several cis-quantitative trait loci (QTL) were shown to affect SERPINA1 expression (eQTL) and splicing (sQTL; Figure 4).

Figure 4 (A) SERPINA1 cis-quantitative trait loci (QTL) view in Genotype-Tissue Expression (GTEx) Locus Browser (https://gtexportal.org/home/). The top section shows the genomic region 1Mb up- and downstream of SERPINA1. In the middle, expression QTLs (eQTL) and splicing QTLs (sQTL) of SERPINA1 are portrayed as bar charts in lung and whole blood with the –log₁₀(p-value) represented in the y-axis. The liver tissue was not available for QTLs analysis. The variant with the most significant p-value – rs2854254 – is highlighted by a red star symbol. The last row shows the functional annotation of the variants according to Ensembl’s Variant Predictor (VEP), Loss-Of-Function Transcript Effect (LOFTEE) and Ensembl Regulatory Build. At the bottom, the plot depicts the pairwise linkage disequilibrium (LD) values (R²) of the QTL variants. TSS, transcription start site; TES, transcription end site. (B) Violin plots of allele-specific cis-sQTL effects of rs2854254 genotypes on SERPINA1 in human lung and whole blood (GTEx database – GTEx Analysis Release v8).

Conclusion

The repertoire of AATD causing mutations has expanded over the last years, and to date, it is acknowledged that along with the major disease risk genotypes, ZZ and SZ, there are many others resulting from different combinations with rare alleles. Indeed, previous studies have shown that rare mutations can be detected in up to 17% of clinically cases^74,75 and like for Z and S alleles, several low-frequency pathogenic variants have been segregating in human populations for hundreds or thousands of years. Therefore, their geographical dispersal could have been facilitated by past human migrations, as well as, by the occurrence of some founders’ effects. Besides these, there are other even rarer variants, which probably arisen recently in time, and thus, are familial or regional confined. Altogether, this data highlights the importance of AATD genetic diagnosis to contemplate more inclusive screenings of SERPINA1 variants, and not only p.Glu264Val and p.Glu342Lys mutations. This goal may be achieved by the traditional AATD laboratorial algorithms combining serum concentration, phenotyping, genotyping and SERPINA1 sequencing, or conversely by the implementation of novel state-of-the-art genotyping/sequencing approaches.

Moreover, it is fundamental to continue unravelling the mechanisms of AATD pathophysiology, for both common and rare pathogenic variants and to define their effects as loss- and/or gain-of-function mutations. Namely, it is imperative to discriminate processes and variants purely associated with an unopposed proteolytic activity, from those stimulating pathways of ER stress and inflammation related to homo- and hetero-polymerization of misfolded alleles. Finally, in a personalized medicine era, it is urgent to solve controversies regarding the disease risk associated with different genotypes (or groups of genotypes for rare variants), and to address the impact in the AATD clinical outcomes of different genetic/genomic modifiers and other factors to which an individual might be exposed (eg smoking, alcoholism, infection, air pollution, obesity).

Acknowledgments

IPATIMUP integrates the i3S Research Unit, which is partially supported by the Portuguese Foundation for Science and Technology (FCT). PIM is supported by the FCT post-doctoral fellowship (SFRH/BPD/120777/2016), financed from the Portuguese State Budget of the Ministry for Science, Technology and High Education and from the European Social Fund, available through the Programa Operacional do Capital Humano.

Disclosure

The authors have no competing interests to declare.